Home
last modified time | relevance | path

Searched refs:ij (Results 1 – 21 of 21) sorted by relevance

/petsc/src/vec/vec/impls/hypre/
H A Dvhyp.c9 PetscErrorCode VecHYPRE_IJVectorCreate(PetscLayout map, VecHYPRE_IJVector *ij) in VecHYPRE_IJVectorCreate() argument
16 …HYPRE(HYPRE_IJVectorCreate(map->comm, (HYPRE_Int)map->rstart, (HYPRE_Int)map->rend - 1, &nij->ij)); in VecHYPRE_IJVectorCreate()
17 PetscCallHYPRE(HYPRE_IJVectorSetObjectType(nij->ij, HYPRE_PARCSR)); in VecHYPRE_IJVectorCreate()
23 PetscCallHYPRE(HYPRE_IJVectorInitialize_v2(nij->ij, memloc)); in VecHYPRE_IJVectorCreate()
26 PetscCallHYPRE(HYPRE_IJVectorInitialize(nij->ij)); in VecHYPRE_IJVectorCreate()
28 PetscCallHYPRE(HYPRE_IJVectorAssemble(nij->ij)); in VecHYPRE_IJVectorCreate()
29 *ij = nij; in VecHYPRE_IJVectorCreate()
33 PetscErrorCode VecHYPRE_IJVectorDestroy(VecHYPRE_IJVector *ij) in VecHYPRE_IJVectorDestroy() argument
36 if (!*ij) PetscFunctionReturn(PETSC_SUCCESS); in VecHYPRE_IJVectorDestroy()
37 …PetscCheck(!(*ij)->pvec, PetscObjectComm((PetscObject)((*ij)->pvec)), PETSC_ERR_ORDER, "Forgot to … in VecHYPRE_IJVectorDestroy()
[all …]
H A Dvhyp.h10 HYPRE_IJVector ij; member
/petsc/src/ksp/pc/impls/sor/
H A Dsor.tex18 …{ij} x_j^{n+1/2} - \sum_{j \ge i} A_{ij} x_j^{n}) = (1 - \omega) x_i^n + \omega A_{ii}^{-1}( b_i -…
22 …ij} x_j^{n+1/2} - \sum_{j > i} A_{ij} x_j^{n+1}) = (1 - \omega) x_i^{n+1/2} + \omega A_{ii}^{-1}(…
28 x_i^{1/2} = \omega A_{ii}^{-1}( b_i - \sum_{j < i} A_{ij} x_j^{1/2})
32 …ega) x_i^{1/2} + \omega A_{ii}^{-1}( b_i - \sum_{j < i} A_{ij} x_j^{1/2} - \sum_{j > i} A_{ij} x_j)
90 x_i = \omega D_{ii}^{-1}(b_i - \sum_{j<i} A_{ij} x_j)
93 x_i = (D_{ii}/\omega)^{-1}(b_i - \sum_{j<i} A_{ij} x_j)
/petsc/src/mat/impls/hypre/
H A Dmhypre.c33 static PetscErrorCode MatHYPRE_IJMatrixPreallocate(Mat A_d, Mat A_o, HYPRE_IJMatrix ij) in MatHYPRE_IJMatrixPreallocate() argument
64 aux_matrix = (hypre_AuxParCSRMatrix *)hypre_IJMatrixTranslator(ij); in MatHYPRE_IJMatrixPreallocate()
66 hypre_IJMatrixTranslator(ij) = NULL; in MatHYPRE_IJMatrixPreallocate()
67 PetscCallHYPRE(HYPRE_IJMatrixSetDiagOffdSizes(ij, nnz_d, nnz_o)); in MatHYPRE_IJMatrixPreallocate()
70 … aux_matrix = (hypre_AuxParCSRMatrix *)hypre_IJMatrixTranslator(ij); in MatHYPRE_IJMatrixPreallocate()
75 PetscCallHYPRE(HYPRE_IJMatrixSetDiagOffdSizes(ij, nnz_d, nnz_o)); in MatHYPRE_IJMatrixPreallocate()
95 if (hA->ij) { in MatHYPRE_CreateFromMat()
96 if (!hA->inner_free) hypre_IJMatrixObject(hA->ij) = NULL; in MatHYPRE_CreateFromMat()
97 PetscCallHYPRE(HYPRE_IJMatrixDestroy(hA->ij)); in MatHYPRE_CreateFromMat()
99 PetscCallHYPRE(HYPRE_IJMatrixCreate(hA->comm, rstart, rend - 1, cstart, cend - 1, &hA->ij)); in MatHYPRE_CreateFromMat()
[all …]
H A Dmhypre.h9 HYPRE_IJMatrix ij; member
/petsc/src/ksp/ksp/tutorials/
H A Dex61f.F9079 integer :: i, j, ij, ij2, ii, jj, nz, ip, dx, dy, icase
131 ij = i + (j - 1)*n
141 is_diag = (ij == ij2)
143 ilist(nz) = ij
H A Dex70.c134 PetscInt ij, r, c, nc; in map_wIwDI_uJuDJ() local
139 ij = r * nc + c; in map_wIwDI_uJuDJ()
140 return ij; in map_wIwDI_uJuDJ()
H A Dex43.c472 PetscInt ij; in ASS_MAP_wIwDI_uJuDJ() local
478 ij = r * nc + c; in ASS_MAP_wIwDI_uJuDJ()
479 return ij; in ASS_MAP_wIwDI_uJuDJ()
H A Dex42.c546 PetscInt ij; in ASS_MAP_wIwDI_uJuDJ() local
555 ij = r * nc + c; in ASS_MAP_wIwDI_uJuDJ()
557 return ij; in ASS_MAP_wIwDI_uJuDJ()
/petsc/src/mat/tests/
H A Dex93.c11 PetscInt ij[] = {0, 1, 2}; in main() local
34 if (rank == 0) PetscCall(MatSetValues(A, 3, ij, 3, ij, a, ADD_VALUES)); in main()
/petsc/src/ksp/pc/impls/hypre/
H A Dhypre.c410 PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->hmnull[i]->ij, (void **)&jac->phmnull[i])); in PCSetUp_HYPRE()
418 … PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->hmnull[nvec]->ij, (void **)&jac->phmnull[nvec])); in PCSetUp_HYPRE()
434 PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->constants[0]->ij, (void **)(&ozz))); in PCSetUp_HYPRE()
435 PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->constants[1]->ij, (void **)(&zoz))); in PCSetUp_HYPRE()
436 …if (jac->constants[2]) PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->constants[2]->ij, (void **)(&zz… in PCSetUp_HYPRE()
444 …if (jac->coords[0]) PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->coords[0]->ij, (void **)(&coords[0… in PCSetUp_HYPRE()
445 …if (jac->coords[1]) PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->coords[1]->ij, (void **)(&coords[1… in PCSetUp_HYPRE()
446 …if (jac->coords[2]) PetscCallHYPRE(HYPRE_IJVectorGetObject(jac->coords[2]->ij, (void **)(&coords[2… in PCSetUp_HYPRE()
451 PetscCallHYPRE(HYPRE_IJMatrixGetObject(hm->ij, (void **)(&parcsr))); in PCSetUp_HYPRE()
455 PetscCallHYPRE(HYPRE_IJMatrixGetObject(hm->ij, (void **)(&parcsr))); in PCSetUp_HYPRE()
[all …]
/petsc/src/dm/impls/plex/
H A Dgmshlex.h276 pair ij[SN2(GMSH_MAX_ORDER)], tmp[SN2(GMSH_MAX_ORDER)]; in GmshLexOrder_PRI() local
284 for (l = 0; l < m; ++l) ij[sub[l]] = tmp[l]; in GmshLexOrder_PRI()
286 i = ij[l].i; in GmshLexOrder_PRI()
287 j = ij[l].j; in GmshLexOrder_PRI()
/petsc/src/mat/impls/aij/seq/
H A Daij.c3223 …PetscInt i, bs = A->rmap->bs, mbs = A->rmap->n / bs, ipvt[5], bs2 = bs * bs, *v_pivots, ij[… in MatInvertBlockDiagonal_SeqAIJ() local
3254 ij[0] = 2 * i; in MatInvertBlockDiagonal_SeqAIJ()
3255 ij[1] = 2 * i + 1; in MatInvertBlockDiagonal_SeqAIJ()
3256 PetscCall(MatGetValues(A, 2, ij, 2, ij, diag)); in MatInvertBlockDiagonal_SeqAIJ()
3265 ij[0] = 3 * i; in MatInvertBlockDiagonal_SeqAIJ()
3266 ij[1] = 3 * i + 1; in MatInvertBlockDiagonal_SeqAIJ()
3267 ij[2] = 3 * i + 2; in MatInvertBlockDiagonal_SeqAIJ()
3268 PetscCall(MatGetValues(A, 3, ij, 3, ij, diag)); in MatInvertBlockDiagonal_SeqAIJ()
3277 ij[0] = 4 * i; in MatInvertBlockDiagonal_SeqAIJ()
3278 ij[1] = 4 * i + 1; in MatInvertBlockDiagonal_SeqAIJ()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij2.c826 const PetscInt *ii, *ij = a->j, *idx; in MatMult_SeqBAIJ_11() local
848 idx = ij + ii[i]; in MatMult_SeqBAIJ_11()
910 const PetscInt *ii, *ij = a->j, *idx; in MatMult_SeqBAIJ_12_ver1() local
932 idx = ij + ii[i]; in MatMult_SeqBAIJ_12_ver1()
994 const PetscInt *ii, *ij = a->j, *idx; in MatMultAdd_SeqBAIJ_12_ver1() local
1016 idx = ij + ii[i]; in MatMultAdd_SeqBAIJ_12_ver1()
1087 const PetscInt *ii, *ij = a->j, *idx, *ridx = NULL; in MatMult_SeqBAIJ_12_ver2() local
1109 idx = ij + ii[i]; in MatMult_SeqBAIJ_12_ver2()
1199 const PetscInt *ii, *ij = a->j, *idx, *ridx = NULL; in MatMultAdd_SeqBAIJ_12_ver2() local
1221 idx = ij + ii[i]; in MatMultAdd_SeqBAIJ_12_ver2()
[all …]
/petsc/src/mat/impls/nest/
H A Dmatnest.c1893 PetscInt ij, rsize = nii[1] - nii[0], ist = ii[ir] + ci[ir]; in MatConvert_Nest_SeqAIJ_fast() local
1895 for (ij = 0; ij < rsize; ij++) { in MatConvert_Nest_SeqAIJ_fast()
1896 jj[ist + ij] = *njj + cst; in MatConvert_Nest_SeqAIJ_fast()
1897 vv[ist + ij] = vscale * *nvv; in MatConvert_Nest_SeqAIJ_fast()
1898 if (PetscUnlikely(vshift != 0.0 && *njj == ir - rst)) vv[ist + ij] += vshift; in MatConvert_Nest_SeqAIJ_fast()
/petsc/doc/manual/
H A Dksp.md1419 Now we can express $P$ as a matrix by looking at the matrix elements $P_{ij} = \phi^F_i P \phi^C_j$…
1424 = &f^F_i - \sum_j P_{ij} f^C_j
1431 \min_{P_{ij}} \| f^F_i - \sum_j P_{ij} f^C_j \|_2
1438 …terpolated accurately, and a single $f$ would not constrain the values $P_{ij}`$ well. Therefore, …
1442 &\min_{P_{ij}} \sum_k w_k \| f^{F,k}_i - \sum_j P_{ij} f^{C,k}_j \|_2 \\
1443 = &\min_{P_{ij}} \sum_k \| \sqrt{w_k} f^{F,k}_i - \sqrt{w_k} \sum_j P_{ij} f^{C,k}_j \|_2 \\
1444 = &\min_{P_{ij}} \| W^{1/2} \mathbf{f}^{F}_i - W^{1/2} \mathbf{f}^{C} p_i \|_2
1466 [p_i]_j &= P_{ij}
1473 \min_{P_{ij}} \| b - A x \|_2
1959 where each $A_{ij}$ is an entire block. The matrices on a parallel computer are not explicitly stor…
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A Dpetscmat.pxi1167 cdef object mat_getitem(Mat self, object ij): argument
1169 rows, cols = ij
1181 cdef PetscErrorCode mat_setitem(Mat self, object ij, object v) except PETSC_ERR_PYTHON: argument
1183 rows, cols = ij
H A DMat.pyx440 def __getitem__(self, ij): argument
441 return mat_getitem(self, ij)
443 def __setitem__(self, ij, v): argument
444 mat_setitem(self, ij, v)
/petsc/src/mat/impls/is/
H A Dmatis.c860 PetscInt l1, l2, lb1, lb2, ij = i * nc + j; in MatConvert_Nest_IS() local
866 PetscCall(PetscObjectTypeCompare((PetscObject)nest[i][j], MATTRANSPOSEVIRTUAL, &istrans[ij])); in MatConvert_Nest_IS()
867 if (istrans[ij]) { in MatConvert_Nest_IS()
874 PetscCall(MatCreateTranspose(lT, &snest[ij])); in MatConvert_Nest_IS()
879 PetscCall(MatISGetLocalMat(nest[i][j], &snest[ij])); in MatConvert_Nest_IS()
885 PetscCall(MatGetSize(snest[ij], &l1, &l2)); in MatConvert_Nest_IS()
886 PetscCall(MatGetBlockSizes(snest[ij], &lb1, &lb2)); in MatConvert_Nest_IS()
894 if (rnest && !rnest[i][j] != !snest[ij]) lreuse = PETSC_FALSE; in MatConvert_Nest_IS()
/petsc/systems/Apple/iOS/PETSc/PETSc.xcodeproj/
H A Dproject.pbxproj668 …D89D1F15207AA1AC00FACC19 /* ij.c in Sources */ = {isa = PBXBuildFile; fileRef = D89D1C2D207AA1AA00…
1420 …89D1C2D207AA1AA00FACC19 /* ij.c */ = {isa = PBXFileReference; fileEncoding = 4; lastKnownFileType …
2191 D89D1C2D207AA1AA00FACC19 /* ij.c */,
2593 D89D1F15207AA1AC00FACC19 /* ij.c in Sources */,
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcprivate.c2355 const PetscInt *ii, *ij; in MatSeqAIJCompress() local
2361 PetscCall(MatGetRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &n, &ii, &ij, &flg_row)); in MatSeqAIJCompress()
2376 if (PetscLikely(PetscAbsScalar(entry) > PETSC_SMALL) || (n == m && ij[j] == i)) { in MatSeqAIJCompress()
2377 bij[nnz] = ij[j]; in MatSeqAIJCompress()
2386 PetscCall(MatRestoreRowIJ(A, 0, PETSC_FALSE, PETSC_FALSE, &n, &ii, &ij, &flg_row)); in MatSeqAIJCompress()