elementwise, where Fp‚q denotes the discrete Fourier transform and Z 1ptq " X ptq ´Φptq Z {ρ. B E is the convolution operator for the derivative with respect to the energy of X-rays, and B J E is its transpose. The role of updating Z is applying the second-order smoothing on spectra, as assumed in (A4).