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