Home
last modified time | relevance | path

Searched refs:ni (Results 1 – 25 of 31) sorted by relevance

12

/petsc/src/mat/tests/
H A Dex152.c29 idx_t ni, isize, *vtxdist, *xadj, *adjncy, *vwgt, *part; in main() local
63 ni = vtxdist[rank + 1] - vtxdist[rank]; in main()
65 PetscCall(PetscMalloc1(ni + 1, &xadj)); in main()
67 red = fread(xadj, sizeof(idx_t), ni + 1, fp); in main()
68 …PetscCheck(red == (size_t)(ni + 1), PETSC_COMM_SELF, PETSC_ERR_SYS, "Unable to read from data file… in main()
70 PetscCall(PetscMalloc1(xadj[ni], &adjncy)); in main()
72 for (PetscInt i = 0; i < ni; i++) { in main()
82 PetscCall(PetscMalloc3(ni * ndims, &xyz, ni, &part, size, &tpwgts)); in main()
83 PetscCall(PetscMalloc1(ni * ndims, &sxyz)); in main()
85 red = fread(xyz, sizeof(PetscReal), ndims * ni, fp); in main()
[all …]
/petsc/src/dm/dt/dualspace/impls/lagrange/
H A Ddspacelagrange.c320 PetscErrorCode PetscLagNodeIndicesGetData_Internal(PetscLagNodeIndices ni, PetscInt *nodeIdxDim, Pe… in PetscLagNodeIndicesGetData_Internal() argument
323 *nodeIdxDim = ni->nodeIdxDim; in PetscLagNodeIndicesGetData_Internal()
324 *nodeVecDim = ni->nodeVecDim; in PetscLagNodeIndicesGetData_Internal()
325 *nNodes = ni->nNodes; in PetscLagNodeIndicesGetData_Internal()
326 *nodeIdx = ni->nodeIdx; in PetscLagNodeIndicesGetData_Internal()
327 *nodeVec = ni->nodeVec; in PetscLagNodeIndicesGetData_Internal()
331 static PetscErrorCode PetscLagNodeIndicesReference(PetscLagNodeIndices ni) in PetscLagNodeIndicesReference() argument
334 if (ni) ni->refct++; in PetscLagNodeIndicesReference()
338 static PetscErrorCode PetscLagNodeIndicesDuplicate(PetscLagNodeIndices ni, PetscLagNodeIndices *niN… in PetscLagNodeIndicesDuplicate() argument
343 (*niNew)->nodeIdxDim = ni->nodeIdxDim; in PetscLagNodeIndicesDuplicate()
[all …]
/petsc/src/ksp/ksp/impls/fetidp/
H A Dfetidp.c591 PetscInt nl, ni, *widxs; in KSPFETIDPSetUpOperators() local
618 PetscCall(ISGetLocalSize(II, &ni)); in KSPFETIDPSetUpOperators()
620 for (i = 0; i < ni; i++) matis->sf_leafdata[idxs[i]] = 1; in KSPFETIDPSetUpOperators()
625 for (i = 0, ni = 0; i < nl; i++) in KSPFETIDPSetUpOperators()
626 if (matis->sf_rootdata[i]) widxs[ni++] = i + rst; in KSPFETIDPSetUpOperators()
627 … PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)ksp), ni, widxs, PETSC_COPY_VALUES, &pII)); in KSPFETIDPSetUpOperators()
699 PetscCall(ISGetLocalSize(Pall, &ni)); in KSPFETIDPSetUpOperators()
701 for (i = 0; i < ni; i++) matis->sf_rootdata[idxs[i] - rst] = 1; in KSPFETIDPSetUpOperators()
705 for (i = 0, ni = 0; i < n; i++) in KSPFETIDPSetUpOperators()
706 if (matis->sf_leafdata[i]) widxs[ni++] = i; in KSPFETIDPSetUpOperators()
[all …]
/petsc/src/dm/partitioner/impls/simple/
H A Dpartsimple.c50 PetscInt Np = 1, Nr, np, nk, nj, ni, pk, pj, pi, ck, cj, ci, i; in PetscPartitionerPartition_Simple_Grid() local
87 for (ni = 0; ni < nodes[0]; ++ni) { in PetscPartitionerPartition_Simple_Grid()
88 const PetscInt nid = (nk * nodes[1] + nj) * nodes[0] + ni; in PetscPartitionerPartition_Simple_Grid()
99 … + ck) * cells[1] + ((nj * procs[1] + pj) * pcells[1] + cj)) * cells[0] + (ni * procs[0] + pi) * p… in PetscPartitionerPartition_Simple_Grid()
/petsc/share/petsc/matlab/
H A DPetscBagRead.m70 ni = strfind(lin,nametag); variable
72 if ni
73 n = str2num(lin(ni+length(nametag):end));
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A Dpetscmat.pxi849 cdef PetscInt ni=0, *i=NULL
853 oi = iarray_i(oi, &ni, &i)
857 if m == PETSC_DECIDE: m = (ni-1)*bs
859 if ((ni-1)*bs != m):
861 (toInt(ni), toInt(m//bs+1)))
865 if (i[ni-1] != nj):
867 (toInt(nj), toInt(i[ni-1])))
929 cdef PetscInt ni=0, *i=NULL
933 oi = iarray_i(oi, &ni, &i)
936 if ni*nj*rbs*cbs != nv: raise ValueError(
[all …]
H A Dpetscvec.pxi375 cdef PetscInt ni=0, nv=0
378 cdef object unused1 = iarray_i(oi, &ni, &i)
380 if ni*bs != nv: raise ValueError(
381 "incompatible array sizes: ni=%d, nv=%d, bs=%d" %
382 (toInt(ni), toInt(nv), toInt(bs)))
387 CHKERR(setvalues(V, ni, i, v, addv))
391 cdef PetscInt ni=0, nv=0
394 cdef object indices = iarray_i(oindices, &ni, &i)
396 values = empty_s(ni)
399 if (ni != nv): raise ValueError(
[all …]
H A DMat.pyx1073 cdef PetscInt ni=0, noi=0, *i=NULL, *oi=NULL
1075 pi = iarray_i(pi, &ni, &i) # Row pointers (diagonal)
1077 if ni != m+1: raise ValueError(
1079 (toInt(m), toInt(m+1), toInt(ni)))
3264 cdef PetscInt ni=0, *i=NULL
3272 rows = iarray_i(rows, &ni, &i)
3273 CHKERR(MatZeroRows(self.mat, ni, i, sval, xvec, bvec))
3297 cdef PetscInt ni=0, *i=NULL
3305 rows = iarray_i(rows, &ni, &i)
3306 CHKERR(MatZeroRowsLocal(self.mat, ni, i, sval, xvec, bvec))
[all …]
H A DPC.pyx2129 cdef PetscInt ni=0, *i=NULL
2132 oi = iarray_i(oi, &ni, &i)
2137 if (i[ni-1] != nj):
2139 (toInt(nj), toInt(i[ni-1])))
2141 CHKERR(PCBDDCSetLocalAdjacencyGraph(self.pc, ni - 1, i, j, PETSC_COPY_VALUES))
/petsc/src/dm/tutorials/
H A Dswarm_ex1.c181 PetscInt is, js, ni, nj, overlap; in ex1_3() local
192 PetscCall(DMDAGetCorners(dmcell, &is, &js, NULL, &ni, &nj, NULL)); in ex1_3()
204 PetscCall(DMSwarmSetLocalSizes(dms, ni * nj * 4, 4)); in ex1_3()
223 for (i = is; i < is + ni; i++) { in ex1_3()
344 PetscInt is, js, ni, nj, overlap, nn; in ex1_4() local
359 PetscCall(DMDAGetCorners(dmcell, &is, &js, NULL, &ni, &nj, NULL)); in ex1_4()
370 PetscCall(DMSwarmSetLocalSizes(dms, ni * nj * 4, 4)); in ex1_4()
389 for (i = is; i < is + ni; i++) { in ex1_4()
H A Dex15.c72 PetscInt i, j, k, l, si, sj, sk, ni, nj, nk, M, N, dof; in DMDAVecGenerateEntries() local
76 PetscCall(DMDAGetCorners(dm, &si, &sj, &sk, &ni, &nj, &nk)); in DMDAVecGenerateEntries()
80 for (i = si; i < si + ni; i++) { in DMDAVecGenerateEntries()
/petsc/src/ksp/pc/impls/telescope/
H A Dtelescope_dmda.c127 PetscInt i, j, si, sj, ni, nj, M, N, Ml, Nl, c, nidx; in PCTelescopeSetUp_dmda_repart_coors2d() local
146 PetscCall(DMDAGetCorners(subdm, &si, &sj, NULL, &ni, &nj, NULL)); in PCTelescopeSetUp_dmda_repart_coors2d()
147 Ml = ni; in PCTelescopeSetUp_dmda_repart_coors2d()
151 ni = nj = 0; in PCTelescopeSetUp_dmda_repart_coors2d()
159 for (i = si; i < si + ni; i++) { in PCTelescopeSetUp_dmda_repart_coors2d()
218 PetscInt i, j, k, si, sj, sk, ni, nj, nk, M, N, P, Ml, Nl, Pl, c, nidx; in PCTelescopeSetUp_dmda_repart_coors3d() local
239 PetscCall(DMDAGetCorners(subdm, &si, &sj, &sk, &ni, &nj, &nk)); in PCTelescopeSetUp_dmda_repart_coors3d()
240 Ml = ni; in PCTelescopeSetUp_dmda_repart_coors3d()
245 ni = nj = nk = 0; in PCTelescopeSetUp_dmda_repart_coors3d()
255 for (i = si; i < si + ni; i++) { in PCTelescopeSetUp_dmda_repart_coors3d()
[all …]
/petsc/src/vec/vec/interface/
H A Drvector.c922 PetscErrorCode VecSetValues(Vec x, PetscInt ni, const PetscInt ix[], const PetscScalar y[], InsertM… in VecSetValues() argument
926 if (!ni) PetscFunctionReturn(PETSC_SUCCESS); in VecSetValues()
932 PetscUseTypeMethod(x, setvalues, ni, ix, y, iora); in VecSetValues()
969 PetscErrorCode VecGetValues(Vec x, PetscInt ni, const PetscInt ix[], PetscScalar y[]) in VecGetValues() argument
973 if (!ni) PetscFunctionReturn(PETSC_SUCCESS); in VecGetValues()
978 PetscUseTypeMethod(x, getvalues, ni, ix, y); in VecGetValues()
1023 PetscErrorCode VecSetValuesBlocked(Vec x, PetscInt ni, const PetscInt ix[], const PetscScalar y[], … in VecSetValuesBlocked() argument
1027 if (!ni) PetscFunctionReturn(PETSC_SUCCESS); in VecSetValuesBlocked()
1033 PetscUseTypeMethod(x, setvaluesblocked, ni, ix, y, iora); in VecSetValuesBlocked()
1075 PetscErrorCode VecSetValuesLocal(Vec x, PetscInt ni, const PetscInt ix[], const PetscScalar y[], In… in VecSetValuesLocal() argument
[all …]
/petsc/src/dm/impls/plex/
H A Dplexdd.c52 PetscInt n, mh, ddovl = 0, pStart, pEnd, ni, no, nl; in DMCreateDomainDecomposition_Plex() local
197 ni = ren - rst; in DMCreateDomainDecomposition_Plex()
212 …PetscCheck(c == ni, PETSC_COMM_SELF, PETSC_ERR_PLIB, "%" PetscInt_FMT " != %" PetscInt_FMT, c, ni); in DMCreateDomainDecomposition_Plex()
213 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)odm), ni, lidxs, PETSC_OWN_POINTER, &li_is)… in DMCreateDomainDecomposition_Plex()
/petsc/src/tao/constrained/tutorials/
H A Dex1.c38 PetscInt ni; /* Global number of inequality constraints */ member
197 user->ni = 2; in InitializeProblem()
199 niloc = (size == 1) ? user->ni : 1; in InitializeProblem()
209 PetscCall(VecSetSizes(user->ci, niloc, user->ni)); in InitializeProblem()
222 PetscCall(MatSetSizes(user->Ai, niloc, nloc, user->ni, user->n)); in InitializeProblem()
/petsc/src/snes/tutorials/
H A Dex30.c96 PetscInt corner, ni, nj, jlid, jfault, inose; member
153 …PetscCall(DMDACreate2d(comm, grid.bx, grid.by, grid.stencil, grid.ni, grid.nj, PETSC_DECIDE, PETSC… in main()
359 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1; in CalcSecInv()
521 PetscInt ilim = grid->ni - 1; in ZMomentumResidual()
588 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, jlid = grid->jlid; in EnergyResidual()
646 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1; in ShearStress()
681 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc; in XNormalStress()
720 PetscInt ilim = grid->ni - 1, jlim = grid->nj - 1, ivisc; in ZNormalStress()
781 grid->ni = 82; in SetParams()
782 PetscCall(PetscOptionsGetInt(NULL, NULL, "-ni", &grid->ni, NULL)); in SetParams()
[all …]
/petsc/src/vec/vec/impls/seq/
H A Dbvec2.c489 PetscErrorCode VecGetValues_Seq(Vec xin, PetscInt ni, const PetscInt ix[], PetscScalar y[]) in VecGetValues_Seq() argument
496 for (PetscInt i = 0; i < ni; ++i) { in VecGetValues_Seq()
508 PetscErrorCode VecSetValues_Seq(Vec xin, PetscInt ni, const PetscInt ix[], const PetscScalar y[], I… in VecSetValues_Seq() argument
518 for (PetscInt i = 0; i < ni; i++) { in VecSetValues_Seq()
535 PetscErrorCode VecSetValuesBlocked_Seq(Vec xin, PetscInt ni, const PetscInt ix[], const PetscScalar… in VecSetValuesBlocked_Seq() argument
544 for (PetscInt i = 0; i < ni; ++i, yin += bs) { in VecSetValuesBlocked_Seq()
/petsc/src/snes/tutorials/output/
H A Dex30_1.out8 [ni,nj] = 18, 18 [dx,dz] = 20., 20. km
/petsc/src/binding/petsc4py/test/
H A Dtest_vec.py52 ni = self.vec.norm(PETSc.NormType.NORM_INFINITY)
56 self.assertAlmostEqual(ni, 1.0)
65 self.assertAlmostEqual(nni, ni)
/petsc/src/vec/vec/impls/mpi/
H A Dpdvec.c710 PetscErrorCode VecGetValues_MPI(Vec xin, PetscInt ni, const PetscInt ix[], PetscScalar y[]) in VecGetValues_MPI() argument
717 for (PetscInt i = 0; i < ni; i++) { in VecGetValues_MPI()
728 PetscErrorCode VecSetValues_MPI(Vec xin, PetscInt ni, const PetscInt ix[], const PetscScalar y[], I… in VecSetValues_MPI() argument
744 for (PetscInt i = 0; i < ni; ++i) { in VecSetValues_MPI()
765 PetscErrorCode VecSetValuesBlocked_MPI(Vec xin, PetscInt ni, const PetscInt ix[], const PetscScalar… in VecSetValuesBlocked_MPI() argument
781 for (i = 0; i < ni; i++) { in VecSetValuesBlocked_MPI()
795 for (i = 0; i < ni; i++) { in VecSetValuesBlocked_MPI()
/petsc/src/dm/impls/da/
H A Dgr2.c730 PetscInt dim, ni, nj, nk, pi, pj, pk, dof, n; in VecView_MPI_DA() local
736 …PetscCall(DMDAGetInfo(da, &dim, &ni, &nj, &nk, &pi, &pj, &pk, &dof, NULL, NULL, NULL, NULL, NULL)); in VecView_MPI_DA()
739 …etscFPrintf(comm, info, "#$$ tmp = reshape(tmp,%" PetscInt_FMT ",%" PetscInt_FMT ");\n", dof, ni)); in VecView_MPI_DA()
740 …, "#$$ tmp = reshape(tmp,%" PetscInt_FMT ",%" PetscInt_FMT ",%" PetscInt_FMT ");\n", dof, ni, nj)); in VecView_MPI_DA()
741 …mp,%" PetscInt_FMT ",%" PetscInt_FMT ",%" PetscInt_FMT ",%" PetscInt_FMT ");\n", dof, ni, nj, nk)); in VecView_MPI_DA()
/petsc/src/ksp/pc/impls/mpi/
H A Dpcmpi.c247 PetscMPIInt ni; in PCMPISetMat() local
249 PetscCall(PetscMPIIntCast(matproperties[7] + 1, &ni)); in PCMPISetMat()
251 PetscCallMPI(MPI_Bcast(cprefix, ni, MPI_CHAR, 0, comm)); in PCMPISetMat()
/petsc/src/mat/impls/nest/
H A Dmatnest.c1103 PetscInt m, n, M, N, mi, ni, Mi, Ni; in MatNestSetSubMat_Nest() local
1113 PetscCall(ISGetLocalSize(bA->isglobal.col[jdxm], &ni)); in MatNestSetSubMat_Nest()
1116 …ni, PetscObjectComm((PetscObject)mat), PETSC_ERR_ARG_INCOMP, "Submatrix local dimension (%" PetscI… in MatNestSetSubMat_Nest()
1718 PetscInt m, n, M, N, mi, ni, Mi, Ni; in MatSetUp_NestIS_Private() local
1726 PetscCall(ISGetLocalSize(vs->isglobal.col[j], &ni)); in MatSetUp_NestIS_Private()
1728 …ni, PetscObjectComm((PetscObject)sub), PETSC_ERR_ARG_INCOMP, "Local sizes (%" PetscInt_FMT ",%" Pe… in MatSetUp_NestIS_Private()
/petsc/src/mat/impls/h2opus/cuda/
H A Dmath2opus.cu956 PetscReal n1, ni, n2; in MatAssemblyEnd_H2OPUS() local
965 PetscCall(MatNorm(E, NORM_INFINITY, &ni)); in MatAssemblyEnd_H2OPUS()
996 …ITY %g, NORM_2 %g (%g %g %g)\n", (double)n1, (double)ni, (double)n2, (double)(n1 / n1A), (double)(… in MatAssemblyEnd_H2OPUS()
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij2.c231 PetscInt n1, n2, ni; in MatCreateSubMatrix_SeqSBAIJ() local
242 PetscCall(ISGetLocalSize(intersect, &ni)); in MatCreateSubMatrix_SeqSBAIJ()
244 if (ni == 0) sym = PETSC_FALSE; in MatCreateSubMatrix_SeqSBAIJ()
248 …PetscCheck(ni == n1 && ni == n2, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Cannot create such a subma… in MatCreateSubMatrix_SeqSBAIJ()

12