parametricDomain.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	import numpy as np
	from rrompy.hfengines.linear_problem import \
	HelmholtzSquareBubbleDomainProblemEngine as HSBDPE
	from rrompy.reduction_methods.centered import RationalPade as Pade
	from rrompy.reduction_methods.centered import RBCentered as RB
	from operator import itemgetter
	def subdict(d, ks):
	return dict(zip(ks, itemgetter(*ks)(d)))

	testNo = 2
	verb = 0

	if testNo == 1:
	mu = 7 ** .5
	solver = HSBDPE(kappa = 12 ** .5, theta = np.pi / 3, n = 40, mu0 = mu,
	degree_threshold = 15, verbosity = verb)

	uh = solver.solve(mu)
	solver.plotmesh()
	print(solver.norm(uh))
	solver.plot(uh)

	############
	if testNo == 2:
	params = {'N':8, 'M':8, 'E':8, 'POD':True}
	mu0 = 7 ** .5
	mutar = (7. + .1j) ** .5
	solver = HSBDPE(kappa = 12 ** .5, theta = np.pi / 3, n = 40, mu0 = mu0,
	degree_threshold = 15, verbosity = verb)
	approxP = Pade(solver, mu0 = mu0, approxParameters = params,
	verbosity = verb)
	paramsRB = subdict(params, ['E', 'POD'])
	approxR = RB(solver, mu0 = mu0, approxParameters = paramsRB,
	verbosity = verb)

	approxP.setupApprox()
	approxR.setupApprox()
	# approxP.plotSamples()
	approxP.plotHF(mutar, name = 'u_HF')
	approxP.plotApprox(mutar, name = 'u_Pade''')
	approxR.plotApprox(mutar, name = 'u_RB')
	approxP.plotErr(mutar, name = 'err_Pade''')
	approxR.plotErr(mutar, name = 'err_RB')

	solNorm = approxP.normHF(mutar)
	appPErr = approxP.normErr(mutar)
	appRErr = approxR.normErr(mutar)
	print(('SolNorm:\t{}\nErrRelP:\t{}\nErrRelR:\t{}').format(solNorm,
	appPErr / solNorm, appRErr / solNorm))
	print('\nPoles Pade'':')
	print(approxP.getPoles())