Understanding Source Localization of bst_inverse_linear[2018]

Hello!
My question may repeat topics discussed here.
But could you explain or give an article with a detailed explanation:
Why do we multiply by raw MEG data when we make an inverse operator MNE or sLORETA?
I know that the notion Regularization plays the crucial role here, but it’s not clear for me how it works.