#199 xsqrtst_pt fails on ARM when tgeqr2 is called

Developer
closed-fixed
5
2012-03-24
2012-03-12
No

If you hack METHOD in atlas_qrmeth.h to always do either PCA copy or no copy, and hardwire nb=4 in ATL_geqrf, then N>10 will get the wrong answer. If tgeqr2 is not called, answer is always correct.

Discussion

  • R. Clint Whaley

    R. Clint Whaley - 2012-03-16

    Not sure if it is related, but xsqrtst_pt also fails sanity checks on 32-bit opt8.

     
  • R. Clint Whaley

    R. Clint Whaley - 2012-03-16

    G5 and G4 xsqrtst_pt sporadically fail as well, but this is almost certainly unrelated since it only occasionally happens.

     
  • R. Clint Whaley

    R. Clint Whaley - 2012-03-24
    • status: open --> closed-fixed
     
  • R. Clint Whaley

    R. Clint Whaley - 2012-03-24

    Both problems caused by out-of-order writes (ARM and PPC both do this). Problem fixed by turning off PCA for these archs.

     

Log in to post a comment.

Get latest updates about Open Source Projects, Conferences and News.

Sign up for the SourceForge newsletter:





No, thanks