* les paramètres minbound et maxbound étaient écrasés et utilisées
  comme des variables locales et donc les paramètres ne servaient à rien
* ajustement des distances d'entrée et de sortie avec minbound et maxbound
* dans le calcul du numéro du voxel d'entrée, il y avait des cast en
  int en plein milieu des formules qui faussaient tout
* dans le calcul initial de Dz, dans le cas "dir.z < 0", il y avait
  un "-1" superflu et faux
* dans la boucle, on en sort lorsque le voxel courant est à une
  distance supérieures à maxbound