MatVec.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef UTL_MATVEC_H
#define UTL_MATVEC_H
 
#include "matrixnd.h"
 
// Some convenience type definitions:
 
typedef utl::vector<Real>   Vector;
typedef utl::matrix<Real>   Matrix;
typedef utl::matrix3d<Real> Matrix3D;
typedef utl::matrix4d<Real> Matrix4D;
 
typedef std::vector<Real>      RealArray;
typedef std::vector<RealArray> Real2DMat;
typedef std::vector<Real2DMat> Real3DMat;
typedef std::vector<Vector>    Vectors;
typedef std::vector<Matrix>    Matrices;
 
 
namespace utl
{
  Vector operator*(const Vector& X, Real c);
  inline Vector operator*(Real c, const Vector& X) { return X*c; }
  inline Vector operator/(const Vector& X, Real d) { return X*(Real(1)/d); }
 
  Vector operator+(const Vector& X, const Vector& Y);
  Vector operator-(const Vector& X, const Vector& Y);
 
  Matrix operator*(const Matrix& A, Real c);
  inline Matrix operator*(Real c, const Matrix& A) { return A*c; }
 
  inline Real operator*(const Vector& X, const Vector& Y) { return X.dot(Y); }
 
  RealArray operator*(const Matrix& A, const Vector& X);
  RealArray operator*(const Vector& X, const Matrix& A);
 
  Matrix operator*(const Matrix& A, const Matrix& B);
 
  bool transform(Matrix& A, const Matrix& Tn);
 
  bool transform(Vector& V, const Matrix& Tn, bool transpose = false);
 
  bool invert(Matrix& A);
  bool solve(Matrix& A, RealArray& b, std::vector<int>* iPivot = nullptr);
 
  void debugPrint(const char* label, const Vector& V, int level = 2);
}
 
#if HAS_CEREAL
namespace cereal
{
  template<class Archive, class T>
  void load(Archive& archive, utl::vector<T>& vec)
  {
    archive(static_cast<std::vector<T>&>(vec));
  }
 
  template<class Archive, class T>
  void save(Archive& archive, const utl::vector<T>& vec)
  {
    archive(static_cast<const std::vector<T>&>(vec));
  }
}
#endif
#endif