Elevating Channel Currents to First-Class States

[huangziwei]

Currently, channel currents are not considered "first-class" states. By this, I mean that channel currents are not included in the .nodes like other states when the channel is inserted into the compartment. Only later, during integration, are the currents added to all_states through get_all_states, which allows us to record them along with other states. Not being first-class states means that a specific channel current cannot be clamped.

Well, to be honest, I don't really have a use case for this, though. It's just fun to think that one could do a current clamp on specific channels instead of the overall current clamp with stimulate().

[huangziwei]

attempted a few times but can't make it works. probably need to refactor _channel_currents() quit a bit but I am not familiar enough with the code to make the change. :D