cookie_single.py
3edebe56625bmaster

cookie_single.py
View Options

			import numpy as np
			import fenics as fen
			from rrompy.hfengines.linear_problem.bidimensional import \
			CookieEngineSingle as CES

			verb = 100

			kappa = 15. ** .5
			theta = - np.pi / 6.
			n = 30
			R = 1.
			L = np.pi
			nX = 3
			nY = 2

			mu0 = [25. ** .5, 1.]
			mutar = [25. ** .5, 1.5]

			solver = CES(kappa = kappa, theta = theta, n = n, R = R, L = L, nX = nX,
			nY = nY, mu0 = mu0, verbosity = verb)
			uh = solver.solve(mutar)[0]
			solver.plotmesh(figsize = (7.5, 4.5))
			fen.plot(fen.project(solver.cookieIn,
			fen.FunctionSpace(solver.V.mesh(), "DG", 0)))
			print(solver.norm(uh))
			solver.plot(uh)
			solver.outParaviewTimeDomain(uh, mutar[0], filename = 'out', folder = True,
			forceNewFile = False)