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
