I made a test case for debugging a potential bug in the Adjoint Refinement Error Estimator. However, I am encountering a somewhat bizarre issue on running it. It compiles fine in both dbg and opt code, but only runs properly in dbg mode.
The output is dbg mode is: Adaptive step 0, we have 256 active elements and 961 active dofs. |F| = 3662.286105743243 |J| = 1137.777777779091 *0 KSP Residual norm 1.317847130831e+01 * * 1 KSP Residual norm 1.260322028533e-14 * |F| = 1.019750353537797e-06 Nonlinear solver converged, step 0, residual reduction 3.93399e-10 < 1e-09 T(H1) = 4.66463 T(L2) = 0.413844 But on running in opt mode, the output is: Adaptive step 0, we have 256 active elements and 961 active dofs. |F| = 3662.286105743243 |J| = 1137.777777779091 * 0 KSP Residual norm -nan * |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 |F| = 3662.286105743243 Warning! Too many iterations used in Brent line search! It looks like the linear solver is just not doing anything. Anyone have an idea as to whats going on ? Thanks. -- Vikram Garg Postdoctoral Associate Predictive Engineering and Computational Science (PECOS) The University of Texas at Austin http://web.mit.edu/vikramvg/www/ http://vikramvgarg.wordpress.com/ http://www.runforindia.org/runners/vikramg ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
