[Ipopt-tickets] [Ipopt] #18: number of nonzeros in inequality Jacobian constraints

[Ipopt]

#18: number of nonzeros in inequality Jacobian constraints
-------------------+--------------------------------------------------------
Reporter:  andanh  |       Owner:  ipopt-team
    Type:  defect  |      Status:  new       
Priority:  normal  |   Component:  Ipopt     
 Version:          |    Severity:  normal    
Keywords:          |  
-------------------+--------------------------------------------------------
 Hi,

 I tested successfully many problems by using IPOPT linking PARDISO.  I can
 now find the optimal solution after few steps.  It is very great!

 I tested another bigger problem.  I met the following difficulty:
 The numbers of nonzeros in Jacobian constraints, which were output by
 IPOPT, is not  equal to NZ which I tranfered to the subroutine EV_JAC_G.
 Although in the previous examples, these values are the same. And IPOPT is
 terminated and produced the following, see the output [1]

 I tried to solve with ma27, the toal number of nonzeros is not equal to
 mine.  But IPOPT can go further (see the output [2])

 What should I do to use Pardiso?

 Thank you so much.

 Cheers.

 Danh.

 [1]
 ******************************************************************************
 This program contains Ipopt, a library for large-scale nonlinear
 optimization.
  Ipopt is released as open source code under the Common Public License
 (CPL).
          For more information visit http://projects.coin-or.org/Ipopt
 ******************************************************************************

 Number of nonzeros in equality constraint Jacobian...:  1107728
 Number of nonzeros in inequality constraint Jacobian.:   767633
 Number of nonzeros in Lagrangian Hessian.............:   486001

 Assertion failed: new_ia[dst->m] == (pint) nnz, file
 /scratch/oschenk/pardiso/src/smat.c, line 597
 IOT/Abort trap (core dumped)


 [2]
 [andanh at janus] /work/janus/andanh/Ipopt_2006Jul26/Ipopt/examples/Johnson1>
 ./pave
  INTYPE: 4  - Done! LM= 5647
  Number of nonzeros in [C]: NZ_CMAT= 1155360
  Successful
  SIZE_PROB:
  NSHAKE= 108001
  NSLACK= 72000
  N= 108001
  M= 77647
  NZ_MAX= 1659360
  NZ= 1659360
  NNZH_MAX= 486001
  NNZH= 486001
  SIZE_PROB
  Subroutine RHS_CONST
  M= 77647
  RHS_CONST
  START_POINT
  Here LM= 5647

 ******************************************************************************
 This program contains Ipopt, a library for large-scale nonlinear
 optimization.
  Ipopt is released as open source code under the Common Public License
 (CPL).
          For more information visit http://projects.coin-or.org/Ipopt
 ******************************************************************************

 Number of nonzeros in equality constraint Jacobian...:  1107728
 Number of nonzeros in inequality constraint Jacobian.:   767633
 Number of nonzeros in Lagrangian Hessian.............:   486001

 Total number of variables............................:   108001
                      variables with only lower bounds:        1
                 variables with lower and upper bounds:        0
                      variables with only upper bounds:        0
 Total number of equality constraints.................:     5647
 Total number of inequality constraints...............:    72000
         inequality constraints with only lower bounds:        0
    inequality constraints with lower and upper bounds:        0
         inequality constraints with only upper bounds:    72000

 iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
 alpha_pr ls
    0 -6.8916711e-02 2.74e+01 2.00e+00   0.0 0.00e+00    -  0.00e+00
 0.00e+000

-- 
Ticket URL: <https://projects.coin-or.org/Ipopt/ticket/18>
Ipopt <http://www.coin-or.org>
Ipopt Tickets