| /petsc/src/ts/interface/sensitivity/ |
| H A D | tssen.c | 295 PetscErrorCode TSSetCostIntegrand(TS ts, PetscInt numcost, Vec costintegral, PetscErrorCode (*rf)(T… in TSSetCostIntegrand() argument 300 …PetscCheck(!ts->numcost || ts->numcost == numcost, PetscObjectComm((PetscObject)ts), PETSC_ERR_USE… in TSSetCostIntegrand() 301 if (!ts->numcost) ts->numcost = numcost; in TSSetCostIntegrand() 309 PetscCall(VecCreateSeq(PETSC_COMM_SELF, numcost, &ts->vec_costintegral)); in TSSetCostIntegrand() 523 for (nadj = 0; nadj < ts->numcost; nadj++) PetscCall(VecScale(VHV[nadj], -1)); in TSComputeIHessianProductFunctionUU() 564 for (nadj = 0; nadj < ts->numcost; nadj++) PetscCall(VecScale(VHV[nadj], -1)); in TSComputeIHessianProductFunctionUP() 605 for (nadj = 0; nadj < ts->numcost; nadj++) PetscCall(VecScale(VHV[nadj], -1)); in TSComputeIHessianProductFunctionPU() 646 for (nadj = 0; nadj < ts->numcost; nadj++) PetscCall(VecScale(VHV[nadj], -1)); in TSComputeIHessianProductFunctionPP() 877 PetscErrorCode TSSetCostGradients(TS ts, PetscInt numcost, Vec lambda[], Vec mu[]) in TSSetCostGradients() argument 884 …PetscCheck(!ts->numcost || ts->numcost == numcost, PetscObjectComm((PetscObject)ts), PETSC_ERR_USE… in TSSetCostGradients() [all …]
|
| /petsc/src/ts/impls/implicit/theta/ |
| H A D | theta.c | 306 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStepBEuler_Private() 324 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStepBEuler_Private() 342 for (nadj = 0; nadj < ts->numcost; nadj++) { /* compute the residual */ in TSAdjointStepBEuler_Private() 349 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStepBEuler_Private() 365 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStepBEuler_Private() 389 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStepBEuler_Private() 461 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_Theta() 483 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_Theta() 505 for (nadj = 0; nadj < ts->numcost; nadj++) { /* compute the residual */ in TSAdjointStep_Theta() 512 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_Theta() [all …]
|
| /petsc/src/ts/impls/explicit/rk/ |
| H A D | rk.c | 856 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam)); in TSAdjointSetUp_RK() 857 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &rk->VecsSensiTemp)); in TSAdjointSetUp_RK() 860 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], s * ts->numcost, &rk->VecsDeltaLam2)); in TSAdjointSetUp_RK() 861 PetscCall(VecDuplicateVecs(ts->vecs_sensi2[0], ts->numcost, &rk->VecsSensi2Temp)); in TSAdjointSetUp_RK() 911 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_RK() 986 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_RK() 1023 for (nadj = 0; nadj < ts->numcost; nadj++) { /* no need to do this for mu's */ in TSAdjointStep_RK() 1037 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &rk->VecsDeltaLam)); in TSAdjointReset_RK() 1038 PetscCall(VecDestroyVecs(ts->numcost, &rk->VecsSensiTemp)); in TSAdjointReset_RK() 1040 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &rk->VecsDeltaLam2)); in TSAdjointReset_RK() [all …]
|
| /petsc/src/ts/tutorials/hybrid/ |
| H A D | ex1fd.c | 50 PetscInt numcost; in ShiftGradients() local 53 PetscCall(TSGetCostGradients(ts, &numcost, &lambda, &mu)); in ShiftGradients()
|
| H A D | ex1adj.c | 50 PetscInt numcost; in ShiftGradients() local 53 PetscCall(TSGetCostGradients(ts, &numcost, &lambda, &mu)); in ShiftGradients()
|
| /petsc/src/ts/impls/arkimex/ |
| H A D | arkimex.c | 1595 for (nadj = 0; nadj < ts->numcost; nadj++) { in TSAdjointStep_ARKIMEX() 1663 for (nadj = 0; nadj < ts->numcost; nadj++) // no need to do this for mu's in TSAdjointStep_ARKIMEX() 1783 PetscCall(VecDestroyVecs(tab->s * ts->numcost, &ark->VecsDeltaLam)); in TSAdjointReset_ARKIMEX() 1784 PetscCall(VecDestroyVecs(ts->numcost, &ark->VecsSensiTemp)); in TSAdjointReset_ARKIMEX() 1785 PetscCall(VecDestroyVecs(ts->numcost, &ark->VecsSensiPTemp)); in TSAdjointReset_ARKIMEX() 2054 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], tab->s * ts->numcost, &ark->VecsDeltaLam)); in TSAdjointSetUp_ARKIMEX() 2055 PetscCall(VecDuplicateVecs(ts->vecs_sensi[0], ts->numcost, &ark->VecsSensiTemp)); in TSAdjointSetUp_ARKIMEX() 2056 …if (ts->vecs_sensip) PetscCall(VecDuplicateVecs(ts->vecs_sensip[0], ts->numcost, &ark->VecsSensiPT… in TSAdjointSetUp_ARKIMEX()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex20td.c | 208 PetscErrorCode AdjointMonitor(TS ts, PetscInt steps, PetscReal time, Vec u, PetscInt numcost, Vec *… in AdjointMonitor() argument
|
| /petsc/include/petsc/private/ |
| H A D | tsimpl.h | 201 PetscInt numcost; /* number of cost functions */ member
|
| /petsc/doc/manual/ |
| H A D | ts.md | 1165 TSSetCostGradients(TS ts,PetscInt numcost, Vec *lambda,Vec *mu); 1168 where `numcost` denotes $n_\text{cost}$. 1275 with respect to the parameter $p_m$. $numcost$ is $1$
|