QMCPACK / QMCPACK/qmcpack

Cray libsci fails on large problems

Open
#4,929 1 comment 0 reactions 0 assignees View on GitHub
stale
Dominant language
C++
Stars
403
Forks
154
Avg merge
1d 12h
Merged PRs (30d)
82

Description

**Describe the bug**
On a new Cray Intel sapphire rapids machine, qmcpack was compiled with gnu compilers + Cray libsci. So far as I can tell, this is like MKL, but Cray. The code compiles and passes 100% of deterministic tests. With batched drivers on a problem containing ~400 electrons, the code runs as expected and gives good results that agree with previous runs on a different machine using legacy drivers. However, on a much larger problem containing > 2800 electrons, VMCBatched segfaults immediately. A backtrace suggests the issue is related to the implementation of dgetrf, which is a LAPACK routine to compute the LU decomposition. This is handled by libsci.

I then recompiled the code using gnu + mkl, and, again, all deterministic tests pass 100%. And it actually works for small and large problems, too.

This isn't really a bug _per se_. But this caused a lot of consternation and time spent debugging. So I mostly wanted to document for posterity.

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.