FEINS
(Finite Element Incompressible Navier-Stokes solver)
- The current development version of the code is always
available via SVN. To get access to the SVN, please contact me: Rene Schneider.
- feins___rev2834_2014-09-09__09_25_35.tgz, SVN snapshot revision 2834, last change 9 September 2014.
Main changes since last release:
- P2 elements in 3D (called E2, using p-hierarchical basis)
- Lame in 3D (even with Neumann BC)
test_lame3D
- OpenMP parallelisation of all assembly routines, matrix multiply
and Gauss-Seidel (smoother for multigrid)
- the most time consuming parts of the code scale well,
tested up to 32 threads
- several parts still serial
(most important: mesh refinement, interpolation/restriction)
- overall speedup is good up to 4 threads,
then Amdahls law becomes significant
- if enough processors are available, there is almost no time
difference between P1 and P2 elements with the same total number
of degrees of freedom (both in 2D and 3D)
this required extensive changes in the flexible sparse matrix
storage format
- two eigenvalue solvers in lin_solver.c
- eigenvalue solver (experimental)
- test_lame3Deig
- test_assem (can be switched on/off)
- mutligrid coarse matrix solves now automatically use UMFPACK if
avialable, the old LAPACK band-matrix solver remains as fallback
- changed boundary meshing (prior to call of "triangle" mesh
generator) to be curvature adaptive for shape segments that are
Bezier-curves, such that the deviation of the line segment to the
Bezier-curve is small in comparison to local mesh-width
- feins___rev2691_2013-04-19__15_51_17.tgz, SVN snapshot revision 2691, last change 19 April 2013.
Main changes since last release:
- mainly added 3D support, E1 tetrahedra only for now
- mesh read and .vtk writing of solution
- solution of Poisson problems with test_assem3d.c
- multigrid for E1 is working
- interface with mesh generator gmsh
- can directly read .msh files
- can use a modified gmsh to project new nodes onto geometry
- adaptivity in 3D:
- still only Poisson equation
- ZZ error estimate (generic)
- Baensch green refinement
- several 3D geometries for tests
- Improvements on Navier-Stokes preconditioner:
- implemented Chebyshev semi iteration in the preconditioner for
Stokes problems. This does indeed allow significant savings. Should
investigate this further.
- new interface to hypre->BoomerAMG:
- algebraic multigrid
- requires hypre to be configured and built without MPI
- Option to use AMG inside Preconditioner for navsto
- select FP_VELOCITY_SOLVER (4) in feins_macros.h to use this
- Other improvements:
- improved OMP implementation on some parts
- convection-diffussion code is now OMP parallel
- feins___rev2579_2011-04-29__13_59_22.tgz, SVN snapshot revision 2579, last change 29 April 2011.
Main changes since last release:
- two new examples for optimisation with linear elasticity equation: pedal crank and support with hole
- generic ZZ error estimate and generalised marking: gen_error_est_zz_tx gen_error_est_marker
- improved the solver_settings section of mesh files (.f1m), now allowing various stopping criteria for the mesh refinement among other options
- backward Euler and ROS3P Rosenbrock ODE solver for time dependent convection diffusion
- loads of changes to make feins compile and run in Windows environment. See README.txt for instructions on building under Windows
- implemented usage of OpenMP parallelisation: for now only matrix vector products and assembly of Navier-Stokes matrices, reasonable speedup, looks promising
- feins___rev1824_2010-05-21__13_00_13.tgz, SVN snapshot revision 1824, last change 21 May 2010.
Main changes since last release:
- Adjoint and shape gradient for Lamé equation (linear elasticity, thanks Andreas Günnel).
- Adaptivity for Lamé equation (Bänsch green refinement,
thanks Andreas Günnel, see his Diplom thesis).
- Redesign of directory layout, see also README.TXT.
- Dropped $ARCHI variable approach and old makefiles, we
now use cmake for automatic configuration, see also README.TXT.
- Added UMFPACK as sparse direct solver, with automatic
detection of availability.
- New visualisation solution: paraview. We can now write
.vtk files for use with paraview.
href="feins_080703.tgz">feins_080703.tgz, released 3 July 2008.
- feins_AD_elem_release_100607.tgz, released 7 June 2010, benchmark for AD community, see also "README_AD.txt". The slides of the presentation leading to this benchmark are here.