increase tolerance when comparing the "velocity" field

this is only relevant for the semi-implicit models.
This commit is contained in:
Andreas Lauser 2012-10-28 11:52:08 +01:00
parent 5c784f59b9
commit bd2f9dda1f

View File

@ -34,15 +34,15 @@ def isFuzzyEqual(vtkFile1, vtkFile2, absTol, relTol):
number2 = curVals2[i]
if curFieldName.startswith("saturation"):
if abs(number1 - number2) > 1e-3:
print 'Difference between %f and %f too large in data field "%s: %s"'%(number1,number2,curFieldName,abs(number1 - number2))
print 'Difference between %f and %f too large in data field "%s": %s'%(number1,number2,curFieldName,abs(number1 - number2))
return False
elif curFieldName == "velocity":
if abs(number1 - number2) > 0.1:
print 'Difference between %f and %f too large in data field "%s: %s"'%(number1,number2,curFieldName,abs(number1 - number2))
if abs(number1 - number2) > 0.2:
print 'Difference between %f and %f too large in data field "%s": %s'%(number1,number2,curFieldName,abs(number1 - number2))
return False
elif curFieldName.startswith("pressure"):
if abs(number1 - number2) > 0.1 and abs(number1 - number2) > 1e-5*abs(number1 + number2):
print 'Difference between %f and %f too large in data field "%s: %s"'%(number1,number2,curFieldName,abs(number1 - number2))
print 'Difference between %f and %f too large in data field "%s": %s'%(number1,number2,curFieldName,abs(number1 - number2))
return False
elif abs(number1 - number2) > absTol and number2 != 0 and abs(number1/number2 - 1) > relTol:
print 'Difference between %f and %f too large (%f%%) in data field "%s"'%(number1,number2,abs(number1/number2 - 1)*100, curFieldName)