WIP: iterative reweighted MxNE with BCD solving MxNE, added depth weighting proposed by Haufe et al.

This pull request adds several functions to the inverse_sparse module:

- BCD code for computing MxNE
- iterative reweighted MxNE for solving L0.5/L2 mixed-norm problem
- depth weighting proposed by Haufe et al. based on source variance estimate
Do you mean only for the fixed orientation case using:
lipschitz_constant = np.diag(np.dot(G.T, G) * 1.1