1766   PetscReal          diagk;  in MatForwardSolve_SeqSBAIJ_1()  local
1785     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1()
1786 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatForwardSolve_SeqSBAIJ_1()
1787     x[k] = xk * PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1()
1802   PetscReal          diagk;  in MatForwardSolve_SeqSBAIJ_1_inplace()  local
1821     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_inplace()
1822 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatForwardSolve_SeqSBAIJ_1_inplace()
1823     x[k] = xk * PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1_inplace()
1838   PetscReal          diagk;  in MatBackwardSolve_SeqSBAIJ_1()  local
1853     diagk = PetscRealPart(aa[adiag[k]]);  in MatBackwardSolve_SeqSBAIJ_1()
1854 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatBackwardSolve_SeqSBAIJ_1()
1855     t[k] = b[k] * PetscSqrtReal(diagk);  in MatBackwardSolve_SeqSBAIJ_1()
1873   PetscReal          diagk;  in MatBackwardSolve_SeqSBAIJ_1_inplace()  local
1888     diagk = PetscRealPart(aa[ai[k]]);  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1889 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatBackwardSolve_SeqSBAIJ_1_inplace()
1890     t[k] = b[k] * PetscSqrtReal(diagk);  in MatBackwardSolve_SeqSBAIJ_1_inplace()
2163   PetscReal          diagk;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2178     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[adiag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2179 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2180     x[k] *= PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering()
2193   PetscReal          diagk;  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2208     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2209 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2210     x[k] *= PetscSqrtReal(diagk);  in MatForwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2223   PetscReal          diagk;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()  local
2236     diagk = PetscRealPart(aa[adiag[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2237 …PetscCheck(!PetscImaginaryPart(aa[adiag[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diago…  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2238     x[k] = PetscSqrtReal(diagk) * b[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering()
2253   PetscReal          diagk;  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()  local
2266     diagk = PetscRealPart(aa[ai[k]]); /* note: aa[diag[k]] = 1/D(k) */  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2267 …PetscCheck(!PetscImaginaryPart(aa[ai[k]]) && diagk >= 0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Diagonal…  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()
2268     x[k] = PetscSqrtReal(diagk) * b[k];  in MatBackwardSolve_SeqSBAIJ_1_NaturalOrdering_inplace()