dmat17.cpp

Go to the documentation of this file.

#include "fvar.hpp"

dmatrix identity_matrix(int min, int max)
{
  dmatrix identity(min, max, min, max);

#ifndef SAFE_INITIALIZE
  identity.initialize();
#endif

  for (int i = min; i <= max; ++i)
  {
    identity.elem(i, i) = 1;
  }

  return identity;
}