Spin Half tests.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "ename": "ImportError",
     "evalue": "cannot import name 'save_data' from 'KitaevHoneycomb' (/home/jerry/KitaevHoneycombTenPy/KitaevHoneycomb.py)",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mImportError\u001b[0m                               Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-12-2b61ca0d4a83>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;32mfrom\u001b[0m \u001b[0mKitaevHoneycomb\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0mrun_atomic\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msave_data\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mImportError\u001b[0m: cannot import name 'save_data' from 'KitaevHoneycomb' (/home/jerry/KitaevHoneycombTenPy/KitaevHoneycomb.py)"
     ]
    }
   ],
   "source": [
    "from KitaevHoneycomb import run_atomic, save_data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  2 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.0 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 24.0 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 30.5 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 37.0 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=2 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  4 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.7 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.4 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 32.2 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 39.0 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=4 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  8 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['P_tol'] for config lanczos_params\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 13.0 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.7 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 32.5 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 39.3 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=8 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  16 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.1 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 31.8 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 38.5 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=16 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  32 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.1 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 31.8 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 38.4 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=32 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  64 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.1 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 31.7 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 38.4 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=64 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  128 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.5 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.1 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 31.7 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 38.6 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=128 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.0 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  256 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 12.9 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1250000000000000, S = nan, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 25.8 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 32.7 s\n",
      "Delta E = 0.0000e+00, Delta S = nan (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1250000000000000, S = 0.0000000000000000, norm_err = 0.0e+00\n",
      "Current memory usage 727.5 MB, time elapsed: 39.5 s\n",
      "Delta E = 0.0000e+00, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 0.0000e+00, max_E_trunc = 0.0000e+00\n",
      "MPS bond dimensions: [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 1\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.0\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=256 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  2 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1252343260471337, S = nan, norm_err = 1.5e-11\n",
      "Current memory usage 727.5 MB, time elapsed: 13.9 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 5.6224e-06, max_E_trunc = 3.9110e-07\n",
      "MPS bond dimensions: [2, 2, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 2]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1252343260351214, S = nan, norm_err = 1.3e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 27.8 s\n",
      "Delta E = 1.2012e-12, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.5599e-06, max_E_trunc = 3.9110e-07\n",
      "MPS bond dimensions: [2, 2, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1252343992457211, S = 0.0019678535461006, norm_err = 1.7e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 34.6 s\n",
      "Delta E = -7.3211e-09, Delta S = nan (per sweep)\n",
      "max_trunc_err = 3.9037e-07, max_E_trunc = 3.9111e-07\n",
      "MPS bond dimensions: [1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1252343829598672, S = 0.0019671137648039, norm_err = 1.4e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 41.5 s\n",
      "Delta E = 1.6286e-09, Delta S = -7.3978e-08 (per sweep)\n",
      "max_trunc_err = 3.9037e-07, max_E_trunc = 3.9111e-07\n",
      "MPS bond dimensions: [1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 2\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  1.000000000000002\n",
      "After the canonicalization:\n",
      "Bond dim =  [1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 1]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=2 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  4 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['P_tol'] for config lanczos_params\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1253516575330078, S = nan, norm_err = 4.1e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 18.5 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 2.4546e-03, max_E_trunc = 1.2278e-03\n",
      "MPS bond dimensions: [4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 3, 4, 3, 3, 4, 4, 4, 4, 3, 4, 4, 3, 4, 4, 4, 4, 4, 3, 4, 3, 4, 4, 4]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1253511279659486, S = nan, norm_err = 1.1e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 37.2 s\n",
      "Delta E = 5.2957e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.2381e-03, max_E_trunc = 6.2032e-04\n",
      "MPS bond dimensions: [4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 3, 4, 4, 4, 3, 4, 2, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1253252288302692, S = 0.0071637786592746, norm_err = 2.6e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 47.0 s\n",
      "Delta E = 2.5899e-06, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.9517e-06, max_E_trunc = 2.3411e-06\n",
      "MPS bond dimensions: [3, 4, 4, 4, 4, 2, 2, 1, 1, 1, 2, 2, 4, 2, 2, 1, 1, 1, 2, 2, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1253252128404891, S = 0.0214747650467492, norm_err = 2.9e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 56.6 s\n",
      "Delta E = 1.5990e-09, Delta S = 1.4311e-03 (per sweep)\n",
      "max_trunc_err = 1.9524e-06, max_E_trunc = 2.3408e-06\n",
      "MPS bond dimensions: [3, 3, 4, 4, 4, 2, 2, 1, 1, 1, 2, 2, 4, 2, 2, 1, 1, 1, 2, 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 8208\n",
      "Energy = -0.1253252128375900, S = 0.0214751908537870, norm_err = 2.7e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 66.2 s\n",
      "Delta E = 2.8990e-13, Delta S = 4.2581e-08 (per sweep)\n",
      "max_trunc_err = 1.9524e-06, max_E_trunc = 2.3408e-06\n",
      "MPS bond dimensions: [3, 3, 4, 4, 4, 2, 2, 1, 1, 1, 2, 2, 4, 4, 2, 1, 1, 1, 2, 2, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4]\n",
      "================================================================================\n",
      "DMRG finished after 50 sweeps.\n",
      "total size = 8208, maximum chi = 4\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/networks/mps.py:2609: UserWarning: canonical_form_infinite: project to smaller bond dimension\n",
      "  Gl, Wr = self._canonical_form_correct_left(j1, Gl, Wr_list[j1 % L])\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['P_tol'] for config lanczos_params\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Before the canonicalization:\n",
      "Bond dim =  [3, 2, 4, 3, 4, 2, 2, 1, 1, 1, 2, 2, 4, 4, 2, 1, 1, 1, 2, 2, 4, 3, 4, 2, 3, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 2, 3, 4, 4, 4, 4, 4, 3]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000073+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [3, 2, 4, 3, 4, 2, 2, 1, 1, 1, 2, 2, 4, 4, 2, 1, 1, 1, 2, 2, 4, 3, 4, 2, 3, 4, 4, 4, 4, 4, 4, 4, 3, 4, 4, 4, 4, 4, 4, 4, 2, 3, 4, 4, 4, 4, 4, 3]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=4 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  8 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256236039874251, S = nan, norm_err = 3.7e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 19.3 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.5776e-06, max_E_trunc = 1.5504e-06\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 7, 8, 8, 8, 8, 8, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256236157926634, S = nan, norm_err = 6.2e-09\n",
      "Current memory usage 727.5 MB, time elapsed: 38.7 s\n",
      "Delta E = -1.1805e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.7677e-06, max_E_trunc = 1.7637e-06\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256236160303388, S = 0.0529087572611060, norm_err = 3.8e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 49.4 s\n",
      "Delta E = -2.3768e-11, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.2373e-06, max_E_trunc = 1.7636e-06\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256236160305638, S = 0.0529087664336673, norm_err = 4.0e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 60.2 s\n",
      "Delta E = -2.2501e-14, Delta S = 9.1726e-10 (per sweep)\n",
      "max_trunc_err = 1.2373e-06, max_E_trunc = 1.7636e-06\n",
      "MPS bond dimensions: [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 8\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.0000000000000067+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=8 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  16 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['P_tol'] for config lanczos_params\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256252689406902, S = nan, norm_err = 2.1e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 27.1 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 4.2875e-08, max_E_trunc = 1.8109e-08\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256253193707006, S = nan, norm_err = 6.4e-09\n",
      "Current memory usage 727.5 MB, time elapsed: 45.0 s\n",
      "Delta E = -5.0430e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.9665e-08, max_E_trunc = 1.2438e-08\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256253196889728, S = 0.0531771490374143, norm_err = 3.9e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 56.2 s\n",
      "Delta E = -3.1827e-11, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.0205e-08, max_E_trunc = 1.4587e-08\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256253196919701, S = 0.0531773738127846, norm_err = 3.2e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 67.5 s\n",
      "Delta E = -2.9973e-13, Delta S = 2.2478e-08 (per sweep)\n",
      "max_trunc_err = 1.0440e-08, max_E_trunc = 1.5312e-08\n",
      "MPS bond dimensions: [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 16\n",
      "================================================================================\n",
      "Before the canonicalization:\n",
      "Bond dim =  [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.000000000000002+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=16 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  32 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256263122784667, S = nan, norm_err = 2.0e-07\n",
      "Current memory usage 727.5 MB, time elapsed: 63.4 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 4.2218e-08, max_E_trunc = 1.5208e-08\n",
      "MPS bond dimensions: [32, 32, 31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256263441226698, S = nan, norm_err = 3.3e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 89.4 s\n",
      "Delta E = -3.1844e-09, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.2183e-08, max_E_trunc = 8.8692e-09\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256263451391035, S = 0.6003613676266510, norm_err = 6.6e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 101.2 s\n",
      "Delta E = -1.0164e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 3.8122e-09, max_E_trunc = 7.9273e-09\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256263452115635, S = 0.6004206852890916, norm_err = 6.6e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 113.0 s\n",
      "Delta E = -7.2460e-12, Delta S = 5.9318e-06 (per sweep)\n",
      "max_trunc_err = 3.9717e-09, max_E_trunc = 7.9912e-09\n",
      "MPS bond dimensions: [32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32]\n",
      "================================================================================\n",
      "DMRG finished after 40 sweeps.\n",
      "total size = 6768, maximum chi = 32\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/networks/mps.py:2588: UserWarning: canonical_form_infinite: project to smaller bond dimension\n",
      "  Gl, Wr = self._canonical_form_correct_left(i1, Gl, Wr)\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/networks/mps.py:2609: UserWarning: canonical_form_infinite: project to smaller bond dimension\n",
      "  Gl, Wr = self._canonical_form_correct_left(j1, Gl, Wr_list[j1 % L])\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Before the canonicalization:\n",
      "Bond dim =  [31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 25, 32, 32, 32, 32, 32, 32, 32, 30, 32, 32, 32, 32, 32, 32, 32, 26, 32, 32, 32, 32, 32, 32, 32, 31, 32, 32, 32, 32, 32, 32, 32]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (0.999999999999999+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 25, 32, 32, 32, 32, 32, 32, 32, 30, 32, 32, 32, 32, 32, 32, 32, 26, 32, 32, 32, 32, 32, 32, 32, 31, 32, 32, 32, 32, 32, 32, 32]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=32 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n",
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  64 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256258211798495, S = nan, norm_err = 1.4e-07\n",
      "Current memory usage 727.5 MB, time elapsed: 76.3 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 8.8472e-08, max_E_trunc = 7.5911e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256264416684090, S = nan, norm_err = 3.0e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 105.4 s\n",
      "Delta E = -6.2049e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 7.1564e-08, max_E_trunc = 6.3181e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256264486272810, S = 0.9042382772105592, norm_err = 1.9e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 123.5 s\n",
      "Delta E = -6.9589e-10, Delta S = nan (per sweep)\n",
      "max_trunc_err = 4.6757e-08, max_E_trunc = 9.1909e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256264494666546, S = 0.9077167151113176, norm_err = 1.5e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 141.4 s\n",
      "Delta E = -8.3937e-11, Delta S = 3.4784e-04 (per sweep)\n",
      "max_trunc_err = 4.6766e-08, max_E_trunc = 9.1777e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 8208\n",
      "Energy = -0.1256264509015589, S = 0.9137153294338171, norm_err = 1.7e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 159.4 s\n",
      "Delta E = -1.4349e-10, Delta S = 5.9986e-04 (per sweep)\n",
      "max_trunc_err = 4.6698e-08, max_E_trunc = 9.1668e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 9648\n",
      "Energy = -0.1256264511174857, S = 0.9142485533576367, norm_err = 1.6e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 177.4 s\n",
      "Delta E = -2.1593e-11, Delta S = 5.3322e-05 (per sweep)\n",
      "max_trunc_err = 4.6658e-08, max_E_trunc = 9.1647e-08\n",
      "MPS bond dimensions: [64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64]\n",
      "================================================================================\n",
      "DMRG finished after 60 sweeps.\n",
      "total size = 9648, maximum chi = 64\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/networks/mps.py:2609: UserWarning: canonical_form_infinite: project to smaller bond dimension\n",
      "  Gl, Wr = self._canonical_form_correct_left(j1, Gl, Wr_list[j1 % L])\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Before the canonicalization:\n",
      "Bond dim =  [64, 64, 64, 64, 64, 64, 64, 64, 58, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 62, 64, 64, 64, 64, 64, 64, 64]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.000000000000003+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [64, 64, 64, 64, 64, 64, 64, 64, 58, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 64, 62, 64, 64, 64, 64, 64, 64, 64]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=64 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  128 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256263068480455, S = nan, norm_err = 6.8e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 108.6 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.0278e-07, max_E_trunc = 8.7779e-08\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256266156579008, S = nan, norm_err = 1.6e-08\n",
      "Current memory usage 727.5 MB, time elapsed: 150.9 s\n",
      "Delta E = -3.0881e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.8915e-07, max_E_trunc = 1.7986e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256266635875841, S = 1.4029175225580308, norm_err = 2.0e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 177.7 s\n",
      "Delta E = -4.7930e-09, Delta S = nan (per sweep)\n",
      "max_trunc_err = 9.5684e-08, max_E_trunc = 1.8738e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256266640987107, S = 1.4039508898996687, norm_err = 2.1e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 204.4 s\n",
      "Delta E = -5.1113e-11, Delta S = 1.0334e-04 (per sweep)\n",
      "max_trunc_err = 9.6443e-08, max_E_trunc = 1.8883e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 8208\n",
      "Energy = -0.1256266645758674, S = 1.4049943335864212, norm_err = 1.7e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 231.0 s\n",
      "Delta E = -4.7716e-11, Delta S = 1.0434e-04 (per sweep)\n",
      "max_trunc_err = 9.6619e-08, max_E_trunc = 1.8912e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 9648\n",
      "Energy = -0.1256266650697266, S = 1.4060391767510678, norm_err = 1.8e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 257.7 s\n",
      "Delta E = -4.9386e-11, Delta S = 1.0448e-04 (per sweep)\n",
      "max_trunc_err = 9.6945e-08, max_E_trunc = 1.8972e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 70, age = 11088\n",
      "Energy = -0.1256266655003770, S = 1.4070851699981410, norm_err = 1.9e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 284.1 s\n",
      "Delta E = -4.3065e-11, Delta S = 1.0460e-04 (per sweep)\n",
      "max_trunc_err = 9.7265e-08, max_E_trunc = 1.9032e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 80, age = 12528\n",
      "Energy = -0.1256266658988826, S = 1.4081290969101892, norm_err = 2.1e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 310.5 s\n",
      "Delta E = -3.9851e-11, Delta S = 1.0439e-04 (per sweep)\n",
      "max_trunc_err = 9.7586e-08, max_E_trunc = 1.9091e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 90, age = 13968\n",
      "Energy = -0.1256266662908037, S = 1.4091771423642159, norm_err = 1.9e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 336.8 s\n",
      "Delta E = -3.9192e-11, Delta S = 1.0480e-04 (per sweep)\n",
      "max_trunc_err = 9.7903e-08, max_E_trunc = 1.9151e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 100, age = 15408\n",
      "Energy = -0.1256266666858631, S = 1.4102250064051678, norm_err = 2.0e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 363.2 s\n",
      "Delta E = -3.9506e-11, Delta S = 1.0479e-04 (per sweep)\n",
      "max_trunc_err = 9.8204e-08, max_E_trunc = 1.9208e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 110, age = 16848\n",
      "Energy = -0.1256266670654848, S = 1.4112707261235162, norm_err = 2.0e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 389.7 s\n",
      "Delta E = -3.7962e-11, Delta S = 1.0457e-04 (per sweep)\n",
      "max_trunc_err = 9.8424e-08, max_E_trunc = 1.9251e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 120, age = 18288\n",
      "Energy = -0.1256266674468994, S = 1.4123141016747913, norm_err = 1.8e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 416.1 s\n",
      "Delta E = -3.8141e-11, Delta S = 1.0434e-04 (per sweep)\n",
      "max_trunc_err = 9.8680e-08, max_E_trunc = 1.9299e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 130, age = 19728\n",
      "Energy = -0.1256266678237381, S = 1.4133574159348405, norm_err = 1.6e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 442.5 s\n",
      "Delta E = -3.7684e-11, Delta S = 1.0433e-04 (per sweep)\n",
      "max_trunc_err = 9.9004e-08, max_E_trunc = 1.9361e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 140, age = 21168\n",
      "Energy = -0.1256266682412578, S = 1.4143981982643561, norm_err = 1.4e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 469.1 s\n",
      "Delta E = -4.1752e-11, Delta S = 1.0408e-04 (per sweep)\n",
      "max_trunc_err = 9.9112e-08, max_E_trunc = 1.9380e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 150, age = 22608\n",
      "Energy = -0.1256266686228523, S = 1.4154406622710205, norm_err = 1.7e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 495.9 s\n",
      "Delta E = -3.8159e-11, Delta S = 1.0425e-04 (per sweep)\n",
      "max_trunc_err = 9.9415e-08, max_E_trunc = 1.9437e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 160, age = 24048\n",
      "Energy = -0.1256266689980947, S = 1.4164838515744067, norm_err = 1.3e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 522.5 s\n",
      "Delta E = -3.7524e-11, Delta S = 1.0432e-04 (per sweep)\n",
      "max_trunc_err = 9.9777e-08, max_E_trunc = 1.9505e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 170, age = 25488\n",
      "Energy = -0.1256266693922091, S = 1.4175278110688767, norm_err = 1.4e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 549.2 s\n",
      "Delta E = -3.9411e-11, Delta S = 1.0440e-04 (per sweep)\n",
      "max_trunc_err = 1.0039e-07, max_E_trunc = 1.9625e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 180, age = 26928\n",
      "Energy = -0.1256266697649693, S = 1.4185693352249784, norm_err = 1.4e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 575.8 s\n",
      "Delta E = -3.7276e-11, Delta S = 1.0415e-04 (per sweep)\n",
      "max_trunc_err = 1.0071e-07, max_E_trunc = 1.9685e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 190, age = 28368\n",
      "Energy = -0.1256266701410880, S = 1.4196134004202141, norm_err = 1.5e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 602.4 s\n",
      "Delta E = -3.7612e-11, Delta S = 1.0441e-04 (per sweep)\n",
      "max_trunc_err = 1.0103e-07, max_E_trunc = 1.9746e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 200, age = 29808\n",
      "Energy = -0.1256266705139903, S = 1.4206577784875130, norm_err = 1.5e-14\n",
      "Current memory usage 727.5 MB, time elapsed: 630.0 s\n",
      "Delta E = -3.7290e-11, Delta S = 1.0444e-04 (per sweep)\n",
      "max_trunc_err = 1.0137e-07, max_E_trunc = 1.9810e-07\n",
      "MPS bond dimensions: [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "================================================================================\n",
      "DMRG finished after 200 sweeps.\n",
      "total size = 29808, maximum chi = 128\n",
      "================================================================================\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/networks/mps.py:2609: UserWarning: canonical_form_infinite: project to smaller bond dimension\n",
      "  Gl, Wr = self._canonical_form_correct_left(j1, Gl, Wr_list[j1 % L])\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Before the canonicalization:\n",
      "Bond dim =  [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 127, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Canonicalizing...\n",
      "The norm is:  1.0\n",
      "The overlap is:  (1.000000000000021+0j)\n",
      "After the canonicalization:\n",
      "Bond dim =  [128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 127, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128, 128]\n",
      "Computing properties\n",
      "Optimizing\n",
      "Loop for chi=128 done.\n",
      "================================================================================\n",
      "============================== END ==============================\n",
      "================================================================================\n",
      "Reading 'Lx'=6 for config KitaevHoneycombModel\n",
      "Reading 'Ly'=4 for config KitaevHoneycombModel\n",
      "Reading 'conserve'='parity' for config KitaevHoneycombModel\n",
      "Reading 'S'=0.5 for config KitaevHoneycombModel\n",
      "Reading 'order'='Cstyle' for config KitaevHoneycombModel\n",
      "Reading 'bc'='periodic' for config KitaevHoneycombModel\n",
      "Reading 'bc_MPS'='infinite' for config KitaevHoneycombModel\n",
      "Reading 'Jx'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jy'=0.1 for config KitaevHoneycombModel\n",
      "Reading 'Jz'=1.0 for config KitaevHoneycombModel\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n",
      "/home/jerry/anaconda3/lib/python3.7/site-packages/tenpy/tools/params.py:225: UserWarning: unused option ['dtype'] for config KitaevHoneycombModel\n",
      "\n",
      "  warnings.warn(msg.format(keys=sorted(unused), name=self.name))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      "================================================================================\n",
      "==============================START==============================\n",
      "================================================================================\n",
      "Chi =  256 \n",
      "\n",
      "Subconfig 'trunc_params'=Config(<6 options>, 'trunc_params') for config TwoSiteDMRGEngine\n",
      "Initial sweeps...\n",
      "Reading 'start_env'=10 for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      "..........\n",
      "Reading 'N_sweeps_check'=10 for config TwoSiteDMRGEngine\n",
      "Reading 'max_sweeps'=200 for config TwoSiteDMRGEngine\n",
      "Reading 'max_E_err'=1e-06 for config TwoSiteDMRGEngine\n",
      "Reading 'max_S_err'=0.0001 for config TwoSiteDMRGEngine\n",
      "Reading 'mixer'=True for config TwoSiteDMRGEngine\n",
      "Subconfig 'mixer_params'=Config(<4 options>, 'mixer_params') for config TwoSiteDMRGEngine\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 10, age = 2448\n",
      "Energy = -0.1256268556299579, S = nan, norm_err = 2.7e-09\n",
      "Current memory usage 727.5 MB, time elapsed: 269.4 s\n",
      "Delta E = nan, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.4763e-09, max_E_trunc = 2.0403e-09\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 20, age = 3888\n",
      "Energy = -0.1256269612940445, S = nan, norm_err = 5.7e-09\n",
      "Current memory usage 727.5 MB, time elapsed: 544.6 s\n",
      "Delta E = -1.0566e-08, Delta S = 0.0000e+00 (per sweep)\n",
      "max_trunc_err = 1.3071e-09, max_E_trunc = 1.1046e-09\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Convergence criterium reached with enabled mixer.\n",
      "disable mixer and continue\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 30, age = 5328\n",
      "Energy = -0.1256270901647246, S = 1.9590296704245505, norm_err = 9.1e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 787.7 s\n",
      "Delta E = -1.2887e-08, Delta S = nan (per sweep)\n",
      "max_trunc_err = 1.9938e-10, max_E_trunc = 1.8395e-10\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 40, age = 6768\n",
      "Energy = -0.1256271202476450, S = 1.9837649237689945, norm_err = 7.0e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 1031.3 s\n",
      "Delta E = -3.0083e-09, Delta S = 2.4735e-03 (per sweep)\n",
      "max_trunc_err = 5.5025e-11, max_E_trunc = 5.5593e-11\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 50, age = 8208\n",
      "Energy = -0.1256271303941799, S = 2.0020232802289195, norm_err = 7.1e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 1278.2 s\n",
      "Delta E = -1.0147e-09, Delta S = 1.8258e-03 (per sweep)\n",
      "max_trunc_err = 2.0364e-11, max_E_trunc = 2.0123e-11\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 60, age = 9648\n",
      "Energy = -0.1256271348880821, S = 2.0106022827079504, norm_err = 7.1e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 1522.6 s\n",
      "Delta E = -4.4939e-10, Delta S = 8.5790e-04 (per sweep)\n",
      "max_trunc_err = 7.3612e-12, max_E_trunc = 9.5497e-12\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n",
      "Updating environment\n",
      ".....\n",
      "================================================================================\n",
      "sweep 70, age = 11088\n",
      "Energy = -0.1256271365717770, S = 2.0134776414528361, norm_err = 7.4e-15\n",
      "Current memory usage 727.5 MB, time elapsed: 1766.5 s\n",
      "Delta E = -1.6837e-10, Delta S = 2.8754e-04 (per sweep)\n",
      "max_trunc_err = 2.8898e-12, max_E_trunc = 7.7307e-12\n",
      "MPS bond dimensions: [256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256, 256]\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-7-36819339b9e7>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m     27\u001b[0m                 \u001b[0mshift\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mshift\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     28\u001b[0m                 \u001b[0mbc\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mbc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbc_MPS\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mbc_MPS\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 29\u001b[0;31m                 \u001b[0minitial_psi\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0minitial_psi\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     30\u001b[0m             )\n\u001b[1;32m     31\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/KitaevHoneycombTenPy/KitaevHoneycomb.py\u001b[0m in \u001b[0;36mrun_atomic\u001b[0;34m(chi, Jx, Jy, Jz, S, Lx, Ly, bc, shift, bc_MPS, conserve, order, dtype, initial_psi, initial, max_E_err, max_S_err, max_sweeps, N_sweeps_check, canonicalized, mixer, mixer_params, trunc_params, verbose, result_verbose)\u001b[0m\n\u001b[1;32m    171\u001b[0m     \u001b[0meng\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mreset_stats\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    172\u001b[0m     \u001b[0meng\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtrunc_params\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'chi_max'\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mchi\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 173\u001b[0;31m     \u001b[0minfo\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0meng\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    174\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    175\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mcanonicalized\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/algorithms/dmrg.py\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    337\u001b[0m             \u001b[0;31m# --------- the main work --------------\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    338\u001b[0m             \u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mN_sweeps_check\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 339\u001b[0;31m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msweep\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmeas_E_trunc\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    340\u001b[0m             \u001b[0mmax_trunc_err\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msweep\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmeas_E_trunc\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    341\u001b[0m             \u001b[0mmax_E_trunc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmax\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mE_trunc_list\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/algorithms/dmrg.py\u001b[0m in \u001b[0;36msweep\u001b[0;34m(self, optimize, meas_E_trunc)\u001b[0m\n\u001b[1;32m    743\u001b[0m         \u001b[0;31m# wrapper around tenpy.algorithms.mps_common.Sweep.sweep()\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    744\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_meas_E_trunc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmeas_E_trunc\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 745\u001b[0;31m         \u001b[0mres\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msuper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msweep\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moptimize\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    746\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0moptimize\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    747\u001b[0m             \u001b[0;31m# update mixer\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/algorithms/mps_common.py\u001b[0m in \u001b[0;36msweep\u001b[0;34m(self, optimize)\u001b[0m\n\u001b[1;32m    298\u001b[0m             \u001b[0;31m# --------- the main work --------------\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    299\u001b[0m             \u001b[0mtheta\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mprepare_update\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 300\u001b[0;31m             \u001b[0mupdate_data\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_local\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtheta\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moptimize\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0moptimize\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    301\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0mupdate_LP\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    302\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mupdate_LP\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mupdate_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'U'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m  \u001b[0;31m# (requires updated B)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/algorithms/dmrg.py\u001b[0m in \u001b[0;36mupdate_local\u001b[0;34m(self, theta, optimize, meas_E_trunc)\u001b[0m\n\u001b[1;32m    884\u001b[0m         \u001b[0mage\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0menv\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_LP_age\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mi0\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0;36m2\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0menv\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_RP_age\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mi0\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    885\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0moptimize\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 886\u001b[0;31m             \u001b[0mE0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mov_change\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdiag\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtheta\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    887\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    888\u001b[0m             \u001b[0mE0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mov_change\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m0.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/algorithms/dmrg.py\u001b[0m in \u001b[0;36mdiag\u001b[0;34m(self, theta_guess)\u001b[0m\n\u001b[1;32m    604\u001b[0m                 \u001b[0mE\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfull_diag_effH\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0meff_H\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta_guess\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeep_sector\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    605\u001b[0m             \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 606\u001b[0;31m                 \u001b[0mE\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlanczos\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0meff_H\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta_guess\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlanczos_params\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    607\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdiag_method\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'lanczos'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    608\u001b[0m             \u001b[0mE\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlanczos\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0meff_H\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtheta_guess\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlanczos_params\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/lanczos.py\u001b[0m in \u001b[0;36mlanczos\u001b[0;34m(H, psi, options, orthogonal_to)\u001b[0m\n\u001b[1;32m    397\u001b[0m         \u001b[0mSee\u001b[0m \u001b[0;34m:\u001b[0m\u001b[0mmeth\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;31m`\u001b[0m\u001b[0mLanczosGroundState\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;31m`\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    398\u001b[0m     \"\"\"\n\u001b[0;32m--> 399\u001b[0;31m     \u001b[0;32mreturn\u001b[0m \u001b[0mLanczosGroundState\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mH\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mpsi\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moptions\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0morthogonal_to\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    400\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    401\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/lanczos.py\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    186\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mN\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    187\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0mE0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpsi0\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcopy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m  \u001b[0;31m# no better estimate available\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 188\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mE0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_calc_result_full\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mN\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    189\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    190\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m_calc_T\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/lanczos.py\u001b[0m in \u001b[0;36m_calc_result_full\u001b[0;34m(self, N)\u001b[0m\n\u001b[1;32m    227\u001b[0m         \u001b[0;31m# and the last len_cache vectors have been cached\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    228\u001b[0m         \u001b[0;32mfor\u001b[0m \u001b[0mk\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmin\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlen_cache\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mN\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 229\u001b[0;31m             \u001b[0mpsif\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0miadd_prefactor_other\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mvf\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mN\u001b[0m \u001b[0;34m-\u001b[0m \u001b[0mk\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_cache\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0mk\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    230\u001b[0m         \u001b[0;31m# other vectors are not cached, so we need to restart the Lanczos iteration.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    231\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_cache\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m  \u001b[0;31m# free memory: we need at least two more vectors\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/_npc_helper.pyx\u001b[0m in \u001b[0;36mtenpy.linalg._npc_helper.Array_iadd_prefactor_other\u001b[0;34m()\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/charges.py\u001b[0m in \u001b[0;36mtest_equal\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m    830\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0moptimize\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mOptimizationFlag\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mskip_arg_checks\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    831\u001b[0m             \u001b[0;32mreturn\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 832\u001b[0;31m         \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchinfo\u001b[0m \u001b[0;34m!=\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchinfo\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    833\u001b[0m             raise ValueError(''.join(\n\u001b[1;32m    834\u001b[0m                 [\"incompatible ChargeInfo\\n\",\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/charges.py\u001b[0m in \u001b[0;36m__ne__\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m    298\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__ne__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    299\u001b[0m         \u001b[0;34mr\"\"\"Define `self != other` as `not (self == other)`\"\"\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 300\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__eq__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mother\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    301\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    302\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/tenpy/linalg/charges.py\u001b[0m in \u001b[0;36m__eq__\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m    289\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    290\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 291\u001b[0;31m         \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmod\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmod\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    292\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0;32mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    293\u001b[0m         \u001b[0;32mfor\u001b[0m \u001b[0ml\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mr\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mzip\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnames\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnames\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m<__array_function__ internals>\u001b[0m in \u001b[0;36mall\u001b[0;34m(*args, **kwargs)\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/site-packages/numpy/core/fromnumeric.py\u001b[0m in \u001b[0;36m_all_dispatcher\u001b[0;34m(a, axis, out, keepdims)\u001b[0m\n\u001b[1;32m   2331\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2332\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2333\u001b[0;31m \u001b[0;32mdef\u001b[0m \u001b[0m_all_dispatcher\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0maxis\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mout\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeepdims\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   2334\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0ma\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mout\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2335\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "Jx = 1.0\n",
    "Jy = 1.0\n",
    "Jz = 1.0\n",
    "\n",
    "# S = 1.5\n",
    "\n",
    "Lx = 6\n",
    "Ly = 4\n",
    "shift = 0\n",
    "bc = 'periodic'\n",
    "bc_MPS = 'infinite'\n",
    "\n",
    "chi_list = [2**i for i in range(1, 9)]\n",
    "\n",
    "result_dict = {}\n",
    "\n",
    "for J in np.linspace(0, 1.0, 11):\n",
    "    J = J.round(decimals=2)\n",
    "    initial_psi = None\n",
    "    for chi in chi_list:\n",
    "        param = (J, J, Jz, chi)\n",
    "        if param not in result_dict:\n",
    "            result = run_atomic(\n",
    "                Jx=J, \n",
    "                Jy=J, \n",
    "                Jz=Jz, \n",
    "                Lx=Lx, Ly=Ly, \n",
    "                chi=chi,\n",
    "                shift=shift,\n",
    "                bc=bc, bc_MPS=bc_MPS,\n",
    "                initial_psi=initial_psi,\n",
    "            )\n",
    "        else:\n",
    "            result = result_dict[param]\n",
    "        initial_psi = result['psi'].copy()\n",
    "        result_dict[param] = result\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(0.0, 0.0, 1.0, 2)\n",
      "(0.0, 0.0, 1.0, 4)\n",
      "(0.0, 0.0, 1.0, 8)\n",
      "(0.0, 0.0, 1.0, 16)\n",
      "(0.0, 0.0, 1.0, 32)\n",
      "(0.0, 0.0, 1.0, 64)\n",
      "(0.0, 0.0, 1.0, 128)\n",
      "(0.0, 0.0, 1.0, 256)\n",
      "(0.1, 0.1, 1.0, 2)\n",
      "(0.1, 0.1, 1.0, 4)\n",
      "(0.1, 0.1, 1.0, 8)\n",
      "(0.1, 0.1, 1.0, 16)\n",
      "(0.1, 0.1, 1.0, 32)\n",
      "(0.1, 0.1, 1.0, 64)\n",
      "(0.1, 0.1, 1.0, 128)\n"
     ]
    }
   ],
   "source": [
    "for key, element in result_dict.items():\n",
    "    Jx, Jy, Jz, chi = key\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "# some api for the file operation\n",
    "import h5py\n",
    "from tenpy.tools import hdf5_io\n",
    "import os.path\n",
    "\n",
    "# functools\n",
    "from functools import wraps\n",
    "\n",
    "# path\n",
    "from pathlib import Path\n",
    "\n",
    "def save_data(\n",
    "    result, \n",
    "    S,\n",
    "    Lx, \n",
    "    Ly, \n",
    "    chi,\n",
    "    Jx, \n",
    "    Jy, \n",
    "    Jz, \n",
    "    prefix='data/', \n",
    "):\n",
    "    file_name = full_path(chi, Jx, Jy, Jz, L, prefix=prefix, suffix='.h5')\n",
    "    with h5py.File(file_name, 'w') as f:\n",
    "        hdf5_io.save_to_hdf5(f, result)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "Jx = 1.0\n",
    "Jy = 1.0\n",
    "Jz = 1.0\n",
    "\n",
    "# isotropic\n",
    "J = 1.0\n",
    "\n",
    "# S = 1.5\n",
    "\n",
    "Lx = 2\n",
    "Ly = 3\n",
    "shift = 0\n",
    "bc = 'periodic'\n",
    "bc_MPS = 'infinite'\n",
    "\n",
    "chi_list = [2**i for i in range(1, 9)]\n",
    "\n",
    "result_dict = {}\n",
    "\n",
    "initial_psi = None\n",
    "for chi in chi_list:\n",
    "    param = (J, J, J, chi)\n",
    "    if param not in result_dict:\n",
    "        result = run_atomic(\n",
    "            Jx=J, Jy=J, Jz=J, \n",
    "            Lx=Lx, Ly=Ly, \n",
    "            chi=chi,\n",
    "            shift=shift,\n",
    "            bc=bc, bc_MPS=bc_MPS,\n",
    "            initial_psi=initial_psi,\n",
    "        )\n",
    "    else:\n",
    "        result = result_dict[param]\n",
    "    initial_psi = result['psi'].copy()\n",
    "    result_dict[param] = result\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}