13   const MatScalar   *aa   = a->a, *v;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()  local
28     v  = aa + 9 * ai[i];  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
40       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
41       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
42       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
43       v += 9;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
51     v   = aa + 9 * diag[i] + 9;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
63       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
64       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
65       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
66       v += 9;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
68     v          = aa + 9 * diag[i];  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
69     x[idt]     = v[0] * s1 + v[3] * s2 + v[6] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
70     x[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
71     x[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace()
86   const MatScalar   *aa = a->a, *v;  in MatSolve_SeqBAIJ_3_NaturalOrdering()  local
100     v   = aa + bs2 * ai[i];  in MatSolve_SeqBAIJ_3_NaturalOrdering()
112       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
113       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
114       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
116       v += bs2;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
126     v   = aa + bs2 * (adiag[i + 1] + 1);  in MatSolve_SeqBAIJ_3_NaturalOrdering()
139       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
140       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
141       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
143       v += bs2;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
146     x[idt]     = v[0] * s1 + v[3] * s2 + v[6] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
147     x[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
148     x[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3;  in MatSolve_SeqBAIJ_3_NaturalOrdering()
163   const MatScalar   *aa = a->a, *v;  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()  local
177     v   = aa + bs2 * ai[i];  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
189       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
190       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
191       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
193       v += bs2;  in MatForwardSolve_SeqBAIJ_3_NaturalOrdering()
213   const MatScalar   *aa = a->a, *v;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()  local
224     v   = aa + bs2 * (adiag[i + 1] + 1);  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
237       s1 -= v[0] * x1 + v[3] * x2 + v[6] * x3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
238       s2 -= v[1] * x1 + v[4] * x2 + v[7] * x3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
239       s3 -= v[2] * x1 + v[5] * x2 + v[8] * x3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
241       v += bs2;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
244     x[idt]     = v[0] * s1 + v[3] * s2 + v[6] * s3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
245     x[1 + idt] = v[1] * s1 + v[4] * s2 + v[7] * s3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()
246     x[2 + idt] = v[2] * s1 + v[5] * s2 + v[8] * s3;  in MatBackwardSolve_SeqBAIJ_3_NaturalOrdering()