linear_sys.H File Reference

Solution of systems of linear equations and eigenvalue decomposition. Some are patterned after the code in Numerical Recipes, but with a bit of the fancy C++ template thing happening. More...

#include <cmath>
#include <algorithm>

Include dependency graph for linear_sys.H:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Namespaces
namespace	mlib
Defines
#define	likely(x)   (x)
#define	unlikely(x)   (x)
Functions
template<class T, int N>
bool	mlib::ludcmp (T a[N][N], int indx[N], T *d=NULL)
	LU decomposition.
template<class T, int N>
void	mlib::lubksb (T a[N][N], int indx[N], T b[N])
	Backsubstitution after ludcmp().
template<class T, int N>
bool	mlib::ldltdc (T A[N][N], T rdiag[N])
	Perform LDL^T decomposition of a symmetric positive definite matrix.
template<class T, int N>
void	mlib::ldltsl (T A[N][N], T rdiag[N], T B[N], T x[N])
	Solve Ax=B after ldltdc().
template<int N, class T>
void	mlib::eigdc (T A[N][N], T d[N])
	Eigenvector decomposition for real, symmetric matrices, a la Bowdler et al. / EISPACK / JAMA.
template<int N, class T>
void	mlib::eigmult (T A[N][N], T d[N], T b[N], T x[N])
	x <- A * d * A' * b

---

Detailed Description

Solution of systems of linear equations and eigenvalue decomposition. Some are patterned after the code in Numerical Recipes, but with a bit of the fancy C++ template thing happening.

Note:

This code is taken almost verbatim from the TriMesh2 library written by Szymon Rusinkiewicz of Princeton University (http://www.cs.princeton.edu/gfx/proj/trimesh2/).

Definition in file linear_sys.H.

---

Define Documentation

#define likely (  x   )     (x)

Definition at line 23 of file linear_sys.H.

#define unlikely (  x   )     (x)

Definition at line 24 of file linear_sys.H. Referenced by mlib::ldltdc().

---