def test_order2(self): dom = RectangularDomain() dom.order = 2 run(dom, sys.stdout)
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.] iob.value_right_bottom_front = [0., 0., 0.]