The ultimate implicit solver for modelling:

• 3D axi-symmetric,
• strongly time-dependent,
• quasi-stationary and steady state,
• compressible,
• weakly compressible,
• radiating,
• magnetized,
• self-gravitating,
• dissipative and diffusive,
• relativistic and
• general-relativistic flows

A non-linear Newton-iterative numerical tool for solving:

• Euler,
• Newtonian Navier-Stokes,
• General Relativistic Navier-Stokes equations,
• the 5D axi-symmetric radiative transfer equation,
• the Poisson equation

Discretization method:

• 3D axi-symmetrical spherical grid
• Conservative formulation (Finite Volume)
• First and second order temporal accuracy
• up to 3rd order spatial accuracy

Properties of the solution method:

• Various degrees of implicitness:
  implicit, semi-implicit, semi-explicit and explicit methods
• Hierarchical Solution Scenario:
  enables gradual coupling of the equations
  

Implicit iterative solution procedure:

• Black-White Line Gauss-Seidel
• Approximate Factorization Method
• Krylov Subspace Iterative Methods:
  GMRES, Bi-CGSTAB, TF-QMR, ...

Basic properties of the simulation code:

• Fortran 90/95
• parallelized using MPI
• PETSc library for Krylov methods
• netCDF data I/O
  (serial and parallel-netCDF)
• visualization scripts using python/matplotlib

Although the program itself is quite complex, it will be clearly structured and easy to use, but still very powerful and flexible.