Add conductivity to implicit voltage solve step for stability. (#735)
Changes implicit solve step from: > solve (c/δt + L) v' = c/δt v - J to > solve (c/δt + g + L)v' = (c/δt +g) v - J where _c_ is capacitance, _g_ is membrane conductance, _J_ is membrane current. * Compute conductivity contribution for mechanisms from symbolic d/dv of current contribution (extracted from linearity test.) * Add new modcc 'source kind' for conductivity; tie to `vec_g_`. * Add conductivity field to fvm shared state. * Include conductivity in matrix assemblies for solution. Fixes #633.
Showing
- arbor/backends/gpu/matrix_assemble.cu 18 additions, 10 deletionsarbor/backends/gpu/matrix_assemble.cu
- arbor/backends/gpu/matrix_fine.cu 7 additions, 4 deletionsarbor/backends/gpu/matrix_fine.cu
- arbor/backends/gpu/matrix_fine.hpp 1 addition, 0 deletionsarbor/backends/gpu/matrix_fine.hpp
- arbor/backends/gpu/matrix_state_fine.hpp 6 additions, 4 deletionsarbor/backends/gpu/matrix_state_fine.hpp
- arbor/backends/gpu/matrix_state_flat.hpp 3 additions, 2 deletionsarbor/backends/gpu/matrix_state_flat.hpp
- arbor/backends/gpu/matrix_state_interleaved.hpp 5 additions, 3 deletionsarbor/backends/gpu/matrix_state_interleaved.hpp
- arbor/backends/gpu/mechanism.cpp 1 addition, 0 deletionsarbor/backends/gpu/mechanism.cpp
- arbor/backends/gpu/mechanism_ppack_base.hpp 1 addition, 0 deletionsarbor/backends/gpu/mechanism_ppack_base.hpp
- arbor/backends/gpu/shared_state.cpp 4 additions, 0 deletionsarbor/backends/gpu/shared_state.cpp
- arbor/backends/gpu/shared_state.hpp 13 additions, 12 deletionsarbor/backends/gpu/shared_state.hpp
- arbor/backends/multicore/matrix_state.hpp 10 additions, 7 deletionsarbor/backends/multicore/matrix_state.hpp
- arbor/backends/multicore/mechanism.cpp 1 addition, 0 deletionsarbor/backends/multicore/mechanism.cpp
- arbor/backends/multicore/mechanism.hpp 1 addition, 0 deletionsarbor/backends/multicore/mechanism.hpp
- arbor/backends/multicore/shared_state.cpp 4 additions, 0 deletionsarbor/backends/multicore/shared_state.cpp
- arbor/backends/multicore/shared_state.hpp 8 additions, 7 deletionsarbor/backends/multicore/shared_state.hpp
- arbor/fvm_lowered_cell_impl.hpp 1 addition, 1 deletionarbor/fvm_lowered_cell_impl.hpp
- arbor/matrix.hpp 2 additions, 2 deletionsarbor/matrix.hpp
- modcc/identifier.hpp 1 addition, 0 deletionsmodcc/identifier.hpp
- modcc/module.cpp 21 additions, 2 deletionsmodcc/module.cpp
- modcc/printer/printerutil.cpp 3 additions, 0 deletionsmodcc/printer/printerutil.cpp
Please register or sign in to comment