31     PetscScalar      s1, s2, s3, s4, x1, x2, x3, x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()  local
47       s1 = b[idx];  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
57         s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
63       x[idx]     = s1;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
75       s1   = x[idt];  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
85         s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
92       x[idt]     = v[0] * s1 + v[4] * s2 + v[8] * s3 + v[12] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
93       x[1 + idt] = v[1] * s1 + v[5] * s2 + v[9] * s3 + v[13] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
94       x[2 + idt] = v[2] * s1 + v[6] * s2 + v[10] * s3 + v[14] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
95       x[3 + idt] = v[3] * s1 + v[7] * s2 + v[11] * s3 + v[15] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace()
116   PetscScalar        s1, s2, s3, s4, x1, x2, x3, x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()  local
132     s1  = b[idx];  in MatSolve_SeqBAIJ_4_NaturalOrdering()
142       s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
150     x[idx]     = s1;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
162     s1  = x[idt];  in MatSolve_SeqBAIJ_4_NaturalOrdering()
173       s1 -= v[0] * x1 + v[4] * x2 + v[8] * x3 + v[12] * x4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
181     x[idt]     = v[0] * s1 + v[4] * s2 + v[8] * s3 + v[12] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
182     x[1 + idt] = v[1] * s1 + v[5] * s2 + v[9] * s3 + v[13] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
183     x[2 + idt] = v[2] * s1 + v[6] * s2 + v[10] * s3 + v[14] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()
184     x[3 + idt] = v[3] * s1 + v[7] * s2 + v[11] * s3 + v[15] * s4;  in MatSolve_SeqBAIJ_4_NaturalOrdering()