vemu07 is a routine to compute a NC-valued function C
at the geometrical nodes of the finite element mesh
using extrapolation from the integration nodes. It works for
isoparametrical meshes and for the element libraries PRFLIB=0 and
PRFLIB=1. C depends on the location and the NK-valued input
solution U and its derivative. U is given at the global nodes of
the FEM mesh. vemu07 can be used to interpolate a solution given
on the global nodes onto the geometrical nodes. The results
can be handed over to a standard postprocessor which
processes nodal results (see veid97, vepa97).
ARGUMENTS
T double precision, scalar, input, global
Real number (e.g. current time).
LCU integer, scalar, input, local
Length of vector CU, LCU>=NC*CU1.
CU double precision, array: CU(LCU), output, local
The values of the function C at the geometrical nodes on the process
MYPROC.
CU(CU1*(j-1)+i) is the value of the j-th component
of C at the geometrical node i, where i refers to the local numbering on the
process. Since the function C is only evaluated for inner
elements it may happen that
a geometrical node gets no contribution from the evaluation
process and CU is undefined at this node. vemu07 sets these entries
to the alternative value 10.D0**IVEM(4).
LU integer, scalar, input, local
Length of solution vector U, LU >=LM.
U double precision, array: U(LU), input, local
The input solution vector at the global nodes. U(i) is the value
of the input solution at the global node i+PTRMBK(MYPROC), see
vemdis.
LIVEM integer, scalar, input, local
Length of the integer information vector,
LIVEM>= MESH+ NINFO.
Length of the real work array. The needed length of LBIG
depends on the given mesh. 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.
USERC external
Name of the subroutine in which the
function C is described, see userc.
The value of the function C at a geometrical node is the
arithmetic mean of all extrapolation of C on the
integration nodes at the elements
containing this geometrical node. The extrapolation is
done by a least square approach. Geometrical nodes contained in
no inner element get the value 10.D0**IVEM(4). To compute the
values of U and its derivatives at a geometrical node, the values of
U at the global nodes are interpolated by a polynomial.
Program by L. Grosz, C. Roll, P. Sternecker, 1989-1996.
Copyrights by Lutz Grosz 1999.
All rights reserved. More details see VECFEM.
by L. Grosz, Auckland , 6. June, 2000.