Computing Relational Fixed Points using Interpolation