search

LutzGross / esys-escript.github.io

Other
29 stars 13 forks source link

Oxley: grid is not refined on both sides of a line. #93

Open LutzGross opened 1 year ago

LutzGross commented 1 year ago

from esys.oxley import Rectangle

domain = Rectangle(n0=10, n1=10, l0=(-5, 5), l1=(-5, 5))

domain.setRefinementLevel(3) boundaryLayerWidth=0
domain.refineRegion(x0=-5, x1=5, y0=-boundaryLayerWidth, y1=boundaryLayerWidth)

the refinement is added symmetrically above and blow line y=0.