Home
last modified time | relevance | path

Searched refs:rr (Results 1 – 25 of 29) sorted by relevance

12

/petsc/src/mat/tests/
H A Dex235.c27 Vec ll, d, rr; in shiftandscale() local
30 PetscCall(MatCreateVecs(A, &ll, &rr)); in shiftandscale()
33 PetscCall(VecSetRandom(rr, NULL)); in shiftandscale()
42 PetscCall(VecSetRandom(rr, NULL)); in shiftandscale()
43 PetscCall(MatDiagonalScale(A, ll, rr)); in shiftandscale()
49 PetscCall(VecSetRandom(rr, NULL)); in shiftandscale()
50 PetscCall(MatDiagonalScale(A, ll, rr)); in shiftandscale()
56 PetscCall(VecDestroy(&rr)); in shiftandscale()
H A Dex23.c305 …PetscInt ri, ci, rr[3] = {0, 1, 0}, cr[4] = {1, 2, 0, 1}, rk[3] = {0, 2, 1}, ck[4] =… in main() local
311 r = (PetscInt *)(ri == 0 ? rr : rk); in main()
/petsc/src/sys/classes/random/interface/
H A Drandomc.c355 PetscRandom rr; in PetscRandomCreate() local
362 …PetscCall(PetscHeaderCreate(rr, PETSC_RANDOM_CLASSID, "PetscRandom", "Random number generator", "S… in PetscRandomCreate()
364 rr->data = NULL; in PetscRandomCreate()
365 rr->low = 0.0; in PetscRandomCreate()
366 rr->width = 1.0; in PetscRandomCreate()
367 rr->iset = PETSC_FALSE; in PetscRandomCreate()
368 rr->seed = 0x12345678 + 76543 * rank; in PetscRandomCreate()
369 PetscCall(PetscRandomSetType(rr, PETSCRANDER48)); in PetscRandomCreate()
370 *r = rr; in PetscRandomCreate()
/petsc/src/ksp/pc/impls/mg/ftn-custom/
H A Dzmgfuncf.c24 MVVVV rr; in pcmgsetresidual_() local
25 if (residual == pcmgresidualdefault_) rr = PCMGResidualDefault; in pcmgsetresidual_()
31 rr = ourresidualfunction; in pcmgsetresidual_()
33 *ierr = PCMGSetResidual(*pc, *l, rr, *mat); in pcmgsetresidual_()
/petsc/src/ksp/ksp/impls/cg/stcg/
H A Dstcg.c20 PetscReal rr, r2, step; in KSPCGSolve_STCG()
66 PetscCall(VecDot(r, r, &rr)); /* rr = r^T r */ in KSPCGSolve_STCG()
67 KSPCheckDot(ksp, rr); in KSPCGSolve_STCG()
87 if (r2 >= rr) { in KSPCGSolve_STCG()
89 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
91 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
122 if (r2 >= rr) { in KSPCGSolve_STCG()
124 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_STCG()
126 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_STCG()
158 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_STCG()
[all …]
/petsc/src/ksp/ksp/impls/cg/nash/
H A Dnash.c21 PetscReal rr, r2, step; in KSPCGSolve_NASH()
69 PetscCall(VecDot(r, r, &rr)); /* rr = r^T r */ in KSPCGSolve_NASH()
70 KSPCheckDot(ksp, rr); in KSPCGSolve_NASH()
90 if (r2 >= rr) { in KSPCGSolve_NASH()
92 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
94 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
125 if (r2 >= rr) { in KSPCGSolve_NASH()
127 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_NASH()
129 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_NASH()
161 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_NASH()
[all …]
/petsc/src/ksp/ksp/impls/cg/gltr/
H A Dgltr.c71 PetscReal rr, r2, piv, step; in KSPCGSolve_GLTR()
122 PetscCall(VecDot(r, r, &rr)); /* rr = r^T r */ in KSPCGSolve_GLTR()
123 KSPCheckDot(ksp, rr); in KSPCGSolve_GLTR()
141 if (r2 >= rr) { in KSPCGSolve_GLTR()
143 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
145 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
172 if (r2 >= rr) { in KSPCGSolve_GLTR()
174 cg->norm_d = PetscSqrtReal(rr); in KSPCGSolve_GLTR()
176 alpha = PetscSqrtReal(r2 / rr); in KSPCGSolve_GLTR()
206 norm_r = PetscSqrtReal(rr); /* norm_r = |r| */ in KSPCGSolve_GLTR()
[all …]
/petsc/src/ts/tutorials/
H A Dex11.h404 …ar *gamr, PetscScalar *rl, PetscScalar *pl, PetscScalar *uxl, PetscScalar *rr, PetscScalar *pr, Pe… in riem1mdt() argument
435 cr = PetscSqrtScalar(*gamr * *pr / *rr); in riem1mdt()
437 wr = *rr * cr; in riem1mdt()
481 *rstarr = *rr * PetscPowScalar(d__1, d__2); in riem1mdt()
504 *rstarr = *rr * PetscPowScalar(d__1, d__2); in riem1mdt()
569 *rstarr = ((*gamr + 1.) * pst + *gamr - 1.) / ((*gamr - 1.) * pst + *gamr + 1.) * *rr; in riem1mdt()
582 …utl, PetscScalar *ubl, PetscScalar *gaml, PetscScalar *rho1l, PetscScalar *rr, PetscScalar *uxr, P… in riemannsolver() argument
592 if (*rl == *rr && *pr == *pl && *uxl == *uxr && *gaml == *gamr) { in riemannsolver()
610 iwave = riem1mdt(gaml, gamr, rl, pl, uxl, rr, pr, uxr, &rstarl, &rstarr, &pstar, &ustar); in riemannsolver()
617 r0 = cvmgm_(rl, rr, &sgn0); in riemannsolver()
[all …]
H A Dex30.c1042 PetscInt cc, rr; in SaveToFile() local
1051 PetscCall(DMGetRefineLevel(dm, &rr)); in SaveToFile()
1055 PetscCall(PetscViewerHDF5WriteAttribute(viewer, NULL, "refinelevel", PETSC_INT, &rr)); in SaveToFile()
1154 PetscInt rr, cc; in LoadFromFile() local
1162 PetscCall(PetscViewerHDF5ReadAttribute(viewer, NULL, "refinelevel", PETSC_INT, NULL, &rr)); in LoadFromFile()
1245 PetscCall(DMSetRefineLevel(dm, rr)); in LoadFromFile()
/petsc/src/ksp/pc/tests/
H A Dex5.c174 PetscErrorCode residual(Mat mat, Vec bb, Vec xx, Vec rr) in residual() argument
184 PetscCall(VecGetArray(rr, &r)); in residual()
191 PetscCall(VecRestoreArray(rr, &r)); in residual()
297 PetscErrorCode restrct(Mat mat, Vec rr, Vec bb) in restrct() argument
304 PetscCall(VecGetSize(rr, &N)); in restrct()
305 PetscCall(VecGetArrayRead(rr, &r)); in restrct()
313 PetscCall(VecRestoreArrayRead(rr, &r)); in restrct()
/petsc/lib/petsc/bin/maint/petsclinter/petsclinter/classes/
H A D_add_line.py43 rl, rr = re_match.group(2).split(',')
/petsc/src/ts/impls/implicit/sundials/
H A Dsundials.c54 Vec rr = cvode->w1, zz = cvode->w2; in TSPSolve_Sundials_Petsc() local
61 PetscCall(VecPlaceArray(rr, r_data)); in TSPSolve_Sundials_Petsc()
66 PetscCall(PCApply(pc, rr, zz)); in TSPSolve_Sundials_Petsc()
67 PetscCall(VecResetArray(rr)); in TSPSolve_Sundials_Petsc()
/petsc/src/mat/impls/baij/seq/baijmkl/
H A Dbaijmkl.c364 static PetscErrorCode MatDiagonalScale_SeqBAIJMKL(Mat A, Vec ll, Vec rr) in MatDiagonalScale_SeqBAIJMKL() argument
367 PetscCall(MatDiagonalScale_SeqBAIJ(A, ll, rr)); in MatDiagonalScale_SeqBAIJMKL()
/petsc/src/mat/impls/sell/mpi/
H A Dmpisell.c797 static PetscErrorCode MatDiagonalScale_MPISELL(Mat mat, Vec ll, Vec rr) in MatDiagonalScale_MPISELL() argument
805 if (rr) { in MatDiagonalScale_MPISELL()
806 PetscCall(VecGetLocalSize(rr, &s1)); in MatDiagonalScale_MPISELL()
809 PetscCall(VecScatterBegin(sell->Mvctx, rr, sell->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPISELL()
817 PetscUseTypeMethod(a, diagonalscale, ll, rr); in MatDiagonalScale_MPISELL()
819 if (rr) { in MatDiagonalScale_MPISELL()
821 PetscCall(VecScatterEnd(sell->Mvctx, rr, sell->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPISELL()
/petsc/src/ksp/pc/impls/gamg/
H A Dgeo.c561 PetscReal rr; in PCGAMGCoarsen_GEO() local
567 PetscCall(PetscRandomGetValueReal(rand, &rr)); in PCGAMGCoarsen_GEO()
568 iSwapIndex = (PetscInt)(rr * nloc); in PCGAMGCoarsen_GEO()
/petsc/src/ksp/ksp/tutorials/
H A Dex42.c1974 PetscMPIInt rr; in solve_stokes_3d_coupled() local
1983 rr = PetscGlobalRank % (px * py); in solve_stokes_3d_coupled()
1984 if (px % 2) jump = (PetscBool)(rr % 2); in solve_stokes_3d_coupled()
1985 else jump = (PetscBool)((rr / px) % 2 ? rr % 2 : !(rr % 2)); in solve_stokes_3d_coupled()
1986 rr = PetscGlobalRank / (px * py); in solve_stokes_3d_coupled()
1987 if (rr % 2) jump = (PetscBool)!jump; in solve_stokes_3d_coupled()
H A Dex70.c1039 PetscReal x_p[2], rr[2]; in SolveTimeDepStokes() local
1042 PetscCall(PetscRandomGetValueReal(r, &rr[0])); in SolveTimeDepStokes()
1043 PetscCall(PetscRandomGetValueReal(r, &rr[1])); in SolveTimeDepStokes()
1044 array_x[2 * p + 0] += rr[0]; in SolveTimeDepStokes()
1045 array_x[2 * p + 1] += rr[1]; in SolveTimeDepStokes()
/petsc/src/mat/impls/aij/seq/kokkos/
H A Daijkok.kokkos.cxx1092 static PetscErrorCode MatDiagonalScale_SeqAIJKokkos(Mat A, Vec ll, Vec rr) in MatDiagonalScale_SeqAIJKokkos() argument
1119 if (rr) { in MatDiagonalScale_SeqAIJKokkos()
1120 PetscCall(VecGetLocalSize(rr, &n)); in MatDiagonalScale_SeqAIJKokkos()
1122 PetscCall(VecGetKokkosView(rr, &rv)); in MatDiagonalScale_SeqAIJKokkos()
1125 PetscCall(VecRestoreKokkosView(rr, &lv)); in MatDiagonalScale_SeqAIJKokkos()
/petsc/src/mat/impls/sell/seq/
H A Dsell.c996 PetscErrorCode MatDiagonalScale_SeqSELL(Mat A, Vec ll, Vec rr) in MatDiagonalScale_SeqSELL() argument
1021 if (rr) { in MatDiagonalScale_SeqSELL()
1022 PetscCall(VecGetLocalSize(rr, &n)); in MatDiagonalScale_SeqSELL()
1024 PetscCall(VecGetArrayRead(rr, &r)); in MatDiagonalScale_SeqSELL()
1034 PetscCall(VecRestoreArrayRead(rr, &r)); in MatDiagonalScale_SeqSELL()
/petsc/src/mat/impls/baij/mpi/
H A Dmpibaij.c1567 static PetscErrorCode MatDiagonalScale_MPIBAIJ(Mat mat, Vec ll, Vec rr) in MatDiagonalScale_MPIBAIJ() argument
1575 if (rr) { in MatDiagonalScale_MPIBAIJ()
1576 PetscCall(VecGetLocalSize(rr, &s1)); in MatDiagonalScale_MPIBAIJ()
1579 PetscCall(VecScatterBegin(baij->Mvctx, rr, baij->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPIBAIJ()
1587 PetscUseTypeMethod(a, diagonalscale, ll, rr); in MatDiagonalScale_MPIBAIJ()
1589 if (rr) { in MatDiagonalScale_MPIBAIJ()
1591 PetscCall(VecScatterEnd(baij->Mvctx, rr, baij->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPIBAIJ()
/petsc/src/mat/impls/sbaij/mpi/
H A Dmpisbaij.c1571 static PetscErrorCode MatDiagonalScale_MPISBAIJ(Mat mat, Vec ll, Vec rr) in MatDiagonalScale_MPISBAIJ() argument
1583 PetscCall(VecGetLocalSize(rr, &nv)); in MatDiagonalScale_MPISBAIJ()
1586 PetscCall(VecScatterBegin(baij->Mvctx, rr, baij->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPISBAIJ()
1592 PetscUseTypeMethod(a, diagonalscale, ll, rr); in MatDiagonalScale_MPISBAIJ()
1595 PetscCall(VecScatterEnd(baij->Mvctx, rr, baij->lvec, INSERT_VALUES, SCATTER_FORWARD)); in MatDiagonalScale_MPISBAIJ()
/petsc/src/mat/impls/is/
H A Dmatis.c1156 Vec ll, rr; in MatDiagonalScale_IS() local
1170 rr = matis->x; in MatDiagonalScale_IS()
1172 PetscCall(VecGetArray(rr, &x)); in MatDiagonalScale_IS()
1175 rr = NULL; in MatDiagonalScale_IS()
1182 if (rr) { in MatDiagonalScale_IS()
1185 PetscCall(VecRestoreArray(rr, &x)); in MatDiagonalScale_IS()
1187 PetscCall(MatDiagonalScale(matis->A, ll, rr)); in MatDiagonalScale_IS()
/petsc/src/mat/impls/dense/mpi/
H A Dmpidense.c973 static PetscErrorCode MatDiagonalScale_MPIDense(Mat A, Vec ll, Vec rr) in MatDiagonalScale_MPIDense() argument
999 if (rr) { in MatDiagonalScale_MPIDense()
1002 PetscCall(VecGetLocalSize(rr, &s3a)); in MatDiagonalScale_MPIDense()
1004 PetscCall(VecGetArrayRead(rr, &ar)); in MatDiagonalScale_MPIDense()
1009 PetscCall(VecRestoreArrayRead(rr, &ar)); in MatDiagonalScale_MPIDense()
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcgraph.c1105 const PetscMPIInt rr = (PetscMPIInt)subs_remote[i].rank; in PCBDDCGraphSetUp() local
1130 iremote[nli].rank = rr; in PCBDDCGraphSetUp()
/petsc/src/mat/impls/dense/seq/
H A Ddense.c1903 static PetscErrorCode MatDiagonalScale_SeqDense(Mat A, Vec ll, Vec rr) in MatDiagonalScale_SeqDense() argument
1927 if (rr) { in MatDiagonalScale_SeqDense()
1928 PetscCall(VecGetSize(rr, &n)); in MatDiagonalScale_SeqDense()
1929 PetscCall(VecGetArrayRead(rr, &r)); in MatDiagonalScale_SeqDense()
1936 PetscCall(VecRestoreArrayRead(rr, &r)); in MatDiagonalScale_SeqDense()

12