[petsc-dev] checking if a vector is ghosted

[mspieg]

Nevermind…I should always think for 10 more minutes before I type…

Just want to check that the following is a kosher use of the return error handler

    bool _is_ghosted = true;
    PetscErrorCode ierr;
    PetscPushErrorHandler(PetscReturnErrorHandler,PETSC_NULL);
    ierr= VecGhostGetLocalForm(*x, x_ghosted.get());
    if (ierr != 0 )
      _is_ghosted = false;
    PetscPopErrorHandler();

cheers
marc

On Nov 8, 2011, at 2:06 PM, mspieg wrote:

> Hi All,
>   is there a clean way to test if a vector is ghosted (i.e created with VecCreateGhost…) without trying to call VecGhostGetLocalForm and throwing an exception?
> 
> In 3.1 I got around this using  PetscExceptonTry1  but apparently in 3.2 and dev  PetscException… is "too cool for school" ;^)
> 
> One possible  simple fix would be something like
> 
> PetscError VecGhostIsGhosted(Vec g, PetscBool isghosted)
> 
> and just return the results of  the logical tests in VecGhostGetLocalForm without throwing an error. 
> 
> (alternatively, I suppose I could mess with the ErrorHandlers but that seems like overkill for this).
> Let me know what might work.  
> cheers
> marc
> 
> 
> 
> 
> 
> ----------------------------------------------------
> Marc Spiegelman
> Lamont-Doherty Earth Observatory
> Dept. of Applied Physics/Applied Math
> Columbia University
> http://www.ldeo.columbia.edu/~mspieg
> tel: 845 704 2323 (SkypeIn)
> ----------------------------------------------------
> 

----------------------------------------------------
Marc Spiegelman
Lamont-Doherty Earth Observatory
Dept. of Applied Physics/Applied Math
Columbia University
http://www.ldeo.columbia.edu/~mspieg
tel: 845 704 2323 (SkypeIn)
----------------------------------------------------