OpenXM/src/asir-contrib/packages/doc/tk_ode_by_mpfr/tk_ode_by_mpfr-ja.texi - diff

Return to tk_ode_by_mpfr-ja.texi CVS log

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

Diff for /OpenXM/src/asir-contrib/packages/doc/tk_ode_by_mpfr/tk_ode_by_mpfr-ja.texi between version 1.3 and 1.4

version 1.3, 2021/03/29 05:08:01

version 1.4, 2021/06/04 06:19:10

Line 1

%% $OpenXM: OpenXM/src/asir-contrib/packages/doc/tk_ode_by_mpfr/tk_ode_by_mpfr-ja.texi,v 1.2 2021/03/17 11:22:31 takayama Exp $

%% %% $OpenXM: OpenXM/src/asir-contrib/packages/doc/tk_ode_by_mpfr/tk_ode_by_mpfr-ja.texi,v 1.3 2021/03/29 05:08:01 takayama Exp $

%% debug に ln -s ../texinfo-ja.tex . をしておく. txi-ja.tex も?

%% xetex tk_ode_by_mpfr-ja.texi (.texi までつける. )

%% 以下コメントは @comment で始める. \input texinfo 以降は普通の tex 命令は使えない.

Line 21

@setfilename asir-contrib-tk_ode_by_mpfr

@comment --- タイトル ---

@settitle ox_pari

@settitle tk_ode_by_mpfr

@comment %**end of header

@comment %@setchapternewpage odd

Line 44

@title tk_ode_by_mpfr

@subtitle MPFR で ODE の数値解析をするC言語のコードを生成

@subtitle 1.0 版

@subtitle 2020 年 9 月 9 日

@subtitle 2021 年 6 月 4 日

@author by Nobuki Takayama

@page

Line 141 ODE dF/dt = P F の係数行列 P(t). t の式を成�

dF/dt = PF を F(T0)=F0 の初期条件の元, 時刻 T1 まで求める.

@item

big float による matrix factorial を用いて計算する.

仮数部(significand)のビットサイズは生成したプログラムに PREC の値として定義されている.

仮数部のサイズを変更するには生成された C プログラムの #define PREC 64 の部分を

変更するかこの関数の option prec を用いる.

@item

defusing heuristics や知られている値を用いて不安定性を

回避コードも含む.

Line 156 defusing heuristics や知られている値を用い�

Line 159 defusing heuristics や知られている値を用い�

@multitable {xxxxxxxxxx} {xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx} {xxxxxxxxxxxxxx}

@item Option @tab @tab default value

@item { verbose} @tab @tab 0

@item { prec} @tab doube size of MPFR @tab 64

@item { prec} @tab significand size of MPFR @tab 64

@item { progname} @tab @tab { tmp-test}

@item { h} @tab step size @tab 0.001

@item { t_noproj}@tab time to apply defusing @tab 0

Line 196 t=8.1 以降は --n_prune で指定した個数の固�

Line 199 t=8.1 以降は --n_prune で指定した個数の固�

を削除する.

--n_defuse 2000 は 2000 個の行列の matrix factorial を計算する.

@sp 3

@noindent

例:

H^k_n(1,t) k=10, n=1, 10000<=t<=10100 の計算とその値fと真の値Hとの relative error (f-H)/Hの計算.

"Hkn10000" は "h2" でもよい(alias).

@smallformat

@example

--> load("tk_ode_assert.rr");

--> tk_ode_assert.usage_assert(); // Usage of assert

--> tk_ode_assert.hkn2()

==> Compile tmp-test.c following the instruction.

==> ./tmp-test --go --output_stepsize 1 | grep '^gnuplot' | awk '{print $2,$3}' > t.txt

// Data is stored in t.txt.

--> tk_ode_assert.output_relative_error("Hkn10000","t.txt"|zoom=10^4301);

// The initial value for tmp-test is multiplied by 10^4301.

// Output tmp-rerror.txt

==> Start the gnuplot and plot "tmp-rerror.txt" w lp to show the relative error

@end example

@end smallformat

@noindent

@sp 3

@noindent

例: Airy や H^k_n(1,t) で 1<=t<=400 で同様.

@smallformat

@example

--> tk_ode_assert.airy1()

==> Run tmp-test

--> tk_ode_assert.output_relative_error("a","t.txt"|zoom=1);

--> tk_ode_assert.hkn1()

==> Run tmp-test

--> tk_ode_assert.output_relative_error("h","t.txt"|zoom=1);

@end example

@end smallformat

@noindent

@comment --- 参照(リンク)を書く ---

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>