Source Codes

Some useful subroutines that I've written in FORTRAN95.


Approximate Riemann Solvers:

  
Flux:
The "flux_" routines compute the convective flux across a face given the left and right conserved states.
(click name for downloading the source code)

Jacobian: (linearization w.r.t state):
The "jac_" routines compute the flux jacobian matrices of the "flux_" routines for use in Newton
solvers. The routines have the ability to produce both exact and inexact Jacobians and
can therefore be used in adjoint codes also.
(click name for downloading the source code)


Geometry and Mesh Generation Tools:

dtmeshgen:
Program that computes the Delaunay Triangulation of the input set of points. Useful for quickly generating
simple unstructured meshes when a set of points is available. (dtmeshgen.tar.bz2)



                    Set of arbitrary points                                                          Delaunay triangulation of set of points
                                                           ( click on images for more detail )




         Movie of triangulation (avi - 27 MB)




Cubic spline routines:
Set of routines that can be used for interpolation using piecewise cubic fit of data.

spline_coefficients.f90 : computes the cubic spline coefficients for the given set of data

solve_tridiag.f90         : This routine is required for 'spline_coefficients.f90'. It solves a tri-diagonal linear
                                    system using the Thomas algorithm.

eval.f90                      : FORTRAN function that interpolates the dependent variable given the
                                    spline coefficients and the independent variable value


© Karthik Mani - 2008