Power law regularization

[parenti1andrea]

I was wondering whether it would be possible to add the option of performing the regularization with a power law fit function, specifying the desired spectral index.