Add functionality of MPO with block diagonal form

[SUSYUSTC]

Currently if one add two MPOs by mpo = mpo1 + mpo2 it creates zeros at off-diagonal blocks. In case the final MPO is composed of many small MPOs the computational efficiency is slow. I'll probably add another PR later for diagonal potential MPO, in which case the final MPO is composed of a dense kinetic block and a diagonal potential MPO block.

This pull request adds a StackedMpo class to improve the efficiency by calculating the Hamiltonians corresponding to each MPO part separately and then adding the Hamiltonians. It works for both the direct and iterative approach.

Example:

// original way: mpo = mpo_k + mpo_p
mpo = StackedMpo([mpo_k, mpo_p])
result = optimize_mps(mps.copy(), mpo)  // expect the same result

[codecov]

Codecov Report

Patch coverage: 98.14% and project coverage change: +0.05% 🎉

Comparison is base (9706a29) 84.76% compared to head (9a985d1) 84.81%.

Additional details and impacted files

@@            Coverage Diff             @@
##           master     #154      +/-   ##
==========================================
+ Coverage   84.76%   84.81%   +0.05%     
==========================================
  Files         105      105              
  Lines       10140    10183      +43     
==========================================
+ Hits         8595     8637      +42     
- Misses       1545     1546       +1     

Files Changed	Coverage Δ
renormalizer/mps/gs.py	95.93% <97.50%> (+0.12%)	⬆️
renormalizer/mps/__init__.py	100.00% <100.00%> (ø)
renormalizer/mps/mpo.py	91.15% <100.00%> (+0.09%)	⬆️
renormalizer/mps/tests/test_gs.py	100.00% <100.00%> (ø)

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

[jjren]

Excellent PR and many thanks to Jiace.
It is better to add a test in mps/tests/test_gs.py. For example, in test_optimization, you can simply add the two mpos together and check if the energy is doubled.

[SUSYUSTC]

Test already added.

[atxy-blip]

The new function looks good to me. Many thanks for your work!

Meanwhile, personally speaking I prefer well-documented code, so maybe you could add new comments and make changes to the original ones accordingly. For example, the following line should contain the new class:

Renormalizer/renormalizer/mps/gs.py

Line 55 in 9706a29

mpo : renormalizer.mps.Mpo

[SUSYUSTC]

Comments & types added