From 9ed646afa65e4ad6df28a2ee946fa64ad292659e Mon Sep 17 00:00:00 2001
From: Saad Jbabdi <saad@fmrib.ox.ac.uk>
Date: Tue, 28 Mar 2006 13:15:28 +0000
Subject: [PATCH] *** empty log message ***
---
miscmaths.cc | 18 ++++++++++++++++++
miscmaths.h | 3 ++-
2 files changed, 20 insertions(+), 1 deletion(-)
diff --git a/miscmaths.cc b/miscmaths.cc
index e9a092b..aba20f3 100644
--- a/miscmaths.cc
+++ b/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)
diff --git a/miscmaths.h b/miscmaths.h
index ddd9c48..84164f5 100644
--- a/miscmaths.h
+++ b/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);
--
GitLab