Test systems never converge

General CULA Sparse support and troubleshooting. Use this forum if you are having a general problem or have encountered a bug.

Test systems never converge

Postby psrgpu » 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.
Attachments
vetor_b.txt
(4.49 KiB) Downloaded 426 times
matrix_coo.txt
(23.53 KiB) Downloaded 463 times
psrgpu
 
Posts: 3
Joined: Mon Mar 12, 2012 3:32 pm

Re: Test systems never converge

Postby john » Mon Aug 12, 2013 5:40 am

Could you be more specific? At minimum we'd need to know your typical restart interval, and typical preconditioners, etc.
john
Administrator
 
Posts: 587
Joined: Thu Jul 23, 2009 2:31 pm

Re: Test systems never converge

Postby psrgpu » Mon Aug 12, 2013 7:17 am

Hi, John.

I have been trying a number of combinations. Solvers, preconditioners, etc.
For example:
Max Iterations: 500
tolerance: 1e-5
Restart: 10
No precond.

In this case, I've got
Residual: aprox. 2.29e-1

Just a little correction. I attached an unsorted matrix (matrix_coo.txt) but before sending it to the solver my program sorts it. I am attaching the correct file now.
Attachments
matrix_coo_sorted.txt
(23.53 KiB) Downloaded 466 times
psrgpu
 
Posts: 3
Joined: Mon Mar 12, 2012 3:32 pm

Re: Test systems never converge

Postby psrgpu » Tue Aug 27, 2013 4:47 am

No light at the end of the tunnel for me?
psrgpu
 
Posts: 3
Joined: Mon Mar 12, 2012 3:32 pm


Return to CULA Sparse Support

Who is online

Users browsing this forum: No registered users and 1 guest

cron