jonhrafe / MCDC_Simulator_public

The Monte Carlo Diffusion and Collision simulator (MC/DC), is a C++ open-source Diffusion-Weighted Magnetic Resonance Imaging (DW-MRI) Monte Carlo Simulator.
https://jonhrafe.github.io/MCDC_Simulator_public/
GNU Lesser General Public License v2.1
19 stars 12 forks source link

Empty file output for the Gamma Sphere Distribution #14

Open mdausort opened 1 year ago

mdausort commented 1 year ago

https://github.com/jonhrafe/MCDC_Simulator_public/blob/a5ffbb0449d36f9015cc1adf9c04d17b149c61e9/src/parallelmcsimulation.cpp#L859

I tried your simulator with Cylinder Gamma Distribution to represent axonal configuration axonal configuration and now I would like to try it for Sphere Gamma Distribution. However, to be able to reproduce my experiment, I would like to save my configuration and load it for each one of my simulations. In the outputs of the simulations, there is a file named "_gamma_distributed_sphere_list.txt" but it does not contain any position or radius of the spheres of the distribution whereas for a simulation of cylinders, the file is not empty. For my test, I used the example (https://github.com/jonhrafe/MCDC_Simulator_public/blob/master/instructions/example_intra-axonal_initialization.md) and just change the parameters to switch from cylinders to spheres. Here is the config file that I used :

N 1000
T 1000
duration 0.0359000000
diffusivity 0.6e-9

exp_prefix instructions/demos/output/sphere_gamma_packing_test

scheme_file docs/scheme_files/PGSE_sample_scheme.scheme
scale_from_stu 1

write_txt 1
write_bin 0
write_traj_file 1

<obstacle>
<sphere_gamma_packing>
alpha 1.5
beta 0.5
icvf 0.70
num_spheres 1000
</sphere_gamma_packing>
</obstacle>

ini_walkers_pos extra
num_process 1
<END>

Version = https://github.com/jonhrafe/MCDC_Simulator_public/releases/tag/1.5

jonhrafe commented 1 year ago

Dear @mdausort

The issue is solved with @RGardier 's last commit.

I updated the current release to this one.

Best.