synthetic_solve.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	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)