| /petsc/src/ts/impls/implicit/discgrad/ |
| H A D | tsdiscgrad.c | 21 Vec X0, X, Xdot; member
29 static PetscErrorCode TSDiscGradGetX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSDiscGradGetX0AndXdot() argument
34 if (X0) { in TSDiscGradGetX0AndXdot()
35 if (dm && dm != ts->dm) PetscCall(DMGetNamedGlobalVector(dm, "TSDiscGrad_X0", X0)); in TSDiscGradGetX0AndXdot()
36 else *X0 = ts->vec_sol; in TSDiscGradGetX0AndXdot()
45 static PetscErrorCode TSDiscGradRestoreX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSDiscGradRestoreX0AndXdot() argument
48 if (X0) { in TSDiscGradRestoreX0AndXdot()
49 if (dm && dm != ts->dm) PetscCall(DMRestoreNamedGlobalVector(dm, "TSDiscGrad_X0", X0)); in TSDiscGradRestoreX0AndXdot()
66 Vec X0, Xdot, X0_c, Xdot_c; in DMRestrictHook_TSDiscGrad() local
69 PetscCall(TSDiscGradGetX0AndXdot(ts, fine, &X0, &Xdot)); in DMRestrictHook_TSDiscGrad()
[all …]
|
| /petsc/src/ts/impls/implicit/alpha/ |
| H A D | alpha1.c | 23 Vec X0, Xa, X1; member
44 PetscCall(TSResizeRegisterVec(ts, "ts:theta:X0", th->X0)); in TSResizeRegister_Alpha()
48 PetscCall(TSResizeRetrieveVec(ts, "ts:theta:X0", &th->X0)); in TSResizeRegister_Alpha()
49 PetscCall(PetscObjectReference((PetscObject)th->X0)); in TSResizeRegister_Alpha()
75 Vec X0 = th->X0, V0 = th->V0; in TSAlpha_StageVecs() local
83 PetscCall(VecWAXPY(V1, -1.0, X0, X1)); in TSAlpha_StageVecs()
86 PetscCall(VecWAXPY(Xa, -1.0, X0, X1)); in TSAlpha_StageVecs()
87 PetscCall(VecAYPX(Xa, Alpha_f, X0)); in TSAlpha_StageVecs()
118 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() local
122 PetscCall(VecDuplicate(X0, &X1)); in TSAlpha_Restart()
[all …]
|
| H A D | alpha2.c | 24 Vec X0, Xa, X1; member
81 if (th->predictor) PetscCall(th->predictor(ts, th->X0, th->V0, th->A0, X1, th->predictor_ctx)); in TSAlpha_ApplyPredictor()
82 else PetscCall(VecCopy(th->X0, X1)); in TSAlpha_ApplyPredictor()
109 Vec X0 = th->X0, V0 = th->V0, A0 = th->A0; in TSAlpha_StageVecs() local
118 PetscCall(VecWAXPY(A1, -1.0, X0, X1)); in TSAlpha_StageVecs()
125 PetscCall(VecWAXPY(Xa, -1.0, X0, X1)); in TSAlpha_StageVecs()
126 PetscCall(VecAYPX(Xa, Alpha_f, X0)); in TSAlpha_StageVecs()
160 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() local
165 PetscCall(VecDuplicate(X0, &X1)); in TSAlpha_Restart()
179 PetscCall(VecCopy(X0, th->X0)); in TSAlpha_Restart()
[all …]
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 13 …Vec X0, X, Xdot; /* Storage for u^n, u^n + dt a_{11} k_1, and time derivative u^{n+1}_t */ member
45 static PetscErrorCode TSThetaGetX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSThetaGetX0AndXdot() argument
50 if (X0) { in TSThetaGetX0AndXdot()
51 if (dm && dm != ts->dm) PetscCall(DMGetNamedGlobalVector(dm, "TSTheta_X0", X0)); in TSThetaGetX0AndXdot()
52 else *X0 = ts->vec_sol; in TSThetaGetX0AndXdot()
61 static PetscErrorCode TSThetaRestoreX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSThetaRestoreX0AndXdot() argument
64 if (X0) { in TSThetaRestoreX0AndXdot()
65 if (dm && dm != ts->dm) PetscCall(DMRestoreNamedGlobalVector(dm, "TSTheta_X0", X0)); in TSThetaRestoreX0AndXdot()
82 Vec X0, Xdot, X0_c, Xdot_c; in DMRestrictHook_TSTheta() local
85 PetscCall(TSThetaGetX0AndXdot(ts, fine, &X0, &Xdot)); in DMRestrictHook_TSTheta()
[all …]
|
| /petsc/src/ts/impls/mimex/ |
| H A D | mimex.c | 15 static PetscErrorCode TSMimexGetX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSMimexGetX0AndXdot() argument
20 if (X0) { in TSMimexGetX0AndXdot()
21 if (dm && dm != ts->dm) PetscCall(DMGetNamedGlobalVector(dm, "TSMimex_X0", X0)); in TSMimexGetX0AndXdot()
22 else *X0 = ts->vec_sol; in TSMimexGetX0AndXdot()
31 static PetscErrorCode TSMimexRestoreX0AndXdot(TS ts, DM dm, Vec *X0, Vec *Xdot) in TSMimexRestoreX0AndXdot() argument
34 if (X0) in TSMimexRestoreX0AndXdot()
35 if (dm && dm != ts->dm) PetscCall(DMRestoreNamedGlobalVector(dm, "TSMimex_X0", X0)); in TSMimexRestoreX0AndXdot()
65 Vec X0, Xdot; in SNESTSFormFunction_Mimex() local
70 PetscCall(TSMimexGetX0AndXdot(ts, dm, &X0, &Xdot)); in SNESTSFormFunction_Mimex()
71 PetscCall(VecAXPBYPCZ(Xdot, -shift, shift, 0, X0, x)); in SNESTSFormFunction_Mimex()
[all …]
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | mrk.c | 25 PetscCall(VecDestroy(&rk->X0)); in TSReset_RK_MultirateNonsplit()
47 PetscCall(VecCopy(rk->X0, X)); in TSInterpolate_RK_MultirateNonsplit()
136 PetscCall(VecCopy(ts->vec_sol, rk->X0)); in TSStep_RK_MultirateNonsplit()
185 PetscCall(VecDuplicate(ts->vec_sol, &rk->X0)); in TSSetUp_RK_MultirateNonsplit()
224 PetscCall(VecDestroy(&rk->X0)); in TSReset_RK_MultirateSplit()
265 PetscCall(VecISCopy(rk->X0, rk->is_slow, SCATTER_REVERSE, Xslow)); in TSInterpolate_RK_MultirateSplit()
358 PetscCall(VecISCopy(rk->X0, rk->is_fast, SCATTER_FORWARD, Xfast)); in TSStepRefine_RK_MultirateSplit()
383 PetscCall(VecCopy(ts->vec_sol, rk->X0)); in TSStep_RK_MultirateSplit()
421 Vec X0; in TSSetUp_RK_MultirateSplit() local
431 PetscCall(VecDuplicate(ts->vec_sol, &X0)); in TSSetUp_RK_MultirateSplit()
[all …]
|
| H A D | rk.h | 23 Vec X0; member
|
| /petsc/src/tao/pde_constrained/impls/lcl/ |
| H A D | lcl.h | 46 Vec X0; member
|
| H A D | lcl.c | 18 PetscCall(VecDestroy(&lclP->X0)); in TaoDestroy_LCL()
104 PetscCall(VecDuplicate(tao->solution, &lclP->X0)); in TaoSetup_LCL()
355 PetscCall(LCLGather(lclP, lclP->U0, lclP->V0, lclP->X0)); in TaoSolve_LCL()
391 …PetscCall(TaoComputeJacobianState(tao, lclP->X0, tao->jacobian_state, tao->jacobian_state_pre, tao… in TaoSolve_LCL()
392 PetscCall(TaoComputeJacobianDesign(tao, lclP->X0, tao->jacobian_design)); in TaoSolve_LCL()
476 …PetscCall(TaoComputeJacobianState(tao, lclP->X0, tao->jacobian_state, tao->jacobian_state_pre, tao… in TaoSolve_LCL()
477 PetscCall(TaoComputeJacobianDesign(tao, lclP->X0, tao->jacobian_design)); in TaoSolve_LCL()
|
| /petsc/src/tao/leastsquares/tutorials/matlab/ |
| H A D | ProblemInitialize.m | 13 X0 = dfoxs(n,nprob,factor^factor_power)'; variable
|
| H A D | ProblemFinalize.m | 35 fminsearch(func,X0,options);
|
| /petsc/src/ts/tutorials/multirate/ |
| H A D | ex5.c | 617 Vec X, X0, R; in main() local
701 PetscCall(VecDuplicate(X, &X0)); in main()
752 PetscCall(FVSample(&ctx, da, 0, X0)); in main()
753 …PetscCall(FVRHSFunction(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, only used to… in main()
754 …PetscCall(VecCopy(X0, X)); /* The function value was not used so we set… in main()
769 PetscCall(VecSum(X0, &mass_initial)); in main()
783 PetscCall(VecDuplicate(X0, &XR)); in main()
796 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
826 PetscCall(VecDestroy(&X0)); in main()
|
| H A D | ex7.c | 665 Vec X, X0, R; in main() local
724 PetscCall(VecDuplicate(X, &X0)); in main()
758 PetscCall(FVSample(&ctx, da, 0, X0)); in main()
759 …PetscCall(FVRHSFunction(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, only used to… in main()
760 …PetscCall(VecCopy(X0, X)); /* The function value was not used so we set… in main()
776 PetscCall(DMDAVecGetArrayRead(da, X0, (void *)&ptr_X0)); in main()
791 PetscCall(DMDAVecRestoreArrayRead(da, X0, (void *)&ptr_X0)); in main()
812 PetscCall(VecDuplicate(X0, &XR)); in main()
829 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
856 PetscCall(VecDestroy(&X0)); in main()
|
| H A D | ex6.c | 746 Vec X, X0, R; in main() local
821 PetscCall(VecDuplicate(X, &X0)); in main()
870 PetscCall(FVSample_2WaySplit(&ctx, da, 0, X0)); in main()
871 …PetscCall(FVRHSFunction_2WaySplit(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, on… in main()
872 …PetscCall(VecCopy(X0, X)); /* The function value was not used… in main()
889 PetscCall(DMDAVecGetArrayRead(da, X0, (void *)&ptr_X0)); in main()
904 PetscCall(DMDAVecRestoreArrayRead(da, X0, (void *)&ptr_X0)); in main()
926 PetscCall(VecDuplicate(X0, &XR)); in main()
945 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
971 PetscCall(VecDestroy(&X0)); in main()
|
| H A D | ex4.c | 1187 Vec X, X0, R; in main() local
1269 PetscCall(VecDuplicate(X, &X0)); in main()
1313 PetscCall(FVSample_2WaySplit(&ctx, da, 0, X0)); in main()
1314 …PetscCall(FVRHSFunction_2WaySplit(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, on… in main()
1315 …PetscCall(VecCopy(X0, X)); /* The function value was not used… in main()
1332 PetscCall(DMDAVecGetArrayRead(da, X0, (void *)&ptr_X0)); in main()
1347 PetscCall(DMDAVecRestoreArrayRead(da, X0, (void *)&ptr_X0)); in main()
1369 PetscCall(VecDuplicate(X0, &XR)); in main()
1388 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
1416 PetscCall(VecDestroy(&X0)); in main()
|
| H A D | ex8.c | 1088 Vec X, X0, R; in main() local
1163 PetscCall(VecDuplicate(X, &X0)); in main()
1232 PetscCall(FVSample_3WaySplit(&ctx, da, 0, X0)); in main()
1233 …PetscCall(FVRHSFunction_3WaySplit(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, on… in main()
1234 …PetscCall(VecCopy(X0, X)); /* The function value was not used… in main()
1252 PetscCall(DMDAVecGetArrayRead(da, X0, (void *)&ptr_X0)); in main()
1272 PetscCall(DMDAVecRestoreArrayRead(da, X0, (void *)&ptr_X0)); in main()
1294 PetscCall(VecDuplicate(X0, &XR)); in main()
1315 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
1341 PetscCall(VecDestroy(&X0)); in main()
|
| /petsc/src/snes/impls/ms/ |
| H A D | ms.c | 270 Vec X0 = snes->work[0]; in SNESMSStep_Basic() local
273 PetscCall(VecCopy(X, X0)); in SNESMSStep_Basic()
277 PetscCall(VecAYPX(X, -alpha[i] * h, X0)); in SNESMSStep_Basic()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex43.c | 53 PetscErrorCode Predictor(TS ts, Vec X0, Vec V0, Vec A0, Vec X1, PetscCtx ctx) in Predictor() argument
61 PetscCall(VecCopy(X0, X1)); in Predictor()
|
| H A D | ex9.c | 1394 Vec X, X0, R; in main() local
1484 PetscCall(VecDuplicate(X, &X0)); in main()
1499 PetscCall(FVSample(&ctx, da, 0, X0)); in main()
1500 …PetscCall(FVRHSFunction(ts, 0, X0, X, (void *)&ctx)); /* Initial function evaluation, only used to… in main()
1501 …PetscCall(VecCopy(X0, X)); /* The function value was not used so we set… in main()
1521 if (draw & 0x1) PetscCall(VecView(X0, PETSC_VIEWER_DRAW_WORLD)); in main()
1547 PetscCall(VecDestroy(&X0)); in main()
|
| /petsc/src/ts/tutorials/hamiltonian/ |
| H A D | ex3.c | 34 X0, enumerator
760 constants[X0] = param->x0; in InitializeConstants()
|
| H A D | ex2.c | 88 X0, enumerator
1445 constants[X0] = param->x0; in InitializeConstants()
|
| H A D | ex4.c | 109 X0, enumerator
2059 constants[X0] = param->x0; in InitializeConstants()
|
| /petsc/include/ |
| H A D | petscts.h | 1577 PETSC_EXTERN_TYPEDEF typedef PetscErrorCode TSAlpha2PredictorFn(TS ts, Vec X0, Vec V0, Vec A0, Vec …
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | TAO.pyx | 61 CONVERGED_GTTOL = TAO_CONVERGED_GTTOL # ||g(X)||/||g(X0)|| < gttol
|