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])
Esempio n. 2
0
 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])