def setUp(self):
"""
Setup mesh and associated field.
"""
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
from spatialdata.units.Nondimensional import Nondimensional
normalizer = Nondimensional()
normalizer._configure()
from pylith.meshio.MeshIOAscii import MeshIOAscii
importer = MeshIOAscii()
importer.inventory.filename = "data/tri3.mesh"
importer.inventory.coordsys = cs
importer._configure()
self.mesh = importer.read(debug=False, interpolate=False)
from pylith.topology.SolutionFields import SolutionFields
fields = SolutionFields(self.mesh)
fields.add("disp t+dt", "displacement")
fields.solutionName("disp t+dt")
solution = fields.solution()
solution.newSection(solution.VERTICES_FIELD, self.mesh.dimension())
solution.allocate()
solution.zero()
self.fields = fields
self.jacobian = Jacobian(solution)
return
def setUp(self):
"""
Setup mesh and associated field.
"""
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
from spatialdata.units.Nondimensional import Nondimensional
normalizer = Nondimensional()
normalizer._configure()
from pylith.meshio.MeshIOAscii import MeshIOAscii
importer = MeshIOAscii()
importer.inventory.filename = "data/tri3.mesh"
importer.inventory.coordsys = cs
importer._configure()
self.mesh = importer.read(debug=False, interpolate=False)
from pylith.topology.SolutionFields import SolutionFields
fields = SolutionFields(self.mesh)
fields.add("disp t+dt", "displacement")
fields.solutionName("disp t+dt")
solution = fields.solution()
solution.newSection(solution.VERTICES_FIELD, self.mesh.dimension())
solution.allocate()
solution.zero()
self.fields = fields
self.jacobian = Jacobian(solution)
return
class TestSolutionFields(unittest.TestCase):
"""
Unit testing of SolutionFields object.
"""
def setUp(self):
"""
Setup mesh and associated field.
"""
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
from spatialdata.units.Nondimensional import Nondimensional
normalizer = Nondimensional()
normalizer._configure()
from pylith.meshio.MeshIOAscii import MeshIOAscii
importer = MeshIOAscii()
importer.inventory.filename = "data/tri3.mesh"
importer.inventory.coordsys = cs
importer._configure()
self.mesh = importer.read(debug=False, interpolate=False)
self.fields = SolutionFields(self.mesh)
return
def test_constructor(self):
"""
Test constructor.
"""
return
def test_solutionName(self):
"""
Test mesh().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
fields.add("field C", "C");
fields.solutionName("field B")
return
def test_solution(self):
"""
Test solution().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
fields.add("field C", "C");
fields.solutionName("field B")
solution = self.fields.solution()
return
def test_fieldAdd(self):
"""
Test fieldAdd().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
helper_fieldAdd(fields)
fieldA = fields.get("field A")
fieldB = fields.get("field B")
fieldA.allocate()
fieldB.allocate()
fieldA.copy(fieldB)
return
class TestSolutionFields(unittest.TestCase):
"""
Unit testing of SolutionFields object.
"""
def setUp(self):
"""
Setup mesh and associated field.
"""
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
from spatialdata.units.Nondimensional import Nondimensional
normalizer = Nondimensional()
normalizer._configure()
from pylith.meshio.MeshIOAscii import MeshIOAscii
importer = MeshIOAscii()
importer.inventory.filename = "data/tri3.mesh"
importer.inventory.coordsys = cs
importer._configure()
self.mesh = importer.read(debug=False, interpolate=False)
self.fields = SolutionFields(self.mesh)
return
def test_constructor(self):
"""
Test constructor.
"""
return
def test_solutionName(self):
"""
Test mesh().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
fields.add("field C", "C");
fields.solutionName("field B")
return
def test_solution(self):
"""
Test solution().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
fields.add("field C", "C");
fields.solutionName("field B")
solution = self.fields.solution()
return
def test_fieldAdd(self):
"""
Test fieldAdd().
"""
fields = self.fields
fields.add("field A", "A");
fields.add("field B", "B");
helper_fieldAdd(fields)
fieldA = fields.get("field A")
fieldB = fields.get("field B")
fieldA.allocate()
fieldB.allocate()
fieldA.copy(fieldB)
return
