252   PetscReal          ai, bi, ci, di, ei;  in MatDFischer()  local
299         bi = PetscRealPart(u[i]) - PetscRealPart(x[i]);  in MatDFischer()
300         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
303         da[i] = bi / ai - 1.0;  in MatDFischer()
313         bi = PetscRealPart(x[i]) - PetscRealPart(l[i]);  in MatDFischer()
314         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
317         da[i] = bi / ai - 1.0;  in MatDFischer()
330         bi = PetscRealPart(x[i]) - PetscRealPart(u[i]);  in MatDFischer()
331         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
334         ci = bi / ai + 1.0;  in MatDFischer()
339         bi = ci + di * PetscRealPart(t2[i]);  in MatDFischer()
340         ai = fischnorm(1.0, bi);  in MatDFischer()
342         bi = bi / ai - 1.0;  in MatDFischer()
349         bi = ei / ai - 1.0;  in MatDFischer()
353       da[i] = ai + bi * ci;  in MatDFischer()
354       db[i] = bi * di;  in MatDFischer()
398   PetscReal          ai, bi, ci, di, ei, fi;  in MatDSFischer()  local
420         bi = PetscRealPart(u[i]) - PetscRealPart(x[i]);  in MatDSFischer()
421         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
424         da[i] = bi / ai - 1.0;  in MatDSFischer()
428         bi = PetscRealPart(x[i]) - PetscRealPart(l[i]);  in MatDSFischer()
429         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
432         da[i] = bi / ai - 1.0;  in MatDSFischer()
440         bi = PetscRealPart(x[i]) - PetscRealPart(u[i]);  in MatDSFischer()
441         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
444         ci = bi / ai + 1.0;  in MatDSFischer()
452         bi = ei / ai - 1.0;  in MatDSFischer()
456         da[i] = ai + bi * ci;  in MatDSFischer()
457         db[i] = bi * di;  in MatDSFischer()
458         dm[i] = ei + bi * fi;  in MatDSFischer()