Skip to content

Prevent compilers from using FMA instructions#1033

Merged
langou merged 2 commits intoReference-LAPACK:masterfrom
ACSimon33:eigenvalue_fma_error
Jul 13, 2024
Merged

Prevent compilers from using FMA instructions#1033
langou merged 2 commits intoReference-LAPACK:masterfrom
ACSimon33:eigenvalue_fma_error

Conversation

@ACSimon33
Copy link
Contributor

Fixes #1031

Description

As described in #1031 when the compiler uses FMA instructions, errors in ?GEEVX can increase up to SQRT(EPS). The different rounding of the FMA instructions causes a change in control flow, which produces the errors. This MR uses explicit parentheses to prevent the compiler from using FMA instructions for those results, which are compared against zero in subsequent if-statements.

Checklist

  • The documentation has been updated.
  • If the PR solves a specific issue, it is set to be closed on merge.

…FMAs where the different rounding can lead to a change in control flow.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Less accuracy due to FMAs

2 participants