*** empty log message ***

miscmaths.cc

@@ -1455,16 +1455,16 @@ ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2)
   res.Release();
   return res;
 }
-ReturnMatrix geqt(const Matrix& mat1,const float a) 
+ReturnMatrix geqt(const Matrix& mat,const float a) 
 {
-  int ncol = mat.Ncols();
-  int nrow = mat.Nrows();
-  Matrix res(nrow,ncol);
+  int ncols = mat.Ncols();
+  int nrows = mat.Nrows();
+  Matrix res(nrows,ncols);
   res=0.0;
 
   for (int ctr1 = 1; ctr1 <= nrows; ctr1++) {
     for (int ctr2 =1; ctr2 <= ncols; ctr2++) {
-      if( mat1(ctr1,ctr2) >= a){
+      if( mat(ctr1,ctr2) >= a){
 	res(ctr1,ctr2) = 1.0;
       }
     }