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)

Most upvoted comments

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.