ALBERTA:  an Adaptive multi Level finite element toolbox using
          Bisectioning refinement and Error control by Residual
          Techniques for scientific Applications. 

This is the demo package for ALBERTA-2.x with the implementation of
the model problems described in Chapter 2 of the book

             ALBERTA: An Adaptive Finite Element Toolbox

Building
========
The programs are generating by cd'ing into src/Xd/ (where
X stands for 1, 2, 3, 4 or 5) and running "make", e.g.

> cd src/2d/
> make ellipt

where PROGRAM stands for the desired program, see below

Run-time configuration
======================
All programs accept a data-file with run-time parameters; the
data-file has the same name as the program but with a ".dat" suffix
and resides in src/Xd/INIT/. The parameters can be modified by using
your favourite editor, e.g.

> cd src/2d/INIT/
> devilsedit ellipt.dat

The macro-triangulations reside in the diretories src/Xd/Macro/ and
can be modified in the same manner

> cd src/2d/Macro/
> hellsedit macro.amc

The conventional suffix is ".amc", but this is really mere convention.

Short description of the available programs
===========================================

1. Programs that work in 1d/2d/3d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt          -- linear Poisson equation with Dirichlet boundary conditions
heat            -- linear heat-equation with theta-splitting scheme
nonlin          -- non-linear elliptic equation
ellipt-periodic -- linear Poisson equation on a periodic domain (torus
                   or Moebius strip/Klein's bottle)

2. Programs that work in 2d and 3d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt-isoparam -- linear Poisson equation with Neumann boundary conditions and
                   iso-parametric curved boundary approximation up to degree 4
stokes          -- simple stationary Stokes example, with Dirichlet boundary
                   conditions
quasi-stokes    -- Quasi-Stokes problem with stress boundary conditions,
                   for testing block matrices, inhomogeneous Neumann-kind
                   boundary conditions

3. Programs that work in 2d/3d/4d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt-sphere,
ellipt-torus    -- embedded sphere or torus (S^1, S^2, S3; T^1, T^2, T^3).
                   interpolated p.w. polynomial approximation up to degree 4;
                   linear Poisson equation as test-case.

4. Programs that work only in 3d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt-moebius  -- Moebius strip embedded into 3d, explicit parameterisation,
                   piece-wise polynomial interpolated surface up to degree 4;
                   linear Poisson equation as test-case.

5. Programs that work only in 4d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt-klein-bottle -- embedded surface with the topology of a Klein's bottle,
                   otherwise same as the other parametric test-programs.

6. Programs that work only in 5d
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ellipt-klein-3-bottle -- embedded surface, something like a 3d Klein's bottle,
                   i.e. a non-orientable 3-manifold, embedded into \R^5.

Graphical output
================

Depends on what is installed on your system. The parametric
test-programs try to pipe their data to "Geomview"
(http://geomview.sourceforge.net/). Others use optionally the
"gltools" package (http://www.wias-berlin.de/software/gltools/).

For other options see the top-level README of the ALBERTA package.

