# Matlab code for Pu and Thomas' Neural Computation paper ## Fast and Accurate Langevin Simulations of Stochas-tic Hodgkin-Huxley Dynamics 1. Figure 1 (HH diagram) is generated by the tex code. 2. Figure 2 (4D and 14D HH models) can be generated from codes in folder "Figure2". It takes roughly 3-5 minutes to run the script Fig2.m on a laptop. 3. Figure 3 (convergence to the multinomial submanifold) can be generated by codes in folder "Figure3". It takes roughly 3-5 minutes to run the script Fig3.m on a laptop. 4. Figure 4 (Edge importance under voltage clamp) is reproduced with permission from Figs. 10 & 13 of [Schmidt and Thmas (2014) paper](https://pubmed.ncbi.nlm.nih.gov/24742077/). Code is reproduced with permission. 5. Figure 5 (edge importance under current clamp) can be reproduced from codes in folder "Figure5". - to generate figure 5, one needs to run hundreds repeated simulations, high performance computing is recommended - one can modify the code gene_data_SS.m to generate as many samples as desired - fig5data.mat can be loaded and viewed by plot_fig5.m 6. Figure 6 (pathwise equivalency) can be generated by the code Fig6.m, which takes roughly 1-2 minutes to run on a laptop. 7. Figures 7, 8 and 9, as well as codes for generating Table 3, are included in the folder "All_models" - cluster computing or high performance computing is recommended for generating the data for these tables and figures - details of the simulation for the paper is specified in Section 5 - simulation efficiency is computed through the same laptop computer; the time might be different using different machines but the ratio should be roughly the same - since the data for all plots in Figs. 7-9 and Table 3 are more than 500 MB, data for the plots are not uploaded here but one should be able to re-generate the data with the provided code