def detrange2(L,a,energies):
res=[0.,0.,0.]
lenergies= len(list(energies))
for i in prange(lenergies):
res[i]= np.linalg.det(H_mat.Hmat(energies[i],L,a,0.,0.,1.0,0.5))
return res
def H2(e, L, a, a2):
return H_mat.Hmat(e, L, a, 0., 0., a2, 0.3)
def H2proj(e, L, a, a2):
return pr.irrep_proj(H_mat.Hmat(e, L, a, 0., 0., a2, 0.5), e, L, "E+")