def __init__(self, name="quadrature"): """ Constructor. """ QuadratureBase.__init__(self, name) ModuleQuadrature.__init__(self) return
def __init__(self, name="quadrature"):
"""
Constructor.
"""
QuadratureBase.__init__(self, name)
ModuleQuadrature.__init__(self)
return
def _initialize(self, cell):
"""
Initialize C++ quadrature object.
"""
import numpy
from pylith.utils.utils import sizeofPylithScalar
size = sizeofPylithScalar()
if 8 == size:
ModuleQuadrature.initialize(self, cell.basis, cell.basisDeriv,
cell.quadPts, cell.quadWts,
cell.geometry.spaceDim())
elif 4 == size:
ModuleQuadrature.initialize(self, numpy.float32(cell.basis),
numpy.float32(cell.basisDeriv),
numpy.float32(cell.quadPts),
numpy.float32(cell.quadWts),
cell.geometry.spaceDim())
else:
raise ValueError("Unknown size for PylithScalar")
return
def _initialize(self, cell):
"""
Initialize C++ quadrature object.
"""
import numpy
from pylith.utils.utils import sizeofPylithScalar
size = sizeofPylithScalar()
if 8 == size:
ModuleQuadrature.initialize(self, cell.basis,
cell.basisDeriv,
cell.quadPts,
cell.quadWts,
cell.geometry.spaceDim())
elif 4 == size:
ModuleQuadrature.initialize(self, numpy.float32(cell.basis),
numpy.float32(cell.basisDeriv),
numpy.float32(cell.quadPts),
numpy.float32(cell.quadWts),
cell.geometry.spaceDim())
else:
raise ValueError("Unknown size for PylithScalar")
return
