Hi K,

These errors aren't particularly informative. Can you do a smaller 
version of the problem so that you only need "-np 1", and run it in 
debug mode, and let me know what happens then?

David


On 12/03/2012 02:21 AM, Kyunghoon Lee wrote:
> Hello all,
>
> I made a simple cantilever beam model parametrized w.r.t tip loads in x & y
> directions with RB methods for steady solutions, then extended the model
> for transient solutions.  For the steady model, I run
>
> $PETSC_DIR/bin/mpirun -np 8 ./beam_2d-opt  -ksp_type preonly -pc_type lu
> -pc_factor_mat_solver_package mumps
>
> and the training was terminated after finding two basis functions.  Now for
> the transient model, if I use the same run options, it crashes while
> computing output dual norms as follows (the outputs are averaged tip
> displacements in x and y):
>
> -----------------------------------------------------------------------------------------------------------------
> $PETSC_DIR/bin/mpirun -np 8 ./beam_2d-opt  -ksp_type preonly -pc_type lu
> -pc_factor_mat_solver_package mumps
> *** Warning, This code is untested, experimental, or likely to see future
> API changes: src/reduced_basis/rb_parametrized.C, line 41, compiled Nov 27
> 2012 at 17:00:18 ***
> Initializing training parameters with random training set...
> Parameter load_Fx: log scaling = 0
> Parameter load_Fy: log scaling = 0
>
>
> RBConstruction parameters:
> system name: RBTransientElasticity
> constrained_problem: 0
> Nmax: 10
> Basis training error tolerance: 1e-06
> Aq operators attached: 1
> Fq functions attached: 2
> n_outputs: 2
> output 0, Q_l = 1
> output 1, Q_l = 1
> Number of parameters: 2
> Parameter load_Fx: Min = -5, Max = 5, value = 0
> Parameter load_Fy: Min = -5, Max = 5, value = 5
> n_training_samples: 100
> single-matrix mode? 0
> reuse preconditioner? 1
> use a relative error bound in greedy? 0
> write out data during basis training? 0
> quiet mode? 1
>
>
> TransientRBConstruction parameters:
> Q_m: 1
> Number of time-steps: 100
> dt: 0.01
> euler_theta (time discretization parameter): 1
> delta_N (number of basis functions to add each POD-Greedy step): 1
> Using zero initial condition
>
> *** Warning, This code is deprecated, and likely to be removed in future
> library versions! src/mesh/boundary_info.C, line 750, compiled Nov 27 2012
> at 16:59:26 ***
> load_Fx: 0.000000e+00
> load_Fy: 1.000000e+00
>
> t_step[0] = 0.0e+00    truth_out[0] = 0.0000e+00 truth_out[1] = 0.0000e+00
> t_step[1] = 1.0e-02    truth_out[0] = 6.4695e-06 truth_out[1] = 2.9555e+01
> t_step[2] = 2.0e-02    truth_out[0] = -1.6811e-05 truth_out[1] = 3.1663e+01
> t_step[3] = 3.0e-02    truth_out[0] = -2.3264e-05 truth_out[1] = 3.3369e+01
> t_step[4] = 4.0e-02    truth_out[0] = -2.3126e-05 truth_out[1] = 3.4889e+01
> t_step[5] = 5.0e-02    truth_out[0] = -2.0980e-05 truth_out[1] = 3.6299e+01
> t_step[6] = 6.0e-02    truth_out[0] = -1.8536e-05 truth_out[1] = 3.7639e+01
> t_step[7] = 7.0e-02    truth_out[0] = -1.6399e-05 truth_out[1] = 3.8929e+01
> t_step[8] = 8.0e-02    truth_out[0] = -1.4733e-05 truth_out[1] = 4.0183e+01
> t_step[9] = 9.0e-02    truth_out[0] = -1.3539e-05 truth_out[1] = 4.1411e+01
> t_step[10] = 1.0e-01    truth_out[0] = -1.2752e-05 truth_out[1] =
> 4.2618e+01
> t_step[11] = 1.1e-01    truth_out[0] = -1.2295e-05 truth_out[1] =
> 4.3811e+01
> t_step[12] = 1.2e-01    truth_out[0] = -1.2092e-05 truth_out[1] =
> 4.4991e+01
> t_step[13] = 1.3e-01    truth_out[0] = -1.2079e-05 truth_out[1] =
> 4.6161e+01
> t_step[14] = 1.4e-01    truth_out[0] = -1.2201e-05 truth_out[1] =
> 4.7324e+01
> t_step[15] = 1.5e-01    truth_out[0] = -1.2414e-05 truth_out[1] =
> 4.8481e+01
> t_step[16] = 1.6e-01    truth_out[0] = -1.2686e-05 truth_out[1] =
> 4.9632e+01
> t_step[17] = 1.7e-01    truth_out[0] = -1.2990e-05 truth_out[1] =
> 5.0780e+01
> t_step[18] = 1.8e-01    truth_out[0] = -1.3308e-05 truth_out[1] =
> 5.1923e+01
> t_step[19] = 1.9e-01    truth_out[0] = -1.3625e-05 truth_out[1] =
> 5.3064e+01
> t_step[20] = 2.0e-01    truth_out[0] = -1.3933e-05 truth_out[1] =
> 5.4202e+01
> t_step[21] = 2.1e-01    truth_out[0] = -1.4224e-05 truth_out[1] =
> 5.5338e+01
> t_step[22] = 2.2e-01    truth_out[0] = -1.4493e-05 truth_out[1] =
> 5.6471e+01
> t_step[23] = 2.3e-01    truth_out[0] = -1.4740e-05 truth_out[1] =
> 5.7602e+01
> t_step[24] = 2.4e-01    truth_out[0] = -1.4962e-05 truth_out[1] =
> 5.8732e+01
> t_step[25] = 2.5e-01    truth_out[0] = -1.5159e-05 truth_out[1] =
> 5.9860e+01
> t_step[26] = 2.6e-01    truth_out[0] = -1.5331e-05 truth_out[1] =
> 6.0986e+01
> t_step[27] = 2.7e-01    truth_out[0] = -1.5480e-05 truth_out[1] =
> 6.2111e+01
> t_step[28] = 2.8e-01    truth_out[0] = -1.5607e-05 truth_out[1] =
> 6.3234e+01
> t_step[29] = 2.9e-01    truth_out[0] = -1.5713e-05 truth_out[1] =
> 6.4356e+01
> t_step[30] = 3.0e-01    truth_out[0] = -1.5800e-05 truth_out[1] =
> 6.5477e+01
> t_step[31] = 3.1e-01    truth_out[0] = -1.5869e-05 truth_out[1] =
> 6.6596e+01
> t_step[32] = 3.2e-01    truth_out[0] = -1.5922e-05 truth_out[1] =
> 6.7713e+01
> t_step[33] = 3.3e-01    truth_out[0] = -1.5960e-05 truth_out[1] =
> 6.8830e+01
> t_step[34] = 3.4e-01    truth_out[0] = -1.5984e-05 truth_out[1] =
> 6.9945e+01
> t_step[35] = 3.5e-01    truth_out[0] = -1.5997e-05 truth_out[1] =
> 7.1059e+01
> t_step[36] = 3.6e-01    truth_out[0] = -1.5999e-05 truth_out[1] =
> 7.2172e+01
> t_step[37] = 3.7e-01    truth_out[0] = -1.5991e-05 truth_out[1] =
> 7.3284e+01
> t_step[38] = 3.8e-01    truth_out[0] = -1.5975e-05 truth_out[1] =
> 7.4394e+01
> t_step[39] = 3.9e-01    truth_out[0] = -1.5951e-05 truth_out[1] =
> 7.5504e+01
> t_step[40] = 4.0e-01    truth_out[0] = -1.5921e-05 truth_out[1] =
> 7.6612e+01
> t_step[41] = 4.1e-01    truth_out[0] = -1.5885e-05 truth_out[1] =
> 7.7719e+01
> t_step[42] = 4.2e-01    truth_out[0] = -1.5843e-05 truth_out[1] =
> 7.8825e+01
> t_step[43] = 4.3e-01    truth_out[0] = -1.5797e-05 truth_out[1] =
> 7.9930e+01
> t_step[44] = 4.4e-01    truth_out[0] = -1.5747e-05 truth_out[1] =
> 8.1034e+01
> t_step[45] = 4.5e-01    truth_out[0] = -1.5694e-05 truth_out[1] =
> 8.2137e+01
> t_step[46] = 4.6e-01    truth_out[0] = -1.5638e-05 truth_out[1] =
> 8.3239e+01
> t_step[47] = 4.7e-01    truth_out[0] = -1.5579e-05 truth_out[1] =
> 8.4339e+01
> t_step[48] = 4.8e-01    truth_out[0] = -1.5518e-05 truth_out[1] =
> 8.5439e+01
> t_step[49] = 4.9e-01    truth_out[0] = -1.5456e-05 truth_out[1] =
> 8.6538e+01
> t_step[50] = 5.0e-01    truth_out[0] = -1.5392e-05 truth_out[1] =
> 8.7635e+01
> t_step[51] = 5.1e-01    truth_out[0] = -1.5327e-05 truth_out[1] =
> 8.8732e+01
> t_step[52] = 5.2e-01    truth_out[0] = -1.5261e-05 truth_out[1] =
> 8.9828e+01
> t_step[53] = 5.3e-01    truth_out[0] = -1.5194e-05 truth_out[1] =
> 9.0922e+01
> t_step[54] = 5.4e-01    truth_out[0] = -1.5127e-05 truth_out[1] =
> 9.2016e+01
> t_step[55] = 5.5e-01    truth_out[0] = -1.5059e-05 truth_out[1] =
> 9.3109e+01
> t_step[56] = 5.6e-01    truth_out[0] = -1.4991e-05 truth_out[1] =
> 9.4200e+01
> t_step[57] = 5.7e-01    truth_out[0] = -1.4923e-05 truth_out[1] =
> 9.5291e+01
> t_step[58] = 5.8e-01    truth_out[0] = -1.4855e-05 truth_out[1] =
> 9.6381e+01
> t_step[59] = 5.9e-01    truth_out[0] = -1.4788e-05 truth_out[1] =
> 9.7470e+01
> t_step[60] = 6.0e-01    truth_out[0] = -1.4720e-05 truth_out[1] =
> 9.8558e+01
> t_step[61] = 6.1e-01    truth_out[0] = -1.4653e-05 truth_out[1] =
> 9.9646e+01
> t_step[62] = 6.2e-01    truth_out[0] = -1.4586e-05 truth_out[1] =
> 1.0073e+02
> t_step[63] = 6.3e-01    truth_out[0] = -1.4519e-05 truth_out[1] =
> 1.0182e+02
> t_step[64] = 6.4e-01    truth_out[0] = -1.4453e-05 truth_out[1] =
> 1.0290e+02
> t_step[65] = 6.5e-01    truth_out[0] = -1.4387e-05 truth_out[1] =
> 1.0399e+02
> t_step[66] = 6.6e-01    truth_out[0] = -1.4322e-05 truth_out[1] =
> 1.0507e+02
> t_step[67] = 6.7e-01    truth_out[0] = -1.4258e-05 truth_out[1] =
> 1.0615e+02
> t_step[68] = 6.8e-01    truth_out[0] = -1.4194e-05 truth_out[1] =
> 1.0723e+02
> t_step[69] = 6.9e-01    truth_out[0] = -1.4130e-05 truth_out[1] =
> 1.0831e+02
> t_step[70] = 7.0e-01    truth_out[0] = -1.4068e-05 truth_out[1] =
> 1.0939e+02
> t_step[71] = 7.1e-01    truth_out[0] = -1.4006e-05 truth_out[1] =
> 1.1047e+02
> t_step[72] = 7.2e-01    truth_out[0] = -1.3944e-05 truth_out[1] =
> 1.1155e+02
> t_step[73] = 7.3e-01    truth_out[0] = -1.3883e-05 truth_out[1] =
> 1.1262e+02
> t_step[74] = 7.4e-01    truth_out[0] = -1.3823e-05 truth_out[1] =
> 1.1370e+02
> t_step[75] = 7.5e-01    truth_out[0] = -1.3764e-05 truth_out[1] =
> 1.1477e+02
> t_step[76] = 7.6e-01    truth_out[0] = -1.3705e-05 truth_out[1] =
> 1.1585e+02
> t_step[77] = 7.7e-01    truth_out[0] = -1.3647e-05 truth_out[1] =
> 1.1692e+02
> t_step[78] = 7.8e-01    truth_out[0] = -1.3589e-05 truth_out[1] =
> 1.1799e+02
> t_step[79] = 7.9e-01    truth_out[0] = -1.3532e-05 truth_out[1] =
> 1.1907e+02
> t_step[80] = 8.0e-01    truth_out[0] = -1.3476e-05 truth_out[1] =
> 1.2014e+02
> t_step[81] = 8.1e-01    truth_out[0] = -1.3421e-05 truth_out[1] =
> 1.2121e+02
> t_step[82] = 8.2e-01    truth_out[0] = -1.3366e-05 truth_out[1] =
> 1.2228e+02
> t_step[83] = 8.3e-01    truth_out[0] = -1.3312e-05 truth_out[1] =
> 1.2335e+02
> t_step[84] = 8.4e-01    truth_out[0] = -1.3258e-05 truth_out[1] =
> 1.2441e+02
> t_step[85] = 8.5e-01    truth_out[0] = -1.3205e-05 truth_out[1] =
> 1.2548e+02
> t_step[86] = 8.6e-01    truth_out[0] = -1.3153e-05 truth_out[1] =
> 1.2655e+02
> t_step[87] = 8.7e-01    truth_out[0] = -1.3101e-05 truth_out[1] =
> 1.2761e+02
> t_step[88] = 8.8e-01    truth_out[0] = -1.3050e-05 truth_out[1] =
> 1.2868e+02
> t_step[89] = 8.9e-01    truth_out[0] = -1.2999e-05 truth_out[1] =
> 1.2974e+02
> t_step[90] = 9.0e-01    truth_out[0] = -1.2949e-05 truth_out[1] =
> 1.3081e+02
> t_step[91] = 9.1e-01    truth_out[0] = -1.2900e-05 truth_out[1] =
> 1.3187e+02
> t_step[92] = 9.2e-01    truth_out[0] = -1.2851e-05 truth_out[1] =
> 1.3293e+02
> t_step[93] = 9.3e-01    truth_out[0] = -1.2803e-05 truth_out[1] =
> 1.3399e+02
> t_step[94] = 9.4e-01    truth_out[0] = -1.2755e-05 truth_out[1] =
> 1.3505e+02
> t_step[95] = 9.5e-01    truth_out[0] = -1.2708e-05 truth_out[1] =
> 1.3611e+02
> t_step[96] = 9.6e-01    truth_out[0] = -1.2661e-05 truth_out[1] =
> 1.3717e+02
> t_step[97] = 9.7e-01    truth_out[0] = -1.2615e-05 truth_out[1] =
> 1.3823e+02
> t_step[98] = 9.8e-01    truth_out[0] = -1.2569e-05 truth_out[1] =
> 1.3929e+02
> t_step[99] = 9.9e-01    truth_out[0] = -1.2524e-05 truth_out[1] =
> 1.4034e+02
> Compute output dual inner products
> [0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [1]PETSC ERROR: Error in external library!
> [1]PETSC ERROR: Error reported by MUMPS in numerical factorization phase:
> INFO(1)=-1, INFO(2)=4
> !
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51
> CDT 2012
> [1]PETSC ERROR: See docs/changes/index.html for recent updates.
> [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [1]PETSC ERROR: See docs/index.html for manual pages.
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: ./beam_2d-opt on a arch-linu named KLee-Linux by aeronova
> Mon Dec  3 14:18:44 2012
> .....
> -----------------------------------------------------------------------------------------------------------------
>
> If I change the run option as
>
> $PETSC_DIR/bin/mpirun -np 8 ./beam_2d-opt -ksp_type preonly
>
> then it looks working, but the output values seem wrong, and somehow the
> output dual norm for the second output is zero:
>
> -----------------------------------------------------------------------------------------------------------------
> *** Warning, This code is untested, experimental, or likely to see future
> API changes: src/reduced_basis/rb_parametrized.C, line 41, compiled Nov 27
> 2012 at 17:00:18 ***
> Initializing training parameters with random training set...
> Parameter load_Fx: log scaling = 0
> Parameter load_Fy: log scaling = 0
>
>
> RBConstruction parameters:
> system name: RBTransientElasticity
> constrained_problem: 0
> Nmax: 10
> Basis training error tolerance: 1e-06
> Aq operators attached: 1
> Fq functions attached: 2
> n_outputs: 2
> output 0, Q_l = 1
> output 1, Q_l = 1
> Number of parameters: 2
> Parameter load_Fx: Min = -5, Max = 5, value = 0
> Parameter load_Fy: Min = -5, Max = 5, value = 5
> n_training_samples: 100
> single-matrix mode? 0
> reuse preconditioner? 1
> use a relative error bound in greedy? 0
> write out data during basis training? 0
> quiet mode? 1
>
>
> TransientRBConstruction parameters:
> Q_m: 1
> Number of time-steps: 100
> dt: 0.01
> euler_theta (time discretization parameter): 1
> delta_N (number of basis functions to add each POD-Greedy step): 1
> Using zero initial condition
>
> *** Warning, This code is deprecated, and likely to be removed in future
> library versions! src/mesh/boundary_info.C, line 750, compiled Nov 27 2012
> at 16:59:26 ***
> load_Fx: 0.000000e+00
> load_Fy: 1.000000e+00
>
> t_step[0] = 0.0e+00    truth_out[0] = 0.0000e+00 truth_out[1] = 0.0000e+00
> t_step[1] = 1.0e-02    truth_out[0] = -1.6962e-04 truth_out[1] = 1.5122e-01
> t_step[2] = 2.0e-02    truth_out[0] = -2.3491e-04 truth_out[1] = 1.5378e-01
> t_step[3] = 3.0e-02    truth_out[0] = -2.9471e-04 truth_out[1] = 1.5459e-01
> t_step[4] = 4.0e-02    truth_out[0] = -3.3419e-04 truth_out[1] = 1.5489e-01
> t_step[5] = 5.0e-02    truth_out[0] = -3.5753e-04 truth_out[1] = 1.5500e-01
> t_step[6] = 6.0e-02    truth_out[0] = -3.7049e-04 truth_out[1] = 1.5505e-01
> t_step[7] = 7.0e-02    truth_out[0] = -3.7738e-04 truth_out[1] = 1.5507e-01
> t_step[8] = 8.0e-02    truth_out[0] = -3.8093e-04 truth_out[1] = 1.5508e-01
> t_step[9] = 9.0e-02    truth_out[0] = -3.8271e-04 truth_out[1] = 1.5508e-01
> t_step[10] = 1.0e-01    truth_out[0] = -3.8359e-04 truth_out[1] =
> 1.5508e-01
> t_step[11] = 1.1e-01    truth_out[0] = -3.8401e-04 truth_out[1] =
> 1.5508e-01
> t_step[12] = 1.2e-01    truth_out[0] = -3.8422e-04 truth_out[1] =
> 1.5508e-01
> t_step[13] = 1.3e-01    truth_out[0] = -3.8431e-04 truth_out[1] =
> 1.5508e-01
> t_step[14] = 1.4e-01    truth_out[0] = -3.8436e-04 truth_out[1] =
> 1.5508e-01
> t_step[15] = 1.5e-01    truth_out[0] = -3.8438e-04 truth_out[1] =
> 1.5508e-01
> t_step[16] = 1.6e-01    truth_out[0] = -3.8439e-04 truth_out[1] =
> 1.5508e-01
> t_step[17] = 1.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[18] = 1.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[19] = 1.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[20] = 2.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[21] = 2.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[22] = 2.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[23] = 2.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[24] = 2.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[25] = 2.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[26] = 2.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[27] = 2.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[28] = 2.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[29] = 2.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[30] = 3.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[31] = 3.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[32] = 3.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[33] = 3.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[34] = 3.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[35] = 3.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[36] = 3.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[37] = 3.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[38] = 3.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[39] = 3.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[40] = 4.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[41] = 4.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[42] = 4.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[43] = 4.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[44] = 4.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[45] = 4.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[46] = 4.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[47] = 4.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[48] = 4.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[49] = 4.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[50] = 5.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[51] = 5.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[52] = 5.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[53] = 5.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[54] = 5.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[55] = 5.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[56] = 5.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[57] = 5.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[58] = 5.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[59] = 5.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[60] = 6.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[61] = 6.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[62] = 6.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[63] = 6.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[64] = 6.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[65] = 6.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[66] = 6.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[67] = 6.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[68] = 6.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[69] = 6.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[70] = 7.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[71] = 7.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[72] = 7.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[73] = 7.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[74] = 7.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[75] = 7.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[76] = 7.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[77] = 7.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[78] = 7.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[79] = 7.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[80] = 8.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[81] = 8.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[82] = 8.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[83] = 8.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[84] = 8.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[85] = 8.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[86] = 8.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[87] = 8.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[88] = 8.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[89] = 8.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[90] = 9.0e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[91] = 9.1e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[92] = 9.2e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[93] = 9.3e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[94] = 9.4e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[95] = 9.5e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[96] = 9.6e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[97] = 9.7e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[98] = 9.8e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> t_step[99] = 9.9e-01    truth_out[0] = -3.8440e-04 truth_out[1] =
> 1.5508e-01
> Compute output dual inner products
> output_dual_innerprods[0][0] = 8.3935e+01
> output_dual_innerprods[1][0] = 0.0000e+00
>
> ---- Performing Greedy basis enrichment ----
>
> ---- Basis dimension: 0 ----
> Performing RB solves on training set
> Maximum (absolute) error bound is 1.8090e+00
>
> Performing truth solve at parameter:
> load_Fx: 4.971986e+00
> load_Fy: 4.072318e+00
> .....
> -----------------------------------------------------------------------------------------------------------------
>
> I'm afraid the problem has something to do with using proper solver
> options, and I'd appreciate any comments/suggestions.
>
> K. Lee.
> ------------------------------------------------------------------------------
> Keep yourself connected to Go Parallel:
> BUILD Helping you discover the best ways to construct your parallel projects.
> http://goparallel.sourceforge.net
> _______________________________________________
> Libmesh-users mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/libmesh-users


------------------------------------------------------------------------------
Keep yourself connected to Go Parallel: 
BUILD Helping you discover the best ways to construct your parallel projects.
http://goparallel.sourceforge.net
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to