Who messed with Identity?

melpca.cc

@@ -35,7 +35,7 @@ namespace Melodic{
 	if(opts.paradigmfname.value().length()>0)
 	{
 		basicGLM tmpglm;
-		tmpglm.olsfit(in,melodat.get_param(),Identity(melodat.get_param().Ncols()));
+		tmpglm.olsfit(in,melodat.get_param(),IdentityMatrix(melodat.get_param().Ncols()));
 		std_pca(tmpglm.get_residu(),weights,Corr,tmpE,tmpD);
 //		std_pca(in,weights,Corr,tmpE,tmpD,melodat.get_param());
 	}