Under some compilers, cmath headers do not include the round() functi…

…on. Implement a local version.
hpc4cmb · Nov 1, 2019 · 7a037d8 · 7a037d8
1 parent 7aa399a
commit 7a037d8
Show file tree

Hide file tree

Showing 3 changed files with 17 additions and 4 deletions.
diff --git a/src/libaatm/include/ATMCommon.h b/src/libaatm/include/ATMCommon.h
@@ -39,4 +39,12 @@
 #define ATM_NAMESPACE atm
 #endif
 
+#include <cmath>
+
+ATM_NAMESPACE_BEGIN
+
+    double atm_round(double number);
+
+ATM_NAMESPACE_END
+
 #endif /*!_ATM_COMMON_H_*/
diff --git a/src/libaatm/src/ATMRefractiveIndex.cpp b/src/libaatm/src/ATMRefractiveIndex.cpp
@@ -36,6 +36,11 @@ using std::cout;
 
 ATM_NAMESPACE_BEGIN
 
+  // round function
+  double atm_round(double number) {
+      return number < 0.0 ? ceil(number - 0.5) : floor(number + 0.5);
+  }
+
   // Constructors
 
 
@@ -176,7 +181,7 @@ ATM_NAMESPACE_BEGIN
     if(nu<1.0){
       vp=0;
     }else{
-      vp = (int) round((nu+1.0)/2.0);
+      vp = (int) atm_round((nu+1.0)/2.0);
       vp=vp-1;
     }
     return vp;

diff --git a/src/libaatm/src/ATMRefractiveIndexProfile.cpp b/src/libaatm/src/ATMRefractiveIndexProfile.cpp
@@ -374,14 +374,14 @@ void RefractiveIndexProfile::mkRefractiveIndexProfile()
     if (v_chanFreq_.size()>1){
       if(nc==0){
 	width = fabs(v_chanFreq_[nc+1]-v_chanFreq_[nc])*1e-9;       // width en GHz para ATM
-	npoints=(size_t)round(width*100);                     // One point every 10 MHz
+	npoints=(size_t)atm_round(width*100);                     // One point every 10 MHz
       }else{
 	if(nc==v_chanFreq_.size()-1){
 	  width = fabs(v_chanFreq_[nc]-v_chanFreq_[nc-1])*1e-9;     // width en GHz para ATM
-	  npoints=(size_t)round(width*100);                   // One point every 10 MHz
+	  npoints=(size_t)atm_round(width*100);                   // One point every 10 MHz
 	}else{
 	  width = fabs((v_chanFreq_[nc+1]-v_chanFreq_[nc-1])/2.0)*1e-9;    // width en GHz para ATM
-	  npoints=(size_t)round(width*100);                          // One point every 10 MHz
+	  npoints=(size_t)atm_round(width*100);                          // One point every 10 MHz
 	}
       }
     }else{