Hi,<div><br></div><div>In FEniCS, we use ghosted vectors for parallel computations, with the functions</div><div><br></div><div>VecCreateGhost</div><div>VecGhostGetLocalForm</div><div><br></div><div>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. </div>

<div><br></div><div>Is it possible to set the type to be mpicusp when using ghost vectors, without changing much of the code? If so, how?</div><div><br></div><div>If not, how would you recommend I proceed to work with mpicusp vectors in this context?</div>

<div><br></div><div>Thanks!</div><div><br></div><div>-- </div><div>Fredrik</div>