def fixme_test_mesh_read_rectangle_from_dudley_file(self):
mydomain1 = Rectangle(n0=8,
n1=10,
order=1,
l0=1.,
l1=1.,
optimize=False)
mydomain2 = ReadMesh(
os.path.join(DUDLEY_TEST_MESH_PATH, "rectangle_8x10.fly"))
self.domainsEqual(mydomain1, mydomain2)
def test_flyTags(self):
dom = ReadMesh(os.path.join(DUDLEY_TEST_MESH_PATH, "tagtest2.fly"))
tags = sorted(dom.showTagNames().split(', '))
self.assertEqual(tags, sorted(['tag1', 'tag2', 'tag3', 'tag4', 'All']))
self.assertEqual(dom.getTag('tag1'), 5, 'error with tag1')
self.assertEqual(dom.getTag('tag2'), 8, 'error with tag2,')
self.assertEqual(dom.getTag('tag3'), 6, 'error with tag3')
self.assertEqual(dom.getTag('tag4'), 7, 'error with tag4')
self.assertEqual(dom.getTag('All'), 10, 'error with All')
self.assertRaises(ValueError, dom.getTag, 'tag6')
def fixme_test_mesh_read_brick_from_dudley_file(self):
mydomain1 = Brick(n0=8,
n1=10,
n2=12,
order=1,
l0=1.,
l1=1.,
l2=1.,
optimize=False)
mydomain2 = ReadMesh(
os.path.join(DUDLEY_TEST_MESH_PATH, "brick_8x10x12.fly"))
self.domainsEqual(mydomain1, mydomain2)
# Parameter # DIM = 2 MESHFILE = "sub.fly" ETA = 1.e22 * U.Pa * U.sec V_MAX = 1. * U.cm / U.yr ALPHA = 30 * U.DEG STRIKE = 10 * U.DEG DIP = 30 * U.DEG N = 1 # boudary layer control g = 9.81 * U.m / U.sec**2 # # derived values # dom = ReadMesh(MESHFILE) DIM = dom.getDim() bb = boundingBox(dom) LX = bb[0][1] - bb[0][0] if DIM == 3: LY = bb[1][1] - bb[1][0] DEPTH = bb[DIM - 1][1] - bb[DIM - 1][0] sc = StokesProblemCartesian(dom) x = dom.getX() # v = Vector(0., Solution(dom)) mask = Vector(0., Solution(dom)) # # in subduction zone: #
# print "p:",inf(p),sup(p)
# print "u:",inf(u),sup(u)
self.__pde_u.setValue(y_contact=-p)
return self.__pde_u.getSolution()
def solve_g(self, u, tol=1.e-8):
dp = Vector(0., FunctionOnContactZero(self.domain))
h = FunctionOnContactZero(self.domain).getSize()
# print jump(u)-self.slip
dp[0] = (self.slip[0] - jump(u[0])) * lam_mu / h
dp[1] = (self.slip[1] - jump(u[1])) * lam_mu / h
dp[2] = (self.slip[2] - jump(u[2])) * lam_mu / h
return dp
dom = ReadMesh("meshfault3D.fly", integrationOrder=-1)
prop = SlippingFault(dom)
d = dom.getDim()
x = dom.getX()[0]
# x=dom.getX()[d-1]
mask = whereZero(x - inf(x)) * numpy.ones((d, ))
x = FunctionOnContactZero(dom).getX()
s = numpy.array([-100000., 1., 1.])
for i in range(3):
d = fend[i] - fstart[i]
if d > 0:
q = (x[i] - fstart[i]) / d
s = q * (1 - q) * 4 * s
elif d < 0:
q = (x[i] - fend[i]) / d
s = q * (1 - q) * 4 * s
def setUp(self):
self.domain = ReadMesh(os.path.join(DUDLEY_TEST_MESH_PATH,"tet_2D_order1.fly"),optimize=False)
self.order = 1
