Use of culaDeviceSgeev

General discussion for CULA. Use this forum for questions, examples, feedback, and feature requests.

Use of culaDeviceSgeev

Postby amdj1 » Thu Mar 19, 2015 3:10 pm

Hi,

I am trying to use the culaDeviceSgeev as following:
culaDeviceSgeev('N', 'V', 3, d_covMatsManagedS, 3, d_eValsManagedS, d_eValsIManagedS, d_eVecsManaged[j], 1, d_eVecsManagedS, 3 );

All the pointers are device float pointers allocated by cudaMalloc. The function runs but when copy the data to host and print it out its not right. I get all zeros and ones in eigen vectors. Its totally different if I use the host function culaSgeev. The input matrix is same in both host and device cases. Any help is welcome.

Here is the sample ouput:
*****Output from host routine
Input Matrix:
Element 0: 507267
Element 1: 369762
Element 2: -902380
Element 3: 369762
Element 4: 269534
Element 5: -657748
Element 6: -902380
Element 7: -657748
Element 8: 1.6054e+06

Eigen Values 0: 2.38214e+06
Eigen Values 1: 0.121296
Eigen Values 2: 65.854

Eigen Vector Matrix 0: -0.461452
Eigen Vector Matrix 1: -0.336359
Eigen Vector Matrix 2: 0.820929
Eigen Vector Matrix 3: 0.704076
Eigen Vector Matrix 4: -0.701832
Eigen Vector Matrix 5: 0.108207
Eigen Vector Matrix 6: -0.541825
Eigen Vector Matrix 7: -0.625864
Eigen Vector Matrix 8: -0.561


*****Output from Device Routine:
Input Matrix:
Element 0: 507267
Element 1: 369762
Element 2: -902380
Element 3: 369762
Element 4: 269534
Element 5: -657748
Element 6: -902380
Element 7: -657748
Element 8: 1.6054e+06

Eigen Values 0: -r 359290
Eigen Values 1: -r 0
Eigen Values 2: -r 0

Eigen Vector Matrix 0: 1
Eigen Vector Matrix 1: 0
Eigen Vector Matrix 2: 0
Eigen Vector Matrix 3: 0
Eigen Vector Matrix 4: 1
Eigen Vector Matrix 5: 0
Eigen Vector Matrix 6: 0
Eigen Vector Matrix 7: 0
Eigen Vector Matrix 8: 1
Press any Key for next Frame.
amdj1
 
Posts: 2
Joined: Wed Mar 18, 2015 3:19 am

Return to General CULA Discussion

Who is online

Users browsing this forum: No registered users and 2 guests

cron