Home
last modified time | relevance | path

Searched refs:r2 (Results 1 – 25 of 32) sorted by relevance

12

/petsc/src/mat/tests/
H A Dex75.c10 PetscReal r1, r2, rnorm, tol = PETSC_SQRT_MACHINE_EPSILON; in main() local
169 PetscCall(MatNorm(sA, NORM_FROBENIUS, &r2)); in main()
170 rnorm = PetscAbsReal(r1 - r2) / r2; in main()
171 …m_FROBENIUS(), Anorm=%16.14e, sAnorm=%16.14e bs=%" PetscInt_FMT "\n", (double)r1, (double)r2, bs)); in main()
173 PetscCall(MatNorm(sA, NORM_INFINITY, &r2)); in main()
174 rnorm = PetscAbsReal(r1 - r2) / r2; in main()
175 …rm_INFINITY(), Anorm=%16.14e, sAnorm=%16.14e bs=%" PetscInt_FMT "\n", (double)r1, (double)r2, bs)); in main()
177 PetscCall(MatNorm(sA, NORM_1, &r2)); in main()
178 rnorm = PetscAbsReal(r1 - r2) / r2; in main()
179 …: MatNorm_1(), Anorm=%16.14e, sAnorm=%16.14e bs=%" PetscInt_FMT "\n", (double)r1, (double)r2, bs)); in main()
[all …]
H A Dex129.c186 PetscScalar *v, *v_neighbor, Hx, Hy, Hz, HxHydHz, HyHzdHx, HxHzdHy, r1, r2; in ComputeMatrix() local
218 PetscCall(PetscRandomGetValue(rand, &r2)); in ComputeMatrix()
221 v_neighbor[k3] = r2; in ComputeMatrix()
/petsc/src/sys/tests/
H A Dex36f.F9010 PetscRandom r, q, r2
17 PetscCallA(PetscRandomCreate(PETSC_COMM_WORLD, r2, ierr))
27 PetscCallA(PetscObjectCompose(r, 'test', r2, ierr))
29 PetscCheckA(q == r2, PETSC_COMM_SELF, PETSC_ERR_PLIB, 'Object compose/query failed')
32 PetscCallA(PetscRandomDestroy(r2, ierr))
/petsc/src/ksp/ksp/tests/
H A Dex10.c86 PetscInt ict, nz, base, r1, r2, N, *rowkeep, nstart; in GetElasticityMatrix() local
119 r2 = base + i2 * shiftx + j2 * shifty + k2 * shiftz; in GetElasticityMatrix()
120 PetscCall(AddElement(mat, r1, r2, K, h1, h2)); in GetElasticityMatrix()
168 PetscErrorCode AddElement(Mat mat, PetscInt r1, PetscInt r2, PetscReal **K, PetscInt h1, PetscInt h… in AddElement() argument
181 col = r2 + l2; in AddElement()
184 row = r2 + l2; in AddElement()
201 PetscReal r2[20] = {-1.0, 0.0, 1.0, -1.0, 1.0, -1.0, 0.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 0.0, 1… variable
372 …N[i][j] = c * (1.0 + r2[i] * rst[0][j]) * (1.0 + s2[i] * rst[1][j]) * (1.0 + t2[i] * rst[2… in paulsetup20()
373 …part_N[0][i][j] = c * r2[i] * (1 + s2[i] * rst[1][j]) * (1 + t2[i] * rst[2][j]) * (-1.0 + 2.0 * r2… in paulsetup20()
374 …part_N[1][i][j] = c * s2[i] * (1 + r2[i] * rst[0][j]) * (1 + t2[i] * rst[2][j]) * (-1.0 + r2[i] * … in paulsetup20()
[all …]
/petsc/src/ksp/ksp/impls/cg/stcg/
H A Dstcg.c20 PetscReal rr, r2, step; in KSPCGSolve_STCG()
37 r2 = cg->radius * cg->radius; in KSPCGSolve_STCG()
87 if (r2 >= rr) { in KSPCGSolve_STCG()
91 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
122 if (r2 >= rr) { in KSPCGSolve_STCG()
126 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
200 if (r2 >= rr) { in KSPCGSolve_STCG()
204 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
259 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_STCG()
275 if (r2 >= rr) { in KSPCGSolve_STCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/nash/
H A Dnash.c21 PetscReal rr, r2, step; in KSPCGSolve_NASH()
40 r2 = cg->radius * cg->radius; in KSPCGSolve_NASH()
90 if (r2 >= rr) { in KSPCGSolve_NASH()
94 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
125 if (r2 >= rr) { in KSPCGSolve_NASH()
129 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
203 if (r2 >= rr) { in KSPCGSolve_NASH()
207 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
262 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_NASH()
278 if (r2 >= rr) { in KSPCGSolve_NASH()
[all …]
/petsc/src/tao/complementarity/impls/asls/
H A Dasils.c60 asls->r2 = NULL; in TaoSetUp_ASILS()
99 PetscCall(VecDestroy(&ssls->r2)); in TaoDestroy_ASILS()
183 PetscCall(TaoVecGetSubVec(asls->da, asls->fixed, tao->subset_type, 1.0, &asls->r2)); in TaoSolve_ASILS()
184 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r2)); in TaoSolve_ASILS()
194 PetscCall(TaoVecGetSubVec(asls->ff, asls->free, tao->subset_type, 0.0, &asls->r2)); in TaoSolve_ASILS()
197 PetscCall(VecPointwiseDivide(asls->r2, asls->r2, asls->r3)); in TaoSolve_ASILS()
209 PetscCall(VecAXPY(asls->r2, -1.0, asls->r3)); in TaoSolve_ASILS()
229 PetscCall(KSPSolve(tao->ksp, asls->r2, asls->dxfree)); in TaoSolve_ASILS()
297 asls->r2 = NULL; in TaoCreate_ASILS()
H A Dasfls.c60 asls->r2 = NULL; in TaoSetUp_ASFLS()
99 PetscCall(VecDestroy(&ssls->r2)); in TaoDestroy_ASFLS()
186 PetscCall(TaoVecGetSubVec(asls->da, asls->fixed, tao->subset_type, 1.0, &asls->r2)); in TaoSolve_ASFLS()
187 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r2)); in TaoSolve_ASFLS()
197 PetscCall(TaoVecGetSubVec(asls->ff, asls->free, tao->subset_type, 0.0, &asls->r2)); in TaoSolve_ASFLS()
200 PetscCall(VecPointwiseDivide(asls->r2, asls->r2, asls->r3)); in TaoSolve_ASFLS()
212 PetscCall(VecAXPY(asls->r2, -1.0, asls->r3)); in TaoSolve_ASFLS()
231 PetscCall(KSPSolve(tao->ksp, asls->r2, asls->dxfree)); in TaoSolve_ASFLS()
302 asls->r2 = NULL; in TaoCreate_ASFLS()
/petsc/src/snes/tutorials/
H A Dex74ftube.h1 double precision p1, p2, p3, p4, r1, r2, r3, r4 double precision e1, e2, e3, e4, v1, v2, v3, v4 dou…
6 …p2, p3, p4, r1, r2, r3, r4, e1, e2, e3, e4, v1, v2, v3, v4, ru1, ru2, ru3, ru4, et1, et2, et3, et4…
H A Dex12.c118 const PetscReal r2 = PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5); in ball_u_2d() local
119 const PetscReal xi = alpha * (radius2 - r2); in ball_u_2d()
143 PetscReal r2, xi; in f0_ball_u() local
145 for (d = 0, r2 = 0.0; d < dim; ++d) r2 += PetscSqr(x[d] - 0.5); in f0_ball_u()
146 xi = alpha * (radius2 - r2); in f0_ball_u()
147 …f0[0] = (-2.0 * dim * alpha - 8.0 * PetscSqr(alpha) * r2 * PetscTanhReal(xi)) * PetscSqr(1.0 / Pet… in f0_ball_u()
372 const PetscReal r2 = PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5) + PetscSqr(x[2] - 0.5); in ball_u_3d() local
373 const PetscReal xi = alpha * (radius2 - r2); in ball_u_3d()
/petsc/include/
H A Dpetscsys.h669 …r2) PetscMallocA(2, PETSC_FALSE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((size_t)((size_t)m1) * … argument
692 …r2) PetscMallocA(2, PETSC_TRUE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((size_t)((size_t)m1) * s… argument
717 #define PetscMalloc3(m1, r1, m2, r2, m3, r3) \ argument
718 …(size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((size_t)m2) * sizeof(**(r2))), (r2), ((size_t)((siz…
743 #define PetscCalloc3(m1, r1, m2, r2, m3, r3) \ argument
744 …(size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((size_t)m2) * sizeof(**(r2))), (r2), ((size_t)((siz…
771 #define PetscMalloc4(m1, r1, m2, r2, m3, r3, m4, r4) \ argument
772 …(size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((size_t)m2) * sizeof(**(r2))), (r2), ((size_t)((siz…
799 #define PetscCalloc4(m1, r1, m2, r2, m3, r3, m4, r4) \ argument
800 …(size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((size_t)m2) * sizeof(**(r2))), (r2), ((size_t)((siz…
[all …]
/petsc/src/ksp/ksp/impls/cg/gltr/
H A Dgltr.c71 PetscReal rr, r2, piv, step; in KSPCGSolve_GLTR()
90 r2 = cg->radius * cg->radius; in KSPCGSolve_GLTR()
141 if (r2 >= rr) { in KSPCGSolve_GLTR()
145 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
172 if (r2 >= rr) { in KSPCGSolve_GLTR()
176 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
241 if (r2 >= rr) { in KSPCGSolve_GLTR()
245 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
298 step = PetscSqrtReal(r2 / norm_p); in KSPCGSolve_GLTR()
308 if (r2 >= rr) { in KSPCGSolve_GLTR()
[all …]
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h73 …PetscReal l, s, r = x[0], z = x[1], i1func, i2func, i3func, ks, es, pi4pow, sqrt_1s, r2, rp2, r2pr… in LandauTensor2D() local
76 r2 = PetscSqr(r); in LandauTensor2D()
80 r2prp2 = r2 + rp2; in LandauTensor2D()
81 l = r2 + rp2 + zmzp2; in LandauTensor2D()
/petsc/src/mat/utils/
H A Dmultequal.c13 PetscReal r1, r2, tol = PETSC_SQRT_MACHINE_EPSILON; in MatMultEqual_Private() local
92 PetscCall(VecNorm(s2, NORM_INFINITY, &r2)); in MatMultEqual_Private()
93 if (r2 < tol) { in MatMultEqual_Private()
98 r1 /= r2; in MatMultEqual_Private()
120 PetscReal r1, r2, tol = PETSC_SQRT_MACHINE_EPSILON; in MatMatMultEqual_Private() local
187 PetscCall(VecNorm(s2, NORM_INFINITY, &r2)); in MatMatMultEqual_Private()
188 if (r2 < tol) { in MatMatMultEqual_Private()
193 r1 /= r2; in MatMatMultEqual_Private()
/petsc/src/dm/field/tutorials/
H A Dex1.c149 PetscReal r2 = 0.; in radiusSquared() local
152 for (i = 0; i < dim; i++) r2 += PetscSqr(x[i]); in radiusSquared()
153 for (i = 0; i < Nf; i++) u[i] = (i + 1) * r2; in radiusSquared()
174 PetscReal r2 = 0.; in TestShellEvaluate() local
176 for (j = 0; j < dim; j++) r2 += PetscSqr(PetscRealPart(x[i * dim + j])); in TestShellEvaluate()
181 ((PetscScalar *)B)[i * Nc + j] = m * r2; in TestShellEvaluate()
183 ((PetscReal *)B)[i * Nc + j] = m * r2; in TestShellEvaluate()
/petsc/src/mat/impls/aij/seq/
H A Dmatrart.c106 PetscScalar r1, r2, r3, r4; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense() local
151 r1 = r2 = r3 = r4 = 0.0; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
157 r2 += (*aa) * b2[*aj]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
162 c[am + i] = r2; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
174 r1 = r2 = r3 = r4 = 0.0; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
180 r2 += (*ra) * c2[*rj]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
185 d[colrm + rm + i] = r2; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
/petsc/src/dm/impls/plex/kokkos/
H A Dplexlocalizationletkf.kokkos.cxx23 const PetscReal r2 = r * r; in GaspariCohn() local
24 const PetscReal r3 = r2 * r; in GaspariCohn()
30 return -0.25 * r5 + 0.5 * r4 + 0.625 * r3 - (5.0 / 3.0) * r2 + 1.0; in GaspariCohn()
33 …return (1.0 / 12.0) * r5 - 0.5 * r4 + 0.625 * r3 + (5.0 / 3.0) * r2 - 5.0 * r + 4.0 - (2.0 / 3.0) … in GaspariCohn()
/petsc/include/petsc/private/valgrind/
H A Dvalgrind.h554 unsigned long int r2; /* what tocptr do we need? */ member
609 _zzq_orig->r2 = __addr; \
631 unsigned long int r2; /* what tocptr do we need? */ member
686 _zzq_orig->r2 = __addr; \
2745 _argvec[1] = (unsigned long)_orig.r2; \
2771 _argvec[1] = (unsigned long)_orig.r2; \
2799 _argvec[1] = (unsigned long)_orig.r2; \
2829 _argvec[1] = (unsigned long)_orig.r2; \
2861 _argvec[1] = (unsigned long)_orig.r2; \
2895 _argvec[1] = (unsigned long)_orig.r2; \
[all …]
/petsc/src/tao/complementarity/impls/ssls/
H A Dssls.h52 Vec r1, r2, r3, w; /* work vectors */ member
/petsc/src/tao/constrained/impls/ipm/
H A Dipm.c677 PetscInt r2, r3, r4; in IPMUpdateAi() local
683 r2 = ipmP->mi; in IPMUpdateAi()
684 r3 = r2 + ipmP->nxlb; in IPMUpdateAi()
700 for (i = r2; i < r4; i++) nonzeros[i] = 1; in IPMUpdateAi()
733 newrow = i + r2; in IPMUpdateAi()
807 PetscInt r1, r2, r3; in IPMUpdateK() local
827 r2 = r1 + ipmP->me; in IPMUpdateK()
829 r3 = c3 = r2 + ipmP->nb; in IPMUpdateK()
844 } else if (i < r2) { in IPMUpdateK()
866 } else if (i < r2) { in IPMUpdateK()
[all …]
/petsc/src/tao/pde_constrained/impls/lcl/
H A Dlcl.c179 PetscReal adec, r2, rGL_U, rWU; in TaoSolve_LCL() local
306 PetscCall(VecNorm(lclP->r, NORM_2, &r2)); in TaoSolve_LCL()
307 r2 = PetscPowScalar(r2, 2.0 + lclP->eps2); in TaoSolve_LCL()
308 adec = lclP->eps1 * r2; in TaoSolve_LCL()
318 PetscCall(VecNorm(lclP->r, NORM_2, &r2)); in TaoSolve_LCL()
319 r2 = PetscPowScalar(r2, 2.0 + lclP->eps2); in TaoSolve_LCL()
321 adec = lclP->eps1 * r2; in TaoSolve_LCL()
/petsc/src/dm/tutorials/
H A Dswarm_ex3.c184 PetscReal px, py, rx, ry, r2; in ex3_1() local
192 r2 = px * px + py * py; in ex3_1()
193 if (r2 < 0.75 * 0.75) { in ex3_1()
H A Dswarm_ex1.c300 PetscReal *array_x, *array_y, r2; in collect_zone() local
310 r2 = zone->radius * zone->radius; in collect_zone()
317 if (sep2 < r2) p2collect++; in collect_zone()
327 if (sep2 < r2) { in collect_zone()
/petsc/src/ksp/ksp/impls/cg/
H A Dcg.c124 PetscReal r2, norm_p, norm_d, dMp; in KSPSolve_CG() local
143 r2 = PetscSqr(cg->radius); in KSPSolve_CG()
267 a = PetscSqrtReal(r2 / norm_p); in KSPSolve_CG()
273 a = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
293 if (norm_dp1 > r2) { in KSPSolve_CG()
297 dp = (PetscSqrtReal(dMp * dMp + norm_p * (r2 - norm_d)) - dMp) / norm_p; in KSPSolve_CG()
/petsc/src/snes/tests/
H A Dex17.c337 Vec x1, x2, r1, r2; in block_system() local
358 PetscCall(VecDuplicate(x2, &r2)); in block_system()
372 bx[1] = r2; in block_system()
375 PetscCall(VecDestroy(&r2)); in block_system()

12