Feature: Building numerical atomic orbitals from spherical Bessel coefficients

[jinzx10]

Reminder

• Have you linked an issue with this 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

#3267

What's changed?

This PR provides a new class "SphbesRadials" derived from the existing class "RadialSet". The new class is responsible for building numerical atomic orbitals from spherical Bessel coefficient file generated by various NAO-generation programs like SIAB or PTG_dpsi (with minor format modification). This is a crucial step towards solving #3267

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

No.

[kirk0830]

do you need to add the sphbes_radials.cpp to CMakeLists.txt in module_basis/module_nao/?

[jinzx10]

do you need to add the sphbes_radials.cpp to CMakeLists.txt in module_basis/module_nao/?

Given that a unified format of spherical Bessel coefficients is still not available, I think the integration of this class into the workflow might be left for a future PR (to keep the present PR minimal)