[petsc-dev] Possible bugs when using TS with ViennaCL

[Karl Rupp]

Hi Mani,

 > Thanks! It worked. However the following cases still don't work:
>
> 1)  Using ComputeResidual (not ComputeResidualViennaCL) but with
> DMSetVecType(da, VECVIENNACL) and DMSetMatType(da, MATAIJVIENNACL). The
> SNES solver converges but gives nonsensical results. I just get a
> flashing blob.
>
> 2) Using ComputeResidualViennaCL (not ComputeResidual) and *without*
> DMSetVecType(da, VECVIENNACL) and DMSetMatType(da, MATAIJVIENNACL). The
> SNES solver converges but gives nonsensical results. I just get a
> flashing blob as 1) above. This case probably should not work anyway
> cause I'm not sure if VecViennaCLGetArray works if the Vecs have not
> been set to ViennaCL.

Both cases worked for me when calling VecView() on the solution vector. 
Maybe I missed something, I'll check again tomorrow.

Best regards,
Karli