19 …PetscCall(MatCreateSubMatrix(tao->hessian, bnk->inactive_idx, bnk->inactive_idx, MAT_INITIAL_MATRI…  in TaoBNKComputeSubHessian()
24 …PetscCall(MatCreateSubMatrix(tao->hessian_pre, bnk->inactive_idx, bnk->inactive_idx, MAT_INITIAL_M…  in TaoBNKComputeSubHessian()
26     if (bnk->bfgs_pre) PetscCall(PCLMVMSetIS(bnk->bfgs_pre, bnk->inactive_idx));  in TaoBNKComputeSubHessian()
167               PetscCall(VecGetSubVector(bnk->W, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKInitialize()
168               PetscCall(VecGetSubVector(bnk->Xwork, bnk->inactive_idx, &bnk->inactive_work));  in TaoBNKInitialize()
176               PetscCall(VecRestoreSubVector(bnk->W, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKInitialize()
177               PetscCall(VecRestoreSubVector(bnk->Xwork, bnk->inactive_idx, &bnk->inactive_work));  in TaoBNKInitialize()
304     PetscCall(ISDestroy(&bnk->inactive_idx));  in TaoBNKEstimateActiveSet()
305 …ctive(tao->XL, tao->solution, bnk->unprojected_gradient, tao->XU, PETSC_TRUE, &bnk->inactive_idx));  in TaoBNKEstimateActiveSet()
307     PetscCall(ISComplementVec(bnk->inactive_idx, tao->solution, &bnk->active_idx));  in TaoBNKEstimateActiveSet()
332 …&bnk->active_lower, &bnk->active_upper, &bnk->active_fixed, &bnk->active_idx, &bnk->inactive_idx));  in TaoBNKEstimateActiveSet()
410   if (!bnk->inactive_idx) {  in TaoBNKComputeStep()
438     PetscCall(VecGetSubVector(bnk->Gwork, bnk->inactive_idx, &bnk->G_inactive));  in TaoBNKComputeStep()
439     PetscCall(VecGetSubVector(tao->stepdirection, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKComputeStep()
483     PetscCall(VecRestoreSubVector(bnk->Gwork, bnk->inactive_idx, &bnk->G_inactive));  in TaoBNKComputeStep()
484     PetscCall(VecRestoreSubVector(tao->stepdirection, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKComputeStep()
530     PetscCall(VecGetSubVector(tao->stepdirection, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKRecomputePred()
531     PetscCall(VecGetSubVector(bnk->Xwork, bnk->inactive_idx, &bnk->inactive_work));  in TaoBNKRecomputePred()
532     PetscCall(VecGetSubVector(bnk->Gwork, bnk->inactive_idx, &bnk->G_inactive));  in TaoBNKRecomputePred()
544     PetscCall(VecRestoreSubVector(tao->stepdirection, bnk->inactive_idx, &bnk->X_inactive));  in TaoBNKRecomputePred()
545     PetscCall(VecRestoreSubVector(bnk->Xwork, bnk->inactive_idx, &bnk->inactive_work));  in TaoBNKRecomputePred()
546     PetscCall(VecRestoreSubVector(bnk->Gwork, bnk->inactive_idx, &bnk->G_inactive));  in TaoBNKRecomputePred()
1029   bnk->inactive_idx  = NULL;  in TaoSetUp_BNK()
1058   PetscCall(ISDestroy(&bnk->inactive_idx));  in TaoDestroy_BNK()