Hello, I am trying to perform a two-dimensional joint inversion of resistivity method and seismic straight-ray tomography imaging using the PGI framework. I used the existing examples in simpeg for their forward simulations, but encountered the following issues during the PGI joint inversion.
pgi_model = inv.run(m0)
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\utils\counter_utils.py”, line 136, in wrapper
out = f(self, *args, **kwargs)
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\inversion.py”, line 111, in run
self.directiveList.call(“initialize”)
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\directives\directives.py”, line 328, in call
getattr(r, ruleType)()
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\directives\directives.py”, line 849, in initialize
dm_eigenvalue_list += [eigenvalue_by_power_iteration(dm, m)]
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\utils\mat_utils.py”, line 196, in eigenvalue_by_power_iteration
fields_list += [obj.simulation.fields(model)]
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\electromagnetics\static\resistivity\simulation.py”, line 106, in fields
A = self.getA()
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\electromagnetics\static\resistivity\simulation.py”, line 522, in getA
MeSigma = self.MeSigma
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\base\pde_simulation.py”, line 112, in Me_prop
prop = getattr(self, arg.lower())
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\props.py”, line 157, in fget
return mapping * self.model
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\maps.py”, line 249, in mul
return self._transform(val)
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\maps.py”, line 471, in _transform
m = map_i * m
File “D:\anaconda\envs\pg\lib\site-packages\SimPEG\maps.py”, line 244, in mul
raise ValueError(
ValueError: Dimension mismatch in InjectActiveCells(560,560) and np.ndarray(1120,).
mesh.nC == 560
During the inversion process, the shape of the input model is incorrect. May I ask how to fix this issue?