We study the efficient numerical solution of infinite matrix equations Au = f for a matrix A in the Jaffard algebra. These matrices appear naturally via frame discretizations in many applications such as Gabor analysis, sampling theory, and quasi-diagonalization of pseudo-differential operators in the weighted Sj