VECFEM3 Reference Manual: vemu13

Type: FORTRAN routine

NAME

vemu13 - prints a solution given at the geometrical nodes

SYNOPSIS

CALL VEMU13(: LU, U, NORMU, LIVEM, IVEM, LNEK, NEK, LRPARM, RPARM, LIPARM, IPARM, LDNOD, DNOD, LRDPRM, RDPARM, LIDPRM, IDPARM, LNODN, NODNUM, LNOD, NOD, LNOPRM, NOPARM, LBIG, RBIG, IBIG)
INTEGER: LU, LIVEM, LNEK, LRPARM, LIPARM, LDNOD, LRDPRM, LIDPRM, LNODN, LNOPRM, LBIG
INTEGER: IVEM(LIVEM), NEK(LNEK), IPARM(LIPARM), DNOD(LDNOD), IDPARM(LIDPRM), NODNUM(LNODN), IBIG(*)
DOUBLE PRECISION: U(LU), NORMU(*), RPARM(LRPARM), RDPARM(LRDPRM), NOD(LDNOD), NOPARM(LNOPRM), RBIG(LBIG)

PURPOSE

vemu13 is a routine to print a NU-valued solution U on the first process. U is given at the geometrical nodes on the processes. Additionally the maximum norms NORMU of the components of U are computed over all processors. Only values lower than 10.D0**IVEM(4) are considered.

ARGUMENTS

LU integer, scalar, input, local

Length of the input solution vector, LU>=NDEG*NU.

U double precision, array: U(LU), input, local

The input solution vector at the geometrical nodes. U(U1*(j-1)+i) is the value of the j-th component at the i-th geometrical node on the process MYPROC. If a geometrical node has different values on two processes, an arbitrarily selected value of the two is printed.

NORMU double precision, array: NORMU(NU), output,global

The maximum norm of the solution vector. NORMU(j) is the norm of the j-th component.

LIVEM integer, scalar, input, local

Length of the integer information vector, LIVEM>= MESH+ NINFO.

IVEM integer, array: IVEM(LIVEM), input/output, local/global

Integer information vector.

(1)=MESH, input, local

Start address of the mesh informations in IVEM, MESH>203+ NPROC.

(2)=ERR, output, global

Error number.

0	program terminated without error.
90	`LBIG` is too small.
95	`IVEM` or `U` is too small.
99	fatal error.

(4), input, local

Only entries lower than 10.D0**IVEM(4) are printed and considered in the norm calculation.

(5)=NIVEM, output, local

Used length of IVEM.

(23)=LOUT, input, local

Unit number of the standard output file, normally 6.

(24)=OUTCNT, input, local

Output control flag.

>0	only error messages are printed.
>1	in addition a protocol and the norm of the inputsolution is printed.
>2	in addition the solution is printed on process 1.

(25)=U1, input, local

Leading dimension of the input solution array, U1>=NDEG.

(26)=NU, input, global

Number of components of the input solution.

(200)=NPROC, input, global

Number of processes, see combgn.

(201)=MYPROC, input, local

Logical process id number, see combgn.

(202)=NMSG, input/output, global

Message counter. The difference of the input and the output values gives the number of communications during the vemu13 call.

(204)=TIDS(1), input, global

Begin of the list TIDS which defines the mapping of the logical process ids to the physical process ids. See combgn.

(MESH), input, local

Start of mesh informations, see mesh.

LNEK integer, scalar, input, local

Length of the element array.

NEK integer, array: NEK(LNEK), input, local

Array of the elements, see mesh.

LRPARM integer, scalar, input, local

Length of the real parameter array.

RPARM double precision, array: RPARM(LRPARM), input, local

Real parameter array, see mesh.

LIPARM integer, scalar, input, local

Length of the integer parameter array.

IPARM integer, array: IPARM(LIPARM), input, local

Integer parameter array, see mesh.

LDNOD integer, scalar, input, local

Length of the array of the Dirichlet nodes.

DNOD integer, array: DNOD(LDNOD), input, local

Array of the Dirichlet nodes, see mesh.

LRDPRM integer, scalar, input, local

Length of the real Dirichlet parameter array.

RDPARM double precision, array: RDPARM(LRDPRM), input, local

Array of the real Dirichlet parameters, see mesh.

LIDPRM integer, scalar, input, local

Length of the integer Dirichlet parameter array.

IDPARM integer, array: IDPARM(LIDPRM), input, local

Array of the integer Dirichlet parameters, see mesh.

LNODN integer, scalar, input, local

Length of the array of the id numbers of the geometrical nodes.

NODNUM integer, array: NODNUM(LNODN), input, local

Array of the id numbers of the geometrical nodes, see mesh.

LNOD integer, scalar, input, local

Length of the array of the coordinates of the geometrical nodes.

NOD double precision, array: NOD(LNOD), input, local

Array of the coordinates of the geometrical nodes, see mesh.

LNOPRM integer, scalar, input, local

Length of the array of the node parameters.

NOPARM double precision, array: NOPARM(LNOPRM), input, local

Array of the node parameters, see mesh.

LBIG integer, scalar, input, local

Length of the real work array. A minimal length of LBIG cannot be given. It should be as large as possible.

RBIG double precision, array: RBIG(LBIG), work array, local

Real work array.

IBIG integer, array: IBIG(*), work array, local

Integer work array, RBIG and IBIG have to be defined by the EQUIVALENCE statement.

EXAMPLE

See vemexamples. If you want to check the quality of an FEM approximation of a problem with a well-known solution, you can use vemu05 to compute the error of the approximation at the geometrical nodes. vemu13 computes the norm of error over all processes.

REFERENCES

[FAQ], [DATAMAN], [DATAMAN2], [P_MPI].

COPYRIGHTS

by L. Grosz, Auckland , 6. June, 2000.