LDL.h

#ifndef MATH_LDL_H
#define MATH_LDL_H

#include "matrix.h"

namespace Math {

template <class T>
struct LDLDecomposition
{
  typedef MatrixTemplate<T> MatrixT;
  typedef VectorTemplate<T> VectorT;

  LDLDecomposition();
  LDLDecomposition(const MatrixT& A);

  void set(const MatrixT& A);
  bool backSub(const VectorT& b, VectorT& x) const;
  void LBackSub(const VectorT& b, VectorT& x) const;
  void LTBackSub(const VectorT& b, VectorT& x) const;
  bool DBackSub(const VectorT& b, VectorT& x) const;
  bool backSub(const MatrixT& B, MatrixT& X) const;
  bool getInverse(MatrixT& Ainv) const;
  void getPseudoInverse(MatrixT& Ainv) const;
  void getL(MatrixT& L) const;
  void getD(VectorT& d) const;
  void getA(MatrixT& A) const;
  void mulL(const Vector& x,Vector& y) const;
  void mulLT(const Vector& x,Vector& y) const;
  void mulD(const Vector& x,Vector& y) const;

  void update(const VectorT& x);
  bool downdate(const VectorT& x);

  MatrixT LDL;
  T zeroTolerance;
  int verbose;
};

} // namespace Math

#endif