Changeset 19573

Timestamp:

09/21/15 16:06:40 (10 years ago)

Author:

schlegel

Message:

BUG: need to add abs to geometry check

Location:

issm/trunk-jpl/src/m/classes

Files:

• geometry.m (modified) (1 diff)
• geometry.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/geometry.m

@@ -70 +70 @@
                                 md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
                                 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
-                                if any((self.thickness-self.surface+self.base)>10^-9),
+                                if any(abs(self.thickness-self.surface+self.base)>10^-9),
                                         md = checkmessage(md,['equality thickness=surface-base violated']);
                                 end 

issm/trunk-jpl/src/m/classes/geometry.py

@@ -49 +49 @@
                 md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'size',[md.mesh.numberofvertices])
                 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices],'>',0,'timeseries',1)
-                if any((self.thickness-self.surface+self.base)>10**-9):
+                if any(abs(self.thickness-self.surface+self.base)>10**-9):
                         md.checkmessage("equality thickness=surface-base violated")
                 if solution==TransientSolutionEnum() and md.transient.isgroundingline: