search

chkwon / Complementarity.jl

provides a modeling interface for mixed complementarity problems (MCP) and math programs with equilibrium problems (MPEC) via JuMP
Other
75 stars 20 forks source link

This is not issue. Just asking for help in specifying and solving a model. #47

Closed Jovansam closed 3 years ago

Jovansam commented 3 years ago

Posted here asking for help but no response so far. I am actually now doubtful I will get any there, so I try here. Essentially, I want to solve a complementarity problem derived from a cake eating model. This link provides for examples solved in GAMS. I would like to do the same in Julia.

I paste example code here that I am working on. If someone can help me make this work, I will appreciate a lot.


using PATH, JuMP, Complementarity

function geological_constraints_MCP()

    Α = 50.0
    Β = 0.1
    δ = 0.05
    R₀ = 0.5
    T = 10
    t = 1:T
    AA = Α*(1+δ)

    gxt =MCPModel()

    @variable(gxt, q[i in 1:T]>=0)
    @variable(gxt, λ[i in 1:T]>=0)
    @variable(gxt, R[i in 1:T]>=0)

    @mapping(gxt, FOCq[i in 1:T], - (Α - Β*q[i] - λ[i]))
    @complementarity(gxt,FOCq,q)
    @mapping(gxt, FOCR[i in 1:T], - ((i<T ? λ[i+1] : AA) - (λ[i]*(1 + δ))))
    @complementarity(gxt,FOCR,λ)
    @mapping(gxt, FOCλ[i in 1:T], - (R[i] + q[i] - (i>1 ? R[i-1] :  R₀)))
    @complementarity(gxt,FOCλ,R)

    print(gxt)
    status = solveMCP(gxt)
    @show result_value.(q)
    @show result_value.(R)
    @show result_value.(λ)

end

geological_constraints_MCP()
chkwon commented 3 years ago

I just pasted your code and it worked in my machine. Try ] up.


Feasibility
Subject to
 q[1] ≥ 0.0
 q[2] ≥ 0.0
 q[3] ≥ 0.0
 q[4] ≥ 0.0
 q[5] ≥ 0.0
 q[6] ≥ 0.0
 q[7] ≥ 0.0
 q[8] ≥ 0.0
 q[9] ≥ 0.0
 q[10] ≥ 0.0
 λ[1] ≥ 0.0
 λ[2] ≥ 0.0
 λ[3] ≥ 0.0
 λ[4] ≥ 0.0
 λ[5] ≥ 0.0
 λ[6] ≥ 0.0
 λ[7] ≥ 0.0
 λ[8] ≥ 0.0
 λ[9] ≥ 0.0
 λ[10] ≥ 0.0
 R[1] ≥ 0.0
 R[2] ≥ 0.0
 R[3] ≥ 0.0
 R[4] ≥ 0.0
 R[5] ≥ 0.0
 R[6] ≥ 0.0
 R[7] ≥ 0.0
 R[8] ≥ 0.0
 R[9] ≥ 0.0
 R[10] ≥ 0.0
