def scalarDomainIntegral(dV,nValueArray,nElements=None):
if useC:
if nElements is None:
partialSum = cfemIntegrals.scalarDomainIntegral(dV,nValueArray,nValueArray.shape[0])
else:
partialSum = cfemIntegrals.scalarDomainIntegral(dV,nValueArray,nElements)
else:
partialSum=0.0
range_nQuadraturePoints_element = range(nValueArray.shape[1])
for eN in range(nValueArray.shape[0]):
for k in range_nQuadraturePoints_element:
partialSum += nValueArray[eN,k]*dV[eN,k]
return globalSum(partialSum)
def scalarDomainIntegral(dV, nValueArray, nElements=None):
if useC:
if nElements == None:
partialSum = cfemIntegrals.scalarDomainIntegral(
dV, nValueArray, nValueArray.shape[0])
else:
partialSum = cfemIntegrals.scalarDomainIntegral(
dV, nValueArray, nElements)
else:
partialSum = 0.0
range_nQuadraturePoints_element = range(nValueArray.shape[1])
for eN in range(nValueArray.shape[0]):
for k in range_nQuadraturePoints_element:
partialSum += nValueArray[eN, k] * dV[eN, k]
return globalSum(partialSum)
def globalScalarDomainIntegral(abs_det_J,quadratureWeightArray,nValueArray):
integral = 0.0
if useC:
integral = cfemIntegrals.scalarDomainIntegral(abs_det_J,quadratureWeightArray,nValueArray)
else:
integral=0.0
range_nQuadraturePoints_element = range(nValueArray.shape[1])
for eN in range(nValueArray.shape[0]):
for k in range_nQuadraturePoints_element:
integral += nValueArray[eN,k]*quadratureWeightArray[k]*abs_det_J[eN,k]
return globalSum(integral)
def globalScalarDomainIntegral(abs_det_J, quadratureWeightArray, nValueArray):
integral = 0.0
if useC:
integral = cfemIntegrals.scalarDomainIntegral(abs_det_J,
quadratureWeightArray,
nValueArray)
else:
integral = 0.0
range_nQuadraturePoints_element = range(nValueArray.shape[1])
for eN in range(nValueArray.shape[0]):
for k in range_nQuadraturePoints_element:
integral += nValueArray[
eN, k] * quadratureWeightArray[k] * abs_det_J[eN, k]
return globalSum(integral)
