How to deal with -1 values in LAT output file

Question

I have been trying to compute the LAT map for a mesh with some infarct regions in it. When trying to construct the LAT map, this leads to -1's being written for these regions because they never depolarize.

When I try to load in the LAT map for prepacing, the -1's in the LAT file cause this error:

L3 : LAT data is not complete. Skipping prepacing.

If I would like to keep the prepacing, what would be the best way to address this issue? I have personally thought of just changing all the -1's to 0's, but am wondering if there is a better workaround.

Answer 1

when dealing with non-depolarizing regions in the LAT map, it is important to consider how they will affect different stages of the simulation, such as construction of the LAT map and prepacing. One possible solution is to replace invalid values with a default value while keeping in mind potential impacts on accuracy.

Answer 2

What's your opinion on how prepacing should be considered in these cells? No prepacing at all?

Comments

This is a bit late but I think I would be fine with no prepacing in non-depolarizing cells. Perhaps there may be slight differences in state variable values but I don't believe they should cause a major issue. The main point for me here would be to preserve the prepacing in the viable cells.

---

I think this is a reasonable suggestion. Please feel free to open a feature request (for example ignoring the -1 values and just using the default initial states in these nodes):
https://opencarp.org/community/contribute#how-to-suggest-a-feature-or-enhancement

If you want to give it a try implementing this yourself, this should be the right place in the code: https://git.opencarp.org/openCARP/openCARP/-/blob/master/physics/electrics.cc#L1798