Open Source Fortran Mathematics Software

Elmer is a finite element software for numerical solution of partial differential equations and multiphysical problems. It includes models of structural mechanics, fluid dynamics, heat transfer, electromagnetics etc. Elmer home is www.csc.fi/elmer

Freemat is an interpreted, matrix-oriented development environment for engineering and scientific applications, similar to the commercial package MATLAB. Freemat provides visualization, image manipulation, and plotting as well as parallel programming.

ATLAS (Automatically Tuned Linear Algebra Software) provides highly optimized Linear Algebra kernels for arbitrary cache-based architectures. ATLAS provides ANSI C and Fortran77 interfaces for the entire BLAS API, and a small portion of the LAPACK AP

The PDL module gives standard perl the ability to COMPACTLY store and SPEEDILY manipulate the large N-dimensional data sets that are the bread and butter of scientific computing.

AUTO is a publicly available software for continuation and bifurcation problems in ordinary differential equations originally written in 1980 and widely used in the dynamical systems community.

librsb is a library for sparse matrix computations featuring the Recursive Sparse Blocks (RSB) matrix format. This format allows cache efficient and multi-threaded (that is, shared memory parallel) operations on large sparse matrices. The most common operations necessary to iterative solvers are available, e.g.: matrix-vector multiplication, triangular solution, rows/columns scaling, diagonal extraction / setting, blocks extraction, norm computation, formats conversion. The RSB format is especially well suited for symmetric and transposed multiplication variants. On these variants, librsb has been found to be faster than Intel MKL's implementation for CSR. Most numerical kernels code is auto generated, and the supported numerical types can be chosen by the user at build time. librsb implements the Sparse BLAS standard, as specified in the BLAS Forum documents.

qrupdate is a Fortran library for fast updates of QR and Cholesky decompositions.

MaxFem is an open software package for electromagnetic simulation by using finite element methods. The package can solve problems in electrostatics, direct current, magnetostatics and eddy-currents. Since version 0.4.0, MaxFEM requires Python 3. We have moved the installers to the MaxFEM website (see below). In order to improve MaxFEM, we will require you to fill out a simple form before downloading them.

Generic data-assimilation toolbox written in java, with native (c and fortran) libraries for high performance computing. Provides tools to couple to your own model and a wide range of algorithms, ranging from parameter calibration to Kalman filters.

Two dimensional (2D) fine mesh finite element (FE) grid editing system. Includes constrained Delaunay triangulation, and automated grid resolution changes based on local attributes. Win32 and Motif GUIs. Mature application, now going open source.

The purpose of this project was to prepare gfortran for integration into gcc. Since this integration was completed in the summer of 2003, this project is now finished. You can find information about gfortran on http://gcc.gnu.org/wiki/GFortran/

Ordinary Differential Equation (ODE) solver for Tcl. Handles both stiff and non-stiff sets of equations; does time-evolved initial value problems and root finding. Based on the popular Lawrence Livermore solver LSODE.

sparse linear algebra library, with emphasis on preconditioned Krylov subspace methods for the solution of linear systems of equations

FFTW++ is a C++ header class for the FFTW Fast Fourier Transform library that automates memory allocation, alignment, planning, wisdom, and communication on both serial and parallel (OpenMP/MPI) architectures. In 2D and 3D, hybrid dealiasing of convolutions substantially reduces memory usage and computation time. Wrappers for C, Python, and Fortran are included.

A Fortran 90 program to model the lattice Boltzmann method

This will be a place to release tools that have not yet been added to SciPy (or will not be added for whatever reason). Currently an interface to ImageMagick is here.

The DUNS (Diagonalized Upwind Navier-Stokes)code is a 2D/3D, structured, multi-block, multi-species,reacting, steady/unsteady, Navier Stokes fluid dynamics code with q-omega turbulence model. It currently uses a diagonalized ADI procedure with upwind diff

High order finite-volume scheme for solving Maxwell's equations on arbitrary 3D domains. This SF project was created mainly for a backup.

Nonlinear Krylov Acceleration (NKA) is a method for accelerating the convergence of fixed-point (Picard) iterations. Many Newton-like and inexact Newton methods are fixed point iterations. The NKA project provides the canonical implementation of the method for several programming languages. The black-box accelerator is simple to integrate into existing code. Placed in the iteration loop, it observes the sequence of solution updates and replaces them with improved updates using information it has gleaned from previous solution iterates. It was only recently recognized (2011) that NLK is essentially equivalent to Anderson Acceleration for a specific choice of mixing parameter. NLK was independently devised in the late 1980's using a very different approach, and though it leads to the same algebraic method, NLK's organization is somewhat different, and arguably superior. The NLK approach also provides clear rationale for the proper choice of Anderson's arbitrary mixing parameter.

Algae is a programming language for numerical analysis. We wrote it because we needed a fast and versatile tool, capable of handling large systems. Algae has been applied to interesting problems in aerospace and related fields for more than a decade.

MatLisp is a package for Common Lisp for handling matrices with real-valued or complex-valued elements. BLAS is used for elementary matrix operations and LAPACK is used for linear algebra routines.

Runge is an Interactive Solver for Systems of Ordinary Differential Equations. It solves initial value problem (aka Cauchy problem). Few advantages Runge provides: It's fast. It utilizes BLAS and LAPACK libraries optimized for modern multi-core processors. It's interactive. It allows you to start a solution by mouse click on a plane. It's precise. It uses Runge Rule to adjust step length to satisfy required precision on each step. It's effective. When it needs to compute derivatives (Jacobian matrix, for example) it does that analytically, i.e. without using numerical methods. It's portable. It works on Windows and Linux 32 and 64 bit and Mac OS X 64 bit. It's open. It allows you to implement and embed your own algorithms (aka "solvers"). It's easy to use. It allows to export results to MS Excel and MATLAB. It's multilingual. It currently supports Englis and Russian languages. It's free. It's distributed under the Boost Software License.

libdspl-2.0 — opensource cross-platform digital signal processing algorithm library, written in C language. Distributed under LGPL v3 license. This allows to use this library in all applications with dynamic linking.

A 2-D inviscid flow and adjoint solver on unstructured triangular grids. It makes use of a vertex-centroid finite volume scheme which is second order accurate. The adjoint solver is developed using the automatic differentiation tool called TAPENADE. Code has been moved to https://github.com/cpraveen/euler2d

Framework for blockstructured adaptive finite volume methods. Provides MPI-parallelized variant of the Berger-Oliger AMR algorithm for Beowulf-clusters. Uses Clawpack. Visualization and conversion tools for HDF4 files included.