impls/kaczmarz/kaczmarz.c

20   const PetscScalar *vals, *xarray;
34     PetscCall(MatGetRow(pc->pmat, i, &ncols, &cols, &vals));
38       if (cols[j] >= ys && cols[j] < ye) r -= yarray[cols[j] - ys] * vals[j];
39       anrm += PetscRealPart(PetscSqr(vals[j]));
43         if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm;
46     PetscCall(MatRestoreRow(pc->pmat, i, &ncols, &cols, &vals));
50       PetscCall(MatGetRow(pc->pmat, i, &ncols, &cols, &vals));
54         if (cols[j] >= ys && cols[j] < ye) r -= yarray[cols[j] - ys] * vals[j];
55         anrm += PetscRealPart(PetscSqr(vals[j]));
59           if (cols[j] >= ys && cols[j] < ye) yarray[cols[j] - ys] += vals[j] * lambda * r / anrm;
62       PetscCall(MatRestoreRow(pc->pmat, i, &ncols, &cols, &vals));