Feature: support k-parallelism diagonalization in lcao basis

[hongriTianqi]

this is a replace of previous PR #4561

Reminder

• Have you linked an issue with this pull request?
• Have you added adequate unit tests and/or case tests for your pull request?
• Have you noticed possible changes of behavior below or in the linked issue?
• Have you explained the changes of codes in core modules of ESolver, HSolver, ElecState, Hamilt, Operator or Psi? (ignore if not applicable)

Linked Issue

Fix #4550

Unit Tests and/or Case Tests for my changes

• Unit tests: source/module_hsolver/test/parallel_k2d_test.cpp
• Case test: tests/integrate/201_NO_15_pseudopots_kpar

What's changed?

• setting kpar > 1 in lcao calculation will trigger k-parallelism diagonallization.

Any changes of core modules? (ignore if not applicable)

• if (KPAR_LCAO > 1) hsolver_lcao will use k-parallelism code