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
