252   PetscReal          ai, bi, ci, di, ei;  in MatDFischer()  local
294         ai = fischnorm(1.0, PetscRealPart(t2[i]));  in MatDFischer()
296         da[i] = -1.0 / ai - 1.0;  in MatDFischer()
297         db[i] = -t2[i] / ai - 1.0;  in MatDFischer()
300         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
301         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDFischer()
303         da[i] = bi / ai - 1.0;  in MatDFischer()
304         db[i] = -f[i] / ai - 1.0;  in MatDFischer()
308         ai = fischnorm(1.0, PetscRealPart(t2[i]));  in MatDFischer()
310         da[i] = 1.0 / ai - 1.0;  in MatDFischer()
311         db[i] = t2[i] / ai - 1.0;  in MatDFischer()
314         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
315         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDFischer()
317         da[i] = bi / ai - 1.0;  in MatDFischer()
318         db[i] = f[i] / ai - 1.0;  in MatDFischer()
325         ai = fischnorm(1.0, PetscRealPart(t2[i]));  in MatDFischer()
327         ci = 1.0 / ai + 1.0;  in MatDFischer()
328         di = PetscRealPart(t2[i]) / ai + 1.0;  in MatDFischer()
331         ai = fischnorm(bi, PetscRealPart(f[i]));  in MatDFischer()
332         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDFischer()
334         ci = bi / ai + 1.0;  in MatDFischer()
335         di = PetscRealPart(f[i]) / ai + 1.0;  in MatDFischer()
340         ai = fischnorm(1.0, bi);  in MatDFischer()
342         bi = bi / ai - 1.0;  in MatDFischer()
343         ai = 1.0 / ai - 1.0;  in MatDFischer()
346         ai = fischnorm(PetscRealPart(x[i]) - PetscRealPart(l[i]), ei);  in MatDFischer()
347         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDFischer()
349         bi = ei / ai - 1.0;  in MatDFischer()
350         ai = (PetscRealPart(x[i]) - PetscRealPart(l[i])) / ai - 1.0;  in MatDFischer()
353       da[i] = ai + bi * ci;  in MatDFischer()
398   PetscReal          ai, bi, ci, di, ei, fi;  in MatDSFischer()  local
421         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
422         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDSFischer()
424         da[i] = bi / ai - 1.0;  in MatDSFischer()
425         db[i] = -PetscRealPart(f[i]) / ai - 1.0;  in MatDSFischer()
426         dm[i] = 2.0 * mu / ai;  in MatDSFischer()
429         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
430         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDSFischer()
432         da[i] = bi / ai - 1.0;  in MatDSFischer()
433         db[i] = PetscRealPart(f[i]) / ai - 1.0;  in MatDSFischer()
434         dm[i] = 2.0 * mu / ai;  in MatDSFischer()
441         ai = fischsnorm(bi, PetscRealPart(f[i]), mu);  in MatDSFischer()
442         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDSFischer()
444         ci = bi / ai + 1.0;  in MatDSFischer()
445         di = PetscRealPart(f[i]) / ai + 1.0;  in MatDSFischer()
446         fi = 2.0 * mu / ai;  in MatDSFischer()
449         ai = fischsnorm(PetscRealPart(x[i]) - PetscRealPart(l[i]), ei, mu);  in MatDSFischer()
450         ai = PetscMax(PETSC_MACHINE_EPSILON, ai);  in MatDSFischer()
452         bi = ei / ai - 1.0;  in MatDSFischer()
454         ai = (PetscRealPart(x[i]) - PetscRealPart(l[i])) / ai - 1.0;  in MatDSFischer()
456         da[i] = ai + bi * ci;  in MatDSFischer()