*** empty log message ***

miscmaths.cc

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

miscmaths.h

@@ -193,7 +193,8 @@ namespace MISCMATHS {
   ReturnMatrix min(const Matrix& mat);
   ReturnMatrix gt(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix lt(const Matrix& mat1,const Matrix& mat2); 
-  ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2); 
+  ReturnMatrix geqt(const Matrix& mat1,const Matrix& mat2);  
+  ReturnMatrix geqt(const Matrix& mat1,const float a); 
   ReturnMatrix leqt(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix eq(const Matrix& mat1,const Matrix& mat2); 
   ReturnMatrix neq(const Matrix& mat1,const Matrix& mat2);