qubits-CUDA/g3.py at master · evalvarez12/qubits-CUDA · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
from pylab import *
import os

#PARAMETROS FIJOS
Js,DJs,Dbs,bx_no,bz=[1.,0.,0.,1.,1.]


#OBTENER DATOS

#for lamb in log10(logspace(0,0.1,100)) :
  #dat_i=".dat"
  #comando = "./abc -o purity_onet -q 17 --x 6 --model modelConexComplete --Jc "+str(lamb)+" --Jp 0.5 --Js 1. --DJs 0. --bx 1. --bz 1. --Dbs 0. --t 1000 --Cseed -1 --Eseed 3462 --PARAMseed 9678  >> grafica3-complete-gamma_0.5"
  #os.system(comando + dat_i)
  #print lamb

#for lamb in log10(logspace(0,0.1,100)) :
  #dat_i=".dat"
  #comando = "./abc -o purity_onet -q 17 --x 6 --model modelConexComplete --Jc "+str(lamb)+" --Jp 0.1 --Js 1. --DJs 0. --bx 1. --bz 1. --Dbs 0. --t 1000 --Cseed -1 --Eseed 3462 --PARAMseed 9678  >> grafica3-complete-gamma_0.1"
  #os.system(comando + dat_i)
  #print lamb

#for lamb in log10(logspace(0,0.1,100)) :
  #dat_i=".dat"
  #comando = "./abc -o purity_onet -q 17 --x 6 --model modelConexComplete --Jc "+str(lamb)+" --Jp 0.01 --Js 1. --DJs 0. --bx 1. --bz 1. --Dbs 0. --t 1000 --Cseed -1 --Eseed 3462 --PARAMseed 9678  >> grafica3-complete-gamma_0.01"
  #os.system(comando + dat_i)
  #print lamb


#HACER grafica3
fig=figure(figsize=(9, 6))
a=loadtxt("grafica3-gamma_0.5.dat")
b=loadtxt("grafica3-gamma_0.1.dat")
c=loadtxt("grafica3-gamma_0.01.dat")
lamb=log10(logspace(0,0.1,100))

plot(lamb,log10(1-a),'k-',label='$\gamma=0.5$')
plot(lamb,log10(1-b),'k--',label='$\gamma=0.1$')
plot(lamb,log10(1-c),'k:',label='$\gamma=0.01$')

xlabel("$\lambda$",fontsize=17)
#xticks(log(arange(0.000001,.1+.001,.05)),('$0$','$0.05$','$0.1$'),fontsize=15)
xscale('log',basex=10)
xticks((.001,0.01,.1),('$10^{-3}$','$10^{-2}$','$10^{-1}$'),fontsize=17)
yticks((0,-1,-2),('$0$','$-1$','$-2$'),fontsize=17)
ylabel("$\log(1-P)$",fontsize=17)

tick_params(axis='both',          # changes apply to the x-axis
    which='both',      # both major and minor ticks are affected
    bottom='on',      # ticks along the bottom edge are off
    top='on',         # ticks along the top edge are off
    labelbottom='on',
    length=10)

axis([.001,.1,-2.7,-0.1])

legend(loc='lower right',fontsize=17)

savefig('art-g3.png')