Add gnuplot script.
This commit is contained in:
parent
7a948bd6d2
commit
49674dbde8
|
@ -0,0 +1,7 @@
|
||||||
|
#!/usr/bin/env gnuplot
|
||||||
|
set terminal pdf enhanced
|
||||||
|
|
||||||
|
set output "mutualinformation.pdf"
|
||||||
|
set xlabel "t"
|
||||||
|
set ylabel "Mutual information"
|
||||||
|
plot "mutualinformation" using 1:2 title "" with lines lw 3, 0 title ""
|
|
@ -31,16 +31,16 @@ int main()
|
||||||
double tPlus = 0; // time on the left boundary
|
double tPlus = 0; // time on the left boundary
|
||||||
double tMinus = 0; // time on the right boundary
|
double tMinus = 0; // time on the right boundary
|
||||||
alpha = 0.4; // encodes the conformal dimention of local operator h_Psi
|
alpha = 0.4; // encodes the conformal dimention of local operator h_Psi
|
||||||
double y = 3; // endpoint of interval A
|
double y = 1; // endpoint of interval A
|
||||||
double L = 15; // length of interval A
|
double L = 5; // length of interval A
|
||||||
double epsilon = 0.01; // smearing parameter
|
double epsilon = 0.01; // smearing parameter
|
||||||
beta = 10; // inverse temperature
|
beta = 10; // inverse temperature
|
||||||
// double tOmega = /*5.4418*/; // thermal state is perturbed by operator inserted at time -tOmega
|
double tOmega = 4; // thermal state is perturbed by operator inserted at time -tOmega
|
||||||
double c = 600; // central charge. Must be large in our approximation
|
double c = 600; // central charge. Must be large in our approximation
|
||||||
// End of parameters
|
// End of parameters
|
||||||
for(unsigned int i = 0; i < 1000; ++i)
|
for(unsigned int i = 0; i < 1000; ++i)
|
||||||
{
|
{
|
||||||
double tOmega = i*2*L/1000;
|
tOmega = i*2*L/1000;
|
||||||
|
|
||||||
// Operator insertion points: Left boundary
|
// Operator insertion points: Left boundary
|
||||||
std::complex<double> x1 (0, -epsilon), x4 (0, epsilon), x1bar, x4bar;
|
std::complex<double> x1 (0, -epsilon), x4 (0, epsilon), x1bar, x4bar;
|
||||||
|
@ -79,7 +79,8 @@ int main()
|
||||||
double S_thermal = 2*c/3 * log(sinh(M_PI*L/beta)/cosh(M_PI/beta*(tMinus-tPlus)));
|
double S_thermal = 2*c/3 * log(sinh(M_PI*L/beta)/cosh(M_PI/beta*(tMinus-tPlus)));
|
||||||
double I = S_A + S_B - S_union + S_thermal;
|
double I = S_A + S_B - S_union + S_thermal;
|
||||||
std::cout << tOmega << "\t" << I << std::endl;
|
std::cout << tOmega << "\t" << I << std::endl;
|
||||||
}
|
}
|
||||||
|
|
||||||
return 0;
|
return 0;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue