My guess is that you haven't attached an "L2" matrix to the system?
If that's the case, refer to reduced_basis_ex3 about how to do that. David On 12/03/2012 07:58 PM, Kyunghoon Lee wrote: > Hi David, > > After running > > gdb ./beam_2d-dbg > > and > > run -ksp_type preonly -pc_type lu -pc_factor_mat_solver_package mumps > > in gdb. I got the following message: > > =================== > ..... > Compute output dual inner products > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Error in external library! > [0]PETSC ERROR: Error reported by MUMPS in numerical factorization phase: > INFO(1)=-9, INFO(2)=1327859 > ! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 4, Fri Oct 26 10:46:51 > CDT 2012 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: > /home/aeronova/Development/RB_models/beam/beam_2d_iso_plane_strain_basic1_transient/offline_model/beam_2d-dbg > on a arch-linu named KLee-Linux by aeronova Tue Dec 4 08:49:23 2012 > [0]PETSC ERROR: Libraries linked from > /home/aeronova/Development/local/lib/petsc/x86_64_linux/3.3_p4_real/lib > [0]PETSC ERROR: Configure run at Fri Nov 23 07:53:58 2012 > [0]PETSC ERROR: Configure options > --prefix=/home/aeronova/Development/local/lib/petsc/x86_64_linux/3.3_p4_real > --download-mpich=/home/aeronova//Development/local/share/petsc/mpich2-1.4.1p1.tar.gz > --download-superlu=1 --download-blacs=1 --download-metis=1 > --download-parmetis=1 --download-scalapack=1 --download-mumps=1 > --download-umfpack=1 --with-clanguage=C++ --with-shared-libraries=1 > --with-debugging=1 > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: MatFactorNumeric_MUMPS() line 670 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/mat/impls/aij/mpi/mumps/mumps.c > [0]PETSC ERROR: MatLUFactorNumeric() line 2803 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/mat/interface/matrix.c > [0]PETSC ERROR: PCSetUp_LU() line 160 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/ksp/pc/impls/factor/lu/lu.c > [0]PETSC ERROR: PCSetUp() line 832 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/ksp/pc/interface/precon.c > [0]PETSC ERROR: KSPSetUp() line 278 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/ksp/ksp/interface/itfunc.c > [0]PETSC ERROR: KSPSolve() line 402 in > /home/aeronova/Development/local/share/petsc/3.3-p4/src/ksp/ksp/interface/itfunc.c > [0]PETSC ERROR: solve() line 678 in > "unknowndirectory/"src/solvers/petsc_linear_solver.C > application called MPI_Abort(comm=0x84000000, 76) - process 0 > [unset]: aborting job: > application called MPI_Abort(comm=0x84000000, 76) - process 0 > [Inferior 1 (process 20258) exited with code 0114] > (gdb) bt > No stack. > > =================== > > K. Lee. > > On Mon, Dec 3, 2012 at 11:16 PM, David Knezevic > <[email protected]>wrote: > >> 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 >> > ------------------------------------------------------------------------------ > LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial > Remotely access PCs and mobile devices and provide instant support > Improve your efficiency, and focus on delivering more value-add services > Discover what IT Professionals Know. Rescue delivers > http://p.sf.net/sfu/logmein_12329d2d > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial Remotely access PCs and mobile devices and provide instant support Improve your efficiency, and focus on delivering more value-add services Discover what IT Professionals Know. Rescue delivers http://p.sf.net/sfu/logmein_12329d2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