With NL expressions
 subexpression[1]: -(((50.0 - 0.1 * q[1]) - λ[1]))
 subexpression[2]: -(((50.0 - 0.1 * q[2]) - λ[2]))
 subexpression[3]: -(((50.0 - 0.1 * q[3]) - λ[3]))
 subexpression[4]: -(((50.0 - 0.1 * q[4]) - λ[4]))
 subexpression[5]: -(((50.0 - 0.1 * q[5]) - λ[5]))
 subexpression[6]: -(((50.0 - 0.1 * q[6]) - λ[6]))
 subexpression[7]: -(((50.0 - 0.1 * q[7]) - λ[7]))
 subexpression[8]: -(((50.0 - 0.1 * q[8]) - λ[8]))
 subexpression[9]: -(((50.0 - 0.1 * q[9]) - λ[9]))
 subexpression[10]: -(((50.0 - 0.1 * q[10]) - λ[10]))
 subexpression[11]: -((λ[2] - λ[1] * (1.0 + 0.05)))
 subexpression[12]: -((λ[3] - λ[2] * (1.0 + 0.05)))
 subexpression[13]: -((λ[4] - λ[3] * (1.0 + 0.05)))
 subexpression[14]: -((λ[5] - λ[4] * (1.0 + 0.05)))
 subexpression[15]: -((λ[6] - λ[5] * (1.0 + 0.05)))
 subexpression[16]: -((λ[7] - λ[6] * (1.0 + 0.05)))
 subexpression[17]: -((λ[8] - λ[7] * (1.0 + 0.05)))
 subexpression[18]: -((λ[9] - λ[8] * (1.0 + 0.05)))
 subexpression[19]: -((λ[10] - λ[9] * (1.0 + 0.05)))
 subexpression[20]: -((52.5 - λ[10] * (1.0 + 0.05)))
 subexpression[21]: -(((R[1] + q[1]) - 0.5))
 subexpression[22]: -(((R[2] + q[2]) - R[1]))
 subexpression[23]: -(((R[3] + q[3]) - R[2]))
 subexpression[24]: -(((R[4] + q[4]) - R[3]))
 subexpression[25]: -(((R[5] + q[5]) - R[4]))
 subexpression[26]: -(((R[6] + q[6]) - R[5]))
 subexpression[27]: -(((R[7] + q[7]) - R[6]))
 subexpression[28]: -(((R[8] + q[8]) - R[7]))
 subexpression[29]: -(((R[9] + q[9]) - R[8]))
 subexpression[30]: -(((R[10] + q[10]) - R[9]))
Could not open options file: path.opt
Using defaults.
30 row/cols, 136 non-zeros, 15.11% dense.

Path 4.7.03 (Thu Jan 24 15:44:12 2013)
Written by Todd Munson, Steven Dirkse, and Michael Ferris

INITIAL POINT STATISTICS
Maximum of X. . . . . . . . . .  0.0000e+00 var: (q[1])
Maximum of F. . . . . . . . . .  5.2500e+01 eqn: (FOCR[(10,)])
Maximum of Grad F . . . . . . .  1.0500e+00 eqn: (FOCR[(1,)])
                                            var: (λ[1])

INITIAL JACOBIAN NORM STATISTICS
Maximum Row Norm. . . . . . . .  3.0000e+00 eqn: (FOCλ[(2,)])
Minimum Row Norm. . . . . . . .  1.0500e+00 eqn: (FOCR[(10,)])
Maximum Column Norm . . . . . .  3.0500e+00 var: (λ[2])
Minimum Column Norm . . . . . .  1.0000e+00 var: (R[10])

Crash Log
major  func  diff  size  residual    step       prox   (label)
    0     0             3.3320e+02             0.0e+00 (FOCR[(10,)])
    1     7    20    11 3.2280e+02  6.9e-02    0.0e+00 (FOCR[(10,)])
pn_search terminated: no progress.

