We present a fast and accurate method for wave propagation through a set of inclined reflecting planes. It is based on the coordinate transformation in reciprocal space leading to a diffraction integral, which can be calculated only by using two 2D Fast Fourier Transforms and one 2D interpolation. The method is numerically tested, and comparisons with standard methods show its superiority in both computational speed and accuracy. The direct application of this method is found in the X-ray phase contrast imaging using the Bragg magnifier-an optics consisting of crystals asymmetrically diffracting in Bragg geometry.