def preinitialize(self, mesh): """ Do pre-initialization setup. """ from pylith.mpi.Communicator import mpi_comm_world comm = mpi_comm_world() if 0 == comm.rank: self._info.log("Pre-initializing fault '%s'." % self.label()) FaultCohesive.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) ModuleFaultCohesiveKin.quadrature(self, self.faultQuadrature) for eqsrc in self.eqsrcs.components(): eqsrc.preinitialize() ModuleFaultCohesiveKin.eqsrcs(self, self.eqsrcs.inventory.facilityNames(), self.eqsrcs.components()) for name in self.eqsrcs.inventory.facilityNames(): self.availableFields['vertex']['info'] += ["final_slip_%s" % name] self.availableFields['vertex']['info'] += ["slip_time_%s" % name] if mesh.dimension() == 2: self.availableFields['vertex']['info'] += ["strike_dir"] elif mesh.dimension() == 3: self.availableFields['vertex']['info'] += ["strike_dir", "dip_dir"] return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ BoundaryCondition.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ from pylith.mpi.Communicator import mpi_comm_world comm = mpi_comm_world() if 0 == comm.rank: self._info.log("Pre-initializing fault '%s'." % self.label()) FaultCohesive.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) ModuleFaultCohesiveDyn.quadrature(self, self.faultQuadrature) if mesh.dimension() == 2: self.availableFields['vertex']['info'] += ["strike_dir"] elif mesh.dimension() == 3: self.availableFields['vertex']['info'] += ["strike_dir", "dip_dir"] if not isinstance(self.tract, NullComponent): self.tract.preinitialize(mesh) self.availableFields['vertex'][ 'info'] += self.tract.availableFields['vertex']['info'] self.availableFields['vertex']['info'] += \ self.friction.availableFields['vertex']['info'] self.availableFields['vertex']['data'] += \ self.friction.availableFields['vertex']['data'] return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ from pylith.mpi.Communicator import mpi_comm_world comm = mpi_comm_world() if 0 == comm.rank: self._info.log("Pre-initializing fault '%s'." % self.label()) FaultCohesive.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) ModuleFaultCohesiveDyn.quadrature(self, self.faultQuadrature) if mesh.dimension() == 2: self.availableFields['vertex']['info'] += ["strike_dir"] elif mesh.dimension() == 3: self.availableFields['vertex']['info'] += ["strike_dir", "dip_dir"] if not isinstance(self.tract, NullComponent): self.tract.preinitialize(mesh) self.availableFields['vertex']['info'] += self.tract.availableFields['vertex']['info'] self.availableFields['vertex']['info'] += \ self.friction.availableFields['vertex']['info'] self.availableFields['vertex']['data'] += \ self.friction.availableFields['vertex']['data'] return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ BoundaryCondition.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) self.bcQuadrature.preinitialize(mesh.coordsys().spaceDim()) self.quadrature(self.bcQuadrature) self.createSubMesh(mesh) return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ self._info.log("Pre-initializing fault '%s'." % self.label()) FaultCohesive.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) ModuleFaultCohesiveTract.quadrature(self, self.faultQuadrature) if mesh.dimension() == 2: self.availableFields['cell']['info'] += ["strike_dir"] elif mesh.dimension() == 3: self.availableFields['cell']['info'] += ["strike_dir", "dip_dir"] return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ BoundaryCondition.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) self.bcQuadrature.preinitialize(mesh.coordsys().spaceDim()) self.quadrature(self.bcQuadrature) self.createSubMesh(mesh) self.output.preinitialize(self) fields = [] if not isinstance(self.inventory.dbInitial, NullComponent): fields += ["initial_value"] if not isinstance(self.inventory.dbRate, NullComponent): fields += ["rate_of_change", "rate_start_time"] if not isinstance(self.inventory.dbChange, NullComponent): fields += ["change_in_value", "change_start_time"] self.availableFields['cell']['info'] += fields return
def preinitialize(self, mesh): """ Do pre-initialization setup. """ from pylith.mpi.Communicator import mpi_comm_world comm = mpi_comm_world() if 0 == comm.rank: self._info.log("Pre-initializing fault '%s'." % self.label()) FaultCohesive.preinitialize(self, mesh) Integrator.preinitialize(self, mesh) ModuleFaultCohesiveImpulses.quadrature(self, self.faultQuadrature) if mesh.dimension() == 2: self.availableFields['vertex']['info'] += ["strike_dir"] elif mesh.dimension() == 3: self.availableFields['vertex']['info'] += ["strike_dir", "dip_dir"] return