search

nakatamaho / sdpa-gmp

SDPA gmp (arbitrary pricision) version
GNU General Public License v2.0
4 stars 6 forks source link

exact rational input? #1

Open dimpase opened 6 years ago

dimpase commented 6 years ago

Is it possible to accept exact rational input?

If not, is there any limit on the size of integers one can have in .dat files?

nakatamaho commented 6 years ago

Hi dimpase! thanks for using SDPA-GMP 1) Unfortunately, not. 2) I've never tried. I tried over 100 decimal digits, and it was okay

dimpase commented 6 years ago

Oh well. In SPDs I have (very badly conditioned) I end up with errors in Cholesky factorisation saying something like "matrix ix not positive semidefinite".

Could it be that some precision is checked using the doubles rather than arbitrary precision?

nakatamaho commented 6 years ago

Hi dimpase, Yes, you can control using 'precision' parameter in param.sdpa .

dimpase commented 6 years ago

No matter how big I make it, it seems not being enough. Perhaps I should really try huge values...

nakatamaho commented 6 years ago

Hi dimpase, could you please send me your logs to maho@riken.jp? thanks Nakata Maho

dimpase commented 6 years ago

I will try to reproduce this on the hardware I have with me (I am on vacation, and do not have access to the machine I used to compute these). It might take a bit of time.