41   PetscInt       *ai,*aj,*r;  in MatReorderingSeqSBAIJ()
54     PetscCall(PetscMalloc2(mbs+1,&ai, 2*a->i[mbs],&aj));  in MatReorderingSeqSBAIJ()
56     ai = a->inew; aj = a->jnew;  in MatReorderingSeqSBAIJ()
58   PetscCall(PetscArraycpy(ai,a->i,mbs+1));  in MatReorderingSeqSBAIJ()
69   PetscCall(PetscMalloc1(ai[mbs],&r));  in MatReorderingSeqSBAIJ()
71   for (i=0; i<ai[mbs]; i++) r[i] = 0;  in MatReorderingSeqSBAIJ()
75     nz = ai[i+1] - ai[i];  in MatReorderingSeqSBAIJ()
76     j  = ai[i];  in MatReorderingSeqSBAIJ()
98     ai[i+1] = ai[i] + nzr[i];  in MatReorderingSeqSBAIJ()
99     nzr[i]  = ai[i+1];  in MatReorderingSeqSBAIJ()
104   jmin = ai[0]; jmax = ai[mbs];  in MatReorderingSeqSBAIJ()
109     if (aj[j] == i) r[j] = ai[i]; /* put diagonal nonzero at beginning of row */  in MatReorderingSeqSBAIJ()
116   a->a2anew = aj + ai[mbs];  in MatReorderingSeqSBAIJ()
117   PetscCall(PetscArraycpy(a->a2anew,r,ai[mbs]));  in MatReorderingSeqSBAIJ()
129   a->inew = ai;  in MatReorderingSeqSBAIJ()