Subspace-based optimization method coupled with multiplicative regularization for edge-preserving inversion