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

Particle initialisation in mesh substrate #8

Open MariamAndersson opened 3 years ago

MariamAndersson commented 3 years ago

When performing simulations on mesh substrates and using a spin initialisation file, is it possible to have some indication of the fraction of spins initialised within the mesh, and the fraction initialised outside? It could also be helpful to see how many unique walker locations are initiated. This would act as a kind of indication that initialisation has occurred as intended, without having to log the particle trajectories and check them after simulation has finished.

I submitted initialisation files in the wrong format (dimension x spin) instead of (spin x dimension) and was happily simulating free diffusion until I realised the error. I guess it appeared to the simulator as if there were thousands of dimensions (and only 3 spins) and it ignored the dimensions above 3? In this case, it could also be useful with an error message if the number of dimensions > 3.