Identity to IdentityMatrix

glim.cc

@@ -81,8 +81,7 @@ namespace FILM {
       pinv_x = pinv(x);
 
       // R = I - x*pinv(x)
-      Matrix I(sizeTS, sizeTS);
-      Identity(I);
+      Matrix I(IdentityMatrix(sizeTS));
 
       RMat = I - x*pinv_x;
       

glimGls.cc

@@ -27,11 +27,7 @@ namespace FILM {
     sigmaSquareds(numTS),
     dof(sizeTS - numParams)
   {
-    I.ReSize(sizeTS);
-    for(int i=1; i<=sizeTS; i++)
-      {
-	I(i,i)=1;
-      }
+    I=IdentityMatrix(sizeTS);
   }
   
   void GlimGls::CleanUp()

glm.cc

@@ -32,8 +32,7 @@ namespace FILM {
       Computeb();
 
       // r = (I - x(x'x)-1x')y
-      Matrix I(sizeTS, sizeTS);
-      Identity(I);      
+      Matrix I(IdentityMatrix(sizeTS));
       r = (I-d*inv_xx*d.t())*(y->Column(1));
 	
       return r;