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)