Thank you for the quick reply Karl, The symmetric matrices are almost certainly positive definite because the system can be solved using the CG solver without preconditioning as well as using the CG solver with the row_scaling preconditioner. It is only the jacobi and chow_patel_icc preconditioners that cause difficulties. I have verified that there are no zeros on the matrix diagonal, but have not performed a positive definite check other than what I stated above.
Thanks again for looking over the code snippet! All the best, Rick
_______________________________________________ ViennaCL-support mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/viennacl-support
