Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

32-site basis fails

See original GitHub issue

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