  We show the flux tubes produced by static quark-antiquark, quark-quark and
quark-gluon charges at finite temperature. The sources are placed in the
lattice with fundamental and adjoint Polyakov loops. We compute the square
densities of the chromomagnetic and chromoelectric fields above and below the
phase transition. Our results are gauge invariant and produced in pure gauge
SU(3). The codes are written in CUDA and the computations are performed with
GPUs.
