Redian新闻
>
解ODE的六阶精度Runge-Kutta法
avatar
解ODE的六阶精度Runge-Kutta法# Computation - 科学计算
m*n
1
请推荐一个常用的格式。
如果你有8阶精度的,更好。
avatar
b*y
2
I don't know if it worthes to do higher order
RK at all.

【在 m****n 的大作中提到】
: 请推荐一个常用的格式。
: 如果你有8阶精度的,更好。

avatar
m*n
3
Surely, it is needed in some scientific calculations. I can give you examples.
J. Chem. Phys. 116, 538 (2002)
Sympletic Integration of Classical Trajectories: A case study
Ch. Schlier and A. Seiter.

【在 b*****y 的大作中提到】
: I don't know if it worthes to do higher order
: RK at all.

avatar
h*o
4
我觉得可以自己推呀, hoho...

【在 m****n 的大作中提到】
: 请推荐一个常用的格式。
: 如果你有8阶精度的,更好。

avatar
w*d
5
sometimes not so easy to get a good high order scheme.
As I remember, the exist highest A-stable(?) RK scheme is
6th order.
If you have some special requirement, e.g. the speed,
the robustness, it is even harder to get
a good one. I once used the 4th order A-stable RK scheme
with only 3 stages. It worked very well sometimes. But for some
problems, the round-off error just killed the solution.
some useful references:
general RK:
Solving ordinary differential equations :
Hairer, E. (Er

【在 h***o 的大作中提到】
: 我觉得可以自己推呀, hoho...
avatar
r*r
6

^^^^^^^^^
it is SSP instead of A-stable.
As for A-stable, RADAU, developed
by Hairer (2002), uses A-stable
implicit RK method of variable
orders (between 5, 9 and 13).
Numerical experiments show that,
for stiff problems

【在 w**d 的大作中提到】
: sometimes not so easy to get a good high order scheme.
: As I remember, the exist highest A-stable(?) RK scheme is
: 6th order.
: If you have some special requirement, e.g. the speed,
: the robustness, it is even harder to get
: a good one. I once used the 4th order A-stable RK scheme
: with only 3 stages. It worked very well sometimes. But for some
: problems, the round-off error just killed the solution.
: some useful references:
: general RK:

avatar
w*d
7
i c. The paper I read was published before 2002.

【在 r*r 的大作中提到】
:
: ^^^^^^^^^
: it is SSP instead of A-stable.
: As for A-stable, RADAU, developed
: by Hairer (2002), uses A-stable
: implicit RK method of variable
: orders (between 5, 9 and 13).
: Numerical experiments show that,
: for stiff problems

avatar
t*o
8
又多学一点
偶最多只用4阶R-K, 一直认为R-K公式是固定的,
u* = u^0 + dt*...
u** = u* + dt*...
想要多精度就是多写几次, u******, 系数不同而已
HOHO

【在 r*r 的大作中提到】
:
: ^^^^^^^^^
: it is SSP instead of A-stable.
: As for A-stable, RADAU, developed
: by Hairer (2002), uses A-stable
: implicit RK method of variable
: orders (between 5, 9 and 13).
: Numerical experiments show that,
: for stiff problems

avatar
d*q
9
taylor的
一种降阶形式
化为代数计算。。
不过用8阶也。。。【 在 twoxiao (老大,一点通) 的大作中提到: 】
相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。