CVS log for OpenXM/src/asir-contrib/packages/src/tk_rk-mfac.rr

Up to [local] / OpenXM / src / asir-contrib / packages / src

Request diff between arbitrary revisions

Default branch: MAIN

Revision 1.7 / (download) - annotate - [select for diffs], Tue Nov 13 02:17:43 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
CVS Tags: HEAD
Changes since 1.6: +10 -3 lines
Diff to previous 1.6 (colored)

Bug fix: truncated step might be 0. In this case, minimal step size is decreased to minimal step size/10.
Bug: print_tex_form(phi) does not return \phi.

Revision 1.6 / (download) - annotate - [select for diffs], Sun Nov 11 01:20:22 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.5: +2 -2 lines
Diff to previous 1.5 (colored)

Bug fix of the previous commit on rounding the step size for RK on rationals.

Revision 1.5 / (download) - annotate - [select for diffs], Sat Nov 10 23:49:14 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.4: +16 -7 lines
Diff to previous 1.4 (colored)

The mesh size h is rounded to avoid the time variable becomes a big rational number
when the Runge-Kutta method is executed over rational numbers.
Option: h_round_size.

Revision 1.4 / (download) - annotate - [select for diffs], Thu Nov 8 10:25:19 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.3: +14 -2 lines
Diff to previous 1.3 (colored)

The option h_rel_err (default value is 10^(-3)) is for the truncation of
the step size H when the rational arithmetic is used (mode=0).

Revision 1.3 / (download) - annotate - [select for diffs], Thu Nov 8 02:53:34 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.2: +4 -11 lines
Diff to previous 1.2 (colored)

rk4_adaptive_step(RK_mat,Y0,X,X0,H,H1)
RK_mat for X and H. Y0 is the value at X=X0. Returns [[NewX0,NewY0],NewH].
Options: mode=0(rational), mode=1(bigfloat), mode=2(double, default).
Example: P=newmat(2,2,[[-2*x,1],[-1,-2*x]]); M=tk_rk.rk4_mat(P,x,x,h); tk_rk.rk4_adaptive_step(M,newvect(2,[1,0]),x,0,h,1/100);"

Revision 1.2 / (download) - annotate - [select for diffs], Tue Nov 6 01:52:39 2018 UTC (5 years, 6 months ago) by takayama
Branch: MAIN
Changes since 1.1: +110 -1 lines
Diff to previous 1.1 (colored)

rk4_adaptive_step(RK_mat,Y0,X,X0,H,H1)
RK_mat for X and H. Y0 is the value at X=X0. Returns [[NewX0,NewY0],NewH]
Options: mode=0 (rational), mode=1 (bigfloat), mode=2 (double)
abs_err, rel_err.

See tk_rk.test2c(Xmax) as for an example call this function.

Revision 1.1 / (download) - annotate - [select for diffs], Sun Mar 12 07:51:24 2017 UTC (7 years, 2 months ago) by takayama
Branch: MAIN

These new modules are for solving approximately linear ODE over Q.
See *test* functions as to examples. The command oxlistusage tk_bs.rr (or tk_rk-mfac.rr)
outputs short usages.

This form allows you to request diff's between any two revisions of a file. You may select a symbolic revision name using the selection box or you may type in a numeric name using the type-in text box.