| /petsc/src/vec/vec/impls/mpi/mpiviennacl/ |
| H A D | mpiviennacl.cxx | 159 static PetscErrorCode VecBindToCPU_MPIViennaCL(Vec vv, PetscBool bind) in VecBindToCPU_MPIViennaCL() argument 162 vv->boundtocpu = bind; in VecBindToCPU_MPIViennaCL() 165 PetscCall(VecViennaCLCopyFromGPU(vv)); in VecBindToCPU_MPIViennaCL() 166 …vv->offloadmask = PETSC_OFFLOAD_CPU; /* since the CPU code will likely change valu… in VecBindToCPU_MPIViennaCL() 167 vv->ops->dotnorm2 = NULL; in VecBindToCPU_MPIViennaCL() 168 vv->ops->waxpy = VecWAXPY_Seq; in VecBindToCPU_MPIViennaCL() 169 vv->ops->dot = VecDot_MPI; in VecBindToCPU_MPIViennaCL() 170 vv->ops->mdot = VecMDot_MPI; in VecBindToCPU_MPIViennaCL() 171 vv->ops->tdot = VecTDot_MPI; in VecBindToCPU_MPIViennaCL() 172 vv->ops->norm = VecNorm_MPI; in VecBindToCPU_MPIViennaCL() [all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex230.c | 43 PetscScalar vv = 0.0; in ex1_nonsquare_bs1() local 47 PetscCall(MatSetValues(preallocator, 1, &ii, 1, &jj, &vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 51 PetscCall(MatSetValues(preallocator, 1, &ii, 1, &jj, &vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 55 PetscCall(MatSetValue(preallocator, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 79 PetscScalar vv; in ex1_nonsquare_bs1() local 83 vv = 0.3; in ex1_nonsquare_bs1() 84 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 88 vv = 3.3; in ex1_nonsquare_bs1() 89 PetscCall(MatSetValues(A, 1, &ii, 1, &jj, &vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 93 vv = 4.3; in ex1_nonsquare_bs1() [all …]
|
| H A D | ex259.c | 38 PetscScalar vv = 22.0; in ex1_nonsquare_bs1() local 42 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 46 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 47 PetscCall(MatSetValue(A, jj, ii, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 51 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 52 PetscCall(MatSetValue(A, jj, ii, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 65 PetscScalar vv; in ex1_nonsquare_bs1() local 69 vv = 0.3; in ex1_nonsquare_bs1() 70 PetscCall(MatSetValue(A, ii, jj, vv, INSERT_VALUES)); in ex1_nonsquare_bs1() 74 vv = 3.3; in ex1_nonsquare_bs1() [all …]
|
| H A D | ex69.c | 36 PetscScalar *vv, *aa; in main() local 94 PetscCall(VecCUDAGetArray(v, &vv)); in main() 99 …PetscCall(MatCreateDenseCUDA(PetscObjectComm((PetscObject)v), nloc, PETSC_DECIDE, n, k - l, vv, &B… in main() 101 …PetscCall(MatDenseCUDAPlaceArray(B, vv + l * nloc)); … in main() 106 PetscCall(MatDenseCUDAPlaceArray(B, vv + l * nloc)); /* set the actual pointer */ in main() 110 …PetscCall(MatCreateDenseCUDA(PetscObjectComm((PetscObject)v), nloc, PETSC_DECIDE, n, k - l, vv + l… in main() 154 PetscCheck(vv == aa - l * nloc, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Wrong array"); in main() 157 PetscCall(VecCUDARestoreArray(v, &vv)); in main()
|
| /petsc/src/vec/vec/impls/shared/ |
| H A D | shvec.c | 36 PETSC_EXTERN PetscErrorCode VecCreate_Shared(Vec vv) in VecCreate_Shared() argument 41 PetscCall(PetscSplitOwnership(PetscObjectComm((PetscObject)vv), &vv->map->n, &vv->map->N)); in VecCreate_Shared() 42 …etscCall(PetscSharedMalloc(PetscObjectComm((PetscObject)vv), vv->map->n * sizeof(PetscScalar), vv-… in VecCreate_Shared() 44 PetscCall(VecCreate_MPI_Private(vv, PETSC_FALSE, 0, array)); in VecCreate_Shared() 45 vv->ops->duplicate = VecDuplicate_Shared; in VecCreate_Shared() 142 PETSC_EXTERN PetscErrorCode VecCreate_Shared(Vec vv) in VecCreate_Shared() argument 147 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)vv), &size)); in VecCreate_Shared() 149 PetscCall(VecCreate_Seq(vv)); in VecCreate_Shared()
|
| /petsc/src/vec/vec/impls/mpi/ |
| H A D | pbvec.c | 650 PetscErrorCode VecCreate_MPI(Vec vv) in VecCreate_MPI() argument 653 PetscCall(VecCreate_MPI_Private(vv, PETSC_TRUE, 0, NULL)); in VecCreate_MPI() 714 …PIWithArray(MPI_Comm comm, PetscInt bs, PetscInt n, PetscInt N, const PetscScalar array[], Vec *vv) in VecCreateMPIWithArray() argument 719 PetscCall(VecCreate(comm, vv)); in VecCreateMPIWithArray() 720 PetscCall(VecSetSizes(*vv, n, N)); in VecCreateMPIWithArray() 721 PetscCall(VecSetBlockSize(*vv, bs)); in VecCreateMPIWithArray() 722 PetscCall(VecCreate_MPI_Private(*vv, PETSC_FALSE, 0, array)); in VecCreateMPIWithArray() 755 …tscInt n, PetscInt N, PetscInt nghost, const PetscInt ghosts[], const PetscScalar array[], Vec *vv) in VecCreateGhostWithArray() argument 762 *vv = NULL; in VecCreateGhostWithArray() 767 PetscCall(VecCreate(comm, vv)); in VecCreateGhostWithArray() [all …]
|
| /petsc/src/mat/utils/ |
| H A D | zerodiag.c | 73 PetscScalar *v, *vv; in MatReorderForNonzeroDiagonal_SeqAIJ() local 114 PetscCall(MatGetRow_SeqAIJ(mat, row[repl], &nnz, &jj, &vv)); in MatReorderForNonzeroDiagonal_SeqAIJ() 116 if (icol[jj[kk]] == prow && PetscAbsScalar(vv[kk]) > abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ() 117 PetscCall(MatRestoreRow_SeqAIJ(mat, row[repl], &nnz, &jj, &vv)); in MatReorderForNonzeroDiagonal_SeqAIJ() 123 PetscCall(MatRestoreRow_SeqAIJ(mat, row[repl], &nnz, &jj, &vv)); in MatReorderForNonzeroDiagonal_SeqAIJ() 131 PetscCall(MatGetRow_SeqAIJ(mat, row[k], &nnz, &jj, &vv)); in MatReorderForNonzeroDiagonal_SeqAIJ() 133 if (icol[jj[kk]] == prow && PetscAbsScalar(vv[kk]) > abstol) { in MatReorderForNonzeroDiagonal_SeqAIJ() 139 PetscCall(MatRestoreRow_SeqAIJ(mat, row[k], &nnz, &jj, &vv)); in MatReorderForNonzeroDiagonal_SeqAIJ()
|
| /petsc/src/mat/tutorials/ |
| H A D | ex12.c | 20 const PetscScalar *vv, *aa; in PadMatrix() local 24 PetscCall(VecGetArrayRead(v, &vv)); in PadMatrix() 32 cnt[i] = nc + (vv[i] != 0.0); in PadMatrix() 36 for (i = 0; i < n; i++) cnt[n] += (vv[i] != 0.0); in PadMatrix() 46 PetscCall(MatSetValues(*B, 1, &n, n, indices, vv, INSERT_VALUES)); in PadMatrix() 47 PetscCall(MatSetValues(*B, n, indices, 1, &n, vv, INSERT_VALUES)); in PadMatrix() 52 PetscCall(VecRestoreArrayRead(v, &vv)); in PadMatrix()
|
| /petsc/src/dm/tests/ |
| H A D | ex21.c | 11 PetscViewer vv; in test1_DAInjection3d() local 80 PetscCall(PetscViewerASCIIOpen(PETSC_COMM_WORLD, "dac_7.vtu", &vv)); in test1_DAInjection3d() 81 PetscCall(VecView(ac, vv)); in test1_DAInjection3d() 82 PetscCall(PetscViewerDestroy(&vv)); in test1_DAInjection3d() 84 PetscCall(PetscViewerASCIIOpen(PETSC_COMM_WORLD, "daf_7.vtu", &vv)); in test1_DAInjection3d() 85 PetscCall(VecView(af, vv)); in test1_DAInjection3d() 86 PetscCall(PetscViewerDestroy(&vv)); in test1_DAInjection3d()
|
| H A D | ex36.c | 312 PetscViewer vv; in da_test_RefineCoords1D() local 372 PetscCall(PetscViewerASCIIOpen(PETSC_COMM_WORLD, "dac_1D.vtr", &vv)); in da_test_RefineCoords1D() 373 PetscCall(VecView(ac, vv)); in da_test_RefineCoords1D() 374 PetscCall(PetscViewerDestroy(&vv)); in da_test_RefineCoords1D() 376 PetscCall(PetscViewerASCIIOpen(PETSC_COMM_WORLD, "daf_1D.vtr", &vv)); in da_test_RefineCoords1D() 377 PetscCall(VecView(af, vv)); in da_test_RefineCoords1D() 378 PetscCall(PetscViewerDestroy(&vv)); in da_test_RefineCoords1D() 392 PetscViewer vv; in da_test_RefineCoords2D() local 458 PetscCall(PetscViewerASCIIOpen(PETSC_COMM_WORLD, "dac_2D.vtr", &vv)); in da_test_RefineCoords2D() 459 PetscCall(VecView(ac, vv)); in da_test_RefineCoords2D() [all …]
|
| /petsc/src/vec/vec/impls/mpi/cupm/ |
| H A D | vecmpicupm_impl.hpp | 276 PETSC_KERNEL_DECL void pack_coo_values(const PetscScalar *PETSC_RESTRICT vv, PetscCount nnz, const … in pack_coo_values() argument 278 …Petsc::device::cupm::kernels::util::grid_stride_1D(nnz, [=](PetscCount i) { buf[i] = vv[perm[i]]; … in pack_coo_values() 282 PETSC_KERNEL_DECL void add_remote_coo_values(const PetscScalar *PETSC_RESTRICT vv, PetscCount nnz2,… in add_remote_coo_values() argument 284 …add_coo_values_impl(vv, nnz2, jmap2, perm2, ADD_VALUES, xv, [=](PetscCount i) { return imap2[i]; }… in add_remote_coo_values() 334 auto vv = const_cast<PetscScalar *>(v); in SetValuesCOO() local 340 PetscCall(PetscDeviceMalloc(dctx, PETSC_MEMTYPE_CUPM(), size, &vv)); in SetValuesCOO() 341 PetscCall(PetscCUPMMemcpyAsync(vv, v, size, cupmMemcpyHostToDevice, stream)); in SetValuesCOO() 346 …PetscCall(PetscCUPMLaunchKernel1D(sendlen, 0, stream, kernels::pack_coo_values, vv, sendlen, vcu->… in SetValuesCOO() 354 …) PetscCall(PetscCUPMLaunchKernel1D(n, 0, stream, kernels::add_coo_values, vv, n, vcu->jmap1_d, vc… in SetValuesCOO() 361 if (PetscMemTypeHost(v_memtype)) PetscCall(PetscDeviceFree(dctx, vv)); in SetValuesCOO()
|
| /petsc/src/dm/impls/network/ |
| H A D | networkmonitor.c | 186 PetscScalar *vv; in DMNetworkMonitorView() local 193 PetscCall(VecGetArray(node->v, &vv)); in DMNetworkMonitorView() 195 for (i = 0; i < node->nodes; i++) vv[i] = xx[start + i * node->blocksize]; in DMNetworkMonitorView() 196 PetscCall(VecRestoreArray(node->v, &vv)); in DMNetworkMonitorView()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex42.c | 31 PetscScalar vv = (PetscReal)Ii + 1; in main() local 32 PetscCall(MatSetValues(A, 1, &Ii, 1, &Ii, &vv, INSERT_VALUES)); in main()
|
| /petsc/src/dm/partitioner/tests/ |
| H A D | ex33.c | 13 PetscInt vv[5] = {0, 2, 4, 6, 8}; in main() local 87 …PetscCall(PetscPartitionerPartition(p, nparts, nv, vv, vadj, vertexSection, NULL, targetSection, p… in main() 113 …PetscCall(PetscPartitionerPartition(p, nparts, nv, vv, pvadj, NULL, NULL, targetSection, partSecti… in main()
|
| /petsc/src/ksp/ksp/impls/bcgs/qmrcgs/ |
| H A D | qmrcgs.c | 20 PetscReal dp = 0.0, final, tau, tau2, theta, theta2, c, F, NV, vv; in KSPSolve_QMRCGS() local 125 PetscCall(VecDotNorm2(S, T, &uu, &vv)); in KSPSolve_QMRCGS() 126 if (vv == 0.0) { in KSPSolve_QMRCGS() 156 omega = uu / vv; /* omega <- uu/vv; */ in KSPSolve_QMRCGS()
|
| /petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/ |
| H A D | perf_profile.m | 76 vv = find(xs==1,1,'last'); variable 77 xs = xs(vv:end); ys = ys(vv:end);
|
| /petsc/src/vec/vec/impls/mpi/kokkos/ |
| H A D | mpikok.kokkos.cxx | 98 static PetscErrorCode VecDuplicate_MPIKokkos(Vec win, Vec *vv) in VecDuplicate_MPIKokkos() argument 119 *vv = v; in VecDuplicate_MPIKokkos() 163 ConstPetscScalarKokkosView vv; in VecSetValuesCOO_MPIKokkos() local 171 …vv = Kokkos::create_mirror_view_and_copy(DefaultMemorySpace(), PetscScalarKokkosViewHost(const_cas… in VecSetValuesCOO_MPIKokkos() 173 vv = ConstPetscScalarKokkosView(v, vecmpi->coo_n); /* Directly use v[]'s memory */ in VecSetValuesCOO_MPIKokkos() 177 …ionSpace(), 0, vecmpi->sendlen), KOKKOS_LAMBDA(const PetscCount i) { sendbuf(i) = vv(Cperm(i)); }); in VecSetValuesCOO_MPIKokkos() 186 for (PetscCount k = jmap1(i); k < jmap1(i + 1); k++) sum += vv(perm1(k)); in VecSetValuesCOO_MPIKokkos()
|
| /petsc/src/sys/classes/random/tutorials/ |
| H A D | ex2.c | 155 double vv, tt; in readData() local 156 …PetscCheck(fscanf(fd, "%s%lf%lf", temp, &vv, &tt) == 3, PETSC_COMM_SELF, PETSC_ERR_FILE_UNEXPECTED… in readData() 157 v[i] = vv; in readData()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex56.c | 95 PetscScalar vv[24], v2[24]; in main() local 163 if (i % 3 == 0) vv[i] = h * h; in main() 164 else if (i % 3 == 1) vv[i] = 2.0 * h * h; in main() 165 else vv[i] = .0; in main() 229 PetscCall(VecSetValuesBlocked(bb, 8, idx, (const PetscScalar *)vv, ADD_VALUES)); in main() 237 PetscCall(VecSetValues(bb, 24, idx3, (const PetscScalar *)vv, ADD_VALUES)); in main()
|
| /petsc/src/binding/petsc4py/src/lib-petsc/compat/ |
| H A D | viennacl.h | 12 …cScalar cpuarray[],PETSC_UNUSED PetscScalar *viennaclvec,PETSC_UNUSED Vec *vv) {PetscViennaCLError… in VecCreateMPIViennaCLWithArrays() argument
|
| /petsc/src/ts/tutorials/ |
| H A D | ex50.c | 324 PetscReal **temp, vv; in MatMult_Laplacian() local 341 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in MatMult_Laplacian() 342 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in MatMult_Laplacian() 416 PetscReal vv; in RHSMatrixLaplaciangllDM() local 435 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in RHSMatrixLaplaciangllDM() 436 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSMatrixLaplaciangllDM()
|
| /petsc/src/dm/impls/plex/adaptors/parmmg/ |
| H A D | parmmgadapt.c | 193 PetscInt vv = vertexNumber[v]; in DMAdaptMetric_ParMmg_Plex() local 194 if (!vv--) continue; in DMAdaptMetric_ParMmg_Plex() 199 if (uniform) metric[Neq * vv + k] = PetscRealPart(met[0]); in DMAdaptMetric_ParMmg_Plex() 200 else metric[Neq * vv + k] = PetscRealPart(met[v]); in DMAdaptMetric_ParMmg_Plex() 201 } else metric[Neq * vv + k] = 0.0; in DMAdaptMetric_ParMmg_Plex() 202 } else metric[Neq * vv + k] = PetscRealPart(met[dim * dim * v + dim * i + j]); in DMAdaptMetric_ParMmg_Plex()
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | spectraladjointassimilation.c | 453 PetscReal vv; in RHSLaplacian() local 466 vv = -appctx->param.mu * 2.0 / appctx->param.Le; in RHSLaplacian() 467 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSLaplacian() 504 PetscReal vv; in RHSAdvection() local 517 vv = -appctx->param.a; in RHSAdvection() 518 for (j = 0; j < appctx->param.N; j++) temp[i][j] = temp[i][j] * vv; in RHSAdvection()
|
| /petsc/src/vec/vec/interface/ |
| H A D | vector.c | 656 PetscErrorCode VecDestroyVecs(PetscInt m, Vec *vv[]) in VecDestroyVecs() argument 659 PetscAssertPointer(vv, 2); in VecDestroyVecs() 661 if (!m || !*vv) { in VecDestroyVecs() 662 *vv = NULL; in VecDestroyVecs() 665 PetscValidHeaderSpecific(**vv, VEC_CLASSID, 2); in VecDestroyVecs() 666 PetscValidType(**vv, 2); in VecDestroyVecs() 667 PetscCall((*(**vv)->ops->destroyvecs)(m, *vv)); in VecDestroyVecs() 668 *vv = NULL; in VecDestroyVecs()
|
| /petsc/src/mat/impls/h2opus/cuda/ |
| H A D | math2opus.cu | 1831 const PetscScalar *u, *v, *uu, *vv; in MatH2OpusLowRankUpdate() local 1867 vv = MatH2OpusGetThrustPointer(*a->yy); in MatH2OpusLowRankUpdate() 1868 PetscCall(PetscSFBcastBegin(vsf, MPIU_SCALAR, v, (PetscScalar *)vv, MPI_REPLACE)); in MatH2OpusLowRankUpdate() 1869 PetscCall(PetscSFBcastEnd(vsf, MPIU_SCALAR, v, (PetscScalar *)vv, MPI_REPLACE)); in MatH2OpusLowRankUpdate() 1874 } else uu = vv; in MatH2OpusLowRankUpdate() 1877 vv = v; in MatH2OpusLowRankUpdate() 1879 hlru_global(*a->hmatrix, uu, ldu, vv, ldv, U->cmap->N, s, handle); in MatH2OpusLowRankUpdate() 1894 vv = MatH2OpusGetThrustPointer(*a->yy_gpu); in MatH2OpusLowRankUpdate() 1895 PetscCall(PetscSFBcastBegin(vsf, MPIU_SCALAR, v, (PetscScalar *)vv, MPI_REPLACE)); in MatH2OpusLowRankUpdate() 1896 PetscCall(PetscSFBcastEnd(vsf, MPIU_SCALAR, v, (PetscScalar *)vv, MPI_REPLACE)); in MatH2OpusLowRankUpdate() [all …]
|