1 #include <../src/tao/bound/impls/bqnk/bqnk.h> 2 3 PETSC_EXTERN PetscErrorCode TaoCreate_BQNLS(Tao); 4 5 static PetscErrorCode TaoSetFromOptions_BLMVM(PetscOptionItems *PetscOptionsObject,Tao tao) 6 { 7 TAO_BNK *bnk = (TAO_BNK *)tao->data; 8 TAO_BQNK *bqnk = (TAO_BQNK*)bnk->ctx; 9 PetscErrorCode ierr; 10 KSPType ksp_type; 11 PetscBool is_spd; 12 13 PetscFunctionBegin; 14 ierr = PetscOptionsHead(PetscOptionsObject,"Quasi-Newton-Krylov method for bound constrained optimization");CHKERRQ(ierr); 15 ierr = PetscOptionsEList("-tao_blmvm_as_type", "active set estimation method", "", BNK_AS, BNK_AS_TYPES, BNK_AS[bnk->as_type], &bnk->as_type, 0);CHKERRQ(ierr); 16 ierr = PetscOptionsReal("-tao_blmvm_epsilon", "(developer) tolerance used when computing actual and predicted reduction", "", bnk->epsilon, &bnk->epsilon,NULL);CHKERRQ(ierr); 17 ierr = PetscOptionsReal("-tao_blmvm_as_tol", "(developer) initial tolerance used when estimating actively bounded variables", "", bnk->as_tol, &bnk->as_tol,NULL);CHKERRQ(ierr); 18 ierr = PetscOptionsReal("-tao_blmvm_as_step", "(developer) step length used when estimating actively bounded variables", "", bnk->as_step, &bnk->as_step,NULL);CHKERRQ(ierr); 19 ierr = PetscOptionsInt("-tao_blmvm_max_cg_its", "number of BNCG iterations to take for each Newton step", "", bnk->max_cg_its, &bnk->max_cg_its,NULL);CHKERRQ(ierr); 20 ierr = PetscOptionsTail();CHKERRQ(ierr); 21 ierr = TaoSetFromOptions(bnk->bncg);CHKERRQ(ierr); 22 ierr = TaoLineSearchSetFromOptions(tao->linesearch);CHKERRQ(ierr); 23 ierr = KSPSetFromOptions(tao->ksp);CHKERRQ(ierr); 24 ierr = KSPGetType(tao->ksp,&ksp_type);CHKERRQ(ierr); 25 bnk->is_nash = bnk->is_gltr = bnk->is_stcg = PETSC_FALSE; 26 ierr = MatSetFromOptions(bqnk->B);CHKERRQ(ierr); 27 ierr = MatGetOption(bqnk->B, MAT_SPD, &is_spd);CHKERRQ(ierr); 28 if (!is_spd) SETERRQ(PetscObjectComm((PetscObject)tao), PETSC_ERR_ARG_INCOMP, "LMVM matrix must be symmetric positive-definite"); 29 PetscFunctionReturn(0); 30 } 31 32 PETSC_EXTERN PetscErrorCode TaoCreate_BLMVM(Tao tao) 33 { 34 PetscErrorCode ierr; 35 36 PetscFunctionBegin; 37 ierr = PetscPrintf(PetscObjectComm((PetscObject)tao), "BLMVM is deprecated and will be removed in a future PETSc/TAO release. Please use BQNLS instead.\n");CHKERRQ(ierr); 38 ierr = TaoCreate_BQNLS(tao);CHKERRQ(ierr); 39 tao->ops->setfromoptions = TaoSetFromOptions_BLMVM; 40 PetscFunctionReturn(0); 41 }