| /petsc/share/petsc/matlab/ |
| H A D | generatePetscTestFiles.m | 17 [At,bt] = PetscBinaryRead('spd-real-int32-float64', 'indices','int32','precision','float64');
18 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0
25 [At,bt] = PetscBinaryRead('spd-real-int64-float64','indices','int64','precision','float64');
26 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0
34 [At,bt] = PetscBinaryRead('spd-real-int32-float32', 'indices','int32','precision','float32');
35 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0
42 [At,bt] = PetscBinaryRead('spd-real-int64-float32','indices','int64','precision','float32');
43 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0
58 [At,bt] = PetscBinaryRead('ns-real-int32-float64','indices','int32','precision','float64');
59 if max(max(max(abs(At-A))),max(abs(bt-b))) ~= 0
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | mattransposematmult.c | 16 PetscCall(VecDestroy(&atb->bt)); in MatProductCtxDestroy_SeqDense_MatTransMatMult()
50 PetscCall(MatCreateVecs(atb->mA, &atb->ct, &atb->bt)); in MatTMatTMultSymbolic_SeqAIJ_SeqDense()
69 Vec bt, ct; in MatTMatTMultNumeric_SeqAIJ_SeqDense() local
76 bt = atb->bt; in MatTMatTMultNumeric_SeqAIJ_SeqDense()
87 PetscCall(VecGetArrayWrite(bt, &btarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
91 PetscCall(VecRestoreArrayWrite(bt, &btarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
94 PetscCall(MatMultTranspose(atb->mA, bt, ct)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
119 PetscCall(MatMult(atb->mA, ct, bt)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
121 PetscCall(VecGetArrayRead(bt, &btarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
125 PetscCall(VecRestoreArrayRead(bt, &btarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpimattransposematmult.c | 16 PetscCall(VecDestroy(&atb->bt)); in MatProductCtxDestroy_MPIDense_MatTransMatMult()
47 PetscCall(MatCreateVecs(atb->mA, &atb->ct, &atb->bt)); in MatTransposeMatMultSymbolic_MPIAIJ_MPIDense()
62 Vec bt, ct; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() local
73 bt = atb->bt; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
79 PetscCall(VecGetArray(bt, &btarray)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
82 PetscCall(VecRestoreArray(bt, &btarray)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
86 PetscCall(MatMultTranspose(atb->mA, bt, ct)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
|
| /petsc/src/snes/linesearch/impls/bt/ |
| H A D | linesearchbt.c | 21 SNESLineSearch_BT *bt = (SNESLineSearch_BT *)linesearch->data; in SNESLineSearchBTSetAlpha() local
27 if (isbt) bt->alpha = alpha; in SNESLineSearchBTSetAlpha()
46 SNESLineSearch_BT *bt = (SNESLineSearch_BT *)linesearch->data; in SNESLineSearchBTGetAlpha() local
53 *alpha = bt->alpha; in SNESLineSearchBTGetAlpha()
59 SNESLineSearch_BT *bt = (SNESLineSearch_BT *)linesearch->data; in SNESLineSearchApply_BT() local
83 alpha = bt->alpha; in SNESLineSearchApply_BT()
352 SNESLineSearch_BT *bt = (SNESLineSearch_BT *)linesearch->data; in SNESLineSearchView_BT() local
362 PetscCall(PetscViewerASCIIPrintf(viewer, " alpha=%e\n", (double)bt->alpha)); in SNESLineSearchView_BT()
376 SNESLineSearch_BT *bt = (SNESLineSearch_BT *)linesearch->data; in SNESLineSearchSetFromOptions_BT() local
380 …l("-snes_linesearch_alpha", "Descent tolerance", "SNESLineSearchBT", bt->alpha, &bt->alpha, NULL)); in SNESLineSearchSetFromOptions_BT()
[all …]
|
| /petsc/include/petsc/private/ |
| H A D | matimpl.h | 882 …LCreate(idx_start, lnk_max, nlnk, lnk, bt) ((PetscErrorCode)(PetscMalloc1(nlnk, &lnk) || PetscBTCr… argument
884 …ew(idx_start, lnk_max, nlnk, lnk, bt, lnk_empty) ((PetscErrorCode)(PetscMalloc1(nlnk, &lnk) || Pet… argument
906 …tart, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt, PetscBool assume_so… in PetscLLAdd_Private() argument
913 …if (!PetscBTLookupSet(bt, entry)) PetscCall(PetscLLInsertLocation_Private(assume_sorted, k, idx_st… in PetscLLAdd_Private()
931 …dices, PetscInt idx_start, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt) in PetscLLAdd() argument
934 PetscCall(PetscLLAdd_Private(nidx, indices, idx_start, nlnk, lnk, bt, PETSC_FALSE)); in PetscLLAdd()
951 …dices, PetscInt idx_start, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt) in PetscLLAddSorted() argument
954 PetscCall(PetscLLAdd_Private(nidx, indices, idx_start, nlnk, lnk, bt, PETSC_TRUE)); in PetscLLAddSorted()
972 … perm, PetscInt idx_start, PetscInt *PETSC_RESTRICT nlnk, PetscInt *PETSC_RESTRICT lnk, PetscBT bt) in PetscLLAddPerm() argument
979 …if (!PetscBTLookupSet(bt, entry)) PetscCall(PetscLLInsertLocation_Private(PETSC_FALSE, k, idx_star… in PetscLLAddPerm()
[all …]
|
| H A D | dmlabelimpl.h | 39 PetscBT bt; /* A bit-wise index */ member
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexreorder.c | 452 PetscBT bt, blst; in DMCreateSectionPermutation_Plex_Cohesive_Old() local
459 PetscCall(PetscBTCreate(pEnd - pStart, &bt)); in DMCreateSectionPermutation_Plex_Cohesive_Old()
469 if (PetscBTLookupSet(bt, p)) continue; in DMCreateSectionPermutation_Plex_Cohesive_Old()
489 if (!PetscBTLookupSet(bt, q)) { in DMCreateSectionPermutation_Plex_Cohesive_Old()
510 if (!PetscBTLookupSet(bt, qsupp[qs])) { in DMCreateSectionPermutation_Plex_Cohesive_Old()
522 if (!PetscBTLookupSet(bt, qq)) { in DMCreateSectionPermutation_Plex_Cohesive_Old()
533 …for (PetscInt p = pStart; p < pEnd; ++p) PetscCheck(PetscBTLookup(bt, p), PETSC_COMM_SELF, PETSC_E… in DMCreateSectionPermutation_Plex_Cohesive_Old()
535 PetscCall(PetscBTDestroy(&bt)); in DMCreateSectionPermutation_Plex_Cohesive_Old()
545 static PetscErrorCode InsertCohesiveBlock_Private(DM dm, PetscBT bt, PetscBT blst, PetscInt p, Pets… in InsertCohesiveBlock_Private() argument
551 if (PetscBTLookupSet(bt, p)) PetscFunctionReturn(PETSC_SUCCESS); in InsertCohesiveBlock_Private()
[all …]
|
| /petsc/src/sys/classes/viewer/utils/ |
| H A D | btview.c | 5 PetscErrorCode PetscBTView(PetscCount m, const PetscBT bt, PetscViewer viewer) in PetscBTView() argument
10 …PetscViewerASCIISynchronizedPrintf(viewer, "%" PetscCount_FMT " %hhu\n", i, PetscBTLookup(bt, i))); in PetscBTView()
|
| /petsc/src/snes/tutorials/network/ |
| H A D | ex1options | 11 # Line search method (bt and secant lead to convergence, basic diverges)
39 # Line search method (bt and secant lead to convergence, basic diverges)
64 # Line search method (bt and secant lead to convergence, basic diverges)
|
| /petsc/src/ts/tests/ |
| H A D | ex21.c | 140 DMBoundaryType bt = DM_BOUNDARY_NONE; in main() local
163 PetscCall(DMDACreate2d(PETSC_COMM_WORLD, bt, bt, st, N, N, n, n, 1, sw, NULL, NULL, &da)); in main()
|
| /petsc/src/snes/tutorials/network/water/ |
| H A D | wateroptions | 7 # Line search method (bt and secant lead to convergence, basic diverges)
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1032 PetscCall(PetscFree6(t->At, t->bt, t->ct, t->A, t->b, t->c)); in TSARKIMEXRegisterDestroy()
1104 …pe name, PetscInt order, PetscInt s, const PetscReal At[], const PetscReal bt[], const PetscReal c… in TSARKIMEXRegister() argument
1124 PetscCall(PetscMalloc6(s * s, &t->At, s, &t->bt, s, &t->ct, s * s, &t->A, s, &t->b, s, &t->c)); in TSARKIMEXRegister()
1131 if (bt) PetscCall(PetscArraycpy(t->bt, bt, s)); in TSARKIMEXRegister()
1133 for (i = 0; i < s; i++) t->bt[i] = At[(s - 1) * s + i]; in TSARKIMEXRegister()
1138 for (i = 0; i < s; i++) t->b[i] = t->bt[i]; in TSARKIMEXRegister()
1155 if (t->At[(s - 1) * s + i] != t->bt[i]) t->stiffly_accurate = PETSC_FALSE; in TSARKIMEXRegister()
1203 …pe name, PetscInt order, PetscInt s, const PetscReal At[], const PetscReal bt[], const PetscReal c… in TSDIRKRegister() argument
1206 …PetscCall(TSARKIMEXRegister(name, order, s, At, bt, ct, NULL, NULL, NULL, bembedt, NULL, pinterp, … in TSDIRKRegister()
1252 for (j = 0; j < s; j++) w[j] = h * tab->bt[j]; in TSEvaluateStep_ARKIMEX()
[all …]
|
| H A D | arkimex.h | 11 PetscReal *At, *bt, *ct; /* Stiff tableau */ member
|
| H A D | fsarkimex.c | 81 PetscScalar *bt = ark->work, *b = ark->work + s; in TSExtrapolate_ARKIMEX_FastSlowSplit() local
90 for (i = 0; i < s; i++) bt[i] = b[i] = 0; in TSExtrapolate_ARKIMEX_FastSlowSplit()
93 bt[i] += h * Bt[i * pinterp + j] * tt; in TSExtrapolate_ARKIMEX_FastSlowSplit()
99 PetscCall(VecMAXPY(X, s, bt, ark->YdotI_prev)); in TSExtrapolate_ARKIMEX_FastSlowSplit()
155 for (j = 0; j < s; j++) w[j] = h * tab->bt[j]; in TSEvaluateStep_ARKIMEX_FastSlowSplit()
190 for (j = 0; j < s; j++) w[j] = h * (tab->bembedt[j] - tab->bt[j]); in TSEvaluateStep_ARKIMEX_FastSlowSplit()
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex12_p4est_solve_fas.out | 96 type: bt
127 type: bt
154 type: bt
|
| H A D | ex12_p4est_exact_q2_nonconformal_serial.out | 56 type: bt
87 type: bt
|
| H A D | ex12_p4est_exact_q2_conformal_serial.out | 56 type: bt
87 type: bt
|
| H A D | ex5_5_aspin.out | 51 type: bt
71 type: bt
|
| H A D | ex12_p4est_exact_q2_nonconformal_parallel_parmetis.out | 56 type: bt
87 type: bt
|
| H A D | ex12_p4est_full_q2_conformal_serial_plex.out | 54 type: bt
85 type: bt
|
| H A D | ex12_p4est_full_q2_nonconformal_serial.out | 57 type: bt
89 type: bt
|
| H A D | ex12_p4est_exact_q2_conformal_parallel.out | 56 type: bt
87 type: bt
|
| /petsc/src/dm/impls/da/ |
| H A D | dacreate.c | 7 DMBoundaryType bt = DM_BOUNDARY_NONE; in DMSetFromOptions_DA() local
38 … every direction", "DMDASetBoundaryType", DMBoundaryTypes, (PetscEnum)bt, (PetscEnum *)&bt, &flg)); in DMSetFromOptions_DA()
39 if (flg) PetscCall(DMDASetBoundaryType(da, bt, bt, bt)); in DMSetFromOptions_DA()
|
| /petsc/src/ts/impls/rosw/ |
| H A D | rosw.c | 37 PetscReal *bt; /* Step completion table in transformed variables */ member
763 PetscCall(PetscFree5(t->At, t->bt, t->GammaInv, t->GammaZeroDiag, t->GammaExplicitCorr)); in TSRosWRegisterDestroy()
851 …PetscCall(PetscMalloc5(s * s, &t->At, s, &t->bt, s * s, &t->GammaInv, s, &t->GammaZeroDiag, s * s,… in TSRosWRegister()
922 t->bt[i] = 0; in TSRosWRegister()
923 for (j = 0; j < s; j++) t->bt[i] += t->b[j] * t->GammaInv[j * s + i]; in TSRosWRegister()
1099 for (i = 0; i < tab->s; i++) w[i] = tab->bt[i]; in TSEvaluateStep_RosW()
1111 for (i = 0; i < tab->s; i++) w[i] = tab->bembedt[i] - tab->bt[i]; in TSEvaluateStep_RosW()
1246 PetscScalar *bt; in TSInterpolate_RosW() local
1268 PetscCall(PetscMalloc1(s, &bt)); in TSInterpolate_RosW()
1269 for (i = 0; i < s; i++) bt[i] = 0; in TSInterpolate_RosW()
[all …]
|
| /petsc/src/dm/label/ |
| H A D | dmlabel.c | 46 (*label)->bt = NULL; in DMLabelCreate()
102 if (label->bt) { in DMLabelMakeValid_Private()
106 PetscCall(PetscBTSet(label->bt, point - label->pStart)); in DMLabelMakeValid_Private()
431 …if (label->bt) PetscCall(PetscViewerASCIIPrintf(viewer, " Index has been calculated in [%" PetscI… in DMLabelView_Concrete_Ascii()
537 PetscCall(PetscBTDestroy(&label->bt)); in DMLabelReset()
621 (*labelnew)->bt = NULL; in DMLabelDuplicate()
799 PetscCall(PetscBTCreate(pEnd - pStart, &label->bt)); in DMLabelCreateIndex()
811 PetscCall(PetscBTSet(label->bt, point - pStart)); in DMLabelCreateIndex()
837 PetscCall(PetscBTDestroy(&label->bt)); in DMLabelDestroyIndex()
933 …*contains = point >= pStart && point < pEnd && (PetscBTLookup(label->bt, point - label->pStart) ? … in DMLabelHasPoint()
[all …]
|