| /petsc/src/sys/tutorials/ |
| H A D | ex9f.F90 | 17 PetscBool, dimension(:), pointer :: barray
40 allocate (barray(nopt), source=PETSC_FALSE)
50 …tionsBoolArray('-boolarray', 'Get an application bool array', 'Man page', barray, nopt, flg, ierr))
52 write (IOBuffer, '("The bool array was set to ",*(L1," "))') barray
99 deallocate (barray)
|
| /petsc/src/mat/impls/aij/mpi/mkl_cpardiso/ |
| H A D | mkl_cpardiso.c | 380 const PetscScalar *barray; in MatSolve_MKL_CPARDISO() local
385 PetscCall(VecGetArrayRead(b, &barray)); in MatSolve_MKL_CPARDISO()
390 …rdiso->nrhs, mat_mkl_cpardiso->iparm, &mat_mkl_cpardiso->msglvl, (void *)barray, (void *)xarray, &… in MatSolve_MKL_CPARDISO()
394 PetscCall(VecRestoreArrayRead(b, &barray)); in MatSolve_MKL_CPARDISO()
403 const PetscScalar *barray; in MatForwardSolve_MKL_CPARDISO() local
408 PetscCall(VecGetArrayRead(b, &barray)); in MatForwardSolve_MKL_CPARDISO()
413 …rdiso->nrhs, mat_mkl_cpardiso->iparm, &mat_mkl_cpardiso->msglvl, (void *)barray, (void *)xarray, &… in MatForwardSolve_MKL_CPARDISO()
417 PetscCall(VecRestoreArrayRead(b, &barray)); in MatForwardSolve_MKL_CPARDISO()
426 const PetscScalar *barray; in MatBackwardSolve_MKL_CPARDISO() local
431 PetscCall(VecGetArrayRead(b, &barray)); in MatBackwardSolve_MKL_CPARDISO()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex45.c | 79 PetscScalar ***barray; in ComputeRHS() local
92 PetscCall(DMDAVecGetArray(dm, b, &barray)); in ComputeRHS()
98 barray[k][j][i] = 2.0 * (HxHydHz + HxHzdHy + HyHzdHx); in ComputeRHS()
100 barray[k][j][i] = Hx * Hy * Hz; in ComputeRHS()
105 PetscCall(DMDAVecRestoreArray(dm, b, &barray)); in ComputeRHS()
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex37.c | 25 PetscScalar *barray, *xarray, *uarray, *array, one = 1.0; in main() local
125 PetscCall(VecGetArray(b, &barray)); in main()
128 PetscCall(VecPlaceArray(subb, barray)); in main()
166 PetscCall(VecRestoreArray(b, &barray)); in main()
|
| /petsc/src/mat/impls/aij/seq/mkl_pardiso/ |
| H A D | mkl_pardiso.c | 380 const PetscScalar *barray; in MatSolve_MKL_PARDISO() local
385 PetscCall(VecGetArrayRead(b, &barray)); in MatSolve_MKL_PARDISO()
390 if (barray == xarray) { /* if the two vectors share the same memory */ in MatSolve_MKL_PARDISO()
404 …pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)barray, (void *)xarray, &… in MatSolve_MKL_PARDISO()
406 PetscCall(VecRestoreArrayRead(b, &barray)); in MatSolve_MKL_PARDISO()
445 const PetscScalar *barray; in MatForwardSolve_MKL_PARDISO() local
452 PetscCall(VecGetArrayRead(b, &barray)); in MatForwardSolve_MKL_PARDISO()
457 …pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)barray, (void *)xarray, &… in MatForwardSolve_MKL_PARDISO()
460 PetscCall(VecRestoreArrayRead(b, &barray)); in MatForwardSolve_MKL_PARDISO()
470 const PetscScalar *barray; in MatBackwardSolve_MKL_PARDISO() local
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex245.c | 10 PetscScalar *v, *barray, rval; in main() local
103 PetscCall(VecGetArray(b, &barray)); in main()
106 barray[j] = rval; in main()
108 PetscCall(VecRestoreArray(b, &barray)); in main()
|
| H A D | ex145.c | 10 PetscScalar *v, *barray, rval; in main() local
104 PetscCall(VecGetArray(b, &barray)); in main()
107 barray[j] = rval; in main()
109 PetscCall(VecRestoreArray(b, &barray)); in main()
|
| /petsc/src/mat/impls/aij/seq/superlu/ |
| H A D | superlu.c | 94 const PetscScalar *barray; in MatSolve_SuperLU_Private() local
114 PetscCall(VecGetArrayRead(b, &barray)); in MatSolve_SuperLU_Private()
115 PetscCall(PetscArraycpy(lu->rhs_dup, barray, n)); in MatSolve_SuperLU_Private()
116 PetscCall(VecRestoreArrayRead(b, &barray)); in MatSolve_SuperLU_Private()
117 barray = lu->rhs_dup; in MatSolve_SuperLU_Private()
119 PetscCall(VecGetArrayRead(b, &barray)); in MatSolve_SuperLU_Private()
125 ((DNformat *)lu->B.Store)->nzval = (singlecomplex *)barray; in MatSolve_SuperLU_Private()
128 ((DNformat *)lu->B.Store)->nzval = (doublecomplex *)barray; in MatSolve_SuperLU_Private()
132 ((DNformat *)lu->B.Store)->nzval = (void *)barray; in MatSolve_SuperLU_Private()
166 if (!lu->options.Equil) PetscCall(VecRestoreArrayRead(b, &barray)); in MatSolve_SuperLU_Private()
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.h | 18 MatScalar *barray; /* Block array of size bs2 */ \
|
| H A D | mpibaij.c | 26 PetscCall(PetscFree(baij->barray)); in MatDestroy_MPIBAIJ()
411 MatScalar *barray = baij->barray; in MatSetValuesBlocked_MPIBAIJ() local
418 if (!barray) { in MatSetValuesBlocked_MPIBAIJ()
419 PetscCall(PetscMalloc1(bs2, &barray)); in MatSetValuesBlocked_MPIBAIJ()
420 baij->barray = barray; in MatSetValuesBlocked_MPIBAIJ()
434 barray = (MatScalar *)v + i * bs2; in MatSetValuesBlocked_MPIBAIJ()
436 barray = (MatScalar *)v + j * bs2; in MatSetValuesBlocked_MPIBAIJ()
444 for (jj = 0; jj < bs; jj++) barray[jj] = value[jj]; in MatSetValuesBlocked_MPIBAIJ()
445 barray += bs; in MatSetValuesBlocked_MPIBAIJ()
447 barray -= bs2; in MatSetValuesBlocked_MPIBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/seq/seqcusparse/ |
| H A D | aijcusparse.cu | 718 const PetscScalar *barray; in MatSolve_SeqAIJCUSPARSE_Cholesky() local
728 PetscCall(VecCUDAGetArrayRead(b, &barray)); in MatSolve_SeqAIJCUSPARSE_Cholesky()
730 bGPU = thrust::device_pointer_cast(barray); in MatSolve_SeqAIJCUSPARSE_Cholesky()
737 PetscCallCUSPARSE(cusparseDnVecSetValues(fs->dnVecDescr_X, (void *)barray)); in MatSolve_SeqAIJCUSPARSE_Cholesky()
766 PetscCall(VecCUDARestoreArrayRead(b, &barray)); in MatSolve_SeqAIJCUSPARSE_Cholesky()
1412 const PetscScalar *barray; in MatSolve_SeqAIJCUSPARSE_LU() local
1425 PetscCall(VecCUDAGetArrayRead(b, &barray)); in MatSolve_SeqAIJCUSPARSE_LU()
1427 bGPU = thrust::device_pointer_cast(barray); in MatSolve_SeqAIJCUSPARSE_LU()
1434 PetscCallCUSPARSE(cusparseDnVecSetValues(fs->dnVecDescr_X, (void *)barray)); in MatSolve_SeqAIJCUSPARSE_LU()
1455 PetscCall(VecCUDARestoreArrayRead(b, &barray)); in MatSolve_SeqAIJCUSPARSE_LU()
[all …]
|
| /petsc/src/mat/impls/aij/seq/seqhipsparse/ |
| H A D | aijhipsparse.hip.cxx | 1088 const PetscScalar *barray; in MatSolveTranspose_SeqAIJHIPSPARSE() local
1107 PetscCall(VecHIPGetArrayRead(bb, &barray)); in MatSolveTranspose_SeqAIJHIPSPARSE()
1109 bGPU = thrust::device_pointer_cast(barray); in MatSolveTranspose_SeqAIJHIPSPARSE()
1130 PetscCall(VecHIPRestoreArrayRead(bb, &barray)); in MatSolveTranspose_SeqAIJHIPSPARSE()
1139 const PetscScalar *barray; in MatSolveTranspose_SeqAIJHIPSPARSE_NaturalOrdering() local
1156 PetscCall(VecHIPGetArrayRead(bb, &barray)); in MatSolveTranspose_SeqAIJHIPSPARSE_NaturalOrdering()
1161 …iFactorT->csrMat->column_indices->data().get(), upTriFactorT->solveInfo, barray, tempGPU->data().g… in MatSolveTranspose_SeqAIJHIPSPARSE_NaturalOrdering()
1168 PetscCall(VecHIPRestoreArrayRead(bb, &barray)); in MatSolveTranspose_SeqAIJHIPSPARSE_NaturalOrdering()
1177 const PetscScalar *barray; in MatSolve_SeqAIJHIPSPARSE() local
1189 PetscCall(VecHIPGetArrayRead(bb, &barray)); in MatSolve_SeqAIJHIPSPARSE()
[all …]
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 31 PetscCall(PetscFree(baij->barray)); in MatDestroy_MPISBAIJ()
568 MatScalar *barray = baij->barray; in MatSetValuesBlocked_MPISBAIJ() local
575 if (!barray) { in MatSetValuesBlocked_MPISBAIJ()
576 PetscCall(PetscMalloc1(bs2, &barray)); in MatSetValuesBlocked_MPISBAIJ()
577 baij->barray = barray; in MatSetValuesBlocked_MPISBAIJ()
597 barray = (MatScalar *)v + i * bs2; in MatSetValuesBlocked_MPISBAIJ()
599 barray = (MatScalar *)v + j * bs2; in MatSetValuesBlocked_MPISBAIJ()
607 for (jj = 0; jj < bs; jj++) *barray++ = *value++; in MatSetValuesBlocked_MPISBAIJ()
609 barray -= bs2; in MatSetValuesBlocked_MPISBAIJ()
614 … PetscCall(MatSetValuesBlocked_SeqSBAIJ_Inlined(baij->A, row, col, barray, addv, im[i], in[j])); in MatSetValuesBlocked_MPISBAIJ()
[all …]
|
| /petsc/src/mat/impls/aij/mpi/mumps/ |
| H A D | mumps.c | 1925 const PetscScalar *barray = NULL; in MatSolve_MUMPS() local
1950 PetscCall(VecGetArrayRead(b, &barray)); in MatSolve_MUMPS()
1951 PetscCall(MatMumpsSetUpDistRHSInfo(A, 1, barray)); in MatSolve_MUMPS()
2020 PetscCall(VecRestoreArrayRead(b, &barray)); in MatSolve_MUMPS()
2053 const PetscScalar *barray; in MatMatSolve_MUMPS() local
2103 PetscCall(MatDenseGetArrayRead(B, &barray)); in MatMatSolve_MUMPS()
2104 PetscCall(PetscArraycpy(array, barray, nrhsM)); in MatMatSolve_MUMPS()
2105 PetscCall(MatDenseRestoreArrayRead(B, &barray)); in MatMatSolve_MUMPS()
2170 PetscCall(MatDenseGetArrayRead(B, &barray)); in MatMatSolve_MUMPS()
2171 PetscCall(MatMumpsSetUpDistRHSInfo(A, nrhs, barray)); // put barray to rhs_loc in MatMatSolve_MUMPS()
[all …]
|
| /petsc/src/mat/impls/scalapack/ |
| H A D | matscalapack.c | 982 PetscScalar *barray; in MatConvert_ScaLAPACK_Dense() local
1018 PetscCall(MatDenseGetArray(Bmpi, &barray)); in MatConvert_ScaLAPACK_Dense()
1019 …PACKgemr2d", SCALAPACKgemr2d_(&a->M, &a->N, a->loc, &one, &one, a->desc, barray, &one, &one, bdesc… in MatConvert_ScaLAPACK_Dense()
1020 PetscCall(MatDenseRestoreArray(Bmpi, &barray)); in MatConvert_ScaLAPACK_Dense()
1057 PetscCall(MatDenseGetArray(Bmpi, &barray)); in MatConvert_ScaLAPACK_Dense()
1058 …PACKgemr2d", SCALAPACKgemr2d_(&a->M, &a->N, a->loc, &one, &one, a->desc, barray, &one, &one, bdesc… in MatConvert_ScaLAPACK_Dense()
1059 PetscCall(MatDenseRestoreArray(Bmpi, &barray)); in MatConvert_ScaLAPACK_Dense()
|
| /petsc/src/mat/impls/nest/ |
| H A D | matnest.c | 98 const PetscScalar *barray; in MatProductNumeric_Nest_Dense() local
121 PetscCall(MatDenseGetArrayRead(B, &barray)); in MatProductNumeric_Nest_Dense()
130 …ts->dn[j], PETSC_DECIDE, M, N, PetscSafePointerPlusOffset((PetscScalar *)barray, contents->dn[j]),… in MatProductNumeric_Nest_Dense()
147 PetscCall(MatDenseRestoreArrayRead(B, &barray)); in MatProductNumeric_Nest_Dense()
172 const PetscScalar *barray; in MatProductSymbolic_Nest_Dense() local
222 PetscCall(MatDenseGetArrayRead(B, &barray)); in MatProductSymbolic_Nest_Dense()
226 …ts->dn[j], PETSC_DECIDE, M, N, PetscSafePointerPlusOffset((PetscScalar *)barray, contents->dn[j]),… in MatProductSymbolic_Nest_Dense()
246 PetscCall(MatDenseRestoreArrayRead(B, &barray)); in MatProductSymbolic_Nest_Dense()
|