be2.c

an introduction to boundary element methods一书源码

#include "cbox2.h"

main()
{
	float X[82],Y[82],F[81],Bc[161],Xi[21],Yi[21],u[21],D;
	int Code[161];
	float G[81][162],H[81][81],q1[21],q2[21];
	int Twice,Dim;
	
	/*[Dim= Maximum dimension of the system AX = F; it is the same as the
	 maximum number of nodes, or maximum number of elements ]*/
		
	Dim = 80;
	Twice = 2*Dim;
	
	
	/*[	Read Data	]*/
	
	Input2(Xi,Yi,X,Y,Code,Bc);
	
	/*[ Compute matrices G and H, and form the system A X = F ]*/
	
	
	Sys2(X,Y,G,H,F,Bc,Code,Dim,Twice);

		
	/*[Solve the system AX = F ]*/
	
	Solve(H,F,&D,N,Dim);	

/*[ Compute the potential and fluxes at interior points ]*/


	Inter2(F,Bc,Code,Xi,Yi,X,Y,u,q1,q2); 
	
	/*[ Output the solution at boundary nodes and interior points ]*/
	
	Out2(X,Y,F,Bc,Xi,Yi,u,q1,q2);
}