[petsc-users] Shell matrix MatCreateVecs

Barry Smith bsmith at mcs.anl.gov
Fri Aug 25 16:56:19 CDT 2017


  Fixed in master, it is MATOP_CREATE_VECS

> On Aug 24, 2017, at 1:51 PM, Greg Meyer <gregory.meyer at gmail.com> wrote:
> 
> Hi,
> 
> I have written a shell matrix for non-standard vectors (CUDA to be specific) that works great. MatMult and MatNorm perform as they should. However, when I try to use it in a SLEPc algorithm, it breaks because MatCreateVecs yields standard vectors even if I set "-vec_type cuda". Looking through the source etc. it seems that MatShellSetOperation can't be used to implement a new MatCreateVecs function (MATOP_CREATEVECS does not exist). Is there a way to change the type of vector that comes from MatCreateVecs for shell matrices?
> 
> Thanks in advance,
> Greg



More information about the petsc-users mailing list