Search found 90 matches

by hartwig anzt
Wed Jan 31, 2018 8:50 am
Forum: User discussion
Topic: Solver parameters for eigenvalue routines
Replies: 7
Views: 2661

Re: Solver parameters for eigenvalue routines

We are looking into this, don't have a solution yet.
Please allow us a few days. Thanks, Hartwig
by hartwig anzt
Sun Nov 12, 2017 11:15 am
Forum: User discussion
Topic: How to use magma_*customspmv in the program?
Replies: 8
Views: 3031

Re: How to use magma_*customspmv in the program?

Mike, sorry for late reply! I think I do know what you mean now: you have two spmv kernels: spmv1 and spmv2, and you want to use them inside CG, time it, quickly switch between them. This would require recompiling magma-sparse. So, yes, unfortunately it is currently not possible to use multiple SpMV...
by hartwig anzt
Tue Nov 07, 2017 5:41 am
Forum: User discussion
Topic: How to use magma_*customspmv in the program?
Replies: 8
Views: 3031

Re: How to use magma_*customspmv in the program?

Mike,
I am not sure I understand correctly. You can use any SpMV inside the CG solver, the solver does not imply any restrictions on the matrix storage format or SpMV kernel.
Please specify.
Thanks!
by hartwig anzt
Tue Nov 07, 2017 3:22 am
Forum: User discussion
Topic: How to use magma_*customspmv in the program?
Replies: 8
Views: 3031

Re: How to use magma_*customspmv in the program?

Dear Mike,
magma_*customspmv.cpp currently does nothing. The idea is that there you call your function.
Assume you have you awesome spmv library. Then, you add in this file:

#include awesome_spmv.h
...
awesome_spmv( alpha, A x, beta, b, );

Does this make sense?
Thanks, Hartwig
by hartwig anzt
Thu Oct 19, 2017 10:44 am
Forum: User discussion
Topic: Matrix export to matrix market file contains only header
Replies: 2
Views: 765

Re: Matrix export to matrix market file contains only header

Klaus,
please always provide the complete code - this makes the analysis much easier.
Thanks, Hartwig
by hartwig anzt
Mon Oct 16, 2017 4:36 am
Forum: User discussion
Topic: How to retrieve final_res and numiter?
Replies: 4
Views: 1060

Re: How to retrieve final_res and numiter?

Klaus,
this is a diagonal matrix. You have the solution after one iteration.
Hartwig
by hartwig anzt
Sat Oct 14, 2017 11:40 am
Forum: User discussion
Topic: How to retrieve final_res and numiter?
Replies: 4
Views: 1060

Re: How to retrieve final_res and numiter?

Klaus,
please post the complete code. From that part, I can not deduce what is happening.
Thanks, Hartwig
by hartwig anzt
Thu Oct 12, 2017 4:20 pm
Forum: User discussion
Topic: How to copy solution vector x back to my application?
Replies: 2
Views: 1659

Re: How to copy solution vector x back to my application?

Klaus, the error code gives you the signature of the function: magma_int_t magma_dvget(magma_d_matrix, magma_int_t*, magma_int_t*, double**, magma_queue_t); so your first argument is of type magma_d_matrix, the second and third are integer pointers ( magma_int_t* ), then a double** and finally the q...
by hartwig anzt
Tue Oct 10, 2017 5:25 am
Forum: User discussion
Topic: Where can I find the linear solver documentation?
Replies: 4
Views: 1800

Re: Where can I find the linear solver documentation?

Dear Klaus, I do admit that this is currently not very well documented. What we have is: http://icl.cs.utk.edu/projectsfiles/magma/doxygen/group__magmasparse__dgesv.html and then click on the specific solvers. For the description of the solver_par object, please look into sparse/include/magmasparse_...
by hartwig anzt
Tue Oct 10, 2017 5:11 am
Forum: User discussion
Topic: Solver parameters for eigenvalue routines
Replies: 7
Views: 2661

Re: Solver parameters for eigenvalue routines

Olafur, 1. Is calling the routine magma_z_solver equivalent to calling the routine corresponding to the solver parameters? > Yes, the magma_z_solver is a wrapper for calling the specific routine. You are free to use the routine directly, if you prefer. 2. zlobpcg doesn't seem to have output paramete...