### Test systems never converge

Posted:

**Sun Aug 11, 2013 1:28 pm**I'm currently testing CULA library, but unfortunately I am having problem with the most part of my test systems. Not only the big ones, but with the small systems too.

I'm using the software R to compare the results between CULA GMRES(for example) and R GMRES implementation (package pracma).

The R version of GMRES has no problem to find the solution (known solution) for 1e-6 residual tolerance, but the CULA version never converges.

I'm attaching the matrix A (coo format) and vector b and I would appreciate if you could take a look at this system.

Thanks.

I'm using the software R to compare the results between CULA GMRES(for example) and R GMRES implementation (package pracma).

The R version of GMRES has no problem to find the solution (known solution) for 1e-6 residual tolerance, but the CULA version never converges.

I'm attaching the matrix A (coo format) and vector b and I would appreciate if you could take a look at this system.

Thanks.