Solution should be f(1,1)=-2 
  1 TAO, Function value: 8.,  Residual: 24496.1 
  1 TAO, Function value: -5.43278,  Residual: 955.461 
  2 TAO, Function value: -1.91991,  Residual: 203.127 
  3 TAO, Function value: -1.89878,  Residual: 2.30255 
  4 TAO, Function value: -1.95114,  Residual: 0.143794 
  5 TAO, Function value: -1.99671,  Residual: 0.0248563 
  6 TAO, Function value: -1.99999,  Residual: 0.00206264 
  7 TAO, Function value: -2.,  Residual: 7.6949e-06 
  8 TAO, Function value: -2.,  Residual: 1.06901e-09 
  9 TAO, Function value: -2.,  Residual: 1.78328e-15 
 10 TAO, Function value: -2.,  Residual: 2.93737e-16 
 11 TAO, Function value: -2.,  Residual: 2.48253e-16 
 12 TAO, Function value: -2.,  Residual: 0. 
TAO solve converged due to CONVERGED_GATOL iterations 13
