*** empty log message ***

glimGls.cc

@@ -36,10 +36,12 @@ namespace FILM {
       // compute r
       Matrix I(sizeTS, sizeTS);
       Identity(I);
-      r = (I-x*inv_xx*x.t())*y;
+      Matrix R = I-x*inv_xx*x.t();
+      r = R*y;
 
       // compute sigma squareds 
-      sigmaSquareds(ind) = (r.t()*r/sizeTS).AsScalar();
+      //      sigmaSquareds(ind) = (r.t()*r/sizeTS).AsScalar();
+      sigmaSquareds(ind) = (r.t()*r).AsScalar()/R.Trace();
 
       // set corrections
       SetCorrection(inv_xx, ind);