[petsc-dev] field split error
Mark F. Adams
mfadams at lbl.gov
Fri Sep 27 21:09:33 CDT 2013
>>
>> No idea. Do you want me to find out? What should a general IS return
>> as 'start'?
>
> ISStrideGetInfo does not dispatch through dynamically registered functionals.
>
> PetscErrorCode ISStrideGetInfo(IS is,PetscInt *first,PetscInt *step)
> {
> IS_Stride *sub;
>
> PetscFunctionBegin;
> PetscValidHeaderSpecific(is,IS_CLASSID,1);
> if (first) PetscValidIntPointer(first,2);
> if (step) PetscValidIntPointer(step,3);
>
> sub = (IS_Stride*)is->data;
> if (first) *first = sub->first;
> if (step) *step = sub->step;
> PetscFunctionReturn(0);
> }
>
As I said earlier sub->first is 1 for ISGeneral, apparently. This is a sort of sensible number. Is there an enum for 'crazy'? or perhaps I can change sub->first to be -1 in ISGeneral and then return an error here if (*first < 0) …
More information about the petsc-dev
mailing list