square.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	import numpy as np
	from rrompy.hfengines.linear_problem.bidimensional import \
	HelmholtzSquareDomainProblemEngine as HSDPE
	verb = 100

	mu0 = [4 .5, 1.5 .5]
	solver = HSDPE(kappa = 2.5, theta = np.pi / 3, mu0 = mu0, n = 50,
	verbosity = verb)
	uh = solver.solve(mu0)[0]
	solver.plotmesh()
	print(solver.norm(uh))
	solver.plot(uh)
	solver.plot(solver.residual(mu0, uh)[0], 'res')