[Ipopt] Unpredictable results with mac OS X Lion vecLib

[Hecht Frédéric]

Thank,
I have make small patch in the IpBlas, because but in not on apple but in 
ipopt due to a miss understanding off the BLAS  limitation. 

the found  blas api in : 

http://www.mathkeisan.com/usersguide/man/ddot.html
http://www.mathkeisan.com/usersguide/man/dcopy.html
http://www.mathkeisan.com/usersguide/man/daxpy.html
remark, They are not correct API of the BLAS lib on net lib 

the  blas  API of DDOT , DCOPY, DAXPY
the  increment in array must be not zero 
otherwise the results will be unpredictable.



My patch are : 



This is a 5 minutes patch. 


the ffassert is just of the check  you can remove it. 

Le 13 mars 2012 à 13:50, Ray Zimmerman a écrit :

> I mentioned that Ipopt does not work with vecLib on Lion at https://projects.coin-or.org/Ipopt/wiki/Ipopt_on_Mac_OS_X. My solution was to use Goto BLAS, which performed almost as well in my tests on Snow Leopard.
> 
> -- 
> Ray Zimmerman
> Senior Research Associate
> 419A Warren Hall, Cornell University, Ithaca, NY 14853
> phone: (607) 255-9645
> 
> 
> 
> 
> On Mar 12, 2012, at 9:01 AM, Sylvain Auliac wrote:
> 
>> Hi,
>> Professor Hecht and I encountered a strange (and bad) behaviour of ipopt while trying to give it an interface in the FEM solver freefem++ developed in our laboratory. These were mainly memory issues. Mr Hecht later found that it was due to some incompatibilities in the file : Ipopt-3.10.2/Ipopt/src/LinAlg/IpBlas.cpp. When used with mac OS X Lion's vecLib BLAS, some of the functions are passed values which lead to unpredictable results (as said in the following page :  http://www.mathkeisan.com/usersguide/man/dcopy.html).
>> I don't know if someone already reported this, but having found nothing about that while doing a quick search, I thought it would be a good thing to do it. 
>> Best regards.
>> 
>> 
>> Sylvain Auliac
>> Laboratoire Jacques-Louis Lions
>> Université Pierre et Marie Curie
>> mail: auliac at ann.jussieu.fr
>> tel : 01 44 27 71 70
>> 
>> 
>> 
>> 
>> _______________________________________________
>> Ipopt mailing list
>> Ipopt at list.coin-or.org
>> http://list.coin-or.org/mailman/listinfo/ipopt
> 



Cordialement, 
Frédéric Hecht       
----------------------------------------------------------------------------------------------
Laboratoire  Jacques-Louis Lions,                        Université P&M Curie
BC187,  4 Place  Jussieu,     75252  PARIS cedex 05,                France 
Campus Jussieu,          Barre 15-25,           3 etage               Bureau 307
tel:  +33 1 44274411,    mob: +33 6 62198986,    fax: +33 1 44277200
mail & web: hecht at ann.jussieu.fr       http://www.ann.jussieu.fr/~hecht
software:  FreeFem++          web site:            http://www.freefem.org/ff++

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20120313/f4347fbe/attachment-0002.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: patch-IpBlas
Type: application/octet-stream
Size: 3012 bytes
Desc: not available
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20120313/f4347fbe/attachment-0001.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20120313/f4347fbe/attachment-0003.html>