2 posts
• Page **1** of **1**

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.

- christianHEL
**Posts:**3**Joined:**Thu Aug 22, 2013 2:20 pm

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

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.

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.

-mark

- 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.

-mark

- mgates3
**Posts:**718**Joined:**Fri Jan 06, 2012 2:13 pm

2 posts
• Page **1** of **1**

Users browsing this forum: No registered users and 1 guest