vemu06 is a routine to compute a NK-valued function C
at the global nodes of the finite element
mesh. It works for isoparametrical meshes and for the
element libraries PRFLIB=0 and PRFLIB=1. C
depends on the location and the NU-valued input
solution U and its derivative. U is given at the geometrical
nodes of the FEM mesh. vemu06 can be used to interpolate
a solution given on the geometrical
nodes onto the global nodes (set NU:=NK), to set the initial
guess for nonlinear solvers (see veme02) or to set the initial
solution for nonsteady problems (see vemp02) (in both cases NU:=0).
ARGUMENTS
T double precision, scalar, input, global
Real number (e.g. current time).
LCU integer, scalar, input, local
Length of vector LCU, LCU >=LM.
CU double precision, array: CU(LCU), output, local
The function C at the global nodes. CU(i) is the value
of the function C at the global node i+PTRMBK(MYPROC), see
vemdis. Since the
function C is only evaluated for inner elements, it may happen that
any global node gets no contribution from the evaluation
process and CU is undefined at this node. vemu06 sets these entries
to the alternative value 0.
LU integer, scalar, input, local
Length of input solution vector U, LU>=NU*U1.
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, one arbitrarily
selected value of the two is used on both processes.
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 input solution is interpolated by the values of the assigned
geometrical nodes. The polynomials are evaluated at the global
nodes to get the values of U and its derivatives. The value of
function C at a global node is the arithmetic mean of all
evaluations of C on all inner elements containing this node. Nodes
contained in no inner element get the value 0.
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.