4.3 _JT_ storage_ EISPACK_ - EISPACK format

• Status: IMPLEMENTED
• Current value: 2
• 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 column of the reduced matrix corresponds to a diagonal of the original matrix. This storage scheme is also used in EISPACK. 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 EISPACK format as the following ia(_JT_ nrows_) x (ia(_JT_ lbw_) + ia(_JT_ lbw_) + 1) array:

  a = $$\left(\vphantom{ \begin{array}{ccccc} 0 & 0 & 0 & 11 & 12 \\ ... ...45 \\ 52 & 0 & 54 & 55 & 56 \\ 63 & 0 & 65 & 66 & 0 \end{array} }\right.$$$$\begin{array}{ccccc} 0 & 0 & 0 & 11 & 12 \\ 0 & 0 & 21 & 22 &... ...3 & 44 & 45 \\ 52 & 0 & 54 & 55 & 56 \\ 63 & 0 & 65 & 66 & 0 \end{array}$$ $$\left.\vphantom{ \begin{array}{ccccc} 0 & 0 & 0 & 11 & 12 \\ ... ...45 \\ 52 & 0 & 54 & 55 & 56 \\ 63 & 0 & 65 & 66 & 0 \end{array} }\right)$$