blob: 0f8a1f32f6e91f2e699eacb3d2c93f6f624d3054 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
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;
}
}
|