ALTree / PLISCuda

PLIS on CUDA
MIT License
1 stars 0 forks source link

algoritmo: calcolo tau per le reazioni critiche #39

Closed ALTree closed 8 years ago

ALTree commented 8 years ago

Nella procedura di calcolo del tau per le reazioni critiche, il valore viene generato come

tau = -logf(rand) / react_rates_sum_cr

dove react_rates_sum_cr è la somma dei tassi di reazione delle reazioni critiche, in acccordo con il passo (5b), sezione 1, in Tao06.

Il tau dovrebbe anche tenere conto in qualche maniera dei tassi di diffusione (non trattati in Tao06, che non si riferisce a simulazioni spaziali). Altrimenti ad un sottovolume con reazioni critiche che non possono però avvenire (i.e. avente react_rates_sum_cr = 0) sarà associato un tau infinito, che non è corretto se possono avvenire diffusioni.

ALTree commented 8 years ago

Fixed.