out1.c

an introduction to boundary element methods一书源码

#include "cbox1.h"
void Out1(Xm,Ym,Bc,F,Xi,Yi,u,q1,q2)
	float Xm[101],Ym[101],Bc[101],F[101],Xi[21],Yi[21];
	float u[21],q1[21],q2[21];
	{
		FILE *outfile;
		int i,j,k;
	
		outfile = fopen(outname,"a");
		
		fprintf(outfile,"Solution: \n");
		fprintf(outfile," \nBoundary Nodes \n\n");
		fprintf(outfile,"    Nodes (X,Y)			u	q\n\n");
		for(i=1;i<=N;i++)
		fprintf(outfile,"(%10.4f, %10.4f) %10.4f %10.4f\n",
		Xm[i],Ym[i],Bc[i],F[i]);
	if(L)
	{
		fprintf(outfile,"\n\nInterior Points (Xi,Yi)  	u	  q_x	    q_y\n\n");
		for(k=1;k<=L;k++)
		fprintf(outfile,"(%10.4f, %10.4f) %10.4f %10.4f %10.4f\n",
				Xi[k],Yi[k],u[k],q1[k],q2[k]);
	} 
	
	fclose(outfile);
}