Search found 32 matches

by hsahasra
Tue Jun 02, 2015 1:14 pm
Forum: User discussion
Topic: What's the role of restart in PGMRES?
Replies: 1
Views: 857

What's the role of restart in PGMRES?

Hi, I have been experimenting solving a linear system with PGMRES+Jacobi preconditioner. The solver diverges with default restart. When I set the restart at 400, it converges in ~700 iterations. The number of iterations increases if I increase or decrease the restart from 400. How should I decide wh...
by hsahasra
Tue Jun 02, 2015 12:41 pm
Forum: User discussion
Topic: Crash in CGESV
Replies: 14
Views: 7987

Re: Crash in CGESV

Which GPU do you have? How big a problem you can solve depends on the GPU and it's memory.
by hsahasra
Tue Jun 02, 2015 12:24 pm
Forum: User discussion
Topic: What's wrong with the ILU sparse preconditioner?
Replies: 1
Views: 1209

What's wrong with the ILU sparse preconditioner?

Hi, I have been trying to use the ILU preconditioner in sparse-iter. I'm using it for iterative linear solvers. I tried the testing_zsolver example for LAPLACE2D of size 200 with PBiCGSTAB (Jacobi and ILU preconditioners). The linear solve finishes in 0.2 seconds with Jacobi precond, and takes 300 s...
by hsahasra
Tue Jun 02, 2015 12:20 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8033

Re: Unable to run zgetri_gpu

Please check this post. viewtopic.php?f=2&t=1062

This was the bug I found in zgetri. I was able to invert matrices with size divisible by 32 in pre 1.5.0 versions. Inversion of matrices of other sizes failed.
by hsahasra
Tue Jun 02, 2015 12:17 pm
Forum: User discussion
Topic: Unable to run zgetri_gpu
Replies: 14
Views: 8033

Re: Unable to run zgetri_gpu

Hi Luigi,

This was a bug with the inversion routine, so there was no way to make this run in the older versions. Apparently it has been fixed now, I haven't checked though. What problems are you facing? Which version of MAGMA are you using?

Harshad
by hsahasra
Tue Oct 07, 2014 12:13 pm
Forum: User discussion
Topic: zlaset_identity
Replies: 2
Views: 1752

Re: zlaset_identity

Thank you!
by hsahasra
Mon Oct 06, 2014 10:10 pm
Forum: User discussion
Topic: zlaset_identity
Replies: 2
Views: 1752

zlaset_identity

Hi,

I was using magmablas_zlaset_identity to create an identity matrix on GPU. I can't find it in MAGMA 1.5.0. Can you suggest an alternative that is present in the newest version?

Thanks,
Harshad
by hsahasra
Mon Jul 14, 2014 4:15 pm
Forum: User discussion
Topic: Possible bug in ztranspose-v2.cu
Replies: 3
Views: 2597

Re: Possible bug in ztranspose-v2.cu

Thank you. Is there a way I can access the code before it is released? Also, when are you planning to release MAGMA 1.5.0?

Thanks,
Harshad
by hsahasra
Mon Jun 30, 2014 2:25 pm
Forum: User discussion
Topic: Error compiling MAGMA 1.5.0 beta 2 in Debian Sid
Replies: 4
Views: 4552

Re: Error compiling MAGMA 1.5.0 beta 2 in Debian Sid

You need to define these preprocessor directives before you include magma.h, or while you compile your code: #ifndef HAVE_CUBLAS #define HAVE_CUBLAS #endif #ifndef ADD_ #define ADD_ #endif #include "magma.h" ---or--- g++-4.8 -DHAVE_CUBLAS -DADD_ This is when you are using CUDA, and BLAS/LAPACK mangl...
by hsahasra
Mon Jun 30, 2014 2:16 pm
Forum: User discussion
Topic: creating libmagma.a in 'lib'.
Replies: 5
Views: 2901

Re: creating libmagma.a in 'lib'.

Did you "make clean" first before compiling with -fPIC?
Also, please add -fPIC to the options, don't copy the options I showed above. I'm using MKL for BLAS/LAPACK and it looks like you are using OpenBLAS so you shouldn't have -DMAGMA_WITH_MKL