Eig,Eigh for Block UniTensor

[j9263178]

#345.

Eig and Eigh now support Block UniTensor, for the behavior I just try a basic 2 site problem S^{+}S^{-} + S^{-}S^{+}:

a = cytnx.Bond(cytnx.BD_IN, [[0],[1]], [1,1]) #Bond for site a
b = cytnx.Bond(cytnx.BD_IN, [[0],[1]], [1,1]) #Bond for site b

H = cytnx.UniTensor([a,b,a.redirect(),b.redirect()],rowrank = 2) # 2 site hamiltonian
H.get_block_(2)[0,0] = 1
H.get_block_(3)[0,0] = 1

res = cytnx.linalg.Eigh(A, True, False)
print(res)

Output:

================================ Eigen vals ================================

-------- start of print ---------
Tensor name: 
braket_form : True
is_diag    : True
[OVERALL] contiguous : True
========================
 *is_diag: True
BLOCK [#0]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [0] U1(0)  -->| 1     1 |-->  [0] U1(0)
                 |         |
                 -----------

Total elem: 1
type  : Double (Float64)
cytnx device: CPU
Shape : (1)
[0.00000e+00 ]

========================
 *is_diag: True
BLOCK [#1]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [1] U1(1)  -->| 2     2 |-->  [1] U1(1)
                 |         |
                 -----------

Total elem: 2
type  : Double (Float64)
cytnx device: CPU
Shape : (2)
[-1.00000e+00 1.00000e+00 ]

========================
 *is_diag: True
BLOCK [#2]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [2] U1(2)  -->| 1     1 |-->  [2] U1(2)
                 |         |
                 -----------

Total elem: 1
type  : Double (Float64)
cytnx device: CPU
Shape : (1)
[0.00000e+00 ]



================================ Eigen vectors ================================

-------- start of print ---------
Tensor name: 
braket_form : True
is_diag    : False
[OVERALL] contiguous : True
========================
BLOCK [#0]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [0] U1(0)  -->| 1     1 |-->  [0] U1(0)
                 |         |
   [0] U1(0)  -->| 1       |
                 |         |
                 -----------

Total elem: 1
type  : Double (Float64)
cytnx device: CPU
Shape : (1,1,1)
[[[1.00000e+00 ]]]

========================
BLOCK [#1]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [0] U1(0)  -->| 1     2 |-->  [1] U1(1)
                 |         |
   [1] U1(1)  -->| 1       |
                 |         |
                 -----------

Total elem: 2
type  : Double (Float64)
cytnx device: CPU
Shape : (1,1,2)
[[[-7.07107e-01 7.07107e-01 ]]]

========================
BLOCK [#2]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [1] U1(1)  -->| 1     2 |-->  [1] U1(1)
                 |         |
   [0] U1(0)  -->| 1       |
                 |         |
                 -----------

Total elem: 2
type  : Double (Float64)
cytnx device: CPU
Shape : (1,1,2)
[[[7.07107e-01 7.07107e-01 ]]]

========================
BLOCK [#3]
 |- []   : Qn index 
 |- Sym(): Qnum of correspond symmetry
                 -----------
                 |         |
   [1] U1(1)  -->| 1     1 |-->  [2] U1(2)
                 |         |
   [1] U1(1)  -->| 1       |
                 |         |
                 -----------

Total elem: 1
type  : Double (Float64)
cytnx device: CPU
Shape : (1,1,1)
[[[1.00000e+00 ]]]

[pcchen]

Can someone fix the failed tests and then approve/merger this pull request. I have tested the basic function and it works ok.

[codecov]

Codecov Report

Attention: Patch coverage is 0.48544% with 205 lines in your changes are missing coverage. Please review.

Project coverage is 17.73%. Comparing base (c93b318) to head (194e558).

Files	Patch %	Lines
src/linalg/Eigh.cpp	0.96%	103 Missing ⚠️
src/linalg/Eig.cpp	0.00%	102 Missing ⚠️

Additional details and impacted files

@@              Coverage Diff               @@
##           dev-master     #372      +/-   ##
==============================================
- Coverage       17.81%   17.73%   -0.08%     
==============================================
  Files             220      220              
  Lines           48761    48965     +204     
  Branches        19955    20113     +158     
==============================================
  Hits             8686     8686              
- Misses          35871    36075     +204     
  Partials         4204     4204              

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