Davidson subspace diagonalization method crashes at basically any CUDA GPU

[Flying-dragon-boxing]

Describe the bug

When using ks_solver=dav_subspace on any CUDA environment I have access to, ABACUS simply can't finish an iter.

CUDA 13.0 on RTX 3090, ABACUS 3.9.0.21

CUDA compatability layer in Hygon DTK 25.04.1 on Hygon's BW500 DCU, ABACUS 3.9.0.18

Expected behavior

Work normally.

To Reproduce

pw_Si2.zip

Environment

No response

Additional Context

No response

Task list for Issue attackers (only for developers)

• Verify the issue is not a duplicate.
• Describe the bug.
• Steps to reproduce.
• Expected behavior.
• Error message.
• Environment details.
• Additional context.
• Assign a priority level (low, medium, high, urgent).
• Assign the issue to a team member.
• Label the issue with relevant tags.
• Identify possible related issues.
• Create a unit test or automated test to reproduce the bug (if applicable).
• Fix the bug.
• Test the fix.
• Update documentation (if necessary).
• Close the issue and inform the reporter (if applicable).

[Cstandardlib]

This bug was introduced in v3.9.0.14, where the PR #6493 was merged to perf GPU dav_subspace.