POSEIDON is a collection of C++ and FORTRAN subroutines for solving Partial Differential Equations using Finite Element Methods. The goal of this collection is to provide with basic finite element tools. The current version of POSEIDON makes available the mesh generation and the numerical linear algebra modules for testing purposes. This modules consists of subroutines for the generation of two-dimensional triangular meshes and various numerical methods for the numerical solution of linear systems with dense and sparse matrices. The methods implemented in these subroutines for the mesh generation are based on Bowyer-Watson algorithm for the initial Delaunay triangulation and a combination of Chew's and Ruppert's refinement algorithms for the generation of good quality triangulations using maximum edge and minimum angle criteria. Other subroutines for global and local mesh refinement are also provided. Screenshots of triangulations and their quality produced with the help of POSEIDON are presented in the pictures.

The library has been tested with gfortran and g95 compilers on OS X and Debian Linux. 

To report bugs, send suggestions and comments please contact me.

The current version is: VERSION 0.0

Download the FORTRAN module for mesh generation and the testing programs here 

Download the C++ module for the solution of sparse and dense linear systems

Download the C++ module for the mesh generation 

Users' Guide