mk3DP1_cube(box,npts);

Parameters:

box - a list of the form [xlo, xhi, ylo, yhi, zlo, zhi], where xlo - zhi are numeric values that specify the range of the three coordinates x,y,z.

npts - a list with three integer entries, specifying the number of nodes in the x,y and z directions.

• Make a uniform 3D grid in a rectangular box. The mesh is saved in the file mesh.scratch.
• After completion the number of nodes and elements are echoed.
• This is included as a demonstration of the 3D element. A coupling to a serious 3D mesh generator is still missing.

Examples:

> # make a 10x10x10 mesh in the unit cube:
> # box:=[xlo, xhi, ylo, yhi, zlo, zhi];
> box:=[0.,1., 0.,1., 0.,1.];
> npts:=[10,10,10];
> mk3DP1_cube(box,npts);
[1331,6000]