synthetic_solve.py
171a647d0193master

synthetic_solve.py
View Options

			import numpy as np
			import fenics as fen
			from rrompy.hfengines.linear_problem.bidimensional import \
			SyntheticBivariateEngine as SBE

			verb = 100

			kappa = 15. ** .5
			theta = - np.pi / 6.
			n = 30
			L = np.pi

			mu0 = [15. .5, 20. .5]
			mutar = [15. .5, 20. .5]

			solver = SBE(kappa = kappa, theta = theta, n = n, L = L,
			mu0 = mu0, verbosity = verb)
			uh = solver.solve(mutar)[0]
			solver.plotmesh(figsize = (7.5, 4.5))
			fen.plot(fen.project(solver._above,
			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)