Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Incorrect spline evaluation index assert error in diamond 1x1x1 (or related bug) #5174

Open
prckent opened this issue Sep 20, 2024 · 2 comments
Labels

Comments

@prckent
Copy link
Contributor

prckent commented Sep 20, 2024

Describe the bug

Amazingly the GCC14-NoMPI-Complex-Mixed-Debug run of short-diamondC_1x1x1_pp-dmc_sdj-1-16 ( https://cdash.qmcpack.org/tests/4312408 ) is giving an assertion failure indicating a bad index, and subsequently aborting. It is the only configuration for this test currently giving an error -- the other real/mixed/complex debug builds do not abort.

qmcpack_complex: /scratch/pk7/QMCPACK_CI_BUILDS/qmcpack/src/einspline/bspline_eval_d_std.cpp:39: void coeff_offset_UBspline_d(const Ugrid*, const BCtype_d*, double, int&, double&): Assertion `ipart <= grid->num -1' failed.

To Reproduce

Nightly run configuration, develop version, ctest -R short-diamondC_1x1x1_pp-dmc_sdj-1-16

Expected behavior

Test passes

System:
nitrogen2, nightly configuration

Additional context

??!!!

@prckent prckent added the bug label Sep 20, 2024
@prckent
Copy link
Contributor Author

prckent commented Sep 20, 2024

@ye-luo Can you please take a look at this? I wonder how this is logically possible. Perhaps we have a <= when we need a < somewhere, or variation in how the pbc is done allowing for numerical error? Happy to instrument this run and do more testing of it.

@ye-luo
Copy link
Contributor

ye-luo commented Sep 20, 2024

It is a MPC bug like #4725. Clearly there is unsafe code. Let me see if I can find the cause.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants