Rheolef 6.2: an efficient FEM C++ finite element library for solving PDE

[Saramito Pierre]

Rheolef: an efficient FEM C++ finite element library for solving PDE

Version : 6.2
Home: http://ljk.imag.fr/membres/Pierre.Saramito/rheolef
Book: http://cel.archives-ouvertes.fr/docs/00/71/37/60/PDF/rheolef.pdf

Distibution: sources and binaries as debian packages

Keywords: finite elements, numerical simulation, partial derivative
equations,
C++, meshes, graphics

News:
* nonlinear solvers improved (see p-laplace example)
* equations on a surface: implements three diferent FEM methods
* improves the high order Pk Lagrange interpolation implementation
* ports on intel c++ 12.0 and gnu c++ 4.7 new compiler versions

Previous features:
Rheolef is a programming environment that serves as a convenient
laboratory for
computations involving finite element methods (FEM) for solving
partial
differential equations (PDE). Rheolef is both a C++ library and a set
of
commands for unix shell programming, providing algorithms and data
structures.

* Algorithms refer to the most up-to-date ones: preconditioned sparse
solvers
for linear systems, incompressible elasticity, Stokes and Navier-
Stokes flows,
characteristic method for convection dominated heat problems, etc.
Also
nonlinear generic algorithms such as fixed point and damped Newton
methods.

* Data structures fit the standard variational formulation concept:
spaces,
discrete fields, bilinear forms are C++ types for variables, that
can be
combined in any expressions, as you write it on the paper.

Combined together, as a Lego game, these bricks allows the user to
solve most
complex nonlinear problems. The concision and readability of codes
written
with Rheolef is certainly a major keypoint of this environment.

Main features

* [NEW] Massively distributed memory finite element environment, based
on MPI.
* [NEW] High-order polynomial approximation.
* Poisson problems in dimension d=1,2,3.
* Stokes problems (d=2,3), with Taylor-Hood or stabilized P1 bubble-P1
elements.
* linear elasticity (d=1,2,3), including the incompressible case.
* characteristic method for time-dependent problems:
transport, convection-difusion, and Navier-Stokes equations.
* input and output in various file format for meshes generators and
numerical
data visualization systems.

Advanced features

* auto-adaptive mesh algorithms.
* axisymetric problems.
* nonlinear problems with either fixed-point algorithms or a provided
generic
damped Newton solver.
* 3d stereo visualization

Both reference manual and users guide are available.

The license is GPL.

Pierre Saramito
--
Pierre....@imag.fr
Directeur de Recherche CNRS
Laboratoire Jean Kuntzmann, Grenoble, France
http://www-ljk.imag.fr/membres/Pierre.Saramito