---- Minimum Surface Area With Plate Problem -----
mx:8, my:6, bmx:3, bmy:3, height:0.2
iter =   0, Function value 1.67039, Residual: 0.47529 
iter =   1, Function value 1.506, Residual: 0.0373927 
iter =   2, Function value 1.50566, Residual: 1.78454e-05 
iter =   3, Function value 1.50566, Residual: < 1.0e-6 
Tao Object: 1 MPI process
  type: tron
  Total PG its: 9,  PG tolerance: 0.001 
  TaoLineSearch Object: 1 MPI process
    type: more-thuente
    maximum function evaluations=30
    tolerances: ftol=0.0001, rtol=1e-10, gtol=0.9
    total number of function evaluations=0
    total number of gradient evaluations=0
    total number of function/gradient evaluations=1
    using variable bounds
    Termination reason: 1
  KSP Object: 1 MPI process
    type: stcg
    maximum iterations=10000, initial guess is zero
    tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.
    left preconditioning
    using UNPRECONDITIONED norm type for convergence test
  PC Object: 1 MPI process
    type: ilu
      out-of-place factorization
      0 levels of fill
      tolerance for zero pivot 2.22045e-14
      matrix ordering: natural
      factor fill ratio given 1., needed 1.
        Factored matrix follows:
          Mat Object: 1 MPI process
            type: seqaij
            rows=40, cols=40
            package used to perform factorization: petsc
            total: nonzeros=228, allocated nonzeros=228
              not using I-node routines
    linear system matrix = precond matrix:
    Mat Object: 1 MPI process
      type: seqaij
      rows=40, cols=40
      total: nonzeros=228, allocated nonzeros=228
      total number of mallocs used during MatSetValues calls=0
        not using I-node routines
  total KSP iterations: 12
  Active Set subset type: subvec
  convergence tolerances: gatol=1e-05,   steptol=0.,   gttol=0.
  Residual in Function/Gradient:=1.4512e-11
  Objective value=1.50566
  total number of iterations=3,                          (max: 50)
  total number of function/gradient evaluations=14,      (max: 4000)
  total number of Hessian evaluations=3
  Solution converged:    ||g(X)|| <= gatol