def test_order1(self): dom = RectangularDomain(esys.dudley) dom.order = 1 run(dom, sys.stdout)
def test_order2(self): dom = RectangularDomain() dom.order = 2 run(dom, sys.stdout)
def test_order2(self): dom=RectangularDomain() dom.order=2 run(dom, sys.stdout)
import os from esys.modellib.mechanics import DruckerPrager from esys.escript.modelframe import Link, Simulation from esys.modellib.geometry import RectangularDomain, VectorConstrainerOverBox, UpdateGeometry from esys.modellib.input import Sequencer, InterpolateOverBox from esys.modellib.visualization import WriteVTK try: WORKDIR = os.environ['MODELLIB_WORKDIR'] except KeyError: WORKDIR = '.' dbg = True dom = RectangularDomain(debug=dbg) dom.dim = 3 dom.l = [0.5, 1., 1.] dom.n = [30, 6, 6] dom.order = 1 dom.integrationOrder = -1 sq = Sequencer(debug=dbg) sq.t = 0 sq.t_end = 0.8 sq.dt_max = 100. iob = InterpolateOverBox(debug=dbg) iob.domain = Link(dom, "domain") iob.value_left_bottom_front = [-1., 0., 0.] iob.value_right_bottom_front = [0., 0., 0.]
def test_order1(self): dom=RectangularDomain(esys.dudley) dom.order=1 run(dom, sys.stdout)
import os from esys.modellib.mechanics import DruckerPrager from esys.escript.modelframe import Link,Simulation from esys.modellib.geometry import RectangularDomain, VectorConstrainerOverBox, UpdateGeometry from esys.modellib.input import Sequencer, InterpolateOverBox from esys.modellib.visualization import WriteVTK try: WORKDIR=os.environ['MODELLIB_WORKDIR'] except KeyError: WORKDIR='.' dbg=True dom=RectangularDomain(debug=dbg) dom.dim=3 dom.l=[0.5,1.,1.] dom.n=[30,6,6] dom.order=1 dom.integrationOrder=-1 sq=Sequencer(debug=dbg) sq.t=0 sq.t_end=0.8 sq.dt_max=100. iob=InterpolateOverBox(debug=dbg) iob.domain=Link(dom,"domain") iob.value_left_bottom_front=[-1.,0.,0.]