def __init__(self, name="explicit"): """ Constructor. """ Formulation.__init__(self, name) ModuleExplicit.__init__(self) self._loggingPrefix = "TSEx " self.dtStable = None return
def _reformJacobian(self, t, dt):
"""
Reform Jacobian matrix for operator.
"""
from pylith.mpi.Communicator import mpi_comm_world
comm = mpi_comm_world()
self._debug.log(resourceUsageString())
if 0 == comm.rank:
self._info.log("Integrating Jacobian operator.")
self._eventLogger.stagePush("Reform Jacobian")
self.updateSettings(self.jacobian, self.fields, t, dt)
ModuleExplicit.reformJacobianLumped(self)
self._eventLogger.stagePop()
if self.viewJacobian:
self.jacobian.view("Jacobian")
self._debug.log(resourceUsageString())
return
def _reformJacobian(self, t, dt):
"""
Reform Jacobian matrix for operator.
"""
from pylith.mpi.Communicator import mpi_comm_world
comm = mpi_comm_world()
self._debug.log(resourceUsageString())
if 0 == comm.rank:
self._info.log("Integrating Jacobian operator.")
self._eventLogger.stagePush("Reform Jacobian")
self.updateSettings(self.jacobian, self.fields, t, dt)
ModuleExplicit.reformJacobianLumped(self)
self._eventLogger.stagePop()
if self.viewJacobian:
self.jacobian.view("Jacobian")
self._debug.log(resourceUsageString())
return