| /petsc/src/ts/impls/explicit/rk/ |
| H A D | mrk.c | 67 PetscCall(VecDuplicate(ts->vec_sol, &vec_fast)); in TSStepRefine_RK_MultirateNonsplit() 76 PetscCall(VecCopy(ts->vec_sol, vec_fast)); in TSStepRefine_RK_MultirateNonsplit() 85 PetscCall(VecCopy(ts->vec_sol, vec_fast)); in TSStepRefine_RK_MultirateNonsplit() 89 PetscCall(VecISCopy(ts->vec_sol, rk->is_fast, SCATTER_FORWARD, subvec_fast)); in TSStepRefine_RK_MultirateNonsplit() 94 PetscCall(VecDuplicateVecs(ts->vec_sol, s, &YdotRHS_copy)); in TSStepRefine_RK_MultirateNonsplit() 134 PetscCall(VecDuplicate(ts->vec_sol, &stage_slow)); in TSStep_RK_MultirateNonsplit() 135 PetscCall(VecDuplicate(ts->vec_sol, &sol_slow)); in TSStep_RK_MultirateNonsplit() 136 PetscCall(VecCopy(ts->vec_sol, rk->X0)); in TSStep_RK_MultirateNonsplit() 140 PetscCall(VecCopy(ts->vec_sol, Y[i])); in TSStep_RK_MultirateNonsplit() 155 PetscCall(VecISCopy(ts->vec_sol, is_slow, SCATTER_FORWARD, subvec_slow)); in TSStep_RK_MultirateNonsplit() [all …]
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 52 else *X0 = ts->vec_sol; in TSThetaGetX0AndXdot() 132 … PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * (1.0 - th->Theta), ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral() 134 PetscCall(TSComputeRHSFunction(quadts, ts->ptime, ts->vec_sol, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral() 135 PetscCall(VecAXPY(quadts->vec_sol, th->time_step0 * th->Theta, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral() 138 PetscCall(VecAXPY(quadts->vec_sol, th->time_step0, ts->vec_costintegrand)); in TSThetaEvaluateCostIntegral() 150 PetscCall(VecCopy(quadts->vec_sol, th->VecCostIntegral0)); in TSForwardCostIntegral_Theta() 207 PetscCall(VecCopy(ts->vec_sol, th->X0)); in TSStep_Theta() 217 if (!th->affine) PetscCall(VecDuplicate(ts->vec_sol, &th->affine)); in TSStep_Theta() 229 PetscCall(VecCopy(th->X, ts->vec_sol)); in TSStep_Theta() 232 …if (th->Theta == 1.0) PetscCall(VecCopy(th->X, ts->vec_sol)); /* BEULER, stage alread… in TSStep_Theta() [all …]
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | fsarkimex.c | 145 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX_FastSlowSplit() 159 } else PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX_FastSlowSplit() 165 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX_FastSlowSplit() 180 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX_FastSlowSplit() 233 PetscCall(VecGetSubVector(ts->vec_sol, ark->is_fast, &Xfast)); in TSStep_ARKIMEX_FastSlowSplit() 237 PetscCall(VecRestoreSubVector(ts->vec_sol, ark->is_fast, &Xfast)); in TSStep_ARKIMEX_FastSlowSplit() 238 PetscCall(VecGetSubVector(ts->vec_sol, ark->is_slow, &Xslow)); in TSStep_ARKIMEX_FastSlowSplit() 239 PetscCall(VecRestoreSubVector(ts->vec_sol, ark->is_fast, &Xslow)); in TSStep_ARKIMEX_FastSlowSplit() 259 PetscCall(TSSetSolution(ts_start, ts->vec_sol)); in TSStep_ARKIMEX_FastSlowSplit() 270 PetscCall(TSSolve(ts_start, ts->vec_sol)); in TSStep_ARKIMEX_FastSlowSplit() [all …]
|
| H A D | arkimex.c | 1251 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX() 1262 } else PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX() 1268 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX() 1279 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_ARKIMEX() 1307 PetscCall(VecDuplicate(ts->vec_sol, &Udot)); in TSARKIMEXTestMassIdentity() 1308 PetscCall(VecDuplicate(ts->vec_sol, &Y1)); in TSARKIMEXTestMassIdentity() 1309 PetscCall(VecDuplicate(ts->vec_sol, &Y2)); in TSARKIMEXTestMassIdentity() 1310 PetscCall(TSComputeIFunction(ts, ts->ptime, ts->vec_sol, Udot, Y1, ark->imex)); in TSARKIMEXTestMassIdentity() 1312 PetscCall(TSComputeIFunction(ts, ts->ptime, ts->vec_sol, Udot, Y2, ark->imex)); in TSARKIMEXTestMassIdentity() 1348 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &ark->Y_prev)); in TSStep_ARKIMEX() [all …]
|
| /petsc/src/ts/impls/implicit/alpha/ |
| H A D | alpha2.c | 160 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() 227 PetscCall(VecCopy(ts->vec_sol, th->X0)); in TSAlpha_Restart() 248 PetscCall(VecCopy(ts->vec_sol, th->X0)); in TSStep_Alpha() 269 PetscCall(VecCopy(th->X1, ts->vec_sol)); in TSStep_Alpha() 274 PetscCall(VecCopy(th->X0, ts->vec_sol)); in TSStep_Alpha() 359 PetscCall(VecCopy(th->X0, ts->vec_sol)); in TSRollBack_Alpha() 448 PetscCall(VecDuplicate(ts->vec_sol, &th->X0)); in TSSetUp_Alpha() 449 PetscCall(VecDuplicate(ts->vec_sol, &th->Xa)); in TSSetUp_Alpha() 450 PetscCall(VecDuplicate(ts->vec_sol, &th->X1)); in TSSetUp_Alpha() 451 PetscCall(VecDuplicate(ts->vec_sol, &th->V0)); in TSSetUp_Alpha() [all …]
|
| H A D | alpha1.c | 118 Vec X0 = ts->vec_sol, X1, X2 = th->X1; in TSAlpha_Restart() 173 PetscCall(VecCopy(ts->vec_sol, th->X0)); in TSAlpha_Restart() 191 PetscCall(VecCopy(ts->vec_sol, th->X0)); in TSStep_Alpha() 211 PetscCall(VecCopy(th->X1, ts->vec_sol)); in TSStep_Alpha() 215 PetscCall(VecCopy(th->X0, ts->vec_sol)); in TSStep_Alpha() 281 PetscCall(VecWAXPY(th->V1, -1.0, th->X0, ts->vec_sol)); in TSInterpolate_Alpha() 283 PetscCall(VecCopy(ts->vec_sol, X)); in TSInterpolate_Alpha() 352 if (!th->X0) PetscCall(VecDuplicate(ts->vec_sol, &th->X0)); in TSSetUp_Alpha() 353 PetscCall(VecDuplicate(ts->vec_sol, &th->Xa)); in TSSetUp_Alpha() 354 PetscCall(VecDuplicate(ts->vec_sol, &th->X1)); in TSSetUp_Alpha() [all …]
|
| /petsc/src/ts/impls/explicit/euler/ |
| H A D | euler.c | 13 Vec solution = ts->vec_sol, update = euler->update; in TSStep_Euler() 51 PetscCall(VecDuplicate(ts->vec_sol, &euler->update)); in TSSetUp_Euler() 93 PetscCall(VecWAXPY(X, -ts->time_step, update, ts->vec_sol)); in TSInterpolate_Euler() 94 PetscCall(VecAXPBY(X, 1.0 - alpha, alpha, ts->vec_sol)); in TSInterpolate_Euler()
|
| /petsc/src/ts/interface/ |
| H A D | tsrhssplit.c | 123 if (!r && ts->vec_sol) { in TSRHSSplitSetRHSFunction() 124 PetscCall(VecGetSubVector(ts->vec_sol, isplit->is, &subvec)); in TSRHSSplitSetRHSFunction() 127 PetscCall(VecRestoreSubVector(ts->vec_sol, isplit->is, &subvec)); in TSRHSSplitSetRHSFunction() 176 if (!r && ts->vec_sol) { in TSRHSSplitSetIFunction() 177 PetscCall(VecGetSubVector(ts->vec_sol, isplit->is, &subvec)); in TSRHSSplitSetIFunction() 180 PetscCall(VecRestoreSubVector(ts->vec_sol, isplit->is, &subvec)); in TSRHSSplitSetIFunction()
|
| H A D | ts.c | 1034 if (!r && !ts->dm && ts->vec_sol) { in TSSetRHSFunction() 1035 PetscCall(VecDuplicate(ts->vec_sol, &ralloc)); in TSSetRHSFunction() 1207 if (!r && !ts->dm && ts->vec_sol) { in TSSetIFunction() 1208 PetscCall(VecDuplicate(ts->vec_sol, &ralloc)); in TSSetIFunction() 1787 if (u) *u = ts->vec_sol; in TS2GetSolution() 1831 PetscCall(DMCreateGlobalVector(ts->dm, &ts->vec_sol)); in TSLoad() 1832 PetscCall(VecLoad(ts->vec_sol, viewer)); in TSLoad() 1970 PetscCall(VecView(ts->vec_sol, viewer)); in TSView() 2279 *v = ts->vec_sol; in TSGetSolution() 2503 if (!ts->vec_sol) { in TSSetUp() [all …]
|
| /petsc/src/ts/impls/eimex/ |
| H A D | eimex.c | 88 PetscCall(VecCopy(ts->vec_sol, ext->VecSolPrev)); in TSStep_EIMEX() 105 PetscCall(TSEvaluateStep(ts, ns, ts->vec_sol, NULL)); /*update ts solution */ in TSStep_EIMEX() 110 …PetscCall(TSErrorWeightedNorm(ts, ts->vec_sol, T[Map(ext->nstages - 1, ext->nstages - 2, ext->nsta… in TSStep_EIMEX() 119 PetscCall(VecDuplicateVecs(ts->vec_sol, (1 + ext->nstages) * ext->nstages / 2, &ext->T)); in TSStep_EIMEX() 135 PetscCall(TSEvaluateStep(ts, ext->nstages, ts->vec_sol, NULL)); in TSStep_EIMEX() 340 …PetscCall(VecDuplicateVecs(ts->vec_sol, (1 + ext->nstages) * ext->nstages / 2, &ext->T)); /* full … in TSSetUp_EIMEX() 341 PetscCall(VecDuplicate(ts->vec_sol, &ext->YdotI)); in TSSetUp_EIMEX() 342 PetscCall(VecDuplicate(ts->vec_sol, &ext->YdotRHS)); in TSSetUp_EIMEX() 343 PetscCall(VecDuplicate(ts->vec_sol, &ext->Ydot)); in TSSetUp_EIMEX() 344 PetscCall(VecDuplicate(ts->vec_sol, &ext->VecSolPrev)); in TSSetUp_EIMEX() [all …]
|
| /petsc/src/ts/impls/implicit/irk/ |
| H A D | irk.c | 291 PetscCall(VecCopy(ts->vec_sol, U)); in TSEvaluateStep_IRK() 302 PetscCall(VecCopy(irk->U0, ts->vec_sol)); in TSRollBack_IRK() 320 if (!ts->steprollback) PetscCall(VecCopy(ts->vec_sol, irk->U0)); in TSStep_IRK() 321 PetscCall(VecGetBlockSize(ts->vec_sol, &bs)); in TSStep_IRK() 322 …for (i = 0; i < nstages; i++) PetscCall(VecStrideScatter(ts->vec_sol, i * bs, irk->Z, INSERT_VALUE… in TSStep_IRK() 326 PetscCall(VecCopy(ts->vec_sol, irk->U)); in TSStep_IRK() 340 PetscCall(TSEvaluateStep_IRK(ts, irk->order, ts->vec_sol, NULL)); in TSStep_IRK() 579 if (!irk->Y) PetscCall(VecDuplicateVecs(ts->vec_sol, irk->nstages, &irk->Y)); in TSSetUp_IRK() 580 if (!irk->YdotI) PetscCall(VecDuplicateVecs(ts->vec_sol, irk->nstages, &irk->YdotI)); in TSSetUp_IRK() 581 if (!irk->Ydot) PetscCall(VecDuplicate(ts->vec_sol, &irk->Ydot)); in TSSetUp_IRK() [all …]
|
| /petsc/src/ts/impls/pseudo/ |
| H A D | posindep.c | 96 PetscCall(SNESSolve(ts->snes, NULL, ts->vec_sol)); in TSStep_Pseudo() 102 PetscCall(TSPostStage(ts, ts->ptime + ts->time_step, 0, &ts->vec_sol)); in TSStep_Pseudo() 104 PetscCall(TSAdaptCheckStage(adapt, ts, ts->ptime + ts->time_step, ts->vec_sol, &accept)); in TSStep_Pseudo() 121 PetscCall(VecCopy(ts->vec_sol0, ts->vec_sol)); in TSStep_Pseudo() 130 PetscCall(TSPseudoComputeFunction(ts, ts->vec_sol, NULL, &fnorm)); in TSStep_Pseudo() 259 PetscCall(VecDuplicate(ts->vec_sol, &pseudo->func)); in TSSetUp_Pseudo() 260 PetscCall(VecDuplicate(ts->vec_sol, &pseudo->xdot)); in TSSetUp_Pseudo() 270 PetscCall(TSPseudoComputeFunction(ts, ts->vec_sol, NULL, NULL)); in TSPseudoMonitorDefault() 340 PetscCall(TSPseudoComputeFunction(ts, ts->vec_sol, NULL, &fnorm)); in TSPseudoTimeStepDefault()
|
| /petsc/src/ksp/ksp/impls/preonly/ |
| H A D | preonly.c | 27 PetscCall(KSP_PCApply(ksp, ksp->vec_rhs, ksp->vec_sol)); in KSPSolve_PREONLY() 30 PetscCall(VecFlag(ksp->vec_sol, pcreason)); in KSPSolve_PREONLY() 44 PetscCall(KSP_MatMult(ksp, A, ksp->vec_sol, v)); in KSPSolve_PREONLY()
|
| /petsc/src/ksp/ksp/interface/ |
| H A D | itfunc.c | 359 if (!ksp->vec_sol) PetscCall(DMCreateGlobalVector(ksp->dm, &ksp->vec_sol)); in KSPSetUp() 371 …"KSP callback initial guess", (*kdm->ops->computeinitialguess)(ksp, ksp->vec_sol, kdm->initialgues… in KSPSetUp() 386 PetscCall(PetscLogEventBegin(KSP_SetUp, ksp, ksp->vec_rhs, ksp->vec_sol, 0)); in KSPSetUp() 426 PetscCall(PetscLogEventEnd(KSP_SetUp, ksp, ksp->vec_rhs, ksp->vec_sol, 0)); in KSPSetUp() 685 PetscCall(VecGetSize(ksp->vec_sol, &n)); in KSPViewEigenvalues_Internal() 759 PetscCall(KSP_MatMult(ksp, A, ksp->vec_sol, t)); in KSPViewFinalResidual_Internal() 826 PetscCall(VecDestroy(&ksp->vec_sol)); in KSPSolve_Private() 827 ksp->vec_sol = x; in KSPSolve_Private() 832 if (ksp->presolve) PetscCall((*ksp->presolve)(ksp, ksp->vec_rhs, ksp->vec_sol, ksp->prectx)); in KSPSolve_Private() 846 PetscCall(PetscObjectStateGet((PetscObject)ksp->vec_sol, &ostate)); in KSPSolve_Private() [all …]
|
| /petsc/src/snes/linesearch/interface/ |
| H A D | linesearch.c | 182 linesearch->vec_sol = NULL; in SNESLineSearchCreate() 233 …if (!linesearch->vec_sol_new) PetscCall(VecDuplicate(linesearch->vec_sol, &linesearch->vec_sol_new… in SNESLineSearchSetUp() 234 …if (!linesearch->vec_func_new) PetscCall(VecDuplicate(linesearch->vec_sol, &linesearch->vec_func_n… in SNESLineSearchSetUp() 633 linesearch->vec_sol = X; in SNESLineSearchApply() 1417 PetscCall(VecNorm(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms() 1419 …PetscCall((*linesearch->ops->vinorm)(snes, linesearch->vec_func, linesearch->vec_sol, &linesearch-… in SNESLineSearchComputeNorms() 1422 PetscCall(VecNormBegin(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms() 1425 PetscCall(VecNormEnd(linesearch->vec_sol, NORM_2, &linesearch->xnorm)); in SNESLineSearchComputeNorms() 1492 *X = linesearch->vec_sol; in SNESLineSearchGetVecs() 1536 linesearch->vec_sol = X; in SNESLineSearchSetVecs() [all …]
|
| /petsc/src/ksp/ksp/impls/bcgs/ |
| H A D | bcgs.c | 27 X = ksp->vec_sol; in KSPSolve_BCGS() 152 PetscCall(KSP_PCApply(ksp, ksp->vec_sol, v)); in KSPBuildSolution_BCGS() 157 PetscCall(VecCopy(ksp->vec_sol, v)); in KSPBuildSolution_BCGS() 159 } else *V = ksp->vec_sol; in KSPBuildSolution_BCGS()
|
| /petsc/src/ts/impls/glee/ |
| H A D | glee.c | 471 } else PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_GLEE() 545 PetscCall(VecCopy(i > 0 ? YStage[i - 1] : ts->vec_sol, YStage[i])); in TSStep_GLEE() 559 PetscCall(TSEvaluateStep(ts, tab->order, ts->vec_sol, NULL)); in TSStep_GLEE() 575 …PetscCall(PetscObjectComposedDataSetReal((PetscObject)ts->vec_sol, explicit_stage_time_id, ts->pti… in TSStep_GLEE() 579 PetscCall(VecMAXPY(ts->vec_sol, r, wr, X)); in TSStep_GLEE() 748 PetscCall(VecDuplicateVecs(ts->vec_sol, r, &glee->Y)); in TSSetUp_GLEE() 749 PetscCall(VecDuplicateVecs(ts->vec_sol, r, &glee->X)); in TSSetUp_GLEE() 750 PetscCall(VecDuplicateVecs(ts->vec_sol, s, &glee->YStage)); in TSSetUp_GLEE() 751 PetscCall(VecDuplicateVecs(ts->vec_sol, s, &glee->YdotStage)); in TSSetUp_GLEE() 752 PetscCall(VecDuplicate(ts->vec_sol, &glee->Ydot)); in TSSetUp_GLEE() [all …]
|
| /petsc/src/ts/impls/mimex/ |
| H A D | mimex.c | 22 else *X0 = ts->vec_sol; in TSMimexGetX0AndXdot() 149 Vec sol = ts->vec_sol, update = mimex->update; in TSStep_Mimex_Split() 216 Vec sol = ts->vec_sol; in TSStep_Mimex_Implicit() 254 PetscCall(VecDuplicate(ts->vec_sol, &mimex->update)); in TSSetUp_Mimex() 255 PetscCall(VecDuplicate(ts->vec_sol, &mimex->Xdot)); in TSSetUp_Mimex() 306 PetscCall(VecAXPBY(ts->vec_sol, 1.0 - alpha, alpha, X)); in TSInterpolate_Mimex()
|
| /petsc/src/ts/impls/bdf/ |
| H A D | bdf.c | 225 PetscCall(TSBDF_Advance(ts, ts->ptime, ts->vec_sol)); in TSBDF_Restart() 260 PetscCall(TSBDF_Advance(ts, ts->ptime, ts->vec_sol)); in TSStep_BDF() 272 else if (!accept) PetscCall(VecCopy(ts->vec_sol, bdf->work[0])); in TSStep_BDF() 289 PetscCall(VecCopy(bdf->work[0], ts->vec_sol)); in TSStep_BDF() 443 if (!bdf->work[i]) PetscCall(VecDuplicate(ts->vec_sol, &bdf->work[i])); in TSSetUp_BDF() 445 …if (i && bdf->transientvar && !bdf->tvwork[i]) PetscCall(VecDuplicate(ts->vec_sol, &bdf->tvwork[i]… in TSSetUp_BDF() 448 PetscCall(VecDuplicate(ts->vec_sol, &bdf->vec_dot)); in TSSetUp_BDF() 449 PetscCall(VecDuplicate(ts->vec_sol, &bdf->vec_wrk)); in TSSetUp_BDF() 450 PetscCall(VecDuplicate(ts->vec_sol, &bdf->vec_lte)); in TSSetUp_BDF()
|
| /petsc/src/ts/impls/multirate/ |
| H A D | mprk.c | 636 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_MPRK() 657 PetscCall(VecCopy(ts->vec_sol, Y[i])); in TSStep_MPRK() 707 PetscCall(TSEvaluateStep(ts, tab->order, ts->vec_sol, NULL)); in TSStep_MPRK() 728 PetscCall(VecCopy(ts->vec_sol, X)); in TSEvaluateStep_MPRKSPLIT() 803 PetscCall(VecCopy(ts->vec_sol, Y[i])); in TSStep_MPRKSPLIT() 894 PetscCall(TSEvaluateStep(ts, tab->order, ts->vec_sol, NULL)); in TSStep_MPRKSPLIT() 970 PetscCall(VecDuplicateVecs(ts->vec_sol, tab->s, &mprk->Y)); in TSMPRKTableauSetUp() 981 PetscCall(VecGetSubVector(ts->vec_sol, mprk->is_slow, &YdotRHS_slow)); in TSMPRKTableauSetUp() 983 PetscCall(VecRestoreSubVector(ts->vec_sol, mprk->is_slow, &YdotRHS_slow)); in TSMPRKTableauSetUp() 985 PetscCall(VecGetSubVector(ts->vec_sol, mprk->is_slowbuffer, &YdotRHS_slowbuffer)); in TSMPRKTableauSetUp() [all …]
|
| /petsc/src/ts/impls/implicit/discgrad/ |
| H A D | tsdiscgrad.c | 36 else *X0 = ts->vec_sol; in TSDiscGradGetX0AndXdot() 112 if (!dg->X) PetscCall(VecDuplicate(ts->vec_sol, &dg->X)); in TSSetUp_DiscGrad() 113 if (!dg->X0) PetscCall(VecDuplicate(ts->vec_sol, &dg->X0)); in TSSetUp_DiscGrad() 114 if (!dg->Xdot) PetscCall(VecDuplicate(ts->vec_sol, &dg->Xdot)); in TSSetUp_DiscGrad() 199 PetscCall(VecCopy(ts->vec_sol, dg->X)); in TSInterpolate_DiscGrad() 230 if (!ts->steprollback) PetscCall(VecCopy(ts->vec_sol, dg->X0)); in TSStep_DiscGrad() 246 PetscCall(VecAXPY(ts->vec_sol, ts->time_step, dg->Xdot)); in TSStep_DiscGrad() 250 PetscCall(VecCopy(dg->X0, ts->vec_sol)); in TSStep_DiscGrad()
|
| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 154 PetscCall(KSPPGMRESBuildSoln(RS(0), ksp->vec_sol, ksp->vec_sol, ksp, it - 2)); in KSPPGMRESCycle() 179 …PetscCall(KSPInitialResidual(ksp, ksp->vec_sol, VEC_TEMP, VEC_TEMP_MATOP, VEC_VV(0), ksp->vec_rhs)… in KSPSolve_PGMRES() 312 if (!pgmres->sol_temp) PetscCall(VecDuplicate(ksp->vec_sol, &pgmres->sol_temp)); in KSPBuildSolution_PGMRES() 320 PetscCall(KSPPGMRESBuildSoln(pgmres->nrs, ksp->vec_sol, ptr, ksp, pgmres->it)); in KSPBuildSolution_PGMRES()
|
| /petsc/src/ts/impls/implicit/sundials/ |
| H A D | sundials.c | 131 PetscCall(VecGetArray(ts->vec_sol, &y_data)); in TSStep_Sundials() 133 PetscCall(VecRestoreArray(ts->vec_sol, NULL)); in TSStep_Sundials() 202 PetscCall(VecCopy(cvode->update, ts->vec_sol)); in TSStep_Sundials() 291 PetscCall(VecGetSize(ts->vec_sol, &glosize)); in TSSetUp_Sundials() 292 PetscCall(VecGetLocalSize(ts->vec_sol, &locsize)); in TSSetUp_Sundials() 307 PetscCall(VecGetArray(ts->vec_sol, &parray)); in TSSetUp_Sundials() 310 PetscCall(VecRestoreArray(ts->vec_sol, NULL)); in TSSetUp_Sundials() 312 PetscCall(VecDuplicate(ts->vec_sol, &cvode->update)); in TSSetUp_Sundials() 313 PetscCall(VecDuplicate(ts->vec_sol, &cvode->ydot)); in TSSetUp_Sundials()
|
| /petsc/src/snes/impls/qn/ |
| H A D | qn.c | 83 X = snes->vec_sol; /* solution vector */ in SNESSolve_QN() 262 if (!snes->vec_sol) { in SNESSetUp_QN() 264 PetscCall(DMCreateGlobalVector(dm, &snes->vec_sol)); in SNESSetUp_QN() 289 PetscCall(VecGetLocalSize(snes->vec_sol, &n)); in SNESSetUp_QN() 290 PetscCall(VecGetSize(snes->vec_sol, &N)); in SNESSetUp_QN() 295 PetscCall(MatLMVMAllocate(qn->B, snes->vec_sol, snes->vec_func)); in SNESSetUp_QN()
|
| /petsc/src/ksp/ksp/impls/gcr/ |
| H A D | gcr.c | 35 x = ksp->vec_sol; in KSPSolve_GCR_cycle() 92 x = ksp->vec_sol; in KSPSolve_GCR() 279 x = ksp->vec_sol; in KSPBuildSolution_GCR() 284 *V = ksp->vec_sol; in KSPBuildSolution_GCR()
|