Home
last modified time | relevance | path

Searched refs:r1 (Results 1 – 19 of 19) sorted by relevance

/petsc/src/mat/tests/
H A Dex75.c10 PetscReal r1, r2, rnorm, tol = PETSC_SQRT_MACHINE_EPSILON; in main() local
168 PetscCall(MatNorm(A, NORM_FROBENIUS, &r1)); 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()
172 PetscCall(MatNorm(A, NORM_INFINITY, &r1)); 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()
176 PetscCall(MatNorm(A, NORM_1, &r1)); 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
217 PetscCall(PetscRandomGetValue(rand, &r1)); in ComputeMatrix()
220 v[k3] = r1; in ComputeMatrix()
/petsc/src/tao/complementarity/impls/asls/
H A Dasils.c59 asls->r1 = NULL; in TaoSetUp_ASILS()
98 PetscCall(VecDestroy(&ssls->r1)); in TaoDestroy_ASILS()
182 PetscCall(TaoVecGetSubVec(asls->ff, asls->fixed, tao->subset_type, 0.0, &asls->r1)); in TaoSolve_ASILS()
184 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r2)); in TaoSolve_ASILS()
186 PetscCall(VecISAXPY(tao->stepdirection, asls->fixed, 1.0, asls->r1)); in TaoSolve_ASILS()
193 PetscCall(TaoVecGetSubVec(asls->da, asls->free, tao->subset_type, 0.0, &asls->r1)); in TaoSolve_ASILS()
196 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r3)); in TaoSolve_ASILS()
221 PetscCall(MatDiagonalSet(asls->J_sub, asls->r1, ADD_VALUES)); in TaoSolve_ASILS()
296 asls->r1 = NULL; in TaoCreate_ASILS()
H A Dasfls.c59 asls->r1 = NULL; in TaoSetUp_ASFLS()
98 PetscCall(VecDestroy(&ssls->r1)); in TaoDestroy_ASFLS()
185 PetscCall(TaoVecGetSubVec(asls->ff, asls->fixed, tao->subset_type, 0.0, &asls->r1)); in TaoSolve_ASFLS()
187 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r2)); in TaoSolve_ASFLS()
189 PetscCall(VecISAXPY(tao->stepdirection, asls->fixed, 1.0, asls->r1)); in TaoSolve_ASFLS()
196 PetscCall(TaoVecGetSubVec(asls->da, asls->free, tao->subset_type, 0.0, &asls->r1)); in TaoSolve_ASFLS()
199 PetscCall(VecPointwiseDivide(asls->r1, asls->r1, asls->r3)); in TaoSolve_ASFLS()
223 PetscCall(MatDiagonalSet(asls->J_sub, asls->r1, ADD_VALUES)); in TaoSolve_ASFLS()
301 asls->r1 = NULL; in TaoCreate_ASFLS()
/petsc/include/
H A Dpetscsys.h622 …PetscMalloc1(m1, r1) PetscMallocA(1, PETSC_FALSE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((size_… argument
646 … PetscCalloc1(m1, r1) PetscMallocA(1, PETSC_TRUE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((size_… argument
669 …cMalloc2(m1, r1, m2, r2) PetscMallocA(2, PETSC_FALSE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((s… argument
692 …scCalloc2(m1, r1, m2, r2) PetscMallocA(2, PETSC_TRUE, __LINE__, PETSC_FUNCTION_NAME, __FILE__, ((s… argument
717 #define PetscMalloc3(m1, r1, m2, r2, m3, r3) \ argument
718 …NE__, PETSC_FUNCTION_NAME, __FILE__, ((size_t)((size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((siz…
743 #define PetscCalloc3(m1, r1, m2, r2, m3, r3) \ argument
744 …NE__, PETSC_FUNCTION_NAME, __FILE__, ((size_t)((size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((siz…
771 #define PetscMalloc4(m1, r1, m2, r2, m3, r3, m4, r4) \ argument
772 …NE__, PETSC_FUNCTION_NAME, __FILE__, ((size_t)((size_t)m1) * sizeof(**(r1))), (r1), ((size_t)((siz…
[all …]
/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()
156 r1 += (*aa) * b1[*aj]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
161 c[i] = r1; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
174 r1 = r2 = r3 = r4 = 0.0; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
179 r1 += (*ra) * c[*rj]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
184 d[colrm + i] = r1; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
192 r1 = 0.0; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
196 for (j = 0; j < n; j++) r1 += (*aa++) * b1[*aj++]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
197 c[i] = r1; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
[all …]
H A Dmatmatmult.c1498 PetscScalar *c, r1, r2, r3, r4, *c1, *c2, *c3, *c4; in MatMatMultNumericAdd_SeqAIJ_SeqDense() local
1530 r1 = r2 = r3 = r4 = 0.0; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1537 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1543 c1[i] += r1; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1548 c1[i] = r1; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1571 r1 = 0.0; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1575 for (j = 0; j < n; j++) r1 += aa[j] * b1[aj[j]]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1576 if (add) c1[i] += r1; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1577 else c1[i] = r1; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1581 r1 = r2 = 0.0; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
[all …]
/petsc/src/mat/utils/
H A Dmultequal.c13 PetscReal r1, r2, tol = PETSC_SQRT_MACHINE_EPSILON; in MatMultEqual_Private() local
94 PetscCall(VecNorm(s1, NORM_INFINITY, &r1)); in MatMultEqual_Private()
97 PetscCall(VecNorm(s2, NORM_INFINITY, &r1)); in MatMultEqual_Private()
98 r1 /= r2; in MatMultEqual_Private()
100 if (r1 > tol) { in MatMultEqual_Private()
102 PetscCall(PetscInfo(A, "Error: %" PetscInt_FMT "-th %s() %g\n", k, sop, (double)r1)); in MatMultEqual_Private()
120 PetscReal r1, r2, tol = PETSC_SQRT_MACHINE_EPSILON; in MatMatMultEqual_Private() local
189 PetscCall(VecNorm(s3, NORM_INFINITY, &r1)); in MatMatMultEqual_Private()
192 PetscCall(VecNorm(s2, NORM_INFINITY, &r1)); in MatMatMultEqual_Private()
193 r1 /= r2; in MatMatMultEqual_Private()
[all …]
/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…
/petsc/src/tao/complementarity/impls/ssls/
H A Dssls.h52 Vec r1, r2, r3, w; /* work vectors */ member
/petsc/src/ksp/ksp/tests/
H A Dex10.c86 PetscInt ict, nz, base, r1, r2, N, *rowkeep, nstart; in GetElasticityMatrix() local
115 r1 = base + i1 * shiftx + j_1 * shifty + k1 * 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
180 row = r1 + l1; in AddElement()
185 col = r1 + l1; in AddElement()
/petsc/src/snes/tests/
H A Dex17.c337 Vec x1, x2, r1, r2; in block_system() local
355 PetscCall(VecDuplicate(x1, &r1)); in block_system()
371 bx[0] = r1; in block_system()
374 PetscCall(VecDestroy(&r1)); in block_system()
/petsc/src/dm/impls/swarm/
H A Ddata_ex.c259 static void _get_tags(PetscInt counter, PetscMPIInt N, PetscMPIInt r0, PetscMPIInt r1, PetscMPIInt … in _get_tags() argument
263 st = (N * r0 + r1 + N * N * counter) % maxtag; in _get_tags()
264 rt = (N * r1 + r0 + N * N * counter) % maxtag; in _get_tags()
360 PetscMPIInt r1 = d->neighbour_procs[n]; in DMSwarmDataExTopologyFinalize() local
362 _get_tags(d->instance, size, r0, r1, *maxtag, &st, &rt); in DMSwarmDataExTopologyFinalize()
/petsc/src/tao/constrained/impls/ipm/
H A Dipm.c807 PetscInt r1, r2, r3; in IPMUpdateK() local
826 r1 = c1 = ipmP->n; in IPMUpdateK()
827 r2 = r1 + ipmP->me; in IPMUpdateK()
839 if (i < r1) { in IPMUpdateK()
862 if (i < r1) { in IPMUpdateK()
902 row = i + r1; in IPMUpdateK()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c793 PetscBool r1, r2, r3; in MatAssemblyEnd_MPISBAIJ() local
819 r1 = baij->roworiented; in MatAssemblyEnd_MPISBAIJ()
844 baij->roworiented = r1; in MatAssemblyEnd_MPISBAIJ()
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c906 PetscBool r1, r2, r3, all_assembled; in MatAssemblyEnd_MPIBAIJ() local
933 r1 = baij->roworiented; in MatAssemblyEnd_MPIBAIJ()
957 baij->roworiented = r1; in MatAssemblyEnd_MPIBAIJ()
/petsc/share/petsc/datafiles/meshes/
H A Dtestcase3D.cas8381 (defaults/pre-2024r1-spalart-allmaras-kappa? #f)
8383 (defaults/pre-2022r1-sbes-bf-include-sub-function1? #f)
8384 (defaults/pre-2022r1-turb-consistent-specified-stress-bc? #f)
8385 (defaults/pre-2022r1-real-gas-species-diffusion-energy? #f)
8390 (defaults/pre-2021r1-node-weight-sum-threshold? #f)
8391 (defaults/pre-2021r1-symmetry-face-grad-weight? #f)
8392 (defaults/pre-2021r1-swirl-axi-mirror-limiting? #f)
8393 (defaults/pre-2021r1-wfgc-jump? #f)
8394 (defaults/pre-2021r1-delta-form? #f)
/petsc/src/dm/impls/forest/p4est/
H A Dpforest.h108 uint32_t r1 = 15; in DMPforestHash() local
121 k = DMPROT32(k, r1); in DMPforestHash()
/petsc/doc/manual/
H A Dksp.md1131 …> - `-pc_gamg_rank_reduction_factors` \<rn,rn-1,...,r1:int> Set a schedule for MPI rank reduction …