The vectorized PR was made by @ChabaneAmaury for accelerating forward computing. It is not implemented on image code yet. But sure, it is very valuable to do so and I will create a Issue on it.

I tried to test forward solver using NY-HEAD 3d model (a very big, refined head mesh) but still failed to load the complete model into RAM. I will working on it for the next releases.

I try to load the NY-HEAD model: https://www.parralab.org/nyhead/ (where the .inp file is used) and placed electrode on the boundary, do the eit forward simulation.

As to the inverse imaging, I think a 3D inverse imaging is not feasible. The Jacobian should be calculated on a 2D and the unkowns are suffice for the inverse problem.

Hy,

Would you give me some insight on what exactly you would like to improve? >> was it for me?

Basically, you add your vector as an *arg in the EitBase constructor, save it and use it when instanciating the solvers jac, bp and greit. Nothing really special. These call your improved solve_eit method. I've already worked out this implementation, I will make a PR after I have tested it!