KMSL95

(Kamitani Mathematical Subroutine Libraries for fortran 95)

補間

1次元スプライン関数　ccspd
2次元スプライン関数　c2csp2
3次元スプライン関数　c3csp
ベクトル場補間用の3次元スプライン関数　c3csv

有限要素法

2次元有限要素法で得られる対称行列の非零要素位置　k_posit_nonzero_2dfem
ICCG法形式を用いて格納された対称行列情報中の非零要素 \(a_{ij}\) の非零要素番号の決定　k_add_nz2_iccg
ICCG法形式を用いて格納された対称行列情報中の非零要素 \(a_{ij}\) の格納位置の追加　k_mod_val_iccg

メッシュレス法

2次元MSL近似（高速）　mlsap2d_95
3次元MSL近似（高速）　mlsap3d_95
CRS法を用いて格納された非正方行列情報中の非零要素 \(a_{ij}\) の非零要素番号の決定　k_add_nz2_crs
CCS法を用いて格納された非正方行列情報中の非零要素 \(a_{ij}\) の非零要素番号の決定　k_add_nz2_ccs
CRS/CCS法を用いて格納された非正方行列情報中の非零要素 \(a_{ij}\) の格納位置の追加　k_mod_val
EFG法で得られる対称行列の非零要素位置　k_posd_nonzero
ICCG法形式とCRS法で圧縮格納されたEFG型連立１次方程式の係数行列をICCG法形式に変換　k_conv_iccg_crs_2_iccg
3つの部分行列がCRS法で圧縮格納されたX-EFG型連立１次方程式の係数行列全体をCRS法で圧縮格納　k_conv_3crs_2_crs

数値積分

