[petsc-dev] TaoTerminationReason -> TaoConvergedReason
Munson, Todd S.
tmunson at anl.gov
Tue Feb 18 10:00:30 CST 2014
I am not sure that VecGetSubVec and MatGetSubMat are needed. These should only be used
for the active-set methods (I hope) and there should be native PETSc equivalents used
by the PETSc active-set variational inequality solvers, for example.
Todd.
On Feb 18, 2014, at 9:22 AM, Jed Brown <jed at jedbrown.org> wrote:
> Lisandro Dalcin <dalcinl at gmail.com> writes:
>
>> Barry, I've pushed a quick fix (missing PETSC_EXTERN in TaoMonitor())
>> in your branch.
>>
>> The following prototypes in petsctao.h are also missing PETSC_EXTERN,
>> and I think they should be more properly namespaced.
>>
>> PetscErrorCode VecGetSubVec(Vec, IS, PetscInt, PetscReal, Vec*);
>> PetscErrorCode MatGetSubMat(Mat, IS, Vec, TaoSubsetType, Mat*);
>>
>> Perhaps you should move these prototypes (as well as TaoSubsetType) to
>> taoimpl.h ? I'm not sure these routines should be part of the public
>> interface.
>
> I would rather not, but if nothing else changes, the PetscInt argument
> should be TaoSubsetType. I'm afraid I don't understand exactly what
> role this function is serving, but the only TAO reference is the enum,
> so it can probably be made generic and moved down to the Vec module.
More information about the petsc-dev
mailing list