Projects

Armadillo C++ Library

Armadillo is a C++ linear algebra library (matrix maths) aiming towards a good balance between speed and ease of use. Integer, floating point and complex numbers are supported, as well as a subset of trigonometric and statistics functions. Various matrix decompositions are provided through optional ... [More] integration with LAPACK and ATLAS libraries. A delayed evaluation approach is employed (during compile time) to combine several operations into one and reduce (or eliminate) the need for temporaries. This is accomplished through recursive templates and template meta-programming. This library is useful if C++ has been decided as the language of choice (due to speed and/or integration capabilities), rather than another language like Matlab or Octave. [Less]

optimal

In my geophysical work, I rely heavily upon a generalized inversion algorithm called Gauss-Newton optimization. The Gauss-Newton method minimizes objective functions that can be iteratively approximated by quadratics. This approach is particularly appropriate for least-squares inversions of ... [More] moderately non-linear transforms. The package also contains code for conjugate-gradient and line-search optimizations. Over two decades I have implemented this package four times, in four different programming languages. Every feature of this package has been motivated by practical application. This particular Java version is only a few years old, but I have already used it for a dozen or more distinct inversion problems. I find no preexisting package so suitable for typical geo [Less]

taucs

TAUCS is a C library of solvers for sparse linear algebra.

gaia-physics

Analyzes accelerations caused by gravity, such as those recorded by GRACE. GAIA inverts these accelerations using a mascon technique to determine the mass distribution that matches the acceleration data. Currently, mass distributions can either be static in time, vary linearly, as well as ... [More] oscillate at tidal and seasonal frequencies. These time dependencies are solved simultaneously. [Less]