def calculateNumericalDiffusion(self, q):
for ci in range(1, self.nc):
if numpy.isnan(q[('pdeResidual', ci)]).any():
import pdb
print "NaN's in res"
pdb.set_trace()
cshockCapturing.calculateNumericalDiffusionResGradQuad(
self.shockCapturingFactor, self.mesh.elementDiametersArray,
q[('pdeResidual', ci)], q[('grad(u)', ci)], self.numDiff[ci])
if numpy.isnan(self.numDiff[ci]).any():
import pdb
print "NaN's in numDiff"
pdb.set_trace()
def calculateNumericalDiffusion(self,q):
for ci in range(1,self.nc):
if numpy.isnan(q[('pdeResidual',ci)]).any():
import pdb
print "NaN's in res"
pdb.set_trace()
cshockCapturing.calculateNumericalDiffusionResGradQuad(self.shockCapturingFactor,
self.mesh.elementDiametersArray,
q[('pdeResidual',ci)],
q[('grad(u)',ci)],
self.numDiff[ci])
if numpy.isnan(self.numDiff[ci]).any():
import pdb
print "NaN's in numDiff"
pdb.set_trace()