L_work and H_work

[clattenb]

Hello,

I'm hoping someone can give me a quick explanation of what the common arguments of L_work and H_work are for. I see these arguments in most of the functions (I am only needing to use gels and gesvd) and I just want to have a better understanding of what they're for beyond what is in the testing .cpp files.

Any help is appreciated.

Thanks,
-Nick

[mgates3]

Most algorithms need some temporary work space. In GELS, that is the array hwork, which has size (dimension) lwork. The exact need varies with the algorithm. For instance, it might be used to store a panel (block column) of A on the CPU, while most of A is stored on the GPU.
-mark

[clattenb]

Hi Mark,

Thanks for that explanation, I'm new to these linear algebra libraries (I have been using my own basic library up till now). For the two functions that I'm using (gels and gesvd) would you recommend that I just use the same formula that's in the code of the respective test application (such as testing_ dgeqrs3_gpu.cpp and testing_dgesvd.cpp)?

Thanks again,
-Nick