[petsc-users] VECMPICUSP with ghosted vector

[Fredrik Heffer Valdmanis]

Hi,

In FEniCS, we use ghosted vectors for parallel computations, with the
functions

VecCreateGhost
VecGhostGetLocalForm

As I am integrating the new GPU-based vectors and matrices in FEniCS, I
want the ghosted vectors to be of type VECMPICUSP. I have tried to do this
by calling VecSetType after creating the vector, but that
makes VecGhostGetLocalForm give an error.

Is it possible to set the type to be mpicusp when using ghost vectors,
without changing much of the code? If so, how?

If not, how would you recommend I proceed to work with mpicusp vectors in
this context?

Thanks!

-- 
Fredrik
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120206/58ad94fb/attachment.htm>