| /petsc/src/mat/impls/sbaij/seq/ ! |
| H A D | sbaijfact12.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering() local
28 PetscCall(PetscMalloc2(49, &dk, 49, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
112 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
113 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
114 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
115 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
116 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
117 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
118 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
120 …uik[7] = -(d[0] * u7 + d[7] * u8 + d[14] * u9 + d[21] * u10 + d[28] * u11 + d[35] * u12 + d[42] *… in MatCholeskyFactorNumeric_SeqSBAIJ_7_NaturalOrdering()
[all …]
|
| H A D | sbaijfact11.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_7() local
27 PetscCall(PetscMalloc2(49, &dk, 49, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
149 …uik[0] = -(d[0] * u0 + d[7] * u1 + d[14] * u2 + d[21] * u3 + d[28] * u4 + d[35] * u5 + d[42] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
150 …uik[1] = -(d[1] * u0 + d[8] * u1 + d[15] * u2 + d[22] * u3 + d[29] * u4 + d[36] * u5 + d[43] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
151 …uik[2] = -(d[2] * u0 + d[9] * u1 + d[16] * u2 + d[23] * u3 + d[30] * u4 + d[37] * u5 + d[44] * u6); in MatCholeskyFactorNumeric_SeqSBAIJ_7()
152 …uik[3] = -(d[3] * u0 + d[10] * u1 + d[17] * u2 + d[24] * u3 + d[31] * u4 + d[38] * u5 + d[45] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
153 …uik[4] = -(d[4] * u0 + d[11] * u1 + d[18] * u2 + d[25] * u3 + d[32] * u4 + d[39] * u5 + d[46] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
154 …uik[5] = -(d[5] * u0 + d[12] * u1 + d[19] * u2 + d[26] * u3 + d[33] * u4 + d[40] * u5 + d[47] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
155 …uik[6] = -(d[6] * u0 + d[13] * u1 + d[20] * u2 + d[27] * u3 + d[34] * u4 + d[41] * u5 + d[48] * u6… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
157 …uik[7] = -(d[0] * u7 + d[7] * u8 + d[14] * u9 + d[21] * u10 + d[28] * u11 + d[35] * u12 + d[42] *… in MatCholeskyFactorNumeric_SeqSBAIJ_7()
[all …]
|
| H A D | sbaijfact8.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
25 PetscCall(PetscMalloc2(25, &dk, 25, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
58 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
59 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
60 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
61 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
62 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
64 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
65 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
66 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
|
| H A D | sbaijfact9.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6() local
26 PetscCall(PetscMalloc2(36, &dk, 36, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
135 uik[0] = -(d[0] * u0 + d[6] * u1 + d[12] * u2 + d[18] * u3 + d[24] * u4 + d[30] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
136 uik[1] = -(d[1] * u0 + d[7] * u1 + d[13] * u2 + d[19] * u3 + d[25] * u4 + d[31] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
137 uik[2] = -(d[2] * u0 + d[8] * u1 + d[14] * u2 + d[20] * u3 + d[26] * u4 + d[32] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
138 uik[3] = -(d[3] * u0 + d[9] * u1 + d[15] * u2 + d[21] * u3 + d[27] * u4 + d[33] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
139 uik[4] = -(d[4] * u0 + d[10] * u1 + d[16] * u2 + d[22] * u3 + d[28] * u4 + d[34] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
140 uik[5] = -(d[5] * u0 + d[11] * u1 + d[17] * u2 + d[23] * u3 + d[29] * u4 + d[35] * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
142 uik[6] = -(d[0] * u6 + d[6] * u7 + d[12] * u8 + d[18] * u9 + d[24] * u10 + d[30] * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
143 uik[7] = -(d[1] * u6 + d[7] * u7 + d[13] * u8 + d[19] * u9 + d[25] * u10 + d[31] * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6()
[all …]
|
| H A D | sbaijfact7.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
24 PetscCall(PetscMalloc2(25, &dk, 25, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
96 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
97 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
98 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
99 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
100 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
102 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
103 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
104 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
[all …]
|
| H A D | sbaijfact5.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
26 PetscCall(PetscMalloc2(16, &dk, 16, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
59 uik[0] = -(diag[0] * u[0] + diag[4] * u[1] + diag[8] * u[2] + diag[12] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
60 uik[1] = -(diag[1] * u[0] + diag[5] * u[1] + diag[9] * u[2] + diag[13] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
61 uik[2] = -(diag[2] * u[0] + diag[6] * u[1] + diag[10] * u[2] + diag[14] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
62 uik[3] = -(diag[3] * u[0] + diag[7] * u[1] + diag[11] * u[2] + diag[15] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
64 uik[4] = -(diag[0] * u[4] + diag[4] * u[5] + diag[8] * u[6] + diag[12] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
65 uik[5] = -(diag[1] * u[4] + diag[5] * u[5] + diag[9] * u[6] + diag[13] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
66 uik[6] = -(diag[2] * u[4] + diag[6] * u[5] + diag[10] * u[6] + diag[14] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
67 uik[7] = -(diag[3] * u[4] + diag[7] * u[5] + diag[11] * u[6] + diag[15] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
[all …]
|
| H A D | sbaijfact6.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
25 PetscCall(PetscMalloc2(16, &dk, 16, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
97 uik[0] = -(diag[0] * u[0] + diag[4] * u[1] + diag[8] * u[2] + diag[12] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
98 uik[1] = -(diag[1] * u[0] + diag[5] * u[1] + diag[9] * u[2] + diag[13] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
99 uik[2] = -(diag[2] * u[0] + diag[6] * u[1] + diag[10] * u[2] + diag[14] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
100 uik[3] = -(diag[3] * u[0] + diag[7] * u[1] + diag[11] * u[2] + diag[15] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
102 uik[4] = -(diag[0] * u[4] + diag[4] * u[5] + diag[8] * u[6] + diag[12] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
103 uik[5] = -(diag[1] * u[4] + diag[5] * u[5] + diag[9] * u[6] + diag[13] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
104 uik[6] = -(diag[2] * u[4] + diag[6] * u[5] + diag[10] * u[6] + diag[14] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
105 uik[7] = -(diag[3] * u[4] + diag[7] * u[5] + diag[11] * u[6] + diag[15] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
[all …]
|
| H A D | sbaijfact4.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
25 PetscCall(PetscMalloc2(9, &dk, 9, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
58 uik[0] = -(diag[0] * u[0] + diag[3] * u[1] + diag[6] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
59 uik[1] = -(diag[1] * u[0] + diag[4] * u[1] + diag[7] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
60 uik[2] = -(diag[2] * u[0] + diag[5] * u[1] + diag[8] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
62 uik[3] = -(diag[0] * u[3] + diag[3] * u[4] + diag[6] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
63 uik[4] = -(diag[1] * u[3] + diag[4] * u[4] + diag[7] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
64 uik[5] = -(diag[2] * u[3] + diag[5] * u[4] + diag[8] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
66 uik[6] = -(diag[0] * u[6] + diag[3] * u[7] + diag[6] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
67 uik[7] = -(diag[1] * u[6] + diag[4] * u[7] + diag[7] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
[all …]
|
| H A D | sbaijfact3.c | 11 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
24 PetscCall(PetscMalloc2(9, &dk, 9, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
96 uik[0] = -(diag[0] * u[0] + diag[3] * u[1] + diag[6] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
97 uik[1] = -(diag[1] * u[0] + diag[4] * u[1] + diag[7] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
98 uik[2] = -(diag[2] * u[0] + diag[5] * u[1] + diag[8] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
100 uik[3] = -(diag[0] * u[3] + diag[3] * u[4] + diag[6] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
101 uik[4] = -(diag[1] * u[3] + diag[4] * u[4] + diag[7] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
102 uik[5] = -(diag[2] * u[3] + diag[5] * u[4] + diag[8] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
104 uik[6] = -(diag[0] * u[6] + diag[3] * u[7] + diag[6] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
105 uik[7] = -(diag[1] * u[6] + diag[4] * u[7] + diag[7] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
[all …]
|
| H A D | sbaijfact10.c | 12 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
33 PetscCall(PetscMalloc2(36, &dk, 36, &uik)); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
140 m0 = uik[0] = -(d0 * u0 + d6 * u1 + d12 * u2 + d18 * u3 + d24 * u4 + d30 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
141 m1 = uik[1] = -(d1 * u0 + d7 * u1 + d13 * u2 + d19 * u3 + d25 * u4 + d31 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
142 m2 = uik[2] = -(d2 * u0 + d8 * u1 + d14 * u2 + d20 * u3 + d26 * u4 + d32 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
143 m3 = uik[3] = -(d3 * u0 + d9 * u1 + d15 * u2 + d21 * u3 + d27 * u4 + d33 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
144 m4 = uik[4] = -(d4 * u0 + d10 * u1 + d16 * u2 + d22 * u3 + d28 * u4 + d34 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
145 m5 = uik[5] = -(d5 * u0 + d11 * u1 + d17 * u2 + d23 * u3 + d29 * u4 + d35 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
147 m6 = uik[6] = -(d0 * u6 + d6 * u7 + d12 * u8 + d18 * u9 + d24 * u10 + d30 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
148 m7 = uik[7] = -(d1 * u6 + d7 * u7 + d13 * u8 + d19 * u9 + d25 * u10 + d31 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
[all …]
|
| H A D | sbaijfact.c | 566 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
581 PetscCall(PetscMalloc3(bs2, &dk, bs2, &uik, bs, &work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
652 PetscCall(PetscArrayzero(uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
653 PetscKernel_A_gets_A_minus_B_times_C(bs, uik, diag, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
656 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, dk, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
660 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
670 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, rtmp_ptr, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
715 PetscCall(PetscFree3(dk, uik, work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
739 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
751 PetscCall(PetscMalloc3(bs2, &dk, bs2, &uik, bs, &work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
[all …]
|