1*e5fecd4eSAlp Dener #include <../src/tao/bound/impls/bqnk/bqnk.h> 2a7e14dcfSSatish Balay 3*e5fecd4eSAlp Dener PETSC_EXTERN PetscErrorCode TaoCreate_BQNLS(Tao); 4a7e14dcfSSatish Balay 52ec5c1acSAlp Dener static PetscErrorCode TaoSetFromOptions_BLMVM(PetscOptionItems *PetscOptionsObject,Tao tao) 6a7e14dcfSSatish Balay { 7*e5fecd4eSAlp Dener TAO_BNK *bnk = (TAO_BNK *)tao->data; 8*e5fecd4eSAlp Dener TAO_BQNK *bqnk = (TAO_BQNK*)bnk->ctx; 9a7e14dcfSSatish Balay PetscErrorCode ierr; 10*e5fecd4eSAlp Dener KSPType ksp_type; 11*e5fecd4eSAlp Dener PetscBool is_spd; 12a7e14dcfSSatish Balay 13a7e14dcfSSatish Balay PetscFunctionBegin; 14*e5fecd4eSAlp Dener ierr = PetscOptionsHead(PetscOptionsObject,"Quasi-Newton-Krylov method for bound constrained optimization");CHKERRQ(ierr); 15*e5fecd4eSAlp Dener 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*e5fecd4eSAlp Dener ierr = PetscOptionsReal("-tao_blmvm_epsilon", "(developer) tolerance used when computing actual and predicted reduction", "", bnk->epsilon, &bnk->epsilon,NULL);CHKERRQ(ierr); 17*e5fecd4eSAlp Dener 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*e5fecd4eSAlp Dener 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*e5fecd4eSAlp Dener 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*e5fecd4eSAlp Dener ierr = PetscOptionsTail();CHKERRQ(ierr); 21*e5fecd4eSAlp Dener ierr = TaoSetFromOptions(bnk->bncg);CHKERRQ(ierr); 22a7e14dcfSSatish Balay ierr = TaoLineSearchSetFromOptions(tao->linesearch);CHKERRQ(ierr); 23*e5fecd4eSAlp Dener ierr = KSPSetFromOptions(tao->ksp);CHKERRQ(ierr); 24*e5fecd4eSAlp Dener ierr = KSPGetType(tao->ksp,&ksp_type);CHKERRQ(ierr); 25*e5fecd4eSAlp Dener bnk->is_nash = bnk->is_gltr = bnk->is_stcg = PETSC_FALSE; 26*e5fecd4eSAlp Dener ierr = MatSetFromOptions(bqnk->B);CHKERRQ(ierr); 27*e5fecd4eSAlp Dener ierr = MatGetOption(bqnk->B, MAT_SPD, &is_spd);CHKERRQ(ierr); 28*e5fecd4eSAlp Dener if (!is_spd) SETERRQ(PetscObjectComm((PetscObject)tao), PETSC_ERR_ARG_INCOMP, "LMVM matrix must be symmetric positive-definite"); 29a7e14dcfSSatish Balay PetscFunctionReturn(0); 30a7e14dcfSSatish Balay } 31a7e14dcfSSatish Balay 32728e0ed0SBarry Smith PETSC_EXTERN PetscErrorCode TaoCreate_BLMVM(Tao tao) 33a7e14dcfSSatish Balay { 34a7e14dcfSSatish Balay PetscErrorCode ierr; 35a7e14dcfSSatish Balay 36a7e14dcfSSatish Balay PetscFunctionBegin; 37*e5fecd4eSAlp Dener 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*e5fecd4eSAlp Dener ierr = TaoCreate_BQNLS(tao);CHKERRQ(ierr); 39a7e14dcfSSatish Balay tao->ops->setfromoptions = TaoSetFromOptions_BLMVM; 40a9603a14SPatrick Farrell PetscFunctionReturn(0); 41a9603a14SPatrick Farrell }