1*c4762a1bSJed Brown 2*c4762a1bSJed Brown---- Elastic-Plastic Torsion Problem ----- 3*c4762a1bSJed Brownmx: 10 my: 10 4*c4762a1bSJed Brown 5*c4762a1bSJed Browniter = 0, Function value -0.330579, Residual: 0.728961 6*c4762a1bSJed Brown ---------- Testing Hessian ------------- 7*c4762a1bSJed Brown Run with -tao_test_hessian_view and optionally -tao_test_hessian <threshold> to show difference 8*c4762a1bSJed Brown of hand-coded and finite difference Hessian entries greater than <threshold>. 9*c4762a1bSJed Brown Testing hand-coded Hessian, if (for double precision runs) ||J - Jfd||_F/||J||_F is 10*c4762a1bSJed Brown O(1.e-8), the hand-coded Hessian is probably correct. 11*c4762a1bSJed Brown ||H - Hfd||_F/||H||_F = 2.97218e-10, ||H - Hfd||_F = 1.31584e-08 12*c4762a1bSJed Browniter = 0, Function value -0.38665, Residual: 0.531327 13*c4762a1bSJed Brown ---------- Testing Hessian ------------- 14*c4762a1bSJed Brown ||H - Hfd||_F/||H||_F = 4.25681e-10, ||H - Hfd||_F = 1.88457e-08 15*c4762a1bSJed Browniter = 1, Function value -0.414886, Residual: 0.204735 16*c4762a1bSJed Brown ---------- Testing Hessian ------------- 17*c4762a1bSJed Brown ||H - Hfd||_F/||H||_F = 1.35328e-10, ||H - Hfd||_F = 5.99124e-09 18*c4762a1bSJed Browniter = 2, Function value -0.427745, Residual: 1.11757e-06 19