Test systems never converge
4 posts
• Page 1 of 1
Test systems never converge
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 1e6 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 1e6 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 1203 times

 matrix_coo.txt
 (23.53 KiB) Downloaded 1217 times
 psrgpu
 Posts: 3
 Joined: Mon Mar 12, 2012 3:32 pm
Re: Test systems never converge
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
Hi, John.
I have been trying a number of combinations. Solvers, preconditioners, etc.
For example:
Max Iterations: 500
tolerance: 1e5
Restart: 10
No precond.
In this case, I've got
Residual: aprox. 2.29e1
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.
I have been trying a number of combinations. Solvers, preconditioners, etc.
For example:
Max Iterations: 500
tolerance: 1e5
Restart: 10
No precond.
In this case, I've got
Residual: aprox. 2.29e1
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 1196 times
 psrgpu
 Posts: 3
 Joined: Mon Mar 12, 2012 3:32 pm
Re: Test systems never converge
No light at the end of the tunnel for me?
 psrgpu
 Posts: 3
 Joined: Mon Mar 12, 2012 3:32 pm
4 posts
• Page 1 of 1
Who is online
Users browsing this forum: No registered users and 1 guest