[Libmesh-users] Sudden explosion of maximum error bound

[Gauvain Wu <cau...@gm...>]

Hi all,
I was running the offline stage of a reduced basis model. At the beginning, 
the maximum error bound decreased at a normal rate and the value appeared 
acceptable (around 20000). However, when the basis dimension increased to 
6, the maximum error bound climbed drastically to e+149 or so and didn't 
decrease any more. What I had modified was the range of the parameters 
which is listed below:
beta_T0 = '1.e3 3.e3' # min, max values
rho_c = '2.4 4.5'
heat_flux = '20 30'
h = '0 1'
h_(Tinf-T0)='0 1'
G = '7.e4 9.e4'
lambda = '2.e5 4.e5'
beta = '5 15'
k = '10 20'
The code was run with option mpirun -np 4 ./Thermoelasticity-opt 
-online_mode 0 -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package 
mumps. Could anyone give a hint of what causes this problem? Thanks in 
advance.
Gauvain

via Newton Mail 
[https://cloudmagic.com/k/d/mailapp?ct=dw&cv=9.8.284&pv=10.0&source=email_footer_2]