/* nag_mesh2d_inc (d06aac) Example Program. * * Copyright 2001 Numerical Algorithms Group. * * Mark 7, 2001. */ #include #include #include #include #define EDGE(I,J) edge[3*((J)-1)+(I)-1] #define CONN(I,J) conn[3*((J)-1)+(I)-1] #define COOR(I,J) coor[2*((J)-1)+(I)-1] int main(void) { const Integer nvmax=250; double coef, power; Integer exit_status, i1, i, itrace, k, nedge, nelt, nv, nvb, reftk; Nag_Boolean smooth; NagError fail; char pmesh[2]; double *bspace=0, *coor=0; Integer *conn=0, *edge=0; INIT_FAIL(fail); exit_status = 0; Vprintf("nag_mesh2d_inc (d06aac) Example Program Results\n\n"); /* Skip heading in data file */ Vscanf("%*[^\n] "); /* Reading of the geometry */ Vscanf("%ld%ld%*[^\n] ", &nvb, &nedge); if (nvb > nvmax) { Vprintf("Problem with the array dimensions\n"); Vprintf(" nvb nvmax %6ld%6ld\n", nvb, nvmax); Vprintf(" Please increase the value of nvmax\n"); exit_status = -1; goto END; } /* Allocate memory */ if ( !(bspace = NAG_ALLOC(nvb, double)) || !(coor = NAG_ALLOC(2*nvmax, double)) || !(conn = NAG_ALLOC(3*(2*nvmax-1), Integer)) || !(edge = NAG_ALLOC(3*nedge, Integer))) { Vprintf("Allocation failure\n"); exit_status = -1; goto END; } /* Coordinates of the boundary mesh vertices and boundary edges */ for (i = 1; i <= nvb; ++i) { Vscanf("%ld", &i1); Vscanf("%lf", &COOR(1,i)); Vscanf("%lf", &COOR(2,i)); Vscanf("%*[^\n] "); } for (i = 1; i <= nedge; ++i) { Vscanf("%ld", &i1); Vscanf("%ld", &EDGE(1,i)); Vscanf("%ld", &EDGE(2,i)); Vscanf("%ld", &EDGE(3,i)); Vscanf("%*[^\n] "); } Vscanf(" ' %1s '%*[^\n]", pmesh); /* Initialise mesh control parameters */ for (i = 0; i < nvb; ++i) bspace[i] = 0.05; smooth = Nag_TRUE; itrace = 0; coef = 0.75; power = 0.25; /* Call to the mesh generator */ /* nag_mesh2d_inc (d06aac). * Generates a two-dimensional mesh using a simple * incremental method */ nag_mesh2d_inc(nvb, nvmax, nedge, edge, &nv, &nelt, coor, conn, bspace, smooth, coef, power, itrace, 0, &fail); if (fail.code == NE_NOERROR) { if (pmesh[0] == 'N') { Vprintf(" nv =%6ld\n", nv); Vprintf(" nelt =%6ld\n", nelt); } else if (pmesh[0] == 'Y') { /* Output the mesh to view it using the NAG Graphics Library */ Vprintf(" %10ld%10ld\n", nv, nelt); for (i = 1; i <= nv; ++i) Vprintf(" %12.6e %12.6e\n", COOR(1,i), COOR(2,i)); reftk = 0; for (k = 1; k <= nelt; ++k) { Vprintf(" %10ld%10ld%10ld%10ld\n", CONN(1,k), CONN(2,k), CONN(3,k), reftk); } } else { Vprintf("Problem with the printing option Y or N\n"); exit_status = -1; goto END; } } else { Vprintf("Error from nag_mesh2d_inc (d06aac).\n%s\n", fail.message); exit_status = 1; goto END; } END: if (bspace) NAG_FREE(bspace); if (coor) NAG_FREE(coor); if (conn) NAG_FREE(conn); if (edge) NAG_FREE(edge); return exit_status; }