vemu03 is a routine to compute a NC-valued function C
at the center points of the inner elements of the finite element
mesh. 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. vemu03 can
be used to compute the heat flow in thermal analysis
problems and the stresses in structural analysis problems. The results
can be handed over to a standard postprocessor which
processes element results (see veid99, vepa99).
ARGUMENTS
T double precision, scalar, input, global
Real number (e.g. current time).
LCU integer, scalar, input, local
Length of vector CU, LCU>=NC*NINE, where NINE is the total
number of inner elements on the process MYPROC.
CU double precision, array: CU(LCU), output, local
The values of the function C at the center points of the inner elements. The
NC components of the function C are stored
separately for each individual group.
CU(ADDCU(G)-1+CU1(G)*(j-1)+i)
is the value of the j-th component of C at the center point of
the i-th element on the process MYPROC in group G. The numbers
ADDCU and CU1 are returned in IVEM.
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.
The center point coordinates are computed by the transformation of
the center point of the reference element to the element by
interpolation of the assigned geometrical nodes. To compute the
values of U and its derivatives on the center point 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 Universitaet Karlsruhe 1989-1996.
Copyrights by Lutz Grosz 1996.
All rights reserved. More details see VECFEM.
by L. Grosz, Auckland , 6. June, 2000.