[SciPy-dev] openopt and algencan

Nils Wagner nwagner at iam.uni-stuttgart.de
Thu Mar 20 15:49:32 EDT 2008


On Thu, 20 Mar 2008 20:32:53 +0200
  dmitrey <dmitrey.kroshko at scipy.org> wrote:
> Nils, try using p.xlabel = 'nf'
> 
> this is temporary name for "number of function 
>evaluations", I will 
> think what is better name for the case, and then commit 
>it into OO doc.
> 
> This is more appropriate vs your "iteration" that you 
>have chosen for 
> your x label.
> Your "iteration" is just a single objfunc evaluation, 
>while users will 
> suppose solvers have different iterations with different 
>number of 
> objFunc evaluations in each one.
> 
> Also, take into account, that ralg can yield infeasible 
>iteration 
> points, so your single-subplot picture doesn't provide 
>all convergence info.
> Ralg is intended to handle lb-ub problems more bette in 
>future.
> Regards, D.
> 
> 
> 
> Nils Wagner wrote:
>> On Thu, 20 Mar 2008 19:03:25 +0200
>>   dmitrey <dmitrey.kroshko at scipy.org> wrote:
>>   
>>> Hi Nils,
>>> I have made some changes to svn, try using ALGENCAN now.
>>>     
>>   
>> Great ! It works for me. Thank you very much !
>>
>> Nils
>> _______________________________________________
>> Scipy-dev mailing list
>> Scipy-dev at scipy.org
>> http://projects.scipy.org/mailman/listinfo/scipy-dev
>>
>>
>>
>>   
> 
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev

Hi Dmitrey,

Thank you for your comments.

BTW, do you have an
optional specification file algencan.dat
that you can send me ?

I just noticed the message

  The optional specification file algencan.dat was not 
found in the current
  directory (this is not a problem nor an error). The 
default values for the
  ALGENCAN parameters will be used.

  ALGENCAN PARAMETERS:
  rhoauto  =                    T
  rhotype  =                    1
  rhomult  =           1.0000E+01
  rhofrac  =           5.0000E-01
  gtype    =                    0
  hptype   =                    6
  intype   =                    0
  precond  =               QNCGNA
  checkder =                    F
  epsfeas  =           1.0000E-06
  epsopt   =           1.0000E-06
  maxoutit =                   50
  maxtotit =                 1000
  maxtotfc =                10000
  iprint   =                    0
  ncomp    =                    5
ALGENCAN has finished solving the problem unnamed
istop:  2 (|| gradient F(X[k]) || < gradtol)
Solver:   Time Elapsed = 0.21   CPU Time Elapsed = 0.2
objFunValue: 443.752171145 (feasible, max constraint =  0)

Nils

  



More information about the SciPy-Dev mailing list