Hello once again.
Due to a known bug on the function find_maximum_on_interval, one cannot use that function to find the global maximum of a function on an interval, Due to the fact that the function solve fails to provide all the solutions in the case of trigonometric functions, one cannot rely on solve(f.diff(x)==0), Due to the fact that I'm not able to use scipy.optimize.brute because Google don't find one single example of that function on the whole web, due to these facts, I'm in a trouble. How does people to find the global maximum of a function ? My aim is to compute bounding boxes of pspictures (I'm using python to produce pstricks code to be embed in LaTeX document), so I'm only interested in numerical answer; a rough precision of 0.05 is much sufficient. What I'm going to do is to take points in the interval separated by dx=0.05 and stupidly compute the value of my function on these points. In my case, it will be efficient, but it is quite frustrating to be reduced to that :( Any better idea ? I suppose that I'm not the only one to miss the fact to numerically compute a global maximum. Have a good night Laurent --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
