int main(int argc, char **argv){ static xyc *Z; static nde *N; FILE *fp; int i; initop(argc, argv); fp = argf(argc, argv); fp2msh(fp,&Z, &N); pltmsh(stdout,Z,N); }
void Mesh::XP2(FILE *pp, vector<Xyc>&Z,vector<Nde>&N) { if ( fork() == 0 ) { long i; { fprintf(pp,"set border 0\n"); fprintf(pp,"set format x \"\"\n"); fprintf(pp,"set format y \"\"\n"); fprintf(pp,"set noxtics \n"); fprintf(pp,"set noytics \n"); } { for ( i=1; i<(long)Z.size(); i++) fprintf(pp,"set label \"%ld\" at %f,%f\n",i,Z[i].x,Z[i].y); } { for ( i=1; i<(long)N.size(); i++) fprintf(pp,"set label \"%ld\" at %f,%f\n", N[i].C, (Z[N[i].a].x + Z[N[i].b].x)/2.0, (Z[N[i].a].y + Z[N[i].b].y)/2.0); for ( i=1; i<(long)N.size(); i++) fprintf(pp,"set label \"%ld\" at %f,%f\n", N[i].A, (Z[N[i].b].x + Z[N[i].c].x)/2.0, (Z[N[i].b].y + Z[N[i].c].y)/2.0); for ( i=1; i<(long)N.size(); i++) fprintf(pp,"set label \"%ld\" at %f,%f\n", N[i].B, (Z[N[i].c].x + Z[N[i].a].x)/2.0, (Z[N[i].c].y + Z[N[i].a].y)/2.0); } { for ( i=1; i<(long)N.size(); i++) fprintf(pp,"set label \"(%ld)\" at %f,%f\n",i, (Z[N[i].a].x + Z[N[i].b].x + Z[N[i].c].x)/3.0, (Z[N[i].a].y + Z[N[i].b].y + Z[N[i].c].y)/3.0); } fprintf(pp,"plot '-' title \"\" with lines\n"); pltmsh(pp,Z,N); fprintf(pp,"e\n"); fflush(pp); } }