Identity to IdentityMatrix

miscmaths.cc

@@ -446,25 +446,6 @@ namespace MISCMATHS {
 	return a*a + b*b + c*c;
     }
 
-  int Identity(Matrix& m)
-    {
-      Tracer tr("Identity");
-      m=0.0;
-      for (int j=1; j<=m.Nrows(); j++)
-	m(j,j)=1.0;
-      return 0;
-    }
-
-  ReturnMatrix Identity(int num)
-    {
-      Tracer tr("Identity");
-      Matrix eye(num,num);
-      Identity(eye);
-      eye.Release();
-      return eye;
-    }
-
-
   int diag(Matrix& m, const float diagvals[])
     {
       Tracer tr("diag");
@@ -954,8 +935,6 @@ namespace MISCMATHS {
       Matrix rotmat(3,3);
       rotmat = aff3 * scales.i() * (skew.SubMatrix(1,3,1,3)).i();
       ColumnVector transl(3);
-      //transl = affmat.SubMatrix(1,3,4,4);
-      //transl = transl - (Identity(3) - rotmat)*centre;
       transl = affmat.SubMatrix(1,3,1,3)*centre + affmat.SubMatrix(1,3,4,4)
 	         - centre;
       for (int i=1; i<=3; i++)  { params(i+3) = transl(i); }

miscmaths.h

@@ -121,7 +121,6 @@ namespace MISCMATHS {
   float norm2sq(float a, float b, float c);
 
   int Identity(Matrix& m);
-  ReturnMatrix Identity(int num);
   ColumnVector seq(const int num);
 
   int diag(Matrix& m, const float diagvals[]);