search

eitcom / pyEIT

Python based toolkit for Electrical Impedance Tomography
Other
177 stars 100 forks source link

The v0 calculated by the forward problem seems to be wrong #95

Closed SimonZhong22 closed 1 year ago

SimonZhong22 commented 1 year ago

When the value of perm defaults to all 1, either through fwd.compute_jac() or fwd.solve_eit(), the solved v0 is clearly wrong. As an example, when the number of electrodes is 16, using the method of adjacent excitation-adjacent measurement, the calculated v0 is not all U-shaped, but appears to plateau at some locations. Figure_1 This is the v0 curve run with pyEIT, but the red circle in the figure below is clearly not correct. 1695655788834 This is not a coincidence! I troubleshot that the problem may be in the solve_vectorizedfunction in pyeit/fem.py, where solving a system of linear equations using the sparse matrices kg and b appears rather obviously off at i=11. Could you please troubleshoot and fix the problem?

liubenyuan commented 1 year ago

@chenmo1515 Hi, any progress on this issue on closing it?

SimonZhong22 commented 1 year ago

Sorry to bother you. When I wrote this comment, I realized that the problem I was talking about was due to a different order of excitation-measurement, that's my mistake.