"Jose E. Roman" <[email protected]> writes: >> El 24 ago 2017, a las 22:51, Greg Meyer <[email protected]> escribió: >> >> 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 > > It is MATOP_GET_VECS. If you implement this operation, most SLEPc > functionality should work with CUDA. > > PETSc developers: shouldn't MATOP_GET_VECS be renamed to MATOP_CREATE_VECS ?
Yes, in their eagerness to rename interfaces, people often forget to complete all parts of the job. The MATOP interface has shown itself to be high maintenance.
signature.asc
Description: PGP signature
