I do have attached L2 as follows:
set_L2_assembly(transient_elasticity_assembly_expansion.M0_assembly);
and M0_assembly is from the below
struct M0: ElemAssembly {
// L2 matrix
virtual void interior_assembly(FEMContext &c) {
const unsigned int u_var = 0;
const std::vector<Real> &JxW =
c.element_fe_var[u_var]->get_JxW();
const std::vector<std::vector<Real> >& phi =
c.element_fe_var[u_var]->get_phi();
// The number of local degrees of freedom in each variable
const unsigned int n_u_dofs = c.dof_indices_var[u_var].size();
// Now we will build the affine operator
unsigned int n_qpoints = c.element_qrule->n_points();
for (unsigned int qp=0; qp<n_qpoints; qp++) {
for (unsigned int i=0; i<n_u_dofs; i++) {
for (unsigned int j=0; j<n_u_dofs; j++) {
c.elem_jacobian(i,j) += JxW[qp] *
phi[j][qp]*phi[i][qp];
}
}
} // end of qp loop
} // end of interior_assembly
};
K. Lee.
On Tue, Dec 4, 2012 at 9:00 AM, David Knezevic
<[email protected]>wrote:
> 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
>
------------------------------------------------------------------------------
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
