def calculateExteriorElementBoundaryQuadrature(self):
log("initalizing ebqe vectors for post-procesing velocity")
self.u[0].femSpace.elementMaps.getValuesGlobalExteriorTrace(
self.elementBoundaryQuadraturePoints, self.ebqe['x'])
self.u[0].femSpace.elementMaps.getJacobianValuesGlobalExteriorTrace(
self.elementBoundaryQuadraturePoints, self.ebqe['inverse(J)'],
self.ebqe['g'], self.ebqe['sqrt(det(g))'], self.ebqe['n'])
cfemIntegrals.calculateIntegrationWeights(
self.ebqe['sqrt(det(g))'],
self.elementBoundaryQuadratureWeights[('u', 0)],
self.ebqe[('dS_u', 0)])
self.u[0].femSpace.elementMaps.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.elementMaps.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
log("setting flux boundary conditions")
self.fluxBoundaryConditionsObjectsDict = dict([
(cj,
FluxBoundaryConditions(
self.mesh,
self.nElementBoundaryQuadraturePoints_elementBoundary,
self.ebqe[('x')],
self.advectiveFluxBoundaryConditionsSetterDict[cj],
self.diffusiveFluxBoundaryConditionsSetterDictDict[cj]))
for cj in list(
self.advectiveFluxBoundaryConditionsSetterDict.keys())
])
log("initializing coefficients ebqe")
self.coefficients.initializeGlobalExteriorElementBoundaryQuadrature(
self.timeIntegration.t, self.ebqe)
def calculateExteriorElementBoundaryQuadrature(self):
log("initalizing ebqe vectors for post-procesing velocity")
self.u[0].femSpace.elementMaps.getValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['x'])
self.u[0].femSpace.elementMaps.getJacobianValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['inverse(J)'],
self.ebqe['g'],
self.ebqe['sqrt(det(g))'],
self.ebqe['n'])
cfemIntegrals.calculateIntegrationWeights(self.ebqe['sqrt(det(g))'],
self.elementBoundaryQuadratureWeights[('u', 0)],
self.ebqe[('dS_u', 0)])
self.u[0].femSpace.elementMaps.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.elementMaps.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
log("setting flux boundary conditions")
self.fluxBoundaryConditionsObjectsDict = dict([(cj, FluxBoundaryConditions(self.mesh,
self.nElementBoundaryQuadraturePoints_elementBoundary,
self.ebqe[('x')],
self.advectiveFluxBoundaryConditionsSetterDict[cj],
self.diffusiveFluxBoundaryConditionsSetterDictDict[cj]))
for cj in self.advectiveFluxBoundaryConditionsSetterDict.keys()])
log("initializing coefficients ebqe")
self.coefficients.initializeGlobalExteriorElementBoundaryQuadrature(self.timeIntegration.t, self.ebqe)
def calculateExteriorElementBoundaryQuadrature(self):
logEvent("initalizing ebqe vectors for post-procesing velocity")
self.u[0].femSpace.elementMaps.getValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['x'])
self.u[0].femSpace.elementMaps.getJacobianValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['inverse(J)'],
self.ebqe['g'],
self.ebqe['sqrt(det(g))'],
self.ebqe['n'])
cfemIntegrals.calculateIntegrationWeights(self.ebqe['sqrt(det(g))'],
self.elementBoundaryQuadratureWeights[('u', 0)],
self.ebqe[('dS_u', 0)])
self.u[0].femSpace.elementMaps.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.elementMaps.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
logEvent("initializing coefficients ebqe")
self.coefficients.initializeGlobalExteriorElementBoundaryQuadrature(self.timeIntegration.t, self.ebqe)
def calculateExteriorElementBoundaryQuadrature(self):
logEvent("initalizing ebqe vectors for post-procesing velocity")
self.u[0].femSpace.elementMaps.getValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['x'])
self.u[0].femSpace.elementMaps.getJacobianValuesGlobalExteriorTrace(self.elementBoundaryQuadraturePoints,
self.ebqe['inverse(J)'],
self.ebqe['g'],
self.ebqe['sqrt(det(g))'],
self.ebqe['n'])
cfemIntegrals.calculateIntegrationWeights(self.ebqe['sqrt(det(g))'],
self.elementBoundaryQuadratureWeights[('u', 0)],
self.ebqe[('dS_u', 0)])
self.u[0].femSpace.elementMaps.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.elementMaps.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisValuesTraceRef(
self.elementBoundaryQuadraturePoints)
self.u[0].femSpace.getBasisGradientValuesTraceRef(
self.elementBoundaryQuadraturePoints)
logEvent("initializing coefficients ebqe")
self.coefficients.initializeGlobalExteriorElementBoundaryQuadrature(self.timeIntegration.t, self.ebqe)