be1.c

an introduction to boundary element methods一书源码

/*[ Program Be1	]*/  
#include"cbox1.h"
main()
{float X[102],Y[102],Xm[101],Ym[101],Bc[101],F[101],Xi[21],Yi[21],u[21];
	float q1[21],q2[21],G[100][100],H[100][100],D;
	int Code[101],Dim;
/*[ Dim is a maximum dimension of the system of equations AX = F; it must be at most equal the dimension of (Xm, Ym)  ]*/
	Dim = 100;
/*[ Read input file ]*/
	Input1(X,Y,Xi,Yi,Code,Bc);
/*[ Compute G and H matrices and form the system A X = F ]*/   
	Sys1(X,Y,Xm,Ym,G,H,Bc,F,Code,Dim);  
/*[ Solve  system AX = F  ]*/
	Solve(G,F,&D,N,Dim);
/*[ Compute the potential and fluxes at interior points  ]*/
	Inter1(Bc,F,Code,Xi,Yi,X,Y,u,q1,q2);
/*[ Output solution at the nodes and interior points  ]*/
	Out1(Xm,Ym,Bc,F,Xi,Yi,u,q1,q2);
}