|
|
74.
|
|
|
Now we will see how to use the lu factors:
|
|
|
|
|
ここで,LU分解を使用する方法を示します:
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
75.
|
|
|
1/ lu factors of a sparse matrix A are obtained through:
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
1/ 疎行列AのLU分解は次のように得られます:
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
76.
|
|
|
lup = umf_lufact(A) 
lup is a pointer to the lu factors (the memory is outside scilab space)
|
|
|
represents a line break.
Start a new line in the equivalent position in the translation.
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
lup = umf_lufact(A)
lup は LU分解へのポインタです (そのメモリはScilab空間の外です)
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
77.
|
|
|
2/ for solving a linear system A*x = b, use:
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
2/ 線形システム A*x = b を解くには, 次のようにします:
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
78.
|
|
|
x = umf_lusolve(lup,b)
or x = umf_lusolve(lup,b,Ax=b,A) to add an iterative refinement step
|
|
|
represents a line break.
Start a new line in the equivalent position in the translation.
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
x = umf_lusolve(lup,b)
または x = umf_lusolve(lup,b,Ax=b,A) (反復改良ステップを追加する場合)
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
79.
|
|
|
3/ to solve A'*x=b you may use:
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
3/ A'*x=b を解く場合, 次のように行います:
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
80.
|
|
|
x = umf_lusolve(lup,b,Ax''=b)
or x = umf_lusolve(lup,b,Ax''=b,A) to add an iterative refinement step
|
|
|
represents a line break.
Start a new line in the equivalent position in the translation.
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
x = umf_lusolve(lup,b,Ax''=b)
または x = umf_lusolve(lup,b,Ax''=b,A) (反復改良ステップを追加する場合)
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
81.
|
|
|
4/ you may also compute the 1-norm condition number quickly with:
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
(no translation yet)
|
|
|
|
|
82.
|
|
|
K1 = condestsp(A,lup)
K1 = condestsp(A) also works but in this case the lu factors are re-computed inside
|
|
|
represents a line break.
Start a new line in the equivalent position in the translation.
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
K1 = condestsp(A,lup)
K1 = condestsp(A) も動作しますが, この場合 LU分解が内部的に再計算されます
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
|
|
83.
|
|
|
5/ if you don't need the lu factors anymore, it is recommended to free them with:
|
|
|
represents a space character.
Enter a space in the equivalent position in the translation.
|
|
|
|
|
5/ LU分解をもう必要としない場合, 次のように解放することを推奨します:
|
|
Translated and reviewed by
Rui Hirokawa
|
|
|
