| Variable | Meaning | Value |
|---|---|---|
| ia(_JT_ nunk_) | number of unknowns | 2 |
| ia(_JT_ nelem_) | number of active elements | 3 |
RSS format is essentially a variation of COO format.
In this format the positions in the full matrix are numbered from 1 to ia(_JT_ nunk_) * ia(_JT_ nunk_) , starting with element (1,1), and proceeding across each row in succession to element (ia(_JT_ nunk_) ,ia(_JT_ nunk_) ). For example, the elements of a 3 x 3 matrix would be numbered as follows:
So again a is a vector of length ia(_JT_ nelem_)
containing the
active elements of A. The map array ja is again a
vector of length ia(_JT_ nelem_)
containing the
sequence number
of the corresponding element of a, that is,
.
The elements of a can be in any order, except that
the main diagonal elements must be in the first
ia(_JT_ nunk_)
elements.
Advantages:
Disadvantages:
Example:
could be represented in row sequential format by the
arrays a and ja as:
a = ( 11, 22, 33, 44, 55, 41, 51, 14, 54, 15, 45 )
ja = ( 1, 7, 13, 19, 25, 16, 21, 4, 24, 5, 20 )