Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21843)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21844)
@@ -10,4 +10,5 @@
 
 #include <stdio.h>
+#include <math.h>
 #include "../cores/cores.h"
 #include "../shared/io/io.h"
@@ -2443,6 +2444,6 @@
 			y = xyz_list[3*i+1];
 			
-			bx=-150.-728.8*std::pow(x/300000.,2)+343.91*std::pow(x/300000.,4)-50.57*std::pow(x/300000.,6);
-			by=500./(1.+std::exp((-2./4000.)*(y-80000./2.-24000.)))+500./(1.+std::exp((2./4000.)*(y-80000./2.+24000.)));
+			bx=-150.-728.8*pow(x/300000.,2)+343.91*pow(x/300000.,4)-50.57*pow(x/300000.,6);
+			by=500./(1.+exp((-2./4000.)*(y-80000./2.-24000.)))+500./(1.+exp((2./4000.)*(y-80000./2.+24000.)));
 			
 			r[i] = -720.;
@@ -3408,9 +3409,10 @@
 				ftxy+=(tauxy[n]-smoothedtauxy[elem_vertices[n]])*basis[n];
 			}
-			error+=Jdet*gauss->weight*( std::pow(ftxx,2)+std::pow(ftyy,2)+std::pow(ftxy,2) ); //e^2 
+			error+=Jdet*gauss->weight*( pow(ftxx,2)+pow(ftyy,2)+pow(ftxy,2) ); //e^2 
 		}
 		/*Set the error in the global vector*/	
       sid=element->Sid();
-		velementerror->SetValue(sid,std::sqrt(error),INS_VAL);//sqrt( e^2 )
+		error = sqrt(error);//sqrt(e^2)
+		velementerror->SetValue(sid,error,INS_VAL);
 		/*Cleanup intermediaries*/
 		xDelete<IssmDouble>(xyz_list);
@@ -3558,9 +3560,10 @@
             fdHdy+=(GradH[1]-smoothed_dHdy[elem_vertices[n]])*basis[n];
          }
-         error+=Jdet*gauss->weight*( std::pow(fdHdx,2)+std::pow(fdHdy,2) ); //e^2
+         error+=Jdet*gauss->weight*(pow(fdHdx,2)+pow(fdHdy,2) ); //e^2
       }
       /*Set the error in the global vector*/
       sid=element->Sid();
-      velementerror->SetValue(sid,std::sqrt(error),INS_VAL);//sqrt( e^2 )
+		error = sqrt(error); //sqrt( e^2 )
+      velementerror->SetValue(sid,error,INS_VAL);
       /*Cleanup intermediaries*/
       xDelete<IssmDouble>(xyz_list);