We are investigating Gap Junctions channels in cells (during cardiac arrhythmia also propagation of excitation in neural networks).
I am trying to implement our model in Julia, but first version ran about 10x slower than C++. Is this normal? because https://aaronang.github.io/2018/stl-benchmark-comparison-cpp-vs-julia/ and https://www.quora.com/Is-it-possible-for-Julia-to-perform-as-fast-as-C++-What-would-be-necessary-for-Julia-to-be-a-suitable-alternative-to-C++ states that Julia is just little slower than C++.
Model is based on Markov chain.
Our paper about model: https://www.sciencedirect.com/science/article/pii/S0006349516001600?via%3Dihub
to run - “runMc36sm.jl” https://pastebin.com/LX5jAFV1
main model - “mc36sm_ss.jl” https://pastebin.com/QyaaW3Ti
commands to run:
using Revise, Plots; plotlyjs(); includet("runMc36sm.jl"); includet("mc36sm_ss.jl"); includet("mc36sm.jl"); (vj, gj) = runGj();
also c++ version: