Help with ERT FOP.response function

[GiuliadePasqualeCEAZA]

Hi all,
I'm working on the cross-gradient joint inversion of (Jordi et al., 2020) "translating" the code to the new API in pygimli and somehow got stuck with the definition of the Forward Operator for the ERT dataset.
In the orgiginal code this was how Claudio defined the ERT fop (invmesh was created with meshtools.createParaMesh from the sensors positios):

fERT = pb.DCSRMultiElectrodeModelling(invmesh, data_ert)
invmesh.createNeighbourInfos()
fERT.setMesh(invmesh) # important --> needed for inversion
fERT.region(1).setBackground(True)
fERT.createRefinedForwardMesh(True)
pd = fERT.regionManager().paraDomain()
pd.createNeighbourInfos()
fERT.createConstraints()

Now what I was trying is this:

mgr = ert.ERTManager(sr=True, useBert=True, verbose=False, debug=False)
mgr.setData(data_ert)
mgr.setMesh(invmesh) # important --> needed for inversion (automatically set region with smallest marker to background)
pd = pg.Mesh(mgr.paraDomain)
pd.createNeighbourInfos()
fERT =mgr.createForwardOperator()
fERT.setMesh(invmesh)  #also tried with pd here but same result
fERT.createConstraints()

butr when I try to use fERT.response as follow:

nM = pd.cellCount()
rhomed = pg.median(data_ert('rhoa'))
startModel_ert = pg.Vector(nM, rhomed)

response_ERT=fERT.response(startModel_ert)

it closes the computation without giving me any result or error...

what am I doing wrong?

thank you in advance

Giulia

[halbmy]

Dear Guilia, thank you for taking this task that had been on the list for a long time.

Maybe we could set up a github repository and work on the same codes. I also have some stuff ready and maybe we start with a simplistic case.

(You should not set the mesh of the forward operator explicitly as this is done by the ERT manager.)

[GiuliadePasqualeCEAZA]

Dear Thomas,
yes sure! I'm not that familiar with github but up for sharing all the little progress: will try and create a repository with what I have at the moment and share it with you