algoritmo: bias nella diffusione leap

[ALTree]

In leap_step, il numero totale di molecole da diffondere nella fase di diffusione leap è determinato in funzione di tau (usando una Poissoniana), e successivamente le molecole vengono suddivise tra i vari sottovolumi adiacenti.

Se il numero di molecole non è divisibile per il numero di vicini, quelle restanti vengono divise in ordine di ID, e per questo motivo le diffusioni sono biased: le molecole si spostano più di frequente verso i vicini con ID più piccolo.

Un esempio del bias. Stato iniziale:

sbv 4: 0 0 256

stato dopo 1000 passi:

sbv 0: 0 0 38 
sbv 1: 0 0 45 
sbv 2: 0 0 40 
sbv 3: 0 0 39 
sbv 4: 0 0 39 
sbv 5: 0 0 38 
sbv 6: 0 0 15 
sbv 7: 0 0 2 
sbv 8: 0 0 0