You could splice the ode-solver to the equation solver, although it might be slow. What sort of problem are you trying to solve? john Nils Wagner wrote:
Hi all,
Has someone implemented the shooting method in scipy ?
Nils
_______________________________________________ SciPy-user mailing list SciPy-user@scipy.org http://projects.scipy.org/mailman/listinfo/scipy-user
On 8/13/06, Stefan van der Walt <stefan@sun.ac.za> wrote:
On Sun, Aug 13, 2006 at 04:24:56PM +0200, Nils Wagner wrote:
Has someone implemented the shooting method in scipy ?
Is this the infamous shooting(self,in_foot) method?
I doubt it, as that method is implemented in every large code-base! Matthew
Matthew Brett wrote:
On 8/13/06, Stefan van der Walt <stefan@sun.ac.za> wrote:
On Sun, Aug 13, 2006 at 04:24:56PM +0200, Nils Wagner wrote:
Has someone implemented the shooting method in scipy ?
Is this the infamous shooting(self,in_foot) method?
I doubt it, as that method is implemented in every large code-base!
Matthew _______________________________________________ SciPy-user mailing list SciPy-user@scipy.org http://projects.scipy.org/mailman/listinfo/scipy-user
The main steps of the shooting method are - Transform the given boundary value problem into an initial value problem with estimated parameters - Adjust the parameters iteratively to reproduce the given boundary values A nice application of the shooting method is given in a recent paper by Coomer et al. "A non-linear eigenvalue problem associated with inextensible whirling strings" Journal of Sound and Vibration, Vol. 239 Issue 5 pp. 969-982 (2001) Nils
Hi! The main steps of the shooting method are
- Transform the given boundary value problem into an initial value problem with estimated parameters - Adjust the parameters iteratively to reproduce the given boundary values
I remember this.. http://www.cambridge.org/catalogue/catalogue.asp?isbn=0521852870&ss=res This is the book.. Numerical Methods in Engineering with Python from the Cambridge Univ Press. This book has a chapter on solving Two point boundary value problems. But since the book is meant to teach Numerical methods, he has developed all the necessary machinery including ODEint and Newton Raphson iteration all by himself using Numpy alone and uses those functions to write the code for solving BVP's. A similar thing will have to be done using Scipy's packages now I guess.. or are we expecting a code fully written in C, C++, FORTRAN ?? I could help if it were the former..not the latter. I have been wanting this for quite some time now !! Too basic a feature to be missing in the Scipy !! Ganesh -- Ganesh V Undergraduate student, Aerospace Engineering, IIT Madras, Chennai-36. My homepage --> http://www.ae.iitm.ac.in/~ae03b007
participants (5)
-
Ganesh V -
John Hassler -
Matthew Brett -
Nils Wagner -
Stefan van der Walt