An automatic algorithm to exploit the symmetries of the system response matrix in PET iterative reconstruction