Rheolef: a C++ finite element library for solving PDE

Version : 6.0 Home page : http://ljk.imag.fr/membres/Pierre.Saramito/rheolef/ User's guide: http://ljk.imag.fr/membres/Pierre.Saramito/rheolef/rheolef.pdf

Distibution: sources and debian packages

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

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. * multi-regions and non-constant coefficients. * 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.Saramito@imag.fr Directeur de Recherche CNRS Laboratoire Jean Kuntzmann, Grenoble, France http://www-ljk.imag.fr/membres/Pierre.Saramito