12   PetscInt    *i, *j; /* offsets of nonzeros by column, non-zero indices by column */  member
18   PetscInt    i, j, *colcnt;  in PCSetUp_CP()  local
32     PetscCall(PetscFree3(cp->a, cp->i, cp->j));  in PCSetUp_CP()
37   if (!cp->a) PetscCall(PetscMalloc3(aij->nz, &cp->a, cp->n + 1, &cp->i, aij->nz, &cp->j));  in PCSetUp_CP()
40   for (i = 0; i < aij->nz; i++) colcnt[aij->j[i]]++;  in PCSetUp_CP()
45     for (j = aij->i[i]; j < aij->i[i + 1]; j++) { /* over columns in row */  in PCSetUp_CP()
46       cp->j[cp->i[aij->j[j]] + colcnt[aij->j[j]]]   = i;  in PCSetUp_CP()
47       cp->a[cp->i[aij->j[j]] + colcnt[aij->j[j]]++] = aij->a[j];  in PCSetUp_CP()
55 …for (j = cp->i[i]; j < cp->i[i + 1]; j++) cp->d[i] += cp->a[j] * cp->a[j]; /* over rows in column …  in PCSetUp_CP()
65   PetscInt     i, j;  in PCApply_CP()  local
74 …  for (j = cp->i[i]; j < cp->i[i + 1]; j++) xt += cp->a[j] * b[cp->j[j]]; /* over rows in column */  in PCApply_CP()
77 …for (j = cp->i[i]; j < cp->i[i + 1]; j++) b[cp->j[j]] -= xt * cp->a[j]; /* over rows in column upd…  in PCApply_CP()
81 …  for (j = cp->i[i]; j < cp->i[i + 1]; j++) xt += cp->a[j] * b[cp->j[j]]; /* over rows in column */  in PCApply_CP()
84 …for (j = cp->i[i]; j < cp->i[i + 1]; j++) b[cp->j[j]] -= xt * cp->a[j]; /* over rows in column upd…  in PCApply_CP()
99   PetscCall(PetscFree3(cp->a, cp->i, cp->j));  in PCReset_CP()
110   PetscCall(PetscFree3(cp->a, cp->i, cp->j));  in PCDestroy_CP()