fracture3.py
3ed8cef19cfcmaster

fracture3.py
View Options

			from rrompy.hfengines.linear_problem.tridimensional import \
			MembraneFractureEngine3 as MFE
			from fracture3_warping import fracture3_warping

			verb = 100

			mu0 = [45. ** .5, .7, .075]

			H = 1.
			L = .75
			delta = .05
			n = 50

			solver = MFE(mu0 = mu0, H = H, L = L, delta = delta,
			n = n, verbosity = verb)
			u0 = solver.liftDirichletData()
			uh = solver.solve(mu0)[0]
			#solver.plotmesh(figsize = (7.5, 4.5))
			#solver.plot(u0, what = 'REAL', figsize = (8, 5))
			print(solver.norm(uh))
			#solver.plot(uh, what = 'REAL', figsize = (8, 5))
			#solver.plot(solver.residual(mu0, uh)[0], name = 'res',
			# what = 'REAL', figsize = (8, 5))
			#solver.outParaviewTimeDomain(uh, mu0[0], filename = 'out', folder = True,
			# forceNewFile = False)

			warps = fracture3_warping(solver.V.mesh(), L, mu0[1], delta, mu0[2])

			#solver.plotmesh(warps, figsize = (7.5, 4.5))
			#solver.plot(u0, warps, what = 'REAL', figsize = (8, 5))
			solver.plot(uh, warps, what = 'REAL', figsize = (8, 5))
			#solver.plot(solver.residual(mu0, uh)[0], warps, name = 'res',
			# what = 'REAL', figsize = (8, 5))
			#solver.outParaviewTimeDomain(uh, mu0[0], warps, filename = 'outW',
			# folder = True, forceNewFile = False)