Major Iteration Log
major minor  func  grad  residual    step  type prox    inorm  (label)
    0     0    20     2 3.2280e+02           I 0.0e+00 9.4e+01 (FOCR[(10,)])
    1    22    21     3 3.5530e+02  1.0e+00 RM 0.0e+00 9.1e+01 (FOCλ[(3,)])
    2    12    22     4 3.5530e+02  1.0e+00 RO 0.0e+00 9.1e+01 (FOCλ[(2,)])
    3    12    23     5 3.5530e+02  1.0e+00 RO 0.0e+00 9.1e+01 (FOCλ[(3,)])
    4    12    24     6 3.5530e+02  1.0e+00 RO 0.0e+00 9.1e+01 (FOCλ[(3,)])
    5    14    25     7 3.5530e+02  1.0e+00 RO 0.0e+00 9.1e+01 (FOCλ[(2,)])
    6    22    30     8 2.5398e+02  4.1e-01 RG 0.0e+00 2.1e+02 (FOCR[(9,)])
    7    10    31     9 2.1635e+02  1.0e+00 RO 0.0e+00 1.8e+02 (FOCR[(9,)])
    8     6    32    10 2.1635e+02  1.0e+00 RO 0.0e+00 1.8e+02 (FOCR[(9,)])
    9     6    33    11 2.1635e+02  1.0e+00 RO 0.0e+00 1.8e+02 (FOCR[(9,)])
   10     6    34    12 2.1635e+02  1.0e+00 RO 0.0e+00 1.8e+02 (FOCR[(9,)])
   11    22    35    13 3.5530e+02  1.0e+00 RL 0.0e+00 9.1e+01 (FOCλ[(3,)])
   12    12    42    14 1.3434e+02  1.7e-01 RG 0.0e+00 1.2e+02 (FOCR[(8,)])
   13     6    43    15 1.3434e+02  1.0e+00 RO 0.0e+00 1.2e+02 (FOCR[(8,)])
   14     6    44    16 1.3434e+02  1.0e+00 RO 0.0e+00 1.2e+02 (FOCR[(8,)])
   15     6    45    17 1.3434e+02  1.0e+00 RO 0.0e+00 1.2e+02 (FOCR[(8,)])
   16    22    46    18 3.5530e+02  1.0e+00 RL 0.0e+00 9.1e+01 (FOCλ[(3,)])
   17    12    47    19 3.5530e+02  1.0e+00 RD 0.0e+00 9.1e+01 (FOCλ[(2,)])
   18    12    53    20 1.3388e+02  2.6e-01 RG 0.0e+00 1.1e+02 (FOCR[(7,)])
   19     7    54    21 1.3388e+02  1.0e+00 RO 0.0e+00 1.1e+02 (FOCR[(7,)])
   20     7    55    22 1.3388e+02  1.0e+00 RO 0.0e+00 1.1e+02 (FOCR[(7,)])
   21    22    56    23 3.5530e+02  1.0e+00 RL 0.0e+00 9.1e+01 (FOCλ[(3,)])
   22    12    57    24 3.5530e+02  1.0e+00 RD 0.0e+00 9.1e+01 (FOCλ[(2,)])
   23    12    58    25 3.5530e+02  1.0e+00 RD 0.0e+00 9.1e+01 (FOCλ[(3,)])
   24    12    64    26 1.1471e+02  2.6e-01 RG 0.0e+00 9.5e+01 (FOCR[(6,)])
   25     8    65    27 1.1471e+02  1.0e+00 RO 0.0e+00 9.5e+01 (FOCR[(6,)])
   26    22    68    28 2.1146e+02  6.4e-01 RB 0.0e+00 6.2e+01 (FOCλ[(7,)])
   27    15    69    29 2.1146e+02  1.0e+00 RD 0.0e+00 6.2e+01 (FOCλ[(7,)])
   28    15    70    30 2.1146e+02  1.0e+00 RD 0.0e+00 6.2e+01 (FOCλ[(7,)])
   29    15    71    31 2.1146e+02  1.0e+00 RD 0.0e+00 6.2e+01 (FOCλ[(7,)])
   30    15    77    32 9.4054e+01  2.6e-01 RG 0.0e+00 7.9e+01 (FOCR[(5,)])
   31    22    82    33 1.1686e+02  3.3e-01 RB 0.0e+00 5.3e+01 (FOCR[(5,)])
   32    12    83    34 1.1686e+02  1.0e+00 RO 0.0e+00 5.3e+01 (FOCR[(5,)])
   33    12    84    35 1.1686e+02  1.0e+00 RO 0.0e+00 5.3e+01 (FOCR[(5,)])
   34    12    85    36 1.1686e+02  1.0e+00 RO 0.0e+00 5.3e+01 (FOCR[(5,)])
   35    12    86    37 1.1686e+02  1.0e+00 RO 0.0e+00 5.3e+01 (FOCR[(5,)])

