def calculateNumericalDiffusion(self,q):
for ci in range(self.nc):
cshockCapturing.calculateNumericalDiffusionHJ(self.shockCapturingFlag,
self.shockCapturingFactor,
self.mesh.elementDiametersArray,
q[('pdeResidual',ci)],
q[('mt',ci)],
q[('H',ci)],
self.numDiff[ci])
def calculateNumericalDiffusion(self,q):
for ci in range(self.nc):
cshockCapturing.calculateNumericalDiffusionHJ(self.shockCapturingFlag,
self.shockCapturingFactor,
self.mesh.elementDiametersArray,
q[('pdeResidual',ci)],
q[('mt',ci)],
q[('H',ci)],
self.numDiff[ci])