Refactor: Replace nlm_save in DeePKS by HContainer object psialpha.

[ErjieWu]

Linked Issue

Fix #5744 , set_size() is defined in BaseMatrix now and will be called when set_size() is called in AtomPair.

Unit Tests and/or Case Tests for my changes

• Change two reference results of integrat test in DeePKS for stress alignment.

What's changed?

• Use Hcontainer to build a new variable psialpha in LCAO_Deepks to replace nlm_save and nlm_save_k, overlap between projected orbitals and basis are now saved in same form for both gamma-only case and multi-k case.
• Add set_size() function in hamilt::BaseMatrix and call it in AtomPair to avoid wrong results.
• Add cal_stress value check for deepks_out_unittest=1.