Restart Log
proximal_perturbation 0
crash_method none
crash_perturb yes
nms_initial_reference_factor 2
lemke_start_type slack
proximal_perturbation 1.0000e-01

Major Iteration Log
major minor  func  grad  residual    step  type prox    inorm  (label)
   35     0    93    38 3.3320e+02           R 1.0e-01 1.0e+02 (FOCR[(10,)])
   36    13    94    39 3.3320e+02  1.0e+00 RO 4.0e-02 1.0e+02 (FOCR[(10,)])
   37    13    95    40 3.3320e+02  1.0e+00 RO 1.6e-02 1.0e+02 (FOCR[(10,)])
   38    13    96    41 3.3320e+02  1.0e+00 RO 6.4e-03 1.0e+02 (FOCR[(10,)])
   39    13    97    42 3.3320e+02  1.0e+00 RO 2.6e-03 1.0e+02 (FOCR[(10,)])
   40    13    98    43 3.3320e+02  1.0e+00 RO 1.0e-03 1.0e+02 (FOCR[(10,)])
   41    22   104    44 2.4694e+02  2.6e-01 RG 4.1e-04 1.3e+02 (FOCR[(9,)])
   42    12   105    45 2.4694e+02  1.0e+00 RO 1.6e-04 1.3e+02 (FOCR[(9,)])
   43    10   106    46 2.4694e+02  1.0e+00 RO 6.6e-05 1.3e+02 (FOCR[(9,)])
   44    12   107    47 2.4694e+02  1.0e+00 RO 2.6e-05 1.3e+02 (FOCR[(9,)])
   45    10   108    48 2.4694e+02  1.0e+00 RO 1.0e-05 1.3e+02 (FOCR[(9,)])
   46    22   110    49 2.9786e+02  8.0e-01 RB 4.2e-06 8.5e+01 (FOCλ[(9,)])
   47    14   115    50 2.4322e+02  4.1e-01 RG 4.1e-04 2.1e+02 (FOCR[(8,)])
   48    13   116    51 2.2607e+02  1.0e+00 RM 1.6e-04 1.9e+02 (FOCR[(8,)])
   49     9   117    52 2.2607e+02  1.0e+00 RO 6.6e-05 1.9e+02 (FOCR[(8,)])
   50     9   118    53 2.2607e+02  1.0e+00 RO 2.6e-05 1.9e+02 (FOCR[(8,)])
   51    22   120    54 2.8794e+02  8.0e-01 RB 1.0e-05 7.9e+01 (FOCq[(10,)])
   52    13   121    55 2.8794e+02  1.0e+00 RD 4.2e-06 7.9e+01 (FOCq[(10,)])
   53    13   128    56 1.4097e+02  1.7e-01 RG 2.6e-05 1.3e+02 (FOCR[(7,)])
   54    11   129    57 1.4097e+02  1.0e+00 RO 1.0e-05 1.3e+02 (FOCR[(7,)])
   55    11   130    58 1.4097e+02  1.0e+00 RO 4.2e-06 1.3e+02 (FOCR[(7,)])
   56    22   133    59 2.2016e+02  6.4e-01 RB 1.7e-06 6.2e+01 (FOCλ[(2,)])
   57    14   134    60 2.2016e+02  1.0e+00 RO 6.7e-07 6.2e+01 (FOCλ[(2,)])
   58    14   135    61 2.2016e+02  1.0e+00 RO 2.7e-07 6.2e+01 (FOCλ[(2,)])
   59    14   142    62 8.2542e+01  1.7e-01 RG 2.6e-05 7.2e+01 (FOCR[(6,)])
   60    11   143    63 8.2542e+01  1.0e+00 RO 1.0e-05 7.2e+01 (FOCR[(6,)])
   61    22   146    64 2.1357e+02  6.4e-01 RB 4.2e-06 6.1e+01 (FOCλ[(9,)])
   62    14   147    65 2.1357e+02  1.0e+00 RO 1.7e-06 6.1e+01 (FOCλ[(9,)])
   63    14   148    66 2.1357e+02  1.0e+00 RO 6.7e-07 6.1e+01 (FOCλ[(9,)])
   64    14   149    67 2.1357e+02  1.0e+00 RO 2.7e-07 6.1e+01 (FOCλ[(9,)])
   65    14   155    68 6.9972e+01  2.6e-01 RG 1.0e-05 5.0e+01 (FOCR[(5,)])
   66    22   158    69 2.1034e+02  6.4e-01 RB 4.2e-06 6.0e+01 (FOCλ[(8,)])
   67    14   159    70 2.1034e+02  1.0e+00 RO 1.7e-06 6.0e+01 (FOCλ[(8,)])
   68    14   160    71 2.1034e+02  1.0e+00 RO 6.7e-07 6.0e+01 (FOCλ[(8,)])
   69    14   161    72 2.1034e+02  1.0e+00 RO 2.7e-07 6.0e+01 (FOCλ[(8,)])
   70    14   162    73 2.1034e+02  1.0e+00 RD 1.1e-07 6.0e+01 (FOCλ[(8,)])

