out5.c

an introduction to boundary element methods一书源码

#include "cbox5.h"

void Out5(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," \nResults \n\n");
	
	fprintf(outfile,"\n            (X,Y)              u         q\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\n\n		Interior Points \n\n");
		fprintf(outfile,"          (Xi,Yi)              u         q_x       q_y\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);
}