3.1 詰め込み方

$ 1\le j\le$   N は当たり前として

   A$\displaystyle (i,j)
=\left\{
\begin{array}{ll}
\mbox{\texttt{AB}}(K_L+K_U+1+...
...min(N,j+K_L)$ のとき)}\\
0 & \mbox{(それ以外)}.
\end{array} \right.
$

(これは LAPACK のドキュメントにも書いてある。)

この逆に ABA で表す公式も 作っておく。

$\displaystyle I:=\{
(p,j);
1\le j\le N, \max(K_L+1,K_L+K_U+2-j)
\le p \le \min(2K_L+K_U+1,N+K_L+K_U+1-j)
\},
$

$\displaystyle J:=\{(p,j); 1\le p\le 2K_L+K_U+1, 1\le j\le N\}\setminus I
$

とするとき、

   AB$\displaystyle (p,j)=
\left\{
\begin{array}{ll}
\mbox{\texttt{A}}(p+j-(K_L+K_...
... のとき)}\\
\mbox{(ゴミ)} & \mbox{($(p,j)\in J$)}.
\end{array} \right.
$

桂田 祐史
2017-10-07