Gauss-Legendre積分(高精度版）　dawg2
Gauss積分（重み関数 \( W(x) = \log \frac{1}{x}\)）　dawlg
適応的自動積分　k_adpint

特殊関数＆関数近似

第1種完全楕円積分　ecei1
第2種完全楕円積分　ecei2
第1種完全楕円積分の多項式近似　pcei1
第2種完全楕円積分の多項式近似　pcei2
Legendre多項式　legen
複素調和多重対数関数　hpl

配列操作

整数型1次元配列の要素をistep個ずつデータファイルに出力　k_1d_int_array_fout
倍精度実数型1次元配列の要素をistep個ずつデータファイルに出力　k_1d_real_array_fout
k_1d_int_array_fout によりデータファイルに書き出された整数型1次元配列の読み込み　k_1d_int_array_finp
k_1d_real_array_fout によりデータファイルに書き出された倍精度実数型1次元配列の読み込み　k_1d_real_array_finp

線形計算と行列操作（圧縮・復元）

ICCG法形式による対称行列の圧縮格納　k_conv_2_iccg
ICCG法形式で圧縮格納された対称行列の復元　k_convm_from_iccg
ICCG法形式に対して，iu(k)=j の解 k 　k__iu_eq_j
ICCG法形式で圧縮格納された対称行列の下三角情報の抽出　k_lower_from_iccg
ICCG法形式で圧縮格納された対称行列とベクトルの積 \(A\boldsymbol{x}\)　prod1d
【最新バージョン】ICCG法形式で圧縮格納された対称行列 \(A\) からの零要素情報の削除　k_r0_iccg_n
ICCG法形式で圧縮格納された対称行列の不完全 Cholesky 分解　k_ic_n
ICCG法形式で圧縮格納された対称行列の不完全 Cholesky 分解結果を用いた行列ベクトル積 \((\hat{L}\hat{D}\hat{L}^T)^{-1}\boldsymbol{x}\)　icsl1d
CRS法による非正方行列の圧縮格納　k_convm_2_crs
CRS法で圧縮格納された非正方行列の復元　k_convm_from_crs
CRS法で圧縮格納された非正方行列のrow_indexと非零行ベクトルの行番号を決定　k_gen_row_ind_crs
CRS法で圧縮格納された非正方行列 \(A\) とベクトル \(\boldsymbol{x}\) との積 \(A\boldsymbol{x}\) 　k_product_ax_crs
CRS法で圧縮格納された非正方行列 \(A\) とベクトル \(\boldsymbol{x}\) との積 \(A\boldsymbol{x}\) を高速計算　k_f_product_ax_crs
CRS法で圧縮格納された非正方行列の転置行列 \(A^T\) とベクトル \(\boldsymbol{x}\) との積 \(A^T\boldsymbol{x}\)　k_product_atx_crs
CRS形式で圧縮格納された非正方行列 \(A\) からの零要素情報の削除　k_r0_crs
CRS法で圧縮格納された非正方行列 \(A\) とベクトル \(\boldsymbol{x}\) から \(A^T A \boldsymbol{x}\) を高速計算　k_product_atax_crs
CRS法で圧縮格納された非正方行列 \(C,D\) とベクトル \(\boldsymbol{x}\) から \(D^T C \boldsymbol{x}\) を高速計算　k_product_dtcx_crs
CRS法で圧縮格納された正方行列の不完全LU分解　k_decmp_ilu_crs
CRS法で圧縮格納された正方行列の不完全LDU分解（格納方式は不完全LU分解と同じ）　k_decmp_ilu_crs_d0
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(L\boldsymbol{x}\)　k_matvec_ilu_crs_el0x
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(L^T \boldsymbol{x}\) 　k_matvec_ilu_crs_e0tx
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(L^{-1} \boldsymbol{x}\)　k_matvec_ilu_crs_el1x
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(L^{-T} \boldsymbol{x}\)　k_matvec_ilu_crs_e1tx
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(U\boldsymbol{x}\)　k_matvec_ilu_crs_uh0x
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(U^T \boldsymbol{x}\)　k_matvec_ilu_crs_u0tx
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(U^{-1} \boldsymbol{x}\)　k_matvec_ilu_crs_uh1x
CRS法で圧縮格納された正方行列の不完全LU分解結果を用いた行列ベクトル積 \(U^{-T} \boldsymbol{x}\)　k_matvec_ilu_crs_u1tx
MRS形式で圧縮格納された正方行列のCRS形式への変換　k_convm_mrs_2_crs
CRS形式で圧縮格納された非正方行列のCCS形式への変換　k_convm_crs_2_ccs
CRS法を用いて格納された非正方行列情報から零要素を削除　k_del_0from_crs
CCS法による非正方行列の圧縮格納　k_convm_2_ccs
CCS法で圧縮格納された非正方行列の復元　k_convm_from_ccs
CCS法で圧縮格納された非正方行列の情報から特定要素の情報を削除　k_rem_el_from_ccs
CCS法で圧縮格納された非正方行列 \(A\) とベクトル \(\boldsymbol{x}\) との積 \(A\boldsymbol{x}\) k_product_ax_ccs
CCS法で圧縮格納された非正方行列の転置行列 \(A^T\) とベクトル \(\boldsymbol{x}\) との積 \(A^T\boldsymbol{x}\)　k_product_atx_ccs
対称行列の1-または\(\infty\)-ノルム：引数改良版　nsmat2
実対称行列の固有値と固有ベクトル：引数改良版　meigen2

連立1次方程式（直接法と定常反復法）

連立1次方程式の解法（LU分解：ludec_p, sollu_pに基づく）　gauss_p
連立1次方程式の解法（LU分解：ludec_n, sollu_nに基づく）　gauss_n
行列のLU分解：行交換実施　ludec_p
行列のLU分解：インデックスベクトル使用　ludec_n
LU分解された係数行列をもつ連立1次方程式の解法：行交換実施　sollu_p
LU分解された係数行列をもつ連立1次方程式の解法：インデックスベクトル使用　sollu_n
連立1次方程式のSOR法による解法：引数改良版　sormt2
3重対角行列を係数行列とする連立1次方程式の解法　sotrd

連立1次方程式（Krylov 部分空間法）

対称係数行列をもつ連立1次方程式のICCG法による解法：引数改良版　iccg1d2
【最新バージョン】対称係数行列をもつ連立1次方程式のCG法による解法　k_cg_hs_n
【最新バージョン】対称係数行列をもつ連立1次方程式の前処理付きCG法による解法　k_pcg_hs_n
【最新バージョン】非対称係数行列をもつ連立１次方程式のGMRES(\(k\))法による解法　k_gmres_n
【最新バージョン】非対称係数行列をもつ連立１次方程式のGCR(\(k\))法による解法　k_gcr_n
【最新バージョン】非対称係数行列をもつ連立１次方程式のOrthomin(\(k\))法による解法　k_orthomin_n
非対称係数行列をもつ連立１次方程式のBiCGSTAB法による解法　k_bicgstab_n
非対称係数行列をもつ連立１次方程式のBiCGSTAB2法による解法　k_bicgstab2_n
【最新バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きGMRES(\(k\))法による解法　k_pgmres_n
【最新バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きGCR(\(k\))法による解法　k_pgcr_n
【最新バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きOrthomin(\(k\))法による解法　k_porthomin_n
【最新バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きBiCGSTAB法による解法　k_pbicgstab_n
非対称係数行列をもつ連立１次方程式の（左）前処理付きGMRES(\(k\))法による解法　lpgmres

非線形方程式

1変数超越方程式の解法（2分法）　zep1b
1変数超越方程式の解法（割線法）　zep1s
1変数超越方程式の解法（Newton法）　zep1n
2次元ベクトル関数の超越方程式の解法（Newton法）　zep2n
1変数超越方程式の解を含む区間の探索（2分法）　ssubi
1変数関数の最小値の探索　sminp

\(\mathcal{H}\)-行列法

実対称密行列の\(\mathcal{H}\)-行列への変換　h_matrix
\(\mathcal{H}\)-行列（近似行列）とベクトルの積　product_hv
\(\mathcal{H}\)-行列とベクトルの近似積　product_hv_app

文字処理

文字列中の最初の半角スペースまでの文字数のカウント ssafw
サブルーチンの実行の終了を知らせる sucmp
エラーメッセージの出力 ermsg

電磁界解析

円柱型永久磁石が生成する磁束密度mflxd_cylind
直方体永久磁石が生成する磁束密度mflxd_rectp
矩形断面コイルが生成する磁束密度 mflxd_coil_rect

その他

machine epsilonの計算　macep
CPU時間の測定　mtime
サブルーチンrandom_numberに再現性をもたせるように設定　k_set_seed_random_rp
Runge-Kutta法による連立１階常微分方程式の解法 k_runge

旧バージョン

2次元スプライン関数　c2csp
【旧バージョン】2次元MSL近似（低速）　mlsap2d
Gauss-Legendre積分　dawgl
CRS法を用いて格納された非正方行列情報中の非零要素 \(a_{ij}\) の格納位置の追加　k_mod_val_crs
連立1次方程式の解法（Gaussの消去法）　gauss
行列のLU分解　ludec
LU分解された係数行列をもつ連立1次方程式の解法　sollu
連立1次方程式の解法（Gaussの消去法：引数改良版）　gauss2
行列のLU分解：引数改良版　ludec2
LU分解された係数行列をもつ連立1次方程式の解法：引数改良版　sollu2
連立1次方程式のSOR法による解法　sormt
対称行列の1-または\(\infty\)-ノルム　nsmat
実対称行列の固有値と固有ベクトル　meigen
対称係数行列をもつ連立1次方程式のICCG法による解法　iccg1d
【旧バージョン】ICCG法形式で圧縮格納された対称行列 \(A\) からの零要素情報の削除　k_r0_iccg
【旧バージョン】対称係数行列をもつ連立1次方程式のCG法による解法　k_cg_hs
【旧バージョン】非対称係数行列をもつ連立１次方程式のGMRES(\(k\))法による解法　gmres
【旧バージョン】非対称係数行列をもつ連立１次方程式のGCR(\(k\))法による解法　k_gcr
【旧バージョン】非対称係数行列をもつ連立１次方程式のOrthomin(\(k\))法による解法　k_orthomin
【旧バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きGMRES(\(k\))法による解法　pgmres
【旧バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きGCR(\(k\))法による解法　k_pgcr
【旧バージョン】非対称係数行列をもつ連立１次方程式の（右）前処理付きOrthomin(\(k\))法による解法　k_porthomin