Home
last modified time | relevance | path

Searched refs:imax (Results 1 – 25 of 50) sorted by relevance

12

/petsc/src/sys/tutorials/
H A Dex3f.F9024 integer imax
26 parameter(imax=10000)
58 PetscCallA(PetscLogFlops(imax*onefp, ierr))
62 PetscCallA(PetscLogFlops(imax*onefp, ierr))
66 PetscCallA(PetscLogFlops(imax*onefp, ierr))
70 PetscCallA(PetscLogFlops(imax*onefp, ierr))
74 PetscCallA(PetscLogFlops(imax*onefp, ierr))
78 PetscCallA(PetscLogFlops(imax*onefp, ierr))
82 PetscCallA(PetscLogFlops(imax*onefp, ierr))
86 PetscCallA(PetscLogFlops(imax*onefp, ierr))
[all …]
H A Dex3.c17 int i, imax = 10000, icount; in main() local
37 for (i = 0; i < imax; i++) icount++; in main()
39 PetscCall(PetscLogFlops(imax)); in main()
/petsc/src/dm/tutorials/
H A Dex13f90aux.F9021 function dfdt_vdp(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f)
27 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n
29 PetscReal, dimension(n, imax, jmax, kmax) :: dfdt_vdp
38 …recursive subroutine forw_euler(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, neq, y, dfd… argument
40 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, neq
46 function dfdt(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n, f)
49 PetscInt, intent(in) :: ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, n
51 PetscReal, dimension(n, imax, jmax, kmax) :: dfdt
56 …y(:, 1:imax, 1:jmax, 1:kmax) = y(:, 1:imax, 1:jmax, 1:kmax) + dt*dfdt(t, dt, ib1, ibn, jb1, jbn, k…
H A Dex13f90.F9037 …PetscInt :: i, j, k, igmax, jgmax, kgmax, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, itime, m…
83 PetscCallA(DMDAGetCorners(SolScal, i, j, k, imax, jmax, kmax, ierr))
87 ibn = imax + stw
104 xmax = grid(1, imax, jmax, kmax)
105 ymax = grid(2, imax, jmax, kmax)
106 zmax = grid(3, imax, jmax, kmax)
140 call forw_euler(t, dt, ib1, ibn, jb1, jbn, kb1, kbn, imax, jmax, kmax, dof, f, dfdt_vdp)
/petsc/src/ts/tutorials/
H A Dex74.c38 PetscInt imax; /* number of grid points */ member
59 ctxt.imax = 40; in main()
64 PetscCall(PetscOptionsInt("-imax", "grid size", "<20>", ctxt.imax, &ctxt.imax, NULL)); in main()
72 PetscCall(VecSetSizes(u, PETSC_DECIDE, ctxt.imax)); in main()
80 PetscCall(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, ctxt.imax, ctxt.imax)); in main()
104 err = PetscSqrtReal(err * err / ((PetscReal)ctxt.imax)); in main()
124 dx = (ctxt->xmax - ctxt->xmin) / ((PetscReal)ctxt->imax); in ExactSolution()
151 dx = (user->xmax - user->xmin) / ((PetscReal)user->imax); in RHSJacobian()
173 col[0] = user->imax - 1; in RHSJacobian()
176 } else if (i == user->imax - 1) { in RHSJacobian()
/petsc/src/sys/tests/
H A Dex26.c7 int i, imax; in main() local
13 imax = 32; in main()
15 imax = 1024; in main()
24 for (i = 0; i < imax; ++i) { in main()
/petsc/src/ksp/ksp/tutorials/
H A Dex74.c72 PetscInt imax; /* number of grid points */ member
103 ctxt.imax = 20; in main()
110 PetscCall(PetscOptionsInt("-imax", "grid size", "<20>", ctxt.imax, &ctxt.imax, NULL)); in main()
123 PetscCall(VecSetSizes(u, PETSC_DECIDE, ctxt.imax)); in main()
175 PetscCall(MatSetSizes(Identity, matie - matis, matie - matis, ctxt.imax, ctxt.imax)); in main()
198 PetscCall(VecSetSizes(z, PETSC_DECIDE, ctxt.imax * nstages)); in main()
230 err = PetscSqrtReal(err * err / ((PetscReal)ctxt.imax)); in main()
259 dx = (ctxt->xmax - ctxt->xmin) / ((PetscReal)ctxt->imax); in ExactSolution()
323 dx = (user->xmax - user->xmin) / ((PetscReal)user->imax); in Assemble_AdvDiff()
327 PetscCall(MatSetSizes(*J, PETSC_DECIDE, PETSC_DECIDE, user->imax, user->imax)); in Assemble_AdvDiff()
[all …]
/petsc/src/vec/is/is/utils/
H A Disdiff.c34 PetscInt i, n1, n2, imin, imax, nout, *iout; in ISDifference() local
54 imax = 0; in ISDifference()
58 imax = PetscMax(imax, i1[i]); in ISDifference()
60 } else imin = imax = 0; in ISDifference()
62 PetscCall(PetscBTCreate(imax - imin, &mask)); in ISDifference()
73 if (i2[i] < imin || i2[i] > imax) continue; in ISDifference()
80 for (i = 0; i < imax - imin + 1; i++) { in ISDifference()
87 for (i = 0; i < imax - imin + 1; i++) { in ISDifference()
270 PetscInt i, n1, n2, imin, imax, nout, *iout; in ISExpand() local
298 imax = 0; in ISExpand()
[all …]
H A Disblock.c23 PetscErrorCode ISCompressIndicesGeneral(PetscInt n, PetscInt nkeys, PetscInt bs, PetscInt imax, con… in ISCompressIndicesGeneral() argument
46 for (i = 0; i < imax; i++) { in ISCompressIndicesGeneral()
126 PetscErrorCode ISExpandIndicesGeneral(PetscInt n, PetscInt nkeys, PetscInt bs, PetscInt imax, const… in ISExpandIndicesGeneral() argument
135 for (i = 0; i < imax; i++) { in ISExpandIndicesGeneral()
141 for (i = 0; i < imax; i++) { in ISExpandIndicesGeneral()
/petsc/src/sys/classes/draw/impls/x/
H A Dxtext.c121 int i, max, imax, tmp; in PetscDrawXiMatchFontSize() local
134 imax = 0; in PetscDrawXiMatchFontSize()
140 imax = i; in PetscDrawXiMatchFontSize()
145 font->font_w = nfonts[imax].w; in PetscDrawXiMatchFontSize()
146 font->font_h = nfonts[imax].h; in PetscDrawXiMatchFontSize()
147 font->font_descent = nfonts[imax].descent; in PetscDrawXiMatchFontSize()
/petsc/src/mat/tests/
H A Dex114.c12 PetscInt m, n, j, imin[M], imax[M], imaxabs[M], indices[N], row, testcase = 0; in main() local
138 PetscCall(MatGetRowMax(A, max, imax)); in main()
141 PetscCall(PetscIntView(n, imax, PETSC_VIEWER_STDOUT_WORLD)); in main()
153 …k(imin[j] == imax[j], PETSC_COMM_SELF, PETSC_ERR_PLIB, "imin[%" PetscInt_FMT "] %" PetscInt_FMT " … in main()
179 PetscCall(MatGetRowMax(Adense, max_d, imax)); in main()
192 …= imax[j], PETSC_COMM_SELF, PETSC_ERR_PLIB, "imin[%" PetscInt_FMT "] %" PetscInt_FMT " != imax %" … in main()
225 …B, "imaxabs[%" PetscInt_FMT "] %" PetscInt_FMT " != imaxabsB %" PetscInt_FMT, j, imin[j], imax[j]); in main()
/petsc/src/binding/petsc4py/demo/legacy/kspsolve/
H A Dpetsc-cg.py1 def cg(A, b, x, imax=50, eps=1e-6): argument
19 while i < imax and \
/petsc/src/tao/leastsquares/tutorials/matlab/
H A Dmatlab_ls_test.c104 int i, i0, imax; in main() local
127 imax = 53; in main()
130 imax = (int)prob_id; in main()
133 for (i = i0; i <= imax; ++i) { in main()
/petsc/src/mat/impls/baij/mpi/
H A Dbaijov.c11 PetscErrorCode MatIncreaseOverlap_MPIBAIJ(Mat C, PetscInt imax, IS is[], PetscInt ov) in MatIncreaseOverlap_MPIBAIJ() argument
18 PetscCall(PetscMalloc1(imax, &is_new)); in MatIncreaseOverlap_MPIBAIJ()
20 PetscCall(ISCompressIndicesGeneral(N, C->rmap->n, bs, imax, is, is_new)); in MatIncreaseOverlap_MPIBAIJ()
22 for (i = 0; i < ov; ++i) PetscCall(MatIncreaseOverlap_MPIBAIJ_Once(C, imax, is_new)); in MatIncreaseOverlap_MPIBAIJ()
23 for (i = 0; i < imax; i++) { in MatIncreaseOverlap_MPIBAIJ()
57 PetscErrorCode MatIncreaseOverlap_MPIBAIJ_Once(Mat C, PetscInt imax, IS is[]) in MatIncreaseOverlap_MPIBAIJ_Once() argument
80 PetscCall(PetscMalloc2(imax, (PetscInt ***)&idx, imax, &n)); in MatIncreaseOverlap_MPIBAIJ_Once()
82 for (PetscInt i = 0; i < imax; i++) { in MatIncreaseOverlap_MPIBAIJ_Once()
90 for (PetscInt i = 0; i < imax; i++) { in MatIncreaseOverlap_MPIBAIJ_Once()
166 …PetscCall(PetscCalloc5(imax, &table, imax, &data, imax, &isz, Mbs * imax, &d_p, (Mbs / PETSC_BITS_… in MatIncreaseOverlap_MPIBAIJ_Once()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaij.c147 if (a->free_imax_ilen) PetscCall(PetscFree2(a->imax, a->ilen)); in MatDestroy_SeqSBAIJ()
579 PetscInt *imax = a->imax, *ai = a->i, *ailen = a->ilen; in MatSetValuesBlocked_SeqSBAIJ() local
594 rmax = imax[row]; in MatSetValuesBlocked_SeqSBAIJ()
648 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqSBAIJ()
677 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqSBAIJ() local
688 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqSBAIJ()
700 fshift += imax[mbs - 1] - ailen[mbs - 1]; in MatAssemblyEnd_SeqSBAIJ()
704 for (i = 0; i < mbs; i++) ailen[i] = imax[i] = ai[i + 1] - ai[i]; in MatAssemblyEnd_SeqSBAIJ()
741 PetscInt *imax = a->imax, *ai = a->i, *ailen = a->ilen, roworiented = a->roworiented; in MatSetValues_SeqSBAIJ() local
754 rmax = imax[brow]; /* maximum space allocated for this row */ in MatSetValues_SeqSBAIJ()
[all …]
/petsc/src/dm/impls/da/hypre/
H A Dmhyp.c132 ex->hbox.imax[0] = hupper[0]; in MatSetUp_HYPREStruct()
133 ex->hbox.imax[1] = hupper[1]; in MatSetUp_HYPREStruct()
134 ex->hbox.imax[2] = hupper[2]; in MatSetUp_HYPREStruct()
509 …size = (ex->hbox.imax[0] - ex->hbox.imin[0] + 1) * (ex->hbox.imax[1] - ex->hbox.imin[1] + 1) * (ex… in MatZeroEntries_HYPRESStruct_3d()
516 iupper[i] = (HYPRE_Int)ex->hbox.imax[i]; in MatZeroEntries_HYPRESStruct_3d()
556 ex->hbox.imax[0] = (HYPRE_Int)iupper[0]; in MatSetUp_HYPRESStruct()
557 ex->hbox.imax[1] = (HYPRE_Int)iupper[1]; in MatSetUp_HYPRESStruct()
558 ex->hbox.imax[2] = (HYPRE_Int)iupper[2]; in MatSetUp_HYPRESStruct()
568 PetscCallHYPRE(HYPRE_SStructGridSetExtents(ex->ss_grid, part, ex->hbox.imin, ex->hbox.imax)); in MatSetUp_HYPRESStruct()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiov.c29 static PetscErrorCode ISAdjustForBlockSize(PetscInt bs, PetscInt imax, IS is[]) in ISAdjustForBlockSize() argument
32 for (PetscInt i = 0; i < imax; i++) { in ISAdjustForBlockSize()
66 PetscErrorCode MatIncreaseOverlap_MPIAIJ(Mat C, PetscInt imax, IS is[], PetscInt ov) in MatIncreaseOverlap_MPIAIJ() argument
72 for (i = 0; i < ov; ++i) PetscCall(MatIncreaseOverlap_MPIAIJ_Once(C, imax, is)); in MatIncreaseOverlap_MPIAIJ()
73 …rmap->bs > 1 && C->rmap->bs == C->cmap->bs) PetscCall(ISAdjustForBlockSize(C->rmap->bs, imax, is)); in MatIncreaseOverlap_MPIAIJ()
77 PetscErrorCode MatIncreaseOverlap_MPIAIJ_Scalable(Mat C, PetscInt imax, IS is[], PetscInt ov) in MatIncreaseOverlap_MPIAIJ_Scalable() argument
83 for (i = 0; i < ov; ++i) PetscCall(MatIncreaseOverlap_MPIAIJ_Once_Scalable(C, imax, is)); in MatIncreaseOverlap_MPIAIJ_Scalable()
84 …rmap->bs > 1 && C->rmap->bs == C->cmap->bs) PetscCall(ISAdjustForBlockSize(C->rmap->bs, imax, is)); in MatIncreaseOverlap_MPIAIJ_Scalable()
492 static PetscErrorCode MatIncreaseOverlap_MPIAIJ_Once(Mat C, PetscInt imax, IS is[]) in MatIncreaseOverlap_MPIAIJ_Once() argument
523 PetscCall(PetscMalloc2(imax, (PetscInt ***)&idx, imax, &n)); in MatIncreaseOverlap_MPIAIJ_Once()
[all …]
/petsc/src/mat/impls/blockmat/seq/
H A Dblockmat.c262 …MatSeqXAIJReallocateAIJ(A, a->mbs, 1, nrow, brow, bcol, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValues_BlockMat()
599 PetscInt fshift = 0, i, j, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_BlockMat() local
609 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_BlockMat()
623 fshift += imax[m - 1] - ailen[m - 1]; in MatAssemblyEnd_BlockMat()
627 for (i = 0; i < m; i++) ailen[i] = imax[i] = ai[i + 1] - ai[i]; in MatAssemblyEnd_BlockMat()
861 if (!bmat->imax) PetscCall(PetscMalloc2(A->rmap->n, &bmat->imax, A->rmap->n, &bmat->ilen)); in MatBlockMatSetPreallocation_BlockMat()
865 bmat->imax[i] = nnz[i]; in MatBlockMatSetPreallocation_BlockMat()
881 for (i = 1; i < bmat->mbs + 1; i++) bmat->i[i] = bmat->i[i - 1] + bmat->imax[i - 1]; in MatBlockMatSetPreallocation_BlockMat()
/petsc/src/tao/unconstrained/impls/nls/
H A Dnlsimpl.h54 PetscReal imax; /* Maximum perturbation added during initialization */ member
H A Dnls.c350 pert = PetscMin(nlsP->imax, PetscMax(nlsP->imin, nlsP->imfac * gnorm)); in TaoSolve_NLS()
412 pert = PetscMin(nlsP->imax, PetscMax(nlsP->imin, nlsP->imfac * gnorm)); in TaoSolve_NLS()
453 pert = PetscMin(nlsP->imax, PetscMax(nlsP->imin, nlsP->imfac * gnorm)); in TaoSolve_NLS()
739 …scOptionsReal("-tao_nls_imax", "maximum initial perturbation", "", nlsP->imax, &nlsP->imax, NULL)); in TaoSetFromOptions_NLS()
890 nlsP->imax = 1.0e+2; in TaoCreate_NLS()
/petsc/src/mat/impls/aij/seq/
H A Daij.c406 PetscInt *imax = a->imax, *ai = a->i, *ailen = a->ilen; in MatSetValues_SeqAIJ() local
420 rmax = imax[row]; in MatSetValues_SeqAIJ()
456 …cateAIJ_structure_only(A, A->rmap->n, 1, nrow, row, col, rmax, ai, aj, rp, imax, nonew, MatScalar); in MatSetValues_SeqAIJ()
458 …MatSeqXAIJReallocateAIJ(A, A->rmap->n, 1, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, M… in MatSetValues_SeqAIJ()
507 a->imax[row] = n; in MatSetValues_SeqAIJ_SortedFullNoPreallocation()
539 if (!a->imax) PetscCall(PetscMalloc1(A->rmap->n, &a->imax)); in MatSeqAIJSetTotalPreallocation()
574 …PetscCheck(n <= a->imax[row], PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Preallocation for row %" Pets… in MatSetValues_SeqAIJ_SortedFull()
1100 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqAIJ() local
1116 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqAIJ()
1128 fshift += imax[m - 1] - ailen[m - 1]; in MatAssemblyEnd_SeqAIJ()
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaij.c1520 if (a->free_imax_ilen) PetscCall(PetscFree2(a->imax, a->ilen)); in MatDestroy_SeqBAIJ()
1664 for (i = 0; i < nbs; i++) at->ilen[i] = at->imax[i] = ati[i + 1] - ati[i]; in MatTranspose_SeqBAIJ()
2109 PetscInt *imax = a->imax, *ai = a->i, *ailen = a->ilen; in MatSetValuesBlocked_SeqBAIJ() local
2127 rmax = imax[row]; in MatSetValuesBlocked_SeqBAIJ()
2184 …locateAIJ_structure_only(A, a->mbs, bs2, nrow, row, col, rmax, ai, aj, rp, imax, nonew, MatScalar); in MatSetValuesBlocked_SeqBAIJ()
2186 …MatSeqXAIJReallocateAIJ(A, a->mbs, bs2, nrow, row, col, rmax, aa, ai, aj, rp, ap, imax, nonew, Mat… in MatSetValuesBlocked_SeqBAIJ()
2217 PetscInt fshift = 0, i, *ai = a->i, *aj = a->j, *imax = a->imax; in MatAssemblyEnd_SeqBAIJ() local
2229 fshift += imax[i - 1] - ailen[i - 1]; in MatAssemblyEnd_SeqBAIJ()
2241 fshift += imax[mbs - 1] - ailen[mbs - 1]; in MatAssemblyEnd_SeqBAIJ()
2248 PetscCall(PetscFree2(a->imax, a->ilen)); in MatAssemblyEnd_SeqBAIJ()
[all …]
/petsc/src/sys/objects/
H A Ddestroy.c12 for (PetscInt i = 0, imax = *count_star; i < imax; ++i) PetscCall(PetscFree(tmp_star[i])); in DestroyComposedData() local
/petsc/src/mat/impls/aij/seq/seqviennacl/
H A Daijviennacl.cxx123 PetscCall(PetscFree(a->imax)); in MatViennaCLCopyFromGPU()
125 PetscCall(PetscMalloc1(m, &a->imax)); in MatViennaCLCopyFromGPU()
137 …a->imax[i] = a->ilen[i] = a->i[i + 1] - a->i[i]; //Set imax[] and ilen[] arrays at the same time a… in MatViennaCLCopyFromGPU()
/petsc/src/tao/bound/impls/bnk/
H A Dbnk.h80 PetscReal imax; /* Maximum perturbation added during initialization */ member

12