[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