A matrix approach to tomographic reconstruction and its implementation on GPUs