Restart Log
proximal_perturbation 0
crash_method none
crash_perturb no
nms_initial_reference_factor 10
nms_memory_size 2
nms_mstep_frequency 1
lemke_search_type slack

Major Iteration Log
major minor  func  grad  residual    step  type prox    inorm  (label)
   70     0   169    74 3.3320e+02           R 0.0e+00 1.0e+02 (FOCR[(10,)])
   71    13   170    75 3.3320e+02  1.0e+00 RM 0.0e+00 1.0e+02 (FOCR[(10,)])
   72    13   171    76 3.3320e+02  1.0e+00 RM 0.0e+00 1.0e+02 (FOCR[(10,)])
   73    13   172    77 3.3320e+02  1.0e+00 RM 0.0e+00 1.0e+02 (FOCR[(10,)])
   74    13   173    78 3.3320e+02  1.0e+00 RM 0.0e+00 1.0e+02 (FOCR[(10,)])
   75    13   174    79 3.3320e+02  1.0e+00 RM 0.0e+00 1.0e+02 (FOCR[(10,)])
   76    22   180    80 2.4694e+02  2.6e-01 RG 0.0e+00 1.3e+02 (FOCR[(9,)])
   77    10   181    81 2.4694e+02  1.0e+00 RM 0.0e+00 1.3e+02 (FOCR[(9,)])
   78    10   194    82 2.4694e+02  4.4e-05 RB 0.0e+00 1.3e+02 (FOCR[(9,)])
   79    10   207    83 2.4694e+02  4.4e-05 RB 0.0e+00 1.3e+02 (FOCR[(9,)])
   80    10   220    84 2.4694e+02  4.4e-05 RB 0.0e+00 1.3e+02 (FOCR[(9,)])
   81    22   225    85 2.4422e+02  3.3e-01 RB 0.0e+00 9.0e+01 (FOCR[(9,)])
   82    12   230    86 1.9635e+02  4.1e-01 RG 0.0e+00 1.4e+02 (FOCR[(8,)])
   83    12   231    87 1.9475e+02  1.0e+00 RM 0.0e+00 1.4e+02 (FOCR[(8,)])
   84    12   232    88 1.8824e+02  1.0e+00 RM 0.0e+00 1.4e+02 (FOCR[(8,)])
   85    12   233    89 1.8824e+02  1.0e+00 RM 0.0e+00 1.4e+02 (FOCR[(8,)])
   86    22   239    90 1.8255e+02  2.1e-01 RB 0.0e+00 1.1e+02 (FOCR[(8,)])
   87    12   240    91 1.8255e+02  1.0e+00 RM 0.0e+00 1.1e+02 (FOCR[(8,)])
   88    12   246    92 1.7507e+02  2.6e-01 RG 0.0e+00 1.1e+02 (FOCR[(7,)])
   89    12   247    93 1.7446e+02  1.0e+00 RM 0.0e+00 1.1e+02 (FOCR[(7,)])
   90    12   248    94 1.7446e+02  1.0e+00 RM 0.0e+00 1.1e+02 (FOCR[(7,)])
   91    22   254    95 1.6662e+02  2.1e-01 RB 0.0e+00 8.8e+01 (FOCR[(7,)])
   92    12   255    96 1.6662e+02  1.0e+00 RM 0.0e+00 8.8e+01 (FOCR[(7,)])
   93    12   270    97 1.6662e+02  3.5e-08 RB 0.0e+00 8.8e+01 (FOCR[(7,)])
   94    12   276    98 1.3934e+02  2.6e-01 RG 0.0e+00 8.0e+01 (FOCR[(6,)])
   95    13   277    99 1.3895e+02  1.0e+00 RM 0.0e+00 8.0e+01 (FOCR[(6,)])
   96    22   285   100 1.3777e+02  5.5e-02 RB 0.0e+00 7.5e+01 (FOCR[(6,)])
   97    12   286   101 1.3772e+02  1.0e+00 RM 0.0e+00 7.5e+01 (FOCR[(6,)])
   98    12   287   102 1.3772e+02  1.0e+00 RM 0.0e+00 7.5e+01 (FOCR[(6,)])
   99    12   288   103 1.3772e+02  1.0e+00 RM 0.0e+00 7.5e+01 (FOCR[(6,)])
  100    13   294   104 1.0859e+02  2.6e-01 RG 0.0e+00 6.0e+01 (FOCR[(5,)])
  101    22   299   105 1.3462e+02  3.3e-01 RB 0.0e+00 4.8e+01 (FOCλ[(8,)])

Restart Log
proximal_perturbation 0
crash_method pnewton
crash_nbchange_limit 10
nms_initial_reference_factor 2
crash_searchtype arc
nms_searchtype arc
gradient_searchtype arc
lemke_search_type slack

Crash Log
major  func  diff  size  residual    step       prox   (label)
    0     0             3.3320e+02             0.0e+00 (FOCR[(10,)])
    1     7    20    11 3.2280e+02  6.9e-02    0.0e+00 (FOCR[(10,)])
pn_search terminated: no progress.

Major Iteration Log
major minor  func  grad  residual    step  type prox    inorm  (label)
  101     0   339   107 3.2280e+02           R 0.0e+00 9.4e+01 (FOCR[(10,)])
  102    14   340   108 3.2280e+02  1.0e+00 RM 0.0e+00 9.4e+01 (FOCR[(10,)])
  103    14   341   109 3.2280e+02  1.0e+00 RM 0.0e+00 9.4e+01 (FOCR[(10,)])
  104    14   342   110 3.2280e+02  1.0e+00 RM 0.0e+00 9.4e+01 (FOCR[(10,)])
  105    14   361   111 2.5398e+02  4.1e-01 RG 0.0e+00 2.1e+02 (FOCR[(9,)])
  106    10   362   112 2.1635e+02  1.0e+00 RM 0.0e+00 1.8e+02 (FOCR[(9,)])
  107    22   369   113 1.3434e+02  1.7e-01 RG 0.0e+00 1.2e+02 (FOCR[(8,)])
  108     6   370   114 1.3434e+02  1.0e+00 RM 0.0e+00 1.2e+02 (FOCR[(8,)])
  109     6   371   115 1.3434e+02  1.0e+00 RM 0.0e+00 1.2e+02 (FOCR[(8,)])
  110     6   372   116 1.3434e+02  1.0e+00 RM 0.0e+00 1.2e+02 (FOCR[(8,)])
  111     6   373   117 1.3434e+02  1.0e+00 RM 0.0e+00 1.2e+02 (FOCR[(8,)])
  112    22   379   118 1.2404e+02  2.1e-01 RB 0.0e+00 9.2e+01 (FOCR[(8,)])
  113    11   385   119 1.0491e+02  2.6e-01 RG 0.0e+00 8.1e+01 (FOCR[(7,)])
  114    10   386   120 1.0448e+02  1.0e+00 RM 0.0e+00 8.1e+01 (FOCR[(7,)])
  115    10   387   121 1.0445e+02  1.0e+00 RM 0.0e+00 8.0e+01 (FOCR[(7,)])
  116    10   388   122 1.0445e+02  1.0e+00 RM 0.0e+00 8.0e+01 (FOCR[(7,)])
  117    22   396   123 1.0297e+02  5.5e-02 RB 0.0e+00 7.6e+01 (FOCR[(7,)])
  118    10   397   124 1.0297e+02  1.0e+00 RM 0.0e+00 7.6e+01 (FOCR[(7,)])
  119    10   403   125 7.6392e+01  2.6e-01 RG 0.0e+00 5.8e+01 (FOCR[(6,)])
  120     9   404   126 7.6428e+01  1.0e+00 RM 0.0e+00 5.8e+01 (FOCR[(6,)])
  121     9   424   127 4.9384e+01  2.6e-01 RG 0.0e+00 3.6e+01 (FOCR[(5,)])

FINAL STATISTICS
Inf-Norm of Complementarity . .  1.7820e+03 eqn: (FOCR[(6,)])
Inf-Norm of Normal Map. . . . .  2.9888e+01 eqn: (FOCq[(6,)])
Inf-Norm of Minimum Map . . . .  2.8624e+01 eqn: (FOCR[(5,)])
Inf-Norm of Fischer Function. .  3.6419e+01 eqn: (FOCR[(5,)])
Inf-Norm of Grad Fischer Fcn. .  7.0381e+01 eqn: (FOCR[(6,)])
Two-Norm of Grad Fischer Fcn. .  1.1243e+02

FINAL POINT STATISTICS
Maximum of X. . . . . . . . . .  7.9715e+01 var: (λ[6])
Maximum of F. . . . . . . . . .  2.9888e+01 eqn: (FOCq[(6,)])
Maximum of Grad F . . . . . . .  1.0500e+00 eqn: (FOCR[(1,)])
                                            var: (λ[1])

 ** EXIT - other error.

Major Iterations. . . . 121
Minor Iterations. . . . 1711
Restarts. . . . . . . . 3
Crash Iterations. . . . 2
Gradient Steps. . . . . 24
Function Evaluations. . 443
Gradient Evaluations. . 127
Basis Time. . . . . . . 0.005547
Total Time. . . . . . . 0.243248
Residual. . . . . . . . 4.938402e+01
result_value.(q) = [0.13314581932310116, 0.0, 3.297604972219143, 0.43501426593664494, 1.4507193921845492, 1.728811012326308, 0.0, 0.0006232525850461908, 0.0, 3.65934977817265]
result_value.(R) = [13.064075308561971, 2.32675119679776, 5.988229839631646, 1.8425076882914881, 2.4685911138888557, 5.302465786515632, 2.0783714858701843, 3.635929527247968, 3.091715461652764, 0.0]
result_value.(λ) = [49.542862590599576, 47.86867958176467, 49.40148711615427, 48.21034559783544, 48.658118572699145, 79.71476550261141, 61.34529397987332, 65.27078901319649, 63.90567638783797, 62.74734431905372]
10-element Array{Float64,1}:
 49.542862590599576
 47.86867958176467
 49.40148711615427
 48.21034559783544
 48.658118572699145
 79.71476550261141
 61.34529397987332
 65.27078901319649
 63.90567638783797
 62.74734431905372