Dependencies: Use ItensorMPS

Author: iagoleal

Project.toml

@@ -4,8 +4,8 @@ authors = ["Pedro Maciel Xavier <pedroxavier@psr-inc.com>", "Iago Leal de Freita
 version = "0.1.0"
 
 [deps]
+ITensorMPS = "0d1a4710-d33b-49a5-8f18-73bdf49b47e2"
+ITensors = "9136182c-28ba-11e9-034c-db9fb085ebd5"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-ITensors      = "9136182c-28ba-11e9-034c-db9fb085ebd5"
-# QUBOVerse
-QUBODrivers   = "a3f166f7-2cd3-47b6-9e1e-6fbfe0449eb0"
-QUBOTools     = "60eb5b62-0a39-4ddc-84c5-97d2adff9319"
+QUBODrivers = "a3f166f7-2cd3-47b6-9e1e-6fbfe0449eb0"
+QUBOTools = "60eb5b62-0a39-4ddc-84c5-97d2adff9319"

src/TenSolver.jl

@@ -1,6 +1,7 @@
 module TenSolver
 
 import ITensors
+import ITensorMPS
 import QUBODrivers
 import QUBOTools
 

src/solver.jl

@@ -1,4 +1,4 @@
-using ITensors
+using ITensors, ITensorMPS
 
 function issquare(A :: AbstractMatrix)
   nrows, ncols = size(A)
@@ -23,7 +23,7 @@ function tensorize_qubo(Q :: AbstractMatrix{T}, sites; cutoff = 1e-8) where {T}
   os = ITensors.OpSum{T}()
   # Construct the Hamiltonian H = Σ Q_ij P_i P_j
   # The less operators in the sum, the fastest we can calculate an MPO.
-  # We use the following simmetries to simplify the construction:
+  # We use the following symmetries to simplify the construction:
   # - For x in Bool, x^2 = x.
   #   Thus, the Hamiltonian can be linear in the diagonal.
   # - P_i commutes with P_j.