FAµST 3.16.0 / Fix of the Hadamard 512×512 factorization in float/single precision

This version fixes the issue documented in the FAQ here without disabling the PALM4MSA normalization of factors as it was first proposed. The float precision leads in some rare cases (e.g. 512×512 Hadamard matrix factorization) to errors in the 2-norm of some L or R Fausts in the PALM4MSA algorithm. By recomputing these norms in double precision it allows to avoid NaN norms to make the algorithm fail.

Note: this workaround is only available on the 2020 implementation of PALM4MSA.

Comments are closed.