def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return for material in ["elastic", "viscoelastic"]: ncells = self.mesh['ncells-%s' % material] self.mesh['ncells'] = ncells filename = "%s-%s_info.h5" % (self.outputRoot, material) from axialdisp_soln import p_mu, p_lambda, p_density propMu = p_mu * numpy.ones((1, ncells, 1), dtype=numpy.float64) propLambda = p_lambda * numpy.ones( (1, ncells, 1), dtype=numpy.float64) propDensity = p_density * numpy.ones( (1, ncells, 1), dtype=numpy.float64) properties = { 'mu': propMu, 'lambda': propLambda, 'density': propDensity } from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return ncells= self.mesh['ncells'] filename = "%s-elastic_info.h5" % self.outputRoot m = self.density*self.vs**2 l = self.density*self.vp**2 - 2*m propMu = m*numpy.ones( (1,ncells, 1), dtype=numpy.float64) propLambda = l*numpy.ones( (1,ncells, 1), dtype=numpy.float64) propDensity = self.density*numpy.ones( (1,ncells, 1), dtype=numpy.float64) properties = {'mu': propMu, 'lambda': propLambda, 'density': propDensity} from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return for material in ["elastic", "viscoelastic"]: ncells= self.mesh['ncells-%s' % material] self.mesh['ncells'] = ncells filename = "%s-%s_info.h5" % (self.outputRoot, material) from axialdisp_soln import p_mu,p_lambda,p_density propMu = p_mu*numpy.ones( (1, ncells, 1), dtype=numpy.float64) propLambda = p_lambda*numpy.ones( (1, ncells, 1), dtype=numpy.float64) propDensity = p_density*numpy.ones( (1, ncells, 1), dtype=numpy.float64) properties = {'mu': propMu, 'lambda': propLambda, 'density': propDensity} from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if self.reader is None: return from pylith.tests.PhysicalProperties import check_properties from rigidbody_soln import p_mu,p_lambda,p_density self.mesh['ncells'] = self.mesh['ncells_elastic'] ncells= self.mesh['ncells'] filename = "%s-elastic_info.vtk" % self.outputRoot propMu = p_mu*numpy.ones( (ncells, 1), dtype=numpy.float64) propLambda = p_lambda*numpy.ones( (ncells, 1), dtype=numpy.float64) propDensity = p_density*numpy.ones( (ncells, 2), dtype=numpy.float64) properties = {'mu': propMu, 'lambda': propLambda, 'density': propDensity} check_properties(self, filename, self.mesh, properties) self.mesh['ncells'] = self.mesh['ncells_viscoelastic'] ncells= self.mesh['ncells'] filename = "%s-viscoelastic_info.vtk" % self.outputRoot propMu = p_mu*numpy.ones( (ncells, 1), dtype=numpy.float64) propLambda = p_lambda*numpy.ones( (ncells, 1), dtype=numpy.float64) propDensity = p_density*numpy.ones( (ncells, 2), dtype=numpy.float64) properties = {'mu': propMu, 'lambda': propLambda, 'density': propDensity} check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return ncells = self.mesh['ncells'] filename = "output/%s-statevars_info.h5" % self.outputRoot from rigidrotate_soln import p_mu, p_lambda, p_density propMu = p_mu * numpy.ones((1, ncells, 1), dtype=numpy.float64) propLambda = p_lambda * numpy.ones((1, ncells, 1), dtype=numpy.float64) propDensity = p_density * numpy.ones( (1, ncells, 1), dtype=numpy.float64) properties = { 'mu': propMu, 'lambda': propLambda, 'density': propDensity } from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return ncells = self.mesh["ncells"] filename = "%s-elastic_info.h5" % self.outputRoot from axialdisp_soln import p_mu, p_lambda, p_density propMu = p_mu * numpy.ones((1, ncells, 1), dtype=numpy.float64) propLambda = p_lambda * numpy.ones((1, ncells, 1), dtype=numpy.float64) propDensity = p_density * numpy.ones((1, ncells, 1), dtype=numpy.float64) properties = {"mu": propMu, "lambda": propLambda, "density": propDensity} from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return
def test_elastic_info(self): """ Check elastic info. """ if not self.checkResults: return ncells= self.mesh['ncells'] filename = "output/%s-statevars_info.h5" % self.outputRoot from rigidrotate_soln import p_mu,p_lambda,p_density propMu = p_mu*numpy.ones( (1, ncells, 1), dtype=numpy.float64) propLambda = p_lambda*numpy.ones( (1, ncells, 1), dtype=numpy.float64) propDensity = p_density*numpy.ones( (1, ncells, 1), dtype=numpy.float64) properties = {'mu': propMu, 'lambda': propLambda, 'density': propDensity} from pylith.tests.PhysicalProperties import check_properties check_properties(self, filename, self.mesh, properties) return