out2.c

an introduction to boundary element methods一书源码

#include "cbox2.h"

void Out2(X,Y,F,Bc,Xi,Yi,u,q1,q2)	

float 	X[82],Y[82],F[81],Bc[161],Xi[21],Yi[21],u[21],q1[21],q2[21];
{
		
		FILE *outfile;
		int i,j,k;
	
		outfile = fopen(outname,"a");
		
		fprintf(outfile," \n\n\nSolution: \n\n");
		fprintf(outfile,"Boundary Nodes \n");
		
		fprintf(outfile,"          (X,Y)               u       Prenodal q  Postnodal q\n");
		fprintf(outfile,"(%10.4f, %10.4f) %10.4f %10.4f %10.4f\n",
				X[1],Y[1],F[1],Bc[2*N],Bc[1]); 
		for(i=2;i<=N;i++)
			fprintf(outfile,"(%10.4f, %10.4f) %10.4f %10.4f %10.4f\n",
				X[i],Y[i],F[i],Bc[2*i-2],Bc[2*i-1]); 


		if(L)
		{
			fprintf(outfile,"\n\nInterior Points \n");
			fprintf(outfile,"          (Xi,Yi)            u          q_x          q_y\n");
			for(i=1;i<=L;i++)
				fprintf(outfile,"(%10.4f, %10.4f) %10.4f %10.4f  %10.4f\n",Xi[i],Yi[i],u[i],q1[i],q2[i]);
		}

		fclose(outfile);
}