MilAMinMillion A Minute
Original milamin.org banner: the 1e6/1 road-sign logo on a finite element mesh MILAMINFast Finite Element Solver in MATLAB

The original MILAMIN

The vintage site, preserved: everything on these pages – the slideshow, the applications gallery, the technical notes and the downloads – is original milamin.org content, as it served from 2008.

MILAMIN, the native MATLAB package, set up, solved and post-processed one million degrees of freedom per minute on a desktop – from vectorized blocked assembly to Cholesky with fill-reducing reordering and Powell–Hestenes iterations for incompressible Stokes. Developed at PGP (Physics of Geological Processes, University of Oslo) by Marcin Dabrowski, Marcin Krotkiewski and Dani Schmid. The method paper: Dabrowski, M., Krotkiewski, M. & Schmid, D. W. (2008), MILAMIN: MATLAB-based finite element method solver for large problems, Geochemistry, Geophysics, Geosystems 9, Q04030, doi:10.1029/2007GC001719.

MILAMIN is a finite element method implementation in native MATLAB that is capable of doing one million degrees of freedom per minute on a modern desktop computer. This includes pre-processing, solving, and post-processing. The MILAMIN strategies and package are applicable to a broad class of problems in Earth science.

Pressure and stress fields around circular inclusions
Pressure- and stress-distributions in materials embedding circular inclusions: application of MILAMINs mechanical solver.
Folded multilayer stack
Folding of a multilayer: application of MILAMINs incompressible Stokes solver in studying problems from structural geology.
Settling elliptical inclusions
Elliptical inclusions in a box subjected to vertical gravity field: application of MILAMINs mechanical solver.
Thermomechanical convection collage
Thermomechanical convection: application of MILAMINs thermal and mechanical solvers to study mantle dynamics.
Mixing of two materials
Mixing of two materials with different densities: application of MILAMINs mechanical solver.
Crystal plasticity collage
Crystal plasticity: application of MILAMINs mechanical solver.
Pressure around elliptical inclusions
Pressure-distribution around elliptical inclusions: application of MILAMINs mechanical solver.
Multilayer folding collage
Folding of multilayers: application of MILAMINs incompressible Stokes solver in studying problems from structural geology.
Settling circular inclusions
Circular inclusions in a box subjected to vertical gravity field: application of MILAMINs mechanical solver.
Mantle plumes
Thermomechanical convection: application of MILAMINs thermal and mechanical solvers to study mantle dynamics.
Porous convection with cracks
Porous Convection: applications of MILAMINs thermal and mechanical solvers.
Stress fields around circular inclusions
Pressure- and stress-distributions in materials embedding circular inclusions: application of MILAMINs mechanical solver.
Folded multilayers
Folding of multilayers: application of MILAMINs incompressible Stokes solver in studying problems from structural geology.
Convection of a heterogeneous fluid with a dense layer
Thermomechanical convection of a heterogeneous fluid: application of MILAMIN to study mantle dynamics.
Sedimentation velocity field from the applications gallery

Applications gallery

What MILAMIN was used for: rigid inclusions, multilayer folding, sedimentation, thermomechanical convection, mixing, suspensions – with the original simulation videos, still streaming after fifteen years.

Sparse assembly timing comparison from the technical notes

Technical notes

The implementation notes behind the speed: sparse assembly, boundary conditions, fill-reducing reordering, Cholesky factorization and substitution, with the original benchmark figures.

The MILAMIN 1e6/1 road-sign logo

Downloads

Getting MILAMIN: the sources and paper, the MUTILS high-performance utilities – sparse_create, quadtree, tsearch2, einterp, SpMV – and the external packages Triangle, SuiteSparse and METIS.

The people behind MILAMIN

The people who contributed to the development of MILAMIN:

The development of MILAMIN was made possible through the support of PGP: Physics of Geological Processes, University of Oslo – today part of the Njord Centre.