DensityMatrix::cal_DMR() with nspin=4 is significantly slower than nspin=1

[PeizeLin]

Details

job_nspin1.log
job_nspin4.log

In this crystal, the average time of DensityMatrix::cal_DMR() with nspin=1 is 0.02s, while with nspin=4 is 829.62s.

Task list for Issue attackers (only for developers)

• Reproduce the performance issue on a similar system or environment.
• Identify the specific section of the code causing the performance issue.
• Investigate the issue and determine the root cause.
• Research best practices and potential solutions for the identified performance issue.
• Implement the chosen solution to address the performance issue.
• Test the implemented solution to ensure it improves performance without introducing new issues.
• Optimize the solution if necessary, considering trade-offs between performance and other factors (e.g., code complexity, readability, maintainability).
• Review and incorporate any relevant feedback from users or developers.
• Merge the improved solution into the main codebase and notify the issue reporter.