---- Minimum Surface Area Problem -----
mx: 10     my: 10   

iter =   0, Function value 1.46076, Residual: 0.207785 
iter =   1, Function value 1.41871, Residual: 0.0223274 
iter =   2, Function value 1.41847, Residual: 0.000373636 
iter =   3, Function value 1.41847, Residual: < 1.0e-6 
Tao Object: 2 MPI processes
  type: snes
  SNES Object: 2 MPI processes
    type: newtonls
    maximum iterations=50, maximum function evaluations=10000
    tolerances: relative=1e-08, absolute=0.0001, solution=1e-08
    total number of linear solver iterations=35
    total number of function evaluations=4
    norm schedule ALWAYS
    SNESLineSearch Object: 2 MPI processes
      type: bt
        interpolation: cubic
        alpha=1.000000e-04
      maxlambda=1.000000e+00, minlambda=1.000000e-12
      tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08
      maximum iterations=40
    KSP Object: 2 MPI processes
      type: gmres
        restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refinement
        happy breakdown tolerance=1e-30
      maximum iterations=15, initial guess is zero
      tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
      left preconditioning
      using PRECONDITIONED norm type for convergence test
    PC Object: 2 MPI processes
      type: bjacobi
        number of blocks = 2
        Local solver information for first block is in the following KSP and PC objects on rank 0:
        Use -ksp_view ::ascii_info_detail to display information for all blocks
        KSP Object: (sub_) 1 MPI process
          type: preonly
          maximum iterations=10000, initial guess is zero
          tolerances: relative=1e-05, absolute=1e-50, divergence=10000.
          left preconditioning
          not checking for convergence
        PC Object: (sub_) 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:
                Mat Object: (sub_) 1 MPI process
                  type: seqaij
                  rows=50, cols=50
                  package used to perform factorization: petsc
                  total: nonzeros=364, allocated nonzeros=364
                    not using I-node routines
          linear system matrix, which is also used to construct the preconditioner:
          Mat Object: (sub_) 1 MPI process
            type: seqaij
            rows=50, cols=50
            total: nonzeros=364, allocated nonzeros=364
            total number of mallocs used during MatSetValues calls=0
              not using I-node routines
      linear system matrix, which is also used to construct the preconditioner:
      Mat Object: 2 MPI processes
        type: mpiaij
        rows=100, cols=100
        total: nonzeros=784, allocated nonzeros=784
        total number of mallocs used during MatSetValues calls=0
  convergence tolerances: gatol=1e-08,   grtol=1e-05,    steptol=0.,   gttol=0.
  Residual in Function/Gradient:=9.23151e-08
  Objective value=1.41847
  total number of iterations=3,                          (max: 10000)
  total number of function evaluations=10,                  (max: unlimited)
  total number of function/gradient evaluations=4,      (max: unlimited)
  total number of Hessian evaluations=3
  Solution converged:    User Terminated