Re: [Libmesh-users] Sudden explosion of maximum error bound

[David K. <dav...@ak...>]

On Mon, Mar 26, 2018 at 1:07 AM, Gauvain Wu <cau...@gm...> wrote:

> 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.
>


I don't know what the problem is, but note that you have 9 parameters here
which is a lot of parameters for an RB model. It may not be practical to
use this many parameters in your case.

I would suggest starting with 1 parameter (set the other parameters to
fixed values) and make sure that works, then go to 2 parameters, etc. That
may help you to figure out what the issue is as well.

David




>
> via Newton Mail [https://cloudmagic.com/k/d/ma
> ilapp?ct=dw&cv=9.8.284&pv=10.0&source=email_footer_2]
> ------------------------------------------------------------
> ------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites, Slashdot.org! http://sdm.link/slashdot
> _______________________________________________
> Libmesh-users mailing list
> Lib...@li...
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>