| /petsc/src/mat/impls/baij/seq/ftn-kernels/ |
| H A D | fsolvebaij.F90 | 20 PetscInt :: idx, ax, jdx 42 jdx = 4*aj(j) 44 …s(0) = s(0) - (a(ax + 0)*x(jdx + 0) + a(ax + 4)*x(jdx + 1) + a(ax + 8)*x(jdx + 2) + a(ax + 12)*x(j… 45 …s(1) = s(1) - (a(ax + 1)*x(jdx + 0) + a(ax + 5)*x(jdx + 1) + a(ax + 9)*x(jdx + 2) + a(ax + 13)*x(j… 46 …s(2) = s(2) - (a(ax + 2)*x(jdx + 0) + a(ax + 6)*x(jdx + 1) + a(ax + 10)*x(jdx + 2) + a(ax + 14)*x(… 47 …s(3) = s(3) - (a(ax + 3)*x(jdx + 0) + a(ax + 7)*x(jdx + 1) + a(ax + 11)*x(jdx + 2) + a(ax + 15)*x(… 62 jdx = 4*aj(j) 63 …s(0) = s(0) - (a(ax + 0)*x(jdx + 0) + a(ax + 4)*x(jdx + 1) + a(ax + 8)*x(jdx + 2) + a(ax + 12)*x(j… 64 …s(1) = s(1) - (a(ax + 1)*x(jdx + 0) + a(ax + 5)*x(jdx + 1) + a(ax + 9)*x(jdx + 2) + a(ax + 13)*x(j… 65 …s(2) = s(2) - (a(ax + 2)*x(jdx + 0) + a(ax + 6)*x(jdx + 1) + a(ax + 10)*x(jdx + 2) + a(ax + 14)*x(… [all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvnat5.c | 8 PetscInt i, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() local 34 jdx = 5 * (*vi++); in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 35 x1 = x[jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 36 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 37 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 38 x4 = x[3 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 39 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering_inplace() 96 PetscInt i, k, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_5_NaturalOrdering() local 122 jdx = 5 * vi[k]; in MatSolve_SeqBAIJ_5_NaturalOrdering() 123 x1 = x[jdx]; in MatSolve_SeqBAIJ_5_NaturalOrdering() [all …]
|
| H A D | baijsolvnat6.c | 7 PetscInt i, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() local 36 jdx = 6 * (*vi++); in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 37 x1 = x[jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 38 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 39 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 40 x4 = x[3 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 41 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 42 x6 = x[5 + jdx]; in MatSolve_SeqBAIJ_6_NaturalOrdering_inplace() 105 PetscInt i, k, nz, idx, jdx, idt; in MatSolve_SeqBAIJ_6_NaturalOrdering() local 135 jdx = bs * vi[k]; in MatSolve_SeqBAIJ_6_NaturalOrdering() [all …]
|
| H A D | baijsolvnat3.c | 16 PetscInt jdx, idt, idx, nz, i; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() local 36 jdx = 3 * (*vi++); in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() 37 x1 = x[jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() 38 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() 39 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering_inplace() 84 PetscInt i, k, nz, idx, jdx, idt; in MatSolve_SeqBAIJ_3_NaturalOrdering() local 108 jdx = bs * vi[k]; in MatSolve_SeqBAIJ_3_NaturalOrdering() 109 x1 = x[jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering() 110 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering() 111 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_3_NaturalOrdering() [all …]
|
| H A D | baijsolvnat2.c | 19 PetscInt jdx, idt, idx, nz, i; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() local 37 jdx = 2 * (*vi++); in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() 38 x1 = x[jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() 39 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering_inplace() 78 PetscInt i, k, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_2_NaturalOrdering() local 100 jdx = 2 * vi[k]; in MatSolve_SeqBAIJ_2_NaturalOrdering() 101 x1 = x[jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering() 102 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_2_NaturalOrdering() 144 PetscInt i, k, nz, idx, jdx; in MatForwardSolve_SeqBAIJ_2_NaturalOrdering() local 166 jdx = 2 * vi[k]; in MatForwardSolve_SeqBAIJ_2_NaturalOrdering() [all …]
|
| H A D | baijsolvnat4.c | 33 PetscInt jdx, idt, idx, nz, i, ai16; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() local 52 jdx = 4 * (*vi++); in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() 53 x1 = x[jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() 54 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() 55 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() 56 x4 = x[3 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering_inplace() 111 PetscInt i, k, nz, idx, jdx, idt; in MatSolve_SeqBAIJ_4_NaturalOrdering() local 137 jdx = bs * vi[k]; in MatSolve_SeqBAIJ_4_NaturalOrdering() 138 x1 = x[jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering() 139 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_4_NaturalOrdering() [all …]
|
| H A D | baijsolvnat7.c | 8 PetscInt i, nz, idx, idt, jdx; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() local 38 jdx = 7 * (*vi++); in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 39 x1 = x[jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 40 x2 = x[1 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 41 x3 = x[2 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 42 x4 = x[3 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 43 x5 = x[4 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 44 x6 = x[5 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 45 x7 = x[6 + jdx]; in MatSolve_SeqBAIJ_7_NaturalOrdering_inplace() 114 PetscInt i, k, nz, idx, jdx, idt; in MatSolve_SeqBAIJ_7_NaturalOrdering() local [all …]
|
| H A D | baijsolvnat1.c | 16 PetscInt jdx, idt, idx, nz, i; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() local 32 jdx = *vi++; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() 33 x1 = x[jdx]; in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
|
| H A D | baijsolv.c | 1146 PetscInt i, nz, idx, jdx, idt, idc, m; in MatSolve_SeqBAIJ_2() local 1174 jdx = 2 * vi[m]; in MatSolve_SeqBAIJ_2() 1175 x1 = t[jdx]; in MatSolve_SeqBAIJ_2() 1176 x2 = t[1 + jdx]; in MatSolve_SeqBAIJ_2()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexpartition.c | 1670 PetscInt ierr, i, j, idx, jdx; in DMPlexRebalanceSharedPoints() local 1806 if (isNonExclusivelyOwned[i]) jdx = globalNumbersOfLocalOwnedVertices[i]; in DMPlexRebalanceSharedPoints() 1807 else if (isLeaf[i]) jdx = leafGlobalNumbers[i]; in DMPlexRebalanceSharedPoints() 1809 PetscCall(MatSetValue(A, idx, jdx, 1, INSERT_VALUES)); in DMPlexRebalanceSharedPoints() 1810 PetscCall(MatSetValue(A, jdx, idx, 1, INSERT_VALUES)); in DMPlexRebalanceSharedPoints()
|