Exemplos de RectangularDomain em Python

Exemplo n.º 1

Arquivo: run_convection.py Projeto: markendr/esys-escript.github.io

 def test_order1(self):
     dom = RectangularDomain(esys.dudley)
     dom.order = 1
     run(dom, sys.stdout)

Exemplo n.º 2

Arquivo: run_convection.py Projeto: markendr/esys-escript.github.io

 def test_order2(self):
     dom = RectangularDomain()
     dom.order = 2
     run(dom, sys.stdout)

Exemplo n.º 3

Arquivo: run_temp.py Projeto: svn2github/Escript

 def test_order2(self):
     dom=RectangularDomain()
     dom.order=2
     run(dom, sys.stdout)

Exemplo n.º 4

Arquivo: drucker_prager.py Projeto: markendr/esys-escript.github.io

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.]

Exemplo n.º 5

Arquivo: run_convection.py Projeto: svn2github/Escript

 def test_order1(self):
     dom=RectangularDomain(esys.dudley)
     dom.order=1
     run(dom, sys.stdout)

Exemplo n.º 6

Arquivo: drucker_prager.py Projeto: svn2github/Escript

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.]