noro_grcrt

noro_grcrt User’s Manual

Edition 1.0

Feb 2019

Copyright © Masayuki Noro 2019. All rights reserved.

1 CRTパッケージ noro_grcrt.rr

このマニュアルでは, asir-contrib パッケージに収録されている, CRTパッケージ ‘noro_grcrt.rr’ について解説する. このパッケージを使うには, まず ‘noro_grcrt.rr’ をロードする.

[1831] load("noro_grcrt.rr");

このパッケージの函数を呼び出すには, 全て noro_grcrt. を先頭につける. このマニュアルでは, 関連する組み込み関数についても解説する.

1.1 CRTによるグレブナー基底計算

1.1.1 f4_chr

noro_grcrt.f4_chr_dec(b,v,ord[|options])

CRT によりグレブナー基底候補を求める.

return

多項式リスト

b

多項式リスト

v

変数リスト

ord

項順序型

options

下の説明参照.

• 有限体上のグレブナー基底を nd_f4 で計算し, それらを CRT で貼り合わせることで, 有理数体上のグレブナー基底候補を求める.
• 正当性のチェックは全く行わない.
• オプション homo=1 を指定すると, 有限体上のグレブナー基底計算を⻫次化付きて計算する.
• オプション weight を指定すると, weight 付きて計算する.
• オプション proc を指定すると, 有限体上のグレブナー基底計算を並列行う. proc て指定するのは ox_asir のストリーム番号のリストで, noro_grcrt.init_pproc により生成する.

[2465] P=noro_grcrt.init_pprocs([[0,4]]|nox=1)$
[2466] B=cyclic(8)$
[2467] V=[c0,c1,c2,c3,c4,c5,c6,c7]$
[2468] G=noro_grcrt.f4_chr(B,V,0|proc=P)$
...
[Template,0,CRT,5.812,IR,11.388,#P,56]
28.74sec(55.12sec)

1.1.2 noro_grcrt.init_pprocs

noro_grcrt.init_pprocs(list[|nox=0|1])

並列計算用サーバを起動する.

return

ストリーム番号のリスト

list

リスト

• ox_asir をまとめて起動する.
• 引数 listは [[host1,n1],[host2,n2],...]の形のリストで, hosti に ni 個の ox_asir を起動することを意味する. ここで hosti=0 の場合, \code{asirが実行されている マシンを意味する.
• オプション nox=1 が指定された場合, サーバ用のメッセージウィンドウは表示されない. デバッグ中は このオプションを指定しない方が安全である.

[2465] P=noro_grcrt.init_pprocs([[0,2],["shinohara",2]]|nox=1);
[0,1,2,3]

1.2 CRTに関する関数

1.2.1 noro_grcrt.dp_chrem2

noro_grcrt.dp_chrem2(g,M,gp,p,stat)

CRTによる貼り合わせを実行する.

return

チェックに失敗した個数

g

分散多項式の配列

M

正整数

gp

分散多項式の配列

p

正整数

stat

配列

• g mod M, gp mod p を CRT で貼り合わせた g1 を求め, g に上書きする. g1 は g1 = g mod M, g1 = gp mod gp を満たす.
• p は小さい整数(31bit 未満) でなければなrない.
• g[i], gp[i] は, サポートが一致した分散多項式でなければならない. gp[i], および stat[i]=0,1 の場合のg[i] はモニックでなければならない.
• statは長さがgの長さが等しい整数配列である. 置き換えられ, その後 CRT か実行される. CRT 後 stat[i]=1 となる. stat[i]=0,1 のとき, 単に CRT が行される. stat[i]=2 のとき, g[i] mod p をモニック化したものと gp[i] が比較され, 等しくない場合に, g[i] は, 有理数変換前の g mod M に 戻され, その後CRTが実行される.
• 外部関数において, 整数-有理数変換が成功した場合にstat[i]=2とすることで, dp_chrem2 を実行したあとに, 不適切な有理数への変換を検出することができる. 返される値は, stat[i]=2 の場合に不適切な有理数変換を検出した個数である.

[2465] G=ltov([<<1,2>>+3*<<0,1>>,<<3,4>>+2*<<1,0>>])$
[2466] GM=ltov([<<1,2>>+5*<<0,1>>,<<3,4>>+7*<<1,0>>])$
[2467] Stat=vector(2)$
[2468] Mod=17$
[2469] P=19$
[2470] noro_grcrt.dp_chrem2(G,Mod,GM,P,Stat);
0
[2471] G;
[ (1)*<<1,2>>+(309)*<<0,1>> (1)*<<3,4>>+(121)*<<1,0>> ] 
[2472] Stat;
[ 1 1 ]

1.2.2 noro_grcrt.intdptoratdp

noro_grcrt.intdptoratdp(f,M,B)

整数-有理数変換を行う.

return

多項式

f

分散多項式

M

B

正整数

• fを法Mで整数-有理数変換を行った結果を返す. 一つでも変換できなかった係数がある場合には 0を返す.
• B は M/2の平方根を超えない最大の整数で, noro_grcrt.calcb で計算できる. 複数の多項式の係数を同一の M, Bで変換するので, この値をあらかじめ与えるようにしてある.

[2495] P=62884891$
[2496] A=<<1,2>>+(25632978)*<<0,1>>$
[2497] B=noro_grcrt.calcb(P);
5607
[2498] noro_grcrt.intdptoratdp(A,P,B);
(1)*<<1,2>>+(-5321/1234)*<<0,1>>

Index

目次

簡略化した目次

この文書について

この文書は4月 19, 2024にtexi2html 5.0を用いて生成されました。

ナビゲーションパネル中のボタンには以下の意味があります。

例では、以下に示す構造を持つ文書の1.2.3項を現在位置に仮定しています。

• 1. 第1項
  • 1.1 第1.1項
    • ...
  • 1.2 第1.2項
    • 1.2.1 第1.2.1項
    • 1.2.2 第1.2.2項
    • 1.2.3 第1.2.3項     <== 現在位置
    • 1.2.4 第1.2.4項
  • 1.3 第1.3項
    • ...
  • 1.4 第1.4項

この文書は4月 19, 2024にtexi2html 5.0を用いて生成されました。