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])