Choice of the FOP, mesh and rhomap arguments for the DCIPMSmoothModelling fct

Hi all,

I’d like to adapt the TDIP class from BERT to my problem (roots imaging using the Mise-à-la-masse approach). I’m using the TDIP manager to generate secondary voltages decays (all gates) for synthetic modeling/inversion. I tested the function mipmodelimg.DCIPMSmoothModelling and it works well only using the FOP defined for the inversion.

The code steps are:

# define survey geometry
# define scheme
# define rhomap and charg_map

ert_fwd = ert.ERTManager(mesh=mesh_fwd)
data_ert = ert_fwd.simulate(mesh_fwd, scheme=scheme, res=rhomap,...)
data_ert.save('ERTfwd.dat')
mgr = ert.ERTManager('ERTfwd.dat')
model = mgr.invert(lam=20, verbose=True)
modelPD = mgr.paraModel(model)  
meshPD = pg.Mesh(mgr.paraDomain) 
fop_mgr = mgr.fop # take the ERT fop
TDmanager = mip.DCIPMSmoothModelling(f=fop_mgr,mesh=meshPD,rho=modelPD,t=t) 
chargeability = TDmanager.response(m=charg_map)

I'd like to be able to skip the inversion, and run it on the forward mesh instead (or another mesh of my choice) such as:

fop_fwd_ERT=ert_fwd.fop # take the ERT fop
fop_fwd_ERT.mesh=mesh_fwd # make sure the mesh exist 
fop_fwd_ERT.data=data_ert # make sure the data exist 
TDmanager = mip.DCIPMSmoothModelling(f=fop_fwd,mesh=mesh_fwd,rho=rhomap,t=t) 

But this way it returns a mesh error.

Is there any physical reason that prevents the use of the ERT manager forward operator to be parsed into the DCIPMModelling? Attached TDIP_gitlab_issue.py the full code to reproduce the error (by uncommenting L152-155).

Thank you for your work and your nice support.

Benjamin