[cig-commits] r7971 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-09-14 16:11:27 -0700 (Fri, 14 Sep 2007)
New Revision: 7971

Modified:
   mc/3D/CitcomS/trunk/lib/Obsolete.c
   mc/3D/CitcomS/trunk/lib/Phase_change.c
Log:
Replace abs() by fabs()

Modified: mc/3D/CitcomS/trunk/lib/Obsolete.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Obsolete.c	2007-09-14 23:10:59 UTC (rev 7970)
+++ mc/3D/CitcomS/trunk/lib/Obsolete.c	2007-09-14 23:11:27 UTC (rev 7971)
@@ -75,7 +75,7 @@
           }
 
       if (cof[3][3]!=0.0)
-        jacobian = sqrt(abs(determinant(cof,E->mesh.nsd)))/cof[3][3];
+        jacobian = sqrt(fabs(determinant(cof,E->mesh.nsd)))/cof[3][3];
 
       dGammax->vpt[k] = jacobian;
 

Modified: mc/3D/CitcomS/trunk/lib/Phase_change.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Phase_change.c	2007-09-14 23:10:59 UTC (rev 7970)
+++ mc/3D/CitcomS/trunk/lib/Phase_change.c	2007-09-14 23:11:27 UTC (rev 7971)
@@ -101,7 +101,7 @@
 
 void phase_change_apply_410(struct All_variables *E, double **buoy)
 {
-  if (abs(E->control.Ra_410) > 1e-10) {
+  if (fabs(E->control.Ra_410) > 1e-10) {
     phase_change_apply(E, buoy, E->Fas410, E->Fas410_b, E->control.Ra_410,
 		       E->control.clapeyron410, E->viscosity.z410,
 		       E->control.transT410, E->control.width410);
@@ -112,7 +112,7 @@
 
 void phase_change_apply_670(struct All_variables *E, double **buoy)
 {
-  if (abs(E->control.Ra_670) > 1e-10)
+  if (fabs(E->control.Ra_670) > 1e-10)
     phase_change_apply(E, buoy, E->Fas670, E->Fas670_b, E->control.Ra_670,
 		       E->control.clapeyron670, E->viscosity.zlm,
 		       E->control.transT670, E->control.width670);
@@ -122,7 +122,7 @@
 
 void phase_change_apply_cmb(struct All_variables *E, double **buoy)
 {
-  if (abs(E->control.Ra_cmb) > 1e-10)
+  if (fabs(E->control.Ra_cmb) > 1e-10)
     phase_change_apply(E, buoy, E->Fascmb, E->Fascmb_b, E->control.Ra_cmb,
 		       E->control.clapeyroncmb, E->viscosity.zcmb,
 		       E->control.transTcmb, E->control.widthcmb);