“…Reconstruction, on the other hand, does not assume any model or template for P R (k), but rather infers its shape from the data. Several methods have been developed for model independent reconstruction of P R (k) based on different inference methods, such as Bayesian inference, linear interpolation methods [46,64], combination of top hat functions [65], wavelet expansion of P R (k) [66,67], smoothing splines [68][69][70], fixed wavenumber knots joined with cubic splines [60,62], the critical filter method [71], different spline techniques [72][73][74] or placement of N free knots in the {k,P R } plane [60,62,[72][73][74]; penalized likelihood, using function space generalization of the Fisher matrix formalism [75] or a P R (k) ansatz [60][61][62]; Principal Component Analysis, using expansion of a orthonormal set of basis functions [76]; and sparsity of the primordial power spectrum, using a sparsity-based linear inversion method [77]. These works have been applied to CMB and/or LSS data to reconstruct P R (k) and test for deviations from the standard power law form.…”