blob: 3b62d917f632dcbe6e647e21912fb5863f591dd0 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
void SubMtx_blockDiagonalInfo (int **);
double *ZV_entries (void);
void SubMtx_fillRowZV (int irow)
{
double *rowvec = ZV_entries ();
double *entries;
int ii, ipivot, jrow, kk, m;
int *pivotsizes;
SubMtx_blockDiagonalInfo (&pivotsizes);
for (jrow = ipivot = kk = 0; jrow <= irow; ipivot++)
{
m = pivotsizes[ipivot];
if (jrow <= irow && irow < jrow + m)
for (ii = jrow; ii < irow; ii++)
{
rowvec[2*ii] = entries[2*kk];
rowvec[2*ii+1] = entries[2*kk+1];
}
jrow += m;
}
}
|