-
-
Notifications
You must be signed in to change notification settings - Fork 5.1k
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Segmentation fault in large scale lbfgs optimization #5102
Comments
Some of the Fortran parts of scipy use the Fortran |
Thanks for the reply! If the problem is integer overflow what is the maximal size of the derivative vector that can be used? Just to note that the optimizer usually performs several steps and only then it crashes. |
Here's the Fortran code: https://github.com/scipy/scipy/blob/master/scipy/optimize/lbfgsb/lbfgsb.f.
The expression In my opinion the best way to fix this would be to find compatibly licensed updates of this code elsewhere and re-introduce the code into scipy from that other source, if this exists. The |
@borislavsm You could try running https://github.com/stephenbeckr/L-BFGS-B-C. A header file has a note about integer sizes https://github.com/stephenbeckr/L-BFGS-B-C/blob/master/src/lbfgsb.h#L10. |
|
I know I'm a bit late to the party, but here goes. I downloaded the original L-BFGS-B code: http://users.iems.northwestern.edu/~nocedal/lbfgsb.html It has a test routine driver1.f which minimizes the generalized Rosenbrock function. The default test size is n = 25. I took @borislavsm's n = 100000000 and ran This is rather tentative, and I need to test it more, but perhaps the issue can be solved simply by changing how the code is compiled? Sorry if I'm missing something and am talking nonsense. |
Hello,
I receive a segmentation fault when I perform a large scale lbfgs optimisation in scipy. When I perform an optimisation with around 100000 parameters everything works OK, but when I do an optimisation with a significantly larger number of parameters (100000000) I get segfault.
The gdb trace is:
OS - CentOS
Python 2.7.5
numpy 1.9.2
scipy 0.12.1
Code that showcases the problem:
The text was updated successfully, but these errors were encountered: