Use random_itensor and random_mps

TODO

@@ -314,7 +314,7 @@ Done: -add Printf for
 Done: -simplify MPOpbc code
 Done: -Get rid of MPOpbs struct.
 Done: -Prevent Documter.jl from showing ITensor docs.
-Done: -Seeing flakey GS tests.  probably need to seed the randomMPS for consistency.
+Done: -Seeing flakey GS tests.  probably need to seed the random_mps for consistency.
 Done: -Add a big lattice test and look at SVs
 Plan:
   90% Document what I have so far

examples/apply_gates.jl

@@ -70,7 +70,7 @@ N,NNN = 10,2
 sites = siteinds("Electron", N;conserve_qns=true)
 H = reg_form_MPO(Hubbard_AutoMPO(ElT,sites, NNN))
 state = [isodd(n) ? "Up" : "Dn" for n in 1:N]
-psi = randomMPS(sites, state)
+psi = random_mps(sites, state)
 # show_directions(H)
 @show get_Dw(H)
 orthogonalize!(H,right) 

examples/directions.jl

@@ -4,7 +4,7 @@ include("../test/hamiltonians/hamiltonians.jl")
 N = 10
 sites = siteinds("S=1/2", N; conserve_qns=true)
 state = [isodd(n) ? "Up" : "Dn" for n in 1:length(sites)]
-psi = randomMPS(sites, state)
+psi = random_mps(sites, state)
 println("Show QN directions.  Arrows should point *away* from ortho centers")
 println("MPS as constructed")
 show_directions(psi)

test/blocking.jl

@@ -149,7 +149,7 @@ end
 @testset "Sub tensor assign for block sparse matrices with compatable QN spaces" begin 
     qns=[QN("Sz",0)=>1,QN("Sz",0)=>3,QN("Sz",0)=>2]
     i,j=Index(qns,"i"),Index(qns,"j")
-    A=randomITensor(i,j)
+    A=random_itensor(i,j)
     nr,nc=dims(A)
     B=copy(A)
     for dr in 0:nr-1
@@ -170,10 +170,10 @@ end
 #     qns=[QN("Sz",0)=>1,QN("Sz",0)=>3,QN("Sz",0)=>2]
 #     qnsC=[QN("Sz",0)=>2,QN("Sz",0)=>2,QN("Sz",0)=>2] #purposely miss allgined.
 #     i,j=Index(qns,"i"),Index(qns,"j")
-#     A=randomITensor(i,j)
+#     A=random_itensor(i,j)
 #     nr,nc=dims(A)
 #     ic,jc=Index(qnsC,"i"),Index(qnsC,"j")
-#     C=randomITensor(ic,jc)
+#     C=random_itensor(ic,jc)
 #     @show dense(A) dense(C)
 #     C[ic=>1:nr,jc=>1:nc]=A[i=>1:nr,j=>1:nc]
 #     @show matrix(A)-matrix(C)

test/gauge_fix.jl

@@ -31,7 +31,7 @@ models = [
   state = [isodd(n) ? "Up" : "Dn" for n in 1:N]
 
   H = MPO(Hrf)
-  psi = randomMPS(elt,sites, state)
+  psi = random_mps(elt,sites, state)
   E0 = inner(psi', H, psi)
 
   @test is_regular_form(Hrf)

test/hamiltonians.jl

@@ -118,7 +118,7 @@ end
     sites = siteinds("SpinHalf", N; conserve_qns=false)
     Hnot = three_body_AutoMPO(sites; cutoff=-1.0) #No truncation inside autoMPO
     H = three_body_AutoMPO(sites; cutoff=1e-15) #Truncated by autoMPO
-    psi = randomMPS(sites)
+    psi = random_mps(sites)
     Enot = inner(psi', Hnot, psi)
     E = inner(psi', H, psi)
     #@show E-Enot get_Dw(Hnot) get_Dw(H)

test/hamiltonians/hamiltonians.jl

@@ -456,7 +456,7 @@ end
 
 function fast_GS(H::MPO, sites, nsweeps::Int64=5)::Tuple{Float64,MPS}
   state = [isodd(n) ? "Up" : "Dn" for n in 1:length(sites)]
-  psi0 = randomMPS(sites, state)
+  psi0 = random_mps(sites, state)
   sweeps = Sweeps(nsweeps)
   setmaxdim!(sweeps, 2, 4, 8, 16, 32)
   setcutoff!(sweeps, 1E-10)

test/orthogonalize.jl

@@ -37,7 +37,7 @@ verbose1 = false #verbose inside orth algos
     sites = siteinds(model[2], N; conserve_qns=qns)
     Hrf = reg_form_MPO(model[1](elt,sites, NNN; ul=ul))
     state = [isodd(n) ? "Up" : "Dn" for n in 1:N]
-    psi = randomMPS(elt,sites, state)
+    psi = random_mps(elt,sites, state)
     E0 = inner(psi', MPO(Hrf), psi)
 
     @test is_regular_form(Hrf)

test/truncate.jl

@@ -12,8 +12,8 @@ using NDTensors: Diag, BlockSparse, tensor
 Base.show(io::IO, f::Float64) = @printf(io, "%1.1e", f)
 
 #
-#  We need consistent output from randomMPS in order to avoid flakey unit testset
-#  So we just pick any seed so that randomMPS (hopefull) gives us the same 
+#  We need consistent output from random_mps in order to avoid flakey unit testset
+#  So we just pick any seed so that random_mps (hopefull) gives us the same 
 #  pseudo-random output for each test run.
 #
 using Random
@@ -46,7 +46,7 @@ verbose1 = false #verbose inside orth algos
     sites = siteinds(model[2], N; conserve_qns=qns)
     Hrf = reg_form_MPO(model[1](elt,sites, NNN; ul=ul))
     state = [isodd(n) ? "Up" : "Dn" for n in 1:N]
-    psi = randomMPS(elt,sites, state)
+    psi = random_mps(elt,sites, state)
     E0 = inner(psi', MPO(Hrf), psi)
 
     bs = truncate!(Hrf, lr)
@@ -135,7 +135,7 @@ verbose1 = false #verbose inside orth algos
   #     H=transIsing_MPO(sites,NNN)
   #     #@show get_Dw(H)
   #     state=[isodd(n) ? "Up" : "Dn" for n=1:N]
-  #     psi=randomMPS(sites,state)
+  #     psi=random_mps(sites,state)
   #     Ea=inner(psi',Ha,psi)
   #     E0=inner(psi',H,psi)
   #     @test E0 ≈ Ea atol = eps
@@ -234,7 +234,7 @@ verbose1 = false #verbose inside orth algos
   #         si = infsiteinds("S=1/2", N; initstate, conserve_szparity=qns)
   #         ψ = InfMPS(si, initstate)
   #         for n in 1:N
-  #             ψ[n] = randomITensor(inds(ψ[n]))
+  #             ψ[n] = random_itensor(inds(ψ[n]))
   #         end
   #         H0=transIsing_iMPO(si,NNN;ul=ul)
   #         H0.llim=-1
@@ -310,7 +310,7 @@ verbose1 = false #verbose inside orth algos
   #             H=three_body_AutoMPO(sites) #Truncated by autoMPO
   #             #@show get_Dw(Hnot)
   #             Dw_auto = get_Dw(H)
-  #             psi=randomMPS(sites)
+  #             psi=random_mps(sites)
   #             Enot=inner(psi',Hnot,psi)
   #             E=inner(psi',H,psi)
   #             @test E ≈ Enot atol = sqrt(eps)