#*******************************************************************
#**
#** v e m b l d e x m 0 1
#**
#** two coupled Poisson equation on the 2-dimensional unit cube.
#** The mesh is read from a vecfem input file (e.g. generated by
#** vemgen2d).
#**
#** by L. Grosz Karlsruhe, Jan. 1995
#**
#*******************************************************************
#**
#** The data set of this examples has two parts (search for
#** 'cut here'). The first part specifies the problem
#** (please copy it to 'vembldexm01.equation') and the second part
#** defines the control parameters (please copy it to
#** 'vembldexm01.resource'). The FORTRAN code for the solution
#** of the problem is generated by entering
#** 'vembuild vembldexm01' into your shell.
#**
#*******************************************************************
####### cut here to get vembldexm01.equation <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#** The problem is a system of two coupled 2-dimensional Poisson
#** equations with Dirichlet and Neuman boundary conditions.
#**
#** The domain is [0,1] x [0,1], where the mesh is givem by
#** rec.vem or is generated by vemgen2d. The mesh uses
#** triangular elements of order 2 and line elements of order 2.
#** Dirichlet conditions are set on all boundary points for the
#** first component, but there is no Dirichlet condition for
#** the second component.
#**
#*******************************************************************
#**
#** The functions b1, b2, f1, f2, g1 and g2 are selected, so that
#**
#** u1 = 1
#** u2 = 2
#**
#** gets the exact solution of the problem.
#**
b1=1
b2=2
f1=10*b1+b2
f2=b1+10*b2
g1=0
g2=0
#**
#*******************************************************************
#**
#** The Dirichlet conditions:
#**
u1=b1
u2=b2
#**
#** 'u2=b2' is not necessary, since no Dirichlet conditions
#** for the second component are specified in the input mesh.
#**
#*******************************************************************
#**
#** the functional equation :
#**
area{v1x1*u1x1+v1x2*u1x2+(10.*u1+ u2)*v1 - f1*v1
+ v2x1*u2x1+v2x2*u2x2+( u1+10.*u2)*v2 - f2*v2}
+ line{v1*g1+v2*g2} =0
#**
#*******************************************************************
>>>>>>> cut here to get vembldexm01.resource <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#** The problem has a two dimensional domain and two solution
#** component:
#**
DIM=2
NK=2
#**
#*******************************************************************
#**
#** One processor with maximal 20 Mbytes are used. Maximal 1000
#** nodes and 1000 elements are allowed:
#**
PROCESS_STORAGE=20
PROCESS_MAXNN=1000
PROCESS_MAXNE=1000
#**
#*******************************************************************
#**
#** the is read from the file rec.vem. other meshes can be
#** created vemgen2d.
#**
MESH_PREP=print
MESH_FILEIN=rec.vem
#**
#*******************************************************************
#**
#** The solution components are printed into the file solution.out
#** and the indicator, which considers the error of all components
#** is printed into file error.out.
#**
MESH_POSTP=print
OUTPUT_INDEX=11
OUTPUT_FILE=solution.out
OUTPUT_ERRINDEX=11
OUTPUT_ERRFILE=error.out
#**
#*******************************************************************