Beispiel #1
0
 def __init__(self, name="explicit"):
   """
   Constructor.
   """
   Formulation.__init__(self, name)
   ModuleExplicit.__init__(self)
   self._loggingPrefix = "TSEx "
   self.dtStable = None
   return
Beispiel #2
0
  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
Beispiel #3
0
    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