=Gubser solution for 2nd order viscous hydrodynamics=
This is a python version, implemented by cffi. One can install it by
pip install gubser
usage in ipython/python:
from gubser import Gubser
help(Gubser)
This will print the help information for Gubser class:
Example 0:
from gubser import Gubser
# define one Gubser object
gub = Gubser(L=5, tau=1, lam1=-10, etaos=0.08,
x=(-10, 10, 100), y=(-10, 10, 100))
# save x, y, ed, vx, vy, non-zero pi^{mu nu} components
gub.save_data()
# save gubser solution at a new propertime time
gub.save_data(tau=1.01)
Example 1:
from gubser import Gubser
# define one Gubser object
gub = Gubser(L=5, tau=1, lam1=-10, etaos=0.08)
# set the lattice for the ed, vx, vy, pi^{mu nu} output
gub.set_range(x0=-10, x1=-10, nx=100, y0=-10, y1=10, ny=100)
vx = gub.velocity(along='x')
print vx
# if you installed matplotlib
import matplotlib.pyplot as plt
plt.imshow(vx.T, origin='lower')
plt.show()
Example 2:
from gubser import Gubser
# define one Gubser object
gub = Gubser()
vx = gub.velocity(along='x')
print vx
# if you installed matplotlib
import matplotlib.pyplot as plt
plt.imshow(vx.T, origin='lower')
plt.show()