The LAPACK forum has moved to https://github.com/Reference-LAPACK/lapack/discussions.

PDGESVD routine gives incorrect results:NB < M

Open discussion regarding features, bugs, issues, vendors, etc.

PDGESVD routine gives incorrect results:NB < M

Postby vikas » Tue Apr 29, 2008 3:02 am

Hi there,

I am using scalapack's PDGESVD routine. I dont seem to understand the affect of the row and column block size on the output results. Why i say this is due to the fact that i am observing that when i process a 3 x 3 matrix with a MB = NB = 3. I get proper S,U and VT matrices.

But when i run the same program with MB=NB = 2. The last entry in the S matrix is wrong.

It would be of great help some one could tell me whats happening here.

Thanks
-Vikas



The sample input and S matrices are given below:

A:
3 3
1
2
3
4
5
6
7
8
9

Output S matrix:

3 1
0.168481033526142063D+02
0.106836951455470852D+01
0.403135067055869220D-15
Code: Select all
Code: Select all
vikas
 
Posts: 4
Joined: Fri Apr 11, 2008 1:24 am

Return to User Discussion

Who is online

Users browsing this forum: No registered users and 6 guests