Home
last modified time | relevance | path

Searched refs:pi (Results 1 – 25 of 40) sorted by relevance

12

/petsc/src/dm/impls/da/tests/output/
H A Dex1_0.out1 [0] p 0 (+0.4619,+0.1913) r=0.50 th=0.125*pi cellid=0
2 [0] p 0 (+0.1913,+0.4619) r=0.50 th=0.375*pi cellid=0
3 [0] p 0 (-0.1913,+0.4619) r=0.50 th=0.625*pi cellid=0
4 [0] p 0 (-0.4619,+0.1913) r=0.50 th=0.875*pi cellid=0
5 [0] p 0 (-0.4619,-0.1913) r=0.50 th=-0.875*pi cellid=0
6 [0] p 0 (-0.1913,-0.4619) r=0.50 th=-0.625*pi cellid=0
7 [0] p 0 (+0.1913,-0.4619) r=0.50 th=-0.375*pi cellid=0
8 [0] p 0 (+0.4619,-0.1913) r=0.50 th=-0.125*pi cellid=0
12 [0] p 0 (+0.4619,+0.1913) r=0.50 th=0.125*pi cellid=5
13 [0] p 0 (+0.1913,+0.4619) r=0.50 th=0.375*pi cellid=7
[all …]
/petsc/src/ksp/pc/impls/tfs/
H A Divec.c229 PetscInt *pi, *pj, temp; in PCTFS_ivec_sort() local
243 pi = ar + 1; in PCTFS_ivec_sort()
247 SWAP(*(ar + (size >> 1)), *pi); in PCTFS_ivec_sort()
251 if (*pi > *pj) SWAP(*pi, *pj); in PCTFS_ivec_sort()
253 else if (*pi > *ar) SWAP(*(ar), *(ar + 1)); in PCTFS_ivec_sort()
259 do pi++; in PCTFS_ivec_sort()
260 while (*pi < *ar); in PCTFS_ivec_sort()
265 if (pj < pi) break; in PCTFS_ivec_sort()
268 SWAP(*pi, *pj); in PCTFS_ivec_sort()
278 if ((*top_s = size - ((PetscInt)(pi - ar)))) { in PCTFS_ivec_sort()
[all …]
/petsc/src/sys/tests/
H A Dex47-options.yaml6 $$pi: &pi 3.14 anchor
11 real: *pi
37 - *pi
/petsc/share/petsc/matlab/
H A Dlaplacian.m230 a1 = pi/2/(u(1)+1);
233 a1 = pi/2/u(1);
237 a1 = pi/4/(u(1)+0.5);
240 a1 = pi/u(1);
248 a2 = pi/2/(u(2)+1);
251 a2 = pi/2/u(2);
255 a2 = pi/4/(u(2)+0.5);
258 a2 = pi/u(2);
268 a3 = pi/2/(u(3)+1);
271 a3 = pi/2/u(3);
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dmatptap.c68 PetscInt *pti, *ptj, *ptJ, *ai = a->i, *aj = a->j, *ajj, *pi = p->i, *pj = p->j, *pjj; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy() local
94 …PetscCall(PetscFreeSpaceGet(PetscRealIntMultTruncate(fill, PetscIntSumTruncate(ai[am], pi[pm])), &… in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
118 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
119 pjj = pj + pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
169 afill = (PetscReal)ci[pn] / (ai[am] + pi[pm] + 1.e-5); in MatPtAPSymbolic_SeqAIJ_SeqAIJ_SparseAxpy()
193 PetscInt *ai = a->i, *aj = a->j, *apj, *apjdense, *pi = p->i, *pj = p->j, *pJ = p->j, *pjj; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy() local
228 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
229 pjj = pj + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
230 paj = pa + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
247 pnzi = pi[i + 1] - pi[i]; in MatPtAPNumeric_SeqAIJ_SeqAIJ_SparseAxpy()
/petsc/src/ts/tutorials/
H A Dex74.c121 PetscReal x, dx, a = ctxt->a, pi = PETSC_PI; in ExactSolution() local
131 uarr[i - is] = PetscExpScalar(-4.0 * pi * pi * a * t) * PetscSinScalar(2 * pi * x); in ExactSolution()
134 uarr[i - is] = PetscSinScalar(2 * pi * (x - a * t)); in ExactSolution()
/petsc/src/dm/tests/
H A Dex50.c57 PetscInt i, j, k, pi, pj, pk; in test_3d() local
111 for (pi = 0; pi < dof; pi++) { in test_3d()
112 xyz[0] = nodes[pi]; in test_3d()
/petsc/src/ksp/ksp/tutorials/
H A Dex66.c67 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
77 pi = 4 * PetscAtanReal(1.0); in ComputeRHS()
84 …s + xm; i++) array[j][i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(t… in ComputeRHS()
H A Dex50.c71 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
81 pi = 4 * atan(1.0); in ComputeRHS()
88 …s + xm; i++) array[j][i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(t… in ComputeRHS()
H A Dex74.c256 PetscReal x, dx, a = ctxt->a, pi = PETSC_PI; in ExactSolution() local
266 uarr[i - is] = PetscExpScalar(-4.0 * pi * pi * a * t) * PetscSinScalar(2 * pi * x); in ExactSolution()
269 uarr[i - is] = PetscSinScalar(2 * pi * (x - a * t)); in ExactSolution()
/petsc/src/tao/unconstrained/tutorials/
H A Dadjointreads.m114 % plot(xg,2.0*mu*pi*sin(pi*x).*exp(-pi^2*t*mu)./(2.0+exp(-pi^2*t*mu).*cos(pi*x)));
/petsc/src/sys/objects/
H A Dfinit.c66 PetscReal pi = PETSC_PI; in PetscInitializeFortran() local
74 petscsetmoduleblocknumeric_(&pi, &maxreal, &minreal, &eps, &seps, &small, &pinf, &pninf); in PetscInitializeFortran()
/petsc/src/snes/tutorials/
H A Dex74fcomd.h1 integer probnum, ihod, mx, my, neq, ientro, gorder parameter(neq = 3) double precision theta, pi, t…
22 / func / pi,
/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
92 for (pi = 0; pi < procs[0]; ++pi) { in PetscPartitionerPartition_Simple_Grid()
93 const PetscInt pid = ((nid * procs[2] + pk) * procs[1] + pj) * procs[0] + pi; in PetscPartitionerPartition_Simple_Grid()
99 …s[1] + ((nj * procs[1] + pj) * pcells[1] + cj)) * cells[0] + (ni * procs[0] + pi) * pcells[0] + ci; in PetscPartitionerPartition_Simple_Grid()
/petsc/doc/tutorials/physics/
H A Dguide_to_stokes.md156 u = \begin{pmatrix} \sin(\pi x) + \sin(\pi y) \\ -\pi \cos(\pi x) y \end{pmatrix} \quad \mathrm{or}…
157 …\begin{pmatrix} 2 \sin(\pi x) + \sin(\pi y) + \sin(\pi z) \\ -\pi \cos(\pi x) y \\ -\pi \cos(\pi x…
167 p = \sin(2 \pi x) + \sin(2 \pi y) \quad \mathrm{or} \quad \sin(2 \pi x) + \sin(2 \pi y) + \sin(2 \p…
177 … = \begin{pmatrix} 2 \pi \cos(2 \pi x) + \mu \pi^2 \sin(\pi x) + \mu \pi^2 \sin(\pi y) \\ 2 \pi \c…
178 …pi \cos(2 \pi x) + 2\mu \pi^2 \sin(\pi x) + \mu \pi^2 \sin(\pi y) + \mu \pi^2 \sin(\pi z) \\ 2 \pi…
/petsc/src/ts/characteristic/interface/
H A Dcharacteristic.c716 PetscMPIInt **procs, pi, pj, pim, pip, pjm, pjp, PIi, PJi; in DMDAGetNeighborsRank() local
733 for (pi = 0; pi < PI; pi++) { in DMDAGetNeighborsRank()
734 procs[pj][pi] = rank; in DMDAGetNeighborsRank()
739 pi = neighbors[0] % PI; in DMDAGetNeighborsRank()
741 pim = pi - 1; in DMDAGetNeighborsRank()
743 pip = (pi + 1) % PIi; in DMDAGetNeighborsRank()
750 neighbors[3] = procs[pjp][pi]; in DMDAGetNeighborsRank()
754 neighbors[7] = procs[pjm][pi]; in DMDAGetNeighborsRank()
758 if (pi == 0) neighbors[1] = neighbors[2] = neighbors[8] = neighbors[0]; in DMDAGetNeighborsRank()
759 if (pi == PI - 1) neighbors[4] = neighbors[5] = neighbors[6] = neighbors[0]; in DMDAGetNeighborsRank()
/petsc/src/ts/event/tests/
H A Dex1sin.c24 PetscReal pi; member
56 ctx.pi = PetscAcosReal(-1.0); in main()
151 if (Ctx->rank == 0) gval[n++] = PetscSinReal(Ctx->pi * t); in EventFunction()
H A Dex4.c26 PetscReal pi; member
61 ctx.pi = PetscAcosReal(-1.0); in main()
198 …if (Ctx->F == 0 || Ctx->F == -1) gval[n++] = PetscSinReal(Ctx->pi * t) / Ctx->pi; // FUNC-0, roots… in EventFunction()
H A Dex5.c28 PetscReal pi; member
64 ctx.pi = PetscAcosReal(-1.0); in main()
200 …if (Ctx->rank == 0) gval[n++] = (t - 5) * (t - 5) * PetscSinReal(Ctx->pi * t); // sin-event -- on … in EventFunction()
201 …if (Ctx->rank == Ctx->size - 1) gval[n++] = 0.5 * PetscCosReal(Ctx->pi * t); // cos-event -- on … in EventFunction()
H A Dex3.c32 PetscReal pi; member
69 ctx.pi = PetscAcosReal(-1.0); in main()
196 if (Ctx->rank == 1 % Ctx->size) gval[n++] = Ctx->V * PetscSinReal(Ctx->pi * t); in EventFunction()
/petsc/src/tao/bound/tutorials/
H A Djbearing2.c176 PetscReal two = 2.0, pi = 4.0 * atan(1.0); in ComputeB() local
184 hx = two * pi / (nx + 1.0); in ComputeB()
213 PetscReal one = 1.0, two = 2.0, six = 6.0, pi = 4.0 * atan(1.0); in FormFunctionGradient() local
225 hx = two * pi / (nx + 1.0); in FormFunctionGradient()
328 PetscReal one = 1.0, two = 2.0, six = 6.0, pi = 4.0 * atan(1.0); in FormHessian() local
338 hx = two * pi / (nx + 1.0); in FormHessian()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiptap.c236 …PetscInt *api, *apj, *Jptr, apnz, *prmap = p->garray, con, j, Crmax, *aj, *ai, *pi,… in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable() local
302 pi = p_loc->i; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
307 pnz = pi[row + 1] - pi[row]; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
308 Jptr = p_loc->j + pi[row]; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
315 pi = p_oth->i; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
320 pnz = pi[row + 1] - pi[row]; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
321 Jptr = p_oth->j + pi[row]; in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
577 PetscInt *ai, nzi, j, *aj, row, col, *pi, *pj, pnz, nzpi, *p_othcols, k; in MatPtAPSymbolicComputeOneRowOfAP_private() local
610 pi = p_oth->i; in MatPtAPSymbolicComputeOneRowOfAP_private()
617 pnz = pi[row + 1] - pi[row]; in MatPtAPSymbolicComputeOneRowOfAP_private()
[all …]
/petsc/src/ksp/ksp/tutorials/amrex/
H A DMyTest.cxx191 const Real pi = 4.0 * std::atan(1.0); in initData() local
199 fab(i, j, k) = std::sin(rx * 2. * pi + 43.5) * std::sin(ry * 2. * pi + 89.); in initData()
/petsc/src/dm/impls/swarm/tests/
H A Dex7.c286 for (int pi = 0, pp = 0; pi < Np2[0]; pi++) { in go() local
288 xx_t[tid][pp] = lo[0] + hp[0] / 2. + pi * hp[0]; in go()
/petsc/src/mat/impls/maij/
H A Dmaij.c591 const PetscInt *ai = a->i, *aj = a->j, *pi = p->i, *pj = p->j, *pJ, *pjj; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() local
614 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
615 pjj = pj + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
616 paj = pa + pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
636 poffset = pi[prow]; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
637 pnzi = pi[prow + 1] - poffset; in MatPtAPNumeric_SeqAIJ_SeqMAIJ()
679 const PetscInt *pi = p->i, *pj = p->j, *pjj, *ai = a->i, *aj = a->j, *ajj; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() local
731 pnzj = pi[prow + 1] - pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
732 pjj = pj + pi[prow]; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()

12