QuSpin: 32-site basis fails
Dear developers of quspin,
I discovered a bug in my simulations.
Ns = 16
basis_full = spinful_fermion_basis_general(Ns, simple_symm = False)
cdag = quantum_operator(dict(H0=[['+', [[1, 0]]]], H1=[]), basis=basis_full, check_herm = False, check_symm=False)
However, the code falls with
OMP_NUM_THREADS=8 python3 test.py
/home/fkp/astronik/spt_newmodel.py:210: UserWarning: Test for particle conservation not implemented for <class 'quspin.basis.basis_general.fermion.spinful_fermion_basis_general'>, to turn off this warning set check_pcon=False in hamiltonian
cdag = quantum_operator(dict(H0=[['+', [[1, 0]]]], H1=[]), basis=basis_full, check_herm = False, check_symm=False)
Traceback (most recent call last):
File "/home/fkp/astronik/spt_newmodel.py", line 333, in <module>
main()
File "/home/fkp/astronik/spt_newmodel.py", line 210, in main
cdag = quantum_operator(dict(H0=[['+', [[1, 0]]]], H1=[]), basis=basis_full, check_herm = False, check_symm=False)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/operators/quantum_operator_core.py", line 172, in __init__
O = make_static(basis,opstr_list,dtype)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/operators/_make_hamiltonian.py", line 109, in make_static
return basis._make_matrix(static_list,dtype)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/basis/base.py", line 226, in _make_matrix
ME,row,col = self.Op(opstr,indx,J,dtype)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/basis/base.py", line 218, in Op
return self._Op(opstr,indx,J,dtype)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/basis/basis_general/fermion.py", line 618, in _Op
return spinless_fermion_basis_general._Op(self,opstr,indx,J,dtype)
File "/home/fkp/astronik/anaconda3/envs/quspin_env/lib/python3.9/site-packages/quspin/basis/basis_general/base_general.py", line 274, in _Op
self._core.op(row,col,ME,opstr,indx,J,self._basis,self._n,self._basis_begin,self._basis_end,self._N_p)
File "quspin/basis/basis_general/_basis_general_core/source/general_basis_core.pyx", line 224, in quspin.basis.basis_general._basis_general_core.spinful_fermion_core.general_basis_core_wrap.__pyx_fused_cpdef
TypeError: Function call with ambiguous argument types
Linux, 756 Gb RAM
Could you please suggest what is going possibly wrong? 32 sites is the biggest size that is yet supported, right?
Thanks a lot!
Best wishes, Nikita
About this issue
- Original URL
- State: closed
- Created 3 years ago
- Comments: 19 (7 by maintainers)
Dear @JiahaoYao , thank you very much for taking care of this. Indeed, for any Ns < 16, the code works. I am only facing problem with Ns = 16.