# inv = Inversion.BaseInversion(invProb, directiveList=[IRLS, saveDict,
# update_Jacobi])
# # Run inversion
# mrec = inv.run(mstart)
# m0 = np.median(ln_sigback) * np.ones(mapping.nP)
# dmis = DataMisfit.l2_DataMisfit(survey)
n = 3
clf = GaussianMixture(
n_components=n, covariance_type='full', reg_covar=5e-3
)
clf.fit(m0.reshape(-1, 1))
Utils.order_clusters_GM_weight(clf)
print(clf.covariances_)
print(clf.means_)
idenMap = Maps.IdentityMap(nP=m0.shape[0])
wires = Maps.Wires(('m', m0.shape[0]))
reg = Regularization.SimplePetroRegularization(
GMmref=clf, mesh=mesh,
wiresmap=wires,
maplist=[idenMap],
mref=m0*0,
indActive=actv,
cell_weights=wr
)
reg.mrefInSmooth = False
reg.approx_gradient = True
