# forum.alglib.net

ALGLIB forum
 It is currently Wed Jun 26, 2019 4:05 pm

 All times are UTC

### Forum rules

1. This forum can be used for discussion of both ALGLIB-related and general numerical analysis questions
2. This forum is English-only - postings in other languages will be removed.

 Page 1 of 1 [ 2 posts ]
 Print view Previous topic | Next topic
Author Message
 Post subject: Cash-Karp ODE solver needs a maxiter and/or min step sizePosted: Sat May 19, 2018 2:41 am

Joined: Sat May 19, 2018 2:14 am
Posts: 2
I've been trying to integrate rocket equation problems, and one thing that comes up is that for any rocket there is a very finite maximum theoretical burntime beyond which you've burned the last proton in the rocket, and accelerated it to infinite speed, infinitely far away. The value is also not much more than typical burntimes of a real rocket stage (usually about ~20% more). The behavior of the alglib ck integrator is that when you try to integrate beyond that timescale (which can happen when a wrapping optimizer is optimizing the time of the burn and takes too long of a step) is to just spin and never return (i suspect it is working very hard on a Xeno's Paradox kind of problem).

I can solve this reasonably easily in my own code by doing the calculation for when the rocket will burn completely and clipping the maximum x value passed to odesolverrkck(), but for other problems I'd imagine that could be computationally expensive or error prone to compute and it would better to have bounds in the routine itself. It also just feels unsafe to have an infinite loop as a possible failure mode.

Top

 Post subject: Re: Cash-Karp ODE solver needs a maxiter and/or min step sizPosted: Sun Jul 15, 2018 4:10 am

Joined: Sat May 19, 2018 2:14 am
Posts: 2
Another thing about the ODE solver is that it disallows integrating backwards in time and throws a rep.terminationtype of -2 "X is not ordered by ascending/descending".

The problem domain I'm interested in (spacecraft trajectory optimization) is fine with negative time integration of initial and terminal coast periods. They just correspond to either starting the maneuver earlier than estimated, or arriving at the target orbit later than estimated.

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 2 posts ]

 All times are UTC

#### Who is online

Users browsing this forum: Bing [Bot] and 1 guest

 You cannot post new topics in this forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum

Search for: