sddspseudoinverse error

[Hwang, Ji-Gwang]

Hello. There is an error in the sddspseudoinverse while I am trying to execute an example of "DATuneScan" with the recent version sdds-toolkit on the Windows machine. It works fine with the old-version sdds-toolkit.

[soliday]

I will take a look. I think I know what happened.

[soliday]

There is a new version posted for Windows. I had linked to the Intel MKL libraries that use 64bit array indexing which is incompatible with some of our other code. For now, I changed it to use the Intel MKL libraries that use 32bit array indexing.

[Hwang, Ji-Gwang]

Thanks! It works well.

[Slava]

Hello. There is an error in the sddspseudoinverse on the Windows 10 with the newest versions of Elegant and SDDS. When I tried to use sddspseudoinverse I've got messages:

Code: Select all

Compiled with package Intel MKL

Intel MKL ERROR: Parameter 8 was incorrect on entry to DGBBRDM.

Intel MKL ERROR: Parameter 5 was incorrect on entry to DGESDD.

Intel MKL ERROR: Parameter 8 was incorrect on entry to DGBBRDM.

Intel MKL ERROR: Parameter 5 was incorrect on entry to DGESDD.
Error (sddspseudoinverse): Error: no non-zero singular values found, unable to find the inverse response matrix.

I supposed that a problem may be in data type. I tried to use different data types (float, long, long64) but it has no effect. Moreover I can not obtained "double" data type. When I tried to use sddsprocess with -cast=column,C1,C,double the new column C1 had float data type.

Do you have any ideas about it?

[soliday]

I think I have found and fixed the problem. Can you download https://ops.aps.anl.gov/downloads/SDDS_ ... 4-TEST.msi
and tell me if this solves your problem? If it does, I will change the official release for Windows.

Thanks,
--Bob Soliday

[Slava]

Dear Bob, thanks a lot for your reply. This SDDS Toolkit release now works well. Results are coincide with that are from MtLab. Unfortunately I did not find information concerns parameter list. What do Tikhonovalpha and Tikhonovbeta parameters mean? Are they Tikhonov_relative_alpha and Tikhonov_n from &correct command? Does special entry for sddpseudoinverse exist that can change NumberOfSingularValuesUsed parameter?

[soliday]

I will direct your question to Mike Borland. He is lead author, I mostly just help handle releases and ports to the various operating systems.

[Slava]

Dear Bob,
when I continued my playing with sddspseudoinverse I faced with the next strange behaviour of the code.
If I use matrix prepared in other program and then make table data in sdds format via csv2sdds command and finally call sddspseudoinverse the result differs from that is obtained with MtLab. I attached input matrix and the result of the sddspseudoinverse as well as the result from MtLab (look at the picture).
Can you explain what is happened?