Page MenuHomec4science
Files F85087596

funnel_output.py
No OneTemporary
Actions

File Metadata

Created
Thu, Sep 26, 16:59

funnel_output.py
View Options

import numpy as np
from funnel_output_engine import FunnelOutputEngine as engine
from rrompy.reduction_methods import (NearestNeighbor as NN,
RationalInterpolant as RI,
RationalInterpolantGreedy as RIG)
from rrompy.parameter.parameter_sampling import (QuadratureSampler as QS,
EmptySampler as ES)
ks = [5., 10.]
k0, n = np.mean(ks), 50
solver = engine(k0, n)
k = 6.5
for method in ["RI", "RI_GREEDY"]:
print("Testing {} method".format(method))
if "GREEDY" not in method:
params = {'S':20, 'POD':True, 'polybasis':"CHEBYSHEV",
'sampler':QS(ks, "CHEBYSHEV")}
algo = RI
if "GREEDY" in method:
params = {'S':2, 'POD':True, 'polybasis':"LEGENDRE", 'greedyTol':1e-1,
'maxIter':25, 'sampler':QS(ks, "UNIFORM"),
'errorEstimatorKind':"LOOK_AHEAD_OUTPUT"}
algo = RIG
approx = algo(solver, mu0 = k0, approxParameters = params, verbosity = 5)
if "GREEDY" not in method:
approx.setupApprox()
else:
approx.setupApprox("LAST")
print("--- Approximant ---")
approx.plotApprox(k, name = 'u_app')
approx.plotHF(k, name = 'u_HF')
approx.plotErr(k, name = 'err_app')
err = approx.getErr(k)[0]
sol = approx.getHF(k)[0]
normErr = np.abs(solver.L2NormMatrix.dot(err).dot(err.conj())) ** .5
normSol = np.abs(solver.L2NormMatrix.dot(sol).dot(sol.conj())) ** .5
print("SolNorm:\t{:.5e}\nErr_app: \t{:.5e}\nErrRel_app:\t{:.5e}".format(
normSol, normErr, normErr / normSol))
print("--- Closest snapshot ---")
approxNN = NN(solver, mu0 = k0, verbosity = 0,
approxParameters = {'S':approx.samplingEngine.nsamples,
'POD':True, 'sampler':ES()})
approxNN.setSamples(approx.samplingEngine)
approxNN.plotApprox(k, name = 'u_close')
approxNN.plotHF(k, name = 'u_HF')
approxNN.plotErr(k, name = 'err_close')
err = approxNN.getErr(k)[0]
sol = approxNN.getHF(k)[0]
normErr = np.abs(solver.L2NormMatrix.dot(err).dot(err.conj())) ** .5
normSol = np.abs(solver.L2NormMatrix.dot(sol).dot(sol.conj())) ** .5
print("SolNorm:\t{:.5e}\nErr_close:\t{:.5e}\nErrRel_close:\t{:.5e}".format(
normSol, normErr, normErr / normSol))
print("Poles:\n{}".format(approx.getPoles()))
print("\n")

Event Timeline

c4science · Help