Home
last modified time | relevance | path

Searched full:pi (Results 1 – 25 of 105) sorted by relevance

12345

/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/snes/tutorials/
H A Dex62.c168 u = sin(pi x) + sin(pi y)
169 v = -pi cos(pi x) y
170 p = sin(2 pi x) + sin(2 pi y)
171 …f = <2pi cos(2 pi x) + mu pi^2 sin(pi x) + mu pi^2 sin(pi y), 2pi cos(2 pi y) - mu pi^3 cos(pi x) …
175 … e(u) = (grad u + grad u^T) = / 2pi cos(pi x) pi cos(pi y) + pi^2 sin(pi x) y \
176 … \ pi cos(pi y) + pi^2 sin(pi x) y -2pi cos(pi x) /
177 …div mu e(u) - \nabla p + f = mu <-pi^2 sin(pi x) - pi^2 sin(pi y), pi^3 cos(pi x) y> - <2pi cos(2 …
178 \nabla \cdot u = pi cos(pi x) - pi cos(pi x) = 0
182 u = 2 sin(pi x) + sin(pi y) + sin(pi z)
183 v = -pi cos(pi x) y
[all …]
H A Dex17.c163 u = sin(2 pi x)
164 v = sin(2 pi y) - 2xy
165 Delta <u,v> - f = <-4 pi^2 u, -4 pi^2 v> - <-4 pi^2 sin(2 pi x), -4 pi^2 sin(2 pi y)>
167 u = sin(2 pi x)
168 v = sin(2 pi y) - 2xy
169 w = sin(2 pi z) - 2yz
170 …Delta <u,v,2> - f = <-4 pi^2 u, -4 pi^2 v, -4 pi^2 w> - <-4 pi^2 sin(2 pi x), -4 pi^2 sin(2 pi y),…
179 u = sin(2 pi x)
180 v = sin(2 pi y) - 2xy
181 \varepsilon = / 2 pi cos(2 pi x) -y \
[all …]
H A Dex36.c268 static PetscErrorCode DestroyCoarseProjection(Mat Pi) in DestroyCoarseProjection() argument
273 PetscCall(MatShellGetContext(Pi, &ctx)); in DestroyCoarseProjection()
282 PetscCall(MatShellSetContext(Pi, NULL)); in DestroyCoarseProjection()
286 static PetscErrorCode CoarseProjection(Mat Pi, Vec x, Vec y) in CoarseProjection() argument
291 PetscCall(MatShellGetContext(Pi, &ctx)); in CoarseProjection()
306 static PetscErrorCode CreateCoarseProjection(DM dmc, DM dmf, Mat *Pi) in CreateCoarseProjection() argument
326 PetscCall(MatCreateShell(PetscObjectComm((PetscObject)dmc), m, n, M, N, ctx, Pi)); in CreateCoarseProjection()
327 PetscCall(MatShellSetOperation(*Pi, MATOP_DESTROY, (PetscErrorCodeFn *)DestroyCoarseProjection)); in CreateCoarseProjection()
328 PetscCall(MatShellSetOperation(*Pi, MATOP_MULT, (PetscErrorCodeFn *)CoarseProjection)); in CreateCoarseProjection()
334 Mat Pi; /* The L_2 stable projection to the DG coarse space */ member
[all …]
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/ts/tutorials/
H A Dex45.c109 u = 2\pi^2 t + cos(\pi x) + cos(\pi y)
110 F(u) = 2\pi^2 - \pi^2 (cos(\pi x) + cos(\pi y)) + \pi^2 (cos(\pi x) + cos(\pi y)) - 2\pi^2 = 0
113 u = 3\pi^2 t + cos(\pi x) + cos(\pi y) + cos(\pi z)
114 …F(u) = 3\pi^2 - \pi^2 (cos(\pi x) + cos(\pi y) + cos(\pi z)) + \pi^2 (cos(\pi x) + cos(\pi y) + co…
140 u = pi^2 cos(t) + cos(\pi x) + cos(\pi y)
141 u_t = -pi^2 sin(t)
142 \Delta u = -\pi^2 (cos(\pi x) + cos(\pi y))
143 f = pi^2 sin(t) - \pi^2 (cos(\pi x) + cos(\pi y))
144 …F(u) = -\pi^2 sin(t) + \pi^2 (cos(\pi x) + cos(\pi y)) - \pi^2 (cos(\pi x) + cos(\pi y)) + \pi^2 s…
147 u = pi^2 cos(t) + cos(\pi x) + cos(\pi y) + cos(\pi z)
[all …]
H A Dex76.c393 u = 1 - cos(\pi(x - t)) sin(\pi(y - t)) exp(-2 \pi^2 \nu t)
394 v = 1 + sin(\pi(x - t)) cos(\pi(y - t)) exp(-2 \pi^2 \nu t)
395 p = -1/4 [cos(2 \pi(x - t)) + cos(2 \pi(y - t))] exp(-4 \pi^2 \nu t)
397 …f = <\nu \pi^2 exp(-2\nu \pi^2 t) cos(\pi(x-t)) sin(\pi(y-t)), -\nu \pi^2 exp(-2\nu \pi^2 t) sin(\…
398 Q = 3 + sin(\pi(x-y)) exp(-2\nu \pi^2 t)
402 …\nabla \cdot u = \pi sin(\pi(x - t)) sin(\pi(y - t)) exp(-2 \pi^2 \nu t) - \pi sin(\pi(x - t)) sin…
405 …= <-\pi (sin(\pi(x - t)) sin(\pi(y - t)) - cos(\pi(x - t)) cos(\pi(y - t)) - 2\pi cos(\pi(x - t)) …
406 …\pi (sin(\pi(x - t)) sin(\pi(y - t)) - cos(\pi(x - t)) cos(\pi(y - t)) - 2\pi sin(\pi(x - t)) cos(…
407 …+ < \pi (1 - cos(\pi(x - t)) sin(\pi(y - t)) exp(-2 \pi^2 \nu t)) sin(\pi(x - t)) sin(\pi(y - t)) …
408 …\pi (1 - cos(\pi(x - t)) sin(\pi(y - t)) exp(-2 \pi^2 \nu t)) cos(\pi(x - t)) cos(\pi(y - t)) exp(…
[all …]
H A Dex53.c259 u = sin(2 pi x)
260 v = sin(2 pi y) - 2xy
261 \varepsilon = / 2 pi cos(2 pi x) -y \
262 \ -y 2 pi cos(2 pi y) - 2x /
263 Tr(\varepsilon) = div u = 2 pi (cos(2 pi x) + cos(2 pi y)) - 2 x
265 …= \lambda \partial_j 2 pi (cos(2 pi x) + cos(2 pi y)) + 2\mu < -4 pi^2 sin(2 pi x) - 1, -4 pi^2 si…
266 …= \lambda < -4 pi^2 sin(2 pi x) - 2, -4 pi^2 sin(2 pi y) > + \mu < -8 pi^2 sin(2 pi x) - 2, -8 pi^…
269 u = sin(2 pi x)
270 v = sin(2 pi y) - 2xy
271 p = (cos(2 pi x) + cos(2 pi y)) t
[all …]
H A Dex18.c374 u = sin(2 pi x)/4 pi^2
375 v = -y cos(2 pi x)/2 pi
377 f_x = -sin(2 pi x)
378 f_y = 2 pi y cos(2 pi x)
382 -\Delta u + f = <sin(2pi x), -2pi y cos(2pi x)> + <-sin(2pi x), 2pi y cos(2pi x)> = 0
386 \nabla \cdot u = cos(2pi x)/2pi - cos(2pi x)/2pi = 0
398 u = sin(2 pi x) cos(2 pi y)/4 pi^2
399 v = -sin(2 pi y) cos(2 pi x)/4 pi^2
401 f_x = -2sin(2 pi x) cos(2 pi y)
402 f_y = 2sin(2 pi y) cos(2 pi x)
[all …]
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/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
18 $$3: &seq-real [*pi, *eu]
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/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/share/petsc/datafiles/meshes/
H A Dgmsh-B3tet.geo56 Rotate { {0,1,0}, {Cx,Cy,Cz}, +Pi/2 } { Duplicata { Volume {1}; } }
57 Rotate { {1,0,0}, {Cx,Cy,Cz}, -Pi/2 } { Duplicata { Volume {1}; } }
58 Rotate { {0,1,0}, {Cx,Cy,Cz}, -Pi/2 } { Duplicata { Volume {1}; } }
59 Rotate { {1,0,0}, {Cx,Cy,Cz}, +Pi/2 } { Duplicata { Volume {1}; } }
60 Rotate { {1,0,0}, {Cx,Cy,Cz}, +Pi } { Duplicata { Volume {1}; } }
/petsc/src/ksp/ksp/tutorials/
H A Dex66.c6 forcing function f = -cos(m*pi*x)*cos(n*pi*y),
67 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.c6 forcing function f = -cos(m*pi*x)*cos(n*pi*y),
71 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()
/petsc/src/ts/characteristic/interface/
H A Dcharacteristic.c716 PetscMPIInt **procs, pi, pj, pim, pip, pjm, pjp, PIi, PJi; in DMDAGetNeighborsRank() local
717 PetscInt PI, PJ; in DMDAGetNeighborsRank() local
722 …PetscCall(DMDAGetInfo(da, NULL, NULL, NULL, NULL, &PI, &PJ, NULL, NULL, NULL, &bx, &by, NULL, NULL… in DMDAGetNeighborsRank()
723 PetscCall(PetscMPIIntCast(PI, &PIi)); in DMDAGetNeighborsRank()
732 PetscCall(PetscMalloc1(PI, &procs[pj])); in DMDAGetNeighborsRank()
733 for (pi = 0; pi < PI; pi++) { in DMDAGetNeighborsRank()
734 procs[pj][pi] = rank; in DMDAGetNeighborsRank()
739 pi = neighbors[0] % PI; in DMDAGetNeighborsRank()
740 pj = neighbors[0] / PI; in DMDAGetNeighborsRank()
741 pim = pi - 1; in DMDAGetNeighborsRank()
[all …]
/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/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/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/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/config/BuildSystem/config/
H A Dheaders.py187 if self.checkCompile('#include <math.h>\n', 'double pi = M_PI;\n(void)pi'):
189 …lf.checkCompile('#define _USE_MATH_DEFINES 1\n#include <math.h>\n', 'double pi = M_PI;\n(void)pi'):
/petsc/src/dm/dt/tests/
H A Dex1.c52 PetscReal *Pi, *Pj; in CheckQuadrature() local
57 PetscCall(PetscMalloc2(npoints, &Pi, npoints, &Pj)); in CheckQuadrature()
59 PetscCall(PetscDTJacobiEval(npoints, alpha, beta, x, 1, &i, Pi, NULL, NULL)); in CheckQuadrature()
88 for (k = 0; k < npoints; k++) I_quad += w[k] * (Pi[k] * Pj[k]); in CheckQuadrature()
94 PetscCall(PetscFree2(Pi, Pj)); in CheckQuadrature()

12345