dmat2.cpp

Go to the documentation of this file.

#include "fvar.hpp"

dmatrix trans(const dmatrix& m1)
{
  dmatrix tmp(m1.colmin(), m1.colmax(), m1.rowmin(), m1.rowmax());

  for (int i = m1.rowmin(); i <= m1.rowmax(); ++i)
  {
    for (int j = m1.colmin(); j <= m1.colmax(); ++j)
    {
      tmp[j][i] = m1[i][j];
    }
  }
  return tmp;
}
double norm(const dmatrix& m1)
{
  double tmp = 0;
  int min = m1.rowmin();
  int max = m1.rowmax();
  const dvector* pm1i = &m1(min);
  for (int i = min; i <= max; ++i)
  {
    tmp += norm2(*pm1i);

    ++pm1i;
  }
  return sqrt(tmp);
}
double norm2(const dmatrix& m1)
{
  double tmp = 0;
  int min = m1.rowmin();
  int max = m1.rowmax();
  const dvector* pm1i = &m1(min);
  for (int i = min; i <= max; ++i)
  {
    tmp += norm2(*pm1i);

    ++pm1i;
  }
  return tmp;
}
double sumsq(const dmatrix& m1) { return norm2(m1); }