4.4 _JT_ storage_ LAPACK_ - LAPACK format

• Status: NOT IMPLEMENTED
• Current value: 3
• Also requires:

  Variable	Meaning	Value
  ia(_JT_ idim_)	leading dimension	3
  ia(_JT_ nrows_)	number of rows	2
  ia(_JT_ lbw_)	lower bandwidth	5
  ia(_JT_ ubw_)	upper bandwidth	6

• Description:

  In this option, each row of the reduced matrix corresponds to a diagonal of the original matrix. This storage scheme is also used in LAPACK.

  Example:

  A = $$\left(\vphantom{ \begin{array}{cccccc} 11 & 12 & 0 & 0 & 0 & 0 \... ... 0 & 52 & 0 & 54 & 55 & 56 \\ 0 & 0 & 63 & 0 & 65 & 66 \end{array} }\right.$$$$\begin{array}{cccccc} 11 & 12 & 0 & 0 & 0 & 0 \\ 21 & 22 & 23... ...& 0 \\ 0 & 52 & 0 & 54 & 55 & 56 \\ 0 & 0 & 63 & 0 & 65 & 66 \end{array}$$ $$\left.\vphantom{ \begin{array}{cccccc} 11 & 12 & 0 & 0 & 0 & 0 \... ... 0 & 52 & 0 & 54 & 55 & 56 \\ 0 & 0 & 63 & 0 & 65 & 66 \end{array} }\right)$$

  would be represented in full-band row storage as the following (ia(_JT_ lbw_) + ia(_JT_ ubw_) + 1) x ia(_JT_ nrows_) array:

  a = $$\left(\vphantom{ \begin{array}{cccccc} 0 & 0 & 0 & 41 & 52 & 63 ... ...& 22 & 33 & 44 & 55 & 66 \\ 12 & 23 & 34 & 45 & 56 & 0 \end{array} }\right.$$$$\begin{array}{cccccc} 0 & 0 & 0 & 41 & 52 & 63 \\ 0 & 0 & 0 &... ...\\ 11 & 22 & 33 & 44 & 55 & 66 \\ 12 & 23 & 34 & 45 & 56 & 0 \end{array}$$ $$\left.\vphantom{ \begin{array}{cccccc} 0 & 0 & 0 & 41 & 52 & 63 ... ...& 22 & 33 & 44 & 55 & 66 \\ 12 & 23 & 34 & 45 & 56 & 0 \end{array} }\right)$$