Home
last modified time | relevance | path

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

/petsc/include/petsc/private/ !
H A Dpetscimpl.h755 PetscReal b1[5]; \
757 b1[4] = 1; \
759 b1[4] = 0; \
761 b1[0] = -PetscRealPart(b0); \
762 b1[1] = PetscRealPart(b0); \
763 b1[2] = -PetscImaginaryPart(b0); \
764 b1[3] = PetscImaginaryPart(b0); \
765 …PetscCallMPI(MPIU_Allreduce(MPI_IN_PLACE, b1, 5, MPIU_REAL, MPIU_MAX, PetscObjectComm((PetscObject…
766 …PetscCheck(b1[4] > 0 || (PetscEqualReal(-b1[0], b1[1]) && PetscEqualReal(-b1[2], b1[3])), PetscObj…
771 PetscReal b0 = (b), b1[3]; \
[all …]
/petsc/src/mat/tests/ !
H A Dex129.c24 Vec x, b, y, b1; in main() local
81 PetscCall(VecDuplicate(y, &b1)); in main()
85 PetscCall(MatMult(A, x, b1)); in main()
86 PetscCall(VecAXPY(b1, -1.0, b)); in main()
87 PetscCall(VecNorm(b1, NORM_2, &norm)); in main()
92 PetscCall(MatMultTranspose(A, x, b1)); in main()
93 PetscCall(VecAXPY(b1, -1.0, b)); in main()
94 PetscCall(VecNorm(b1, NORM_2, &norm)); in main()
99 PetscCall(MatMult(A, y, b1)); in main()
100 PetscCall(VecScale(b1, -1.0)); in main()
[all …]
H A Dex221.c185 Vec x = NULL, b1 = NULL, b2 = NULL; in main() local
201 PetscCall(MatCreateVecs(A, &x, &b1)); in main()
204 PetscCall(VecSetRandom(b1, NULL)); in main()
207 PetscCall(VecSet(b1, -14.2)); in main()
209 PetscCall(VecDuplicate(b1, &b2)); in main()
210 PetscCall(VecCopy(b1, b2)); in main()
211 PetscCall(PetscObjectSetName((PetscObject)b1, "A_B1")); in main()
214 PetscCall(VecDestroy(&b1)); in main()
217 PetscCall(MatZeroRowsIS(A, zr, diag, x, b1)); in main()
220 PetscCall(MatZeroRowsColumnsIS(A, zr, diag, x, b1)); in main()
[all …]
H A Dex23.c793 PetscInt b1, b2; in main() local
795 for (b1 = 0; b1 < bs; b1++) { in main()
796 …for (b2 = 0; b2 < bs; b2++) vals[b1 * bs + b2] = i * bs * bs + b1 * bs + b2 + 1 + (b1 == b2 ? 1.0 … in main()
808 PetscInt b1, b2; in main() local
810 for (b1 = 0; b1 < bs; b1++) { in main()
812 …if (PetscAbsScalar(isbd[i * bs * bs + b1 * bs + b2] - aijbd[i * bs * bs + b1 * bs + b2]) > PETSC_S… in main()
814 …: %g %g\n", rank, i, b1, b2, (double)PetscAbsScalar(isbd[i * bs * bs + b1 * bs + b2]), (double)Pet… in main()
/petsc/src/binding/petsc4py/demo/legacy/taosolve/ !
H A Dchwirut.py44 b1, b2, b3 = X.array
45 F.array = y - np.exp(-b1*x)/(b2 + b3*x)
52 b1, b2, b3 = X.array
55 v = np.exp(-b1*u)/(b2+b3*u)
/petsc/src/vec/is/tests/ !
H A Dex8.c8 PetscInt n1, n2, b1, b2; in TestEqual() local
14 PetscCall(ISLocalToGlobalMappingGetBlockSize(m1, &b1)); in TestEqual()
19 flg = (PetscBool)(b1 == b2); in TestEqual()
20 …ETSC_COMM_SELF, "%s: different block sizes %" PetscInt_FMT " %" PetscInt_FMT "\n", tname, b1, b2)); in TestEqual()
41 PetscInt n, n1, n2, b1, b2; in main() local
98 PetscCall(ISLocalToGlobalMappingGetBlockSize(lg1v, &b1)); in main()
99 n1 /= b1; in main()
104 PetscCall(ISLocalToGlobalMappingCreate(PETSC_COMM_SELF, b1, n1, idx, PETSC_OWN_POINTER, &lg1lh)); in main()
/petsc/src/ksp/ksp/tests/ !
H A Dex32.c93 Vec b1; in main() local
96 PetscCall(VecDuplicate(b, &b1)); in main()
97 PetscCall(MatMult(A, x, b1)); in main()
98 PetscCall(VecAXPY(b1, -1.0, b)); in main()
99 PetscCall(VecNorm(b1, NORM_2, &norm)); in main()
101 PetscCall(VecDestroy(&b1)); in main()
H A Dex11.c43 PetscErrorCode LSCLoadTestOperators(Mat *A11, Mat *A12, Mat *A21, Mat *A22, Vec *b1, Vec *b2) in LSCLoadTestOperators() argument
58 PetscCall(VecCreate(PETSC_COMM_WORLD, b1)); in LSCLoadTestOperators()
68 PetscCall(VecLoad(*b1, viewer)); in LSCLoadTestOperators()
/petsc/src/ksp/ksp/tutorials/ !
H A Dex9.c24 Vec x1, b1, x2, b2; /* solution and RHS vectors for systems #1 and #2 */ in main() local
84 PetscCall(VecDuplicate(u, &b1)); in main()
231 PetscCall(MatMult(C1, u, b1)); in main()
255 PetscCall(KSPSolve(ksp1, b1, x1)); in main()
261 PetscCall(CheckError(u, x1, b1, its, 1.e-4, CHECK_ERROR)); in main()
374 PetscCall(VecDestroy(&b1)); in main()
/petsc/src/dm/dt/tests/ !
H A Dex1.c207 PetscReal a1, b1, a2, b2; in main() local
211 b1 = 1.1; in main()
216 b1 = 1.; in main()
221 PetscCall(CheckJacobiQuadrature(i, a1, b1, PetscDTGaussJacobiQuadrature, 2 * i - 1)); in main()
225 PetscCall(CheckJacobiQuadrature(i, a1, b1, PetscDTGaussLobattoJacobiQuadrature, 2 * i - 3)); in main()
/petsc/src/snes/impls/al/ !
H A Dal.c408 PetscReal a0, b0, b1, c0, c1, c2; in SNESSolve_NEWTONAL() local
422 b1 = 2 * PetscRealPart(coefs2[0]); in SNESSolve_NEWTONAL()
427 as = b1 * b1 - 4 * a0 * c2; in SNESSolve_NEWTONAL()
428 bs = 2 * (b1 * b0 - 2 * a0 * c1); in SNESSolve_NEWTONAL()
440 deltaLambda = -0.5 * (b0 + b1 * deltaS) / a0; in SNESSolve_NEWTONAL()
445 PetscQuadraticRoots(a0, b0 + b1, c0 + c1 + c2, &dlambda1, &dlambda2); in SNESSolve_NEWTONAL()
/petsc/src/mat/impls/aij/seq/ !
H A Dmatrart.c107 const PetscScalar *b, *b1, *b2, *b3, *b4; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense() local
141 b1 = b; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
142 b2 = b1 + bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
156 r1 += (*aa) * b1[*aj]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
166 b1 += bm4; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
196 for (j = 0; j < n; j++) r1 += (*aa++) * b1[*aj++]; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
199 b1 += bm; in MatMatMatMultNumeric_SeqAIJ_SeqAIJ_SeqDense()
H A Dmatmatmult.c1499 const PetscScalar *aa, *b, *b1, *b2, *b3, *b4, *av; in MatMatMultNumericAdd_SeqAIJ_SeqDense() local
1519 b1 = b; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1520 b2 = b1 + bm; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1523 } else b1 = b2 = b3 = b4 = NULL; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1537 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1555 b1 += bm4; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1575 for (j = 0; j < n; j++) r1 += aa[j] * b1[aj[j]]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1588 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
1608 r1 += aatmp * b1[ajtmp]; in MatMatMultNumericAdd_SeqAIJ_SeqDense()
/petsc/src/snes/tutorials/ !
H A Dex70.c222 Vec b0, b1; in StokesRhs() local
240 PetscCall(VecGetSubVector(s->b, s->isg[1], &b1)); in StokesRhs()
241 PetscCall(VecGetOwnershipRange(b1, &start, &end)); in StokesRhs()
246 PetscCall(VecSetValue(b1, row, val, INSERT_VALUES)); in StokesRhs()
248 PetscCall(VecRestoreSubVector(s->b, s->isg[1], &b1)); in StokesRhs()
657 Vec b0, b1; in StokesCalcResidual() local
671 PetscCall(VecGetSubVector(s->b, s->isg[1], &b1)); in StokesCalcResidual()
672 PetscCall(VecNorm(b1, NORM_2, &val)); in StokesCalcResidual()
674 PetscCall(VecRestoreSubVector(s->b, s->isg[1], &b1)); in StokesCalcResidual()
/petsc/src/ts/adapt/impls/dsp/ !
H A Dadaptdsp.c121 PetscReal b1 = dsp->kBeta[0]; in TSAdaptChoose_DSP() local
137 rho0 = PetscPowReal(ctr0, b1); in TSAdaptChoose_DSP()
188 double b1 = (double)dsp->kBeta[0], b2 = (double)dsp->kBeta[1], b3 = (double)dsp->kBeta[2]; in TSAdaptView_DSP() local
189 …ewerASCIIPrintf(viewer, "filter parameters kBeta=[%g,%g,%g] Alpha=[%g,%g]\n", b1, b2, b3, a2, a3)); in TSAdaptView_DSP()
/petsc/src/sys/utils/ !
H A Dmpits.c31 PetscMPIInt b1[2], b2[2]; in PetscCommBuildTwoSidedSetType() local
32 b1[0] = -(PetscMPIInt)twosided; in PetscCommBuildTwoSidedSetType()
33 b1[1] = (PetscMPIInt)twosided; in PetscCommBuildTwoSidedSetType()
34 PetscCallMPI(MPIU_Allreduce(b1, b2, 2, MPI_INT, MPI_MAX, comm)); in PetscCommBuildTwoSidedSetType()
/petsc/src/ksp/pc/impls/parms/ !
H A Dparms.c302 const PetscScalar *b1; in PCApply_PARMS() local
306 PetscCall(VecGetArrayRead(b, &b1)); in PCApply_PARMS()
308 parms_VecPermAux((PetscScalar *)b1, parms->lvec0, parms->map); in PCApply_PARMS()
311 PetscCall(VecRestoreArrayRead(b, &b1)); in PCApply_PARMS()
/petsc/src/ts/impls/rosw/ !
H A Drosw.c365 const PetscReal b1[2] = {1.0, 0.0}; in TSRosWRegisterAll() local
372 PetscCall(TSRosWRegister(TSROSW2P, 2, 2, &A[0][0], &Gamma[0][0], b, b1, 2, &binterpt[0][0])); in TSRosWRegisterAll()
385 const PetscReal b1[2] = {1.0, 0.0}; in TSRosWRegisterAll() local
392 PetscCall(TSRosWRegister(TSROSW2M, 2, 2, &A[0][0], &Gamma[0][0], b, b1, 2, &binterpt[0][0])); in TSRosWRegisterAll()
965 …PetscReal a4, a32, a42, a43, b1, b2, b4, beta2p, beta3p, beta4p, beta32, beta42, beta43, beta32b… in TSRosWRegisterRos4() local
989 b1 = PetscRealPart(M[0][0] * rhs[0] + M[0][1] * rhs[1] + M[0][2] * rhs[2]); in TSRosWRegisterRos4()
1055 b[0] = b1; in TSRosWRegisterRos4()
/petsc/lib/petsc/bin/maint/ !
H A Dtoclapack.sh2316 static doublereal b1, b2, c1, c2, d1, d2, one, zero, rbase;
2360 b1 = dlamc3_(&d__1, &zero);
2370 a = b1;
2372 b1 = dlamc3_(&d__1, &zero);
2373 d__1 = b1 * *base;
2378 d1 += b1;
4307 static halfreal b1, b2, c1, c2, d1, d2, one, zero, rbase;
4351 b1 = hlamc3_(&d__1, &zero);
4361 a = b1;
4363 b1 = hlamc3_(&d__1, &zero);
[all …]
/petsc/src/sys/logging/handler/impls/default/ !
H A Dlogdefault.c559 PetscInt64 b1[2], b2[2]; in PetscLogHandlerEventBegin_Default() local
561 b1[0] = -o1->cidx; in PetscLogHandlerEventBegin_Default()
562 b1[1] = o1->cidx; in PetscLogHandlerEventBegin_Default()
563 PetscCallMPI(MPIU_Allreduce(b1, b2, 2, MPIU_INT64, MPI_MAX, PetscObjectComm(o1))); in PetscLogHandlerEventBegin_Default()
624 PetscInt64 b1[2], b2[2]; in PetscLogHandlerEventEnd_Default() local
626 b1[0] = -o1->cidx; in PetscLogHandlerEventEnd_Default()
627 b1[1] = o1->cidx; in PetscLogHandlerEventEnd_Default()
628 PetscCallMPI(MPIU_Allreduce(b1, b2, 2, MPIU_INT64, MPI_MAX, PetscObjectComm(o1))); in PetscLogHandlerEventEnd_Default()
/petsc/src/vec/vec/interface/ !
H A Drvector.c214 PetscMPIInt b0 = (PetscMPIInt)flg, b1[2], b2[2]; in VecNorm() local
215 b1[0] = -b0; in VecNorm()
216 b1[1] = b0; in VecNorm()
217 PetscCallMPI(MPIU_Allreduce(b1, b2, 2, MPI_INT, MPI_MAX, PetscObjectComm((PetscObject)x))); in VecNorm()
220 PetscReal b1[2], b2[2]; in VecNorm() local
221 b1[0] = -(*val); in VecNorm()
222 b1[1] = *val; in VecNorm()
223 PetscCallMPI(MPIU_Allreduce(b1, b2, 2, MPIU_REAL, MPIU_MAX, PetscObjectComm((PetscObject)x))); in VecNorm()
/petsc/src/ts/impls/multirate/ !
H A Dmprk.c69 …onst PetscReal Abase[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], PetscReal A2[], P… in TSMPRKGenerateTableau2() argument
86 b1[k * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
93 …onst PetscReal Abase[], const PetscReal bbase[], PetscReal A1[], PetscReal b1[], PetscReal A2[], P… in TSMPRKGenerateTableau3() argument
120 b1[(k * ratio + n) * s + j] = bbase[j] / ratio / ratio; in TSMPRKGenerateTableau3()
/petsc/src/mat/impls/aij/mpi/ !
H A Dmpiaij.c6084 PetscCount t, t1, t2, b1, e1, b2, e2; in MatMergeEntries_Internal() local
6091 b1 = rowBegin1[r]; in MatMergeEntries_Internal()
6095 while (b1 < e1 && b2 < e2) { in MatMergeEntries_Internal()
6096 if (j1[b1] == j2[b2]) { /* Same column index and hence same nonzero */ in MatMergeEntries_Internal()
6097 j[t] = j1[b1]; in MatMergeEntries_Internal()
6100 b1 += jmap1[t1 + 1] - jmap1[t1]; /* Jump to next unique local nonzero */ in MatMergeEntries_Internal()
6105 } else if (j1[b1] < j2[b2]) { in MatMergeEntries_Internal()
6106 j[t] = j1[b1]; in MatMergeEntries_Internal()
6108 b1 += jmap1[t1 + 1] - jmap1[t1]; in MatMergeEntries_Internal()
6120 while (b1 < e1) { in MatMergeEntries_Internal()
[all …]
/petsc/share/petsc/datafiles/meshes/ !
H A Dtestcase3D.cas352 (premix/peters-coeff-b1 2.)
9038 a b b1 aa 22 2a
9043 b b1 b2 e 2a 2b
9050 b c b8 b1 23 2b
9054 b 20 b0 b1 22 23
9204 aa b1 b0 a9 22 0
9205 aa ab b2 b1 2a 0
9211 b1 b8 b7 b0 23 0
9212 b1 b2 b9 b8 2b 0
/petsc/doc/ !
H A Dpetsc.bib5945 @InBook{ b1, inbook