magma_dgeqrf vs. magma_dgeqrf4

Open discussion for MAGMA library (Matrix Algebra on GPU and Multicore Architectures)

magma_dgeqrf vs. magma_dgeqrf4

Postby christianHEL » Mon Sep 16, 2013 5:39 am

I am using MAGMA v1.4 with MKL, CUDA v5.5 and 2 Tesla M2070 GPUs. I am getting wrong results for matrices larger than 6000x6000 when using the magma_dgeqrf() function (with 1 GPU). I am comparing the solution to the corresponding dgeqrf() function in MKL. Interestingly, the results are correct when using the magma_dgeqrf4() with either 1 or 2 GPUs.
Posts: 3
Joined: Thu Aug 22, 2013 2:20 pm

Re: magma_dgeqrf vs. magma_dgeqrf4

Postby mgates3 » Mon Sep 23, 2013 12:32 pm

When you say you're comparing results to MKL, do you mean you're doing something like this?

Code: Select all
  A1 = your matrix
  A2 = A1
  normA = norm( A1 )
  magma_dgeqrf( ... A1 ... )
  mkl_dgeqrf( ... A2 ... )
  R = A1 - A2   // where A1 and A2 are the result of geqrf
  error = norm( R ) / normA

What error do you get? Can you reproduce this problem using the MAGMA tester, testing_dgeqrf? This actually implements two different tests. The -c flag will do || R - Q^T A || and ||I - Q^T Q||. The -c2 flag will do ||A1 - A2||, as above.
Code: Select all
./testing_sgeqrf -c
./testing_sgeqrf -c2

Note that the QR factorization is not unique -- it's possible to have two different, correct QR factorizations of the same matrix. So directly comparing the QR factorizations may indicate an error when none is present. However, our experience is that MAGMA and MKL produce the same QR factorization, since we use MKL for factoring panels.

Posts: 712
Joined: Fri Jan 06, 2012 2:13 pm

Return to User discussion

Who is online

Users browsing this forum: No registered users and 2 guests