Procedure filtreg_hemisph is always called with sdd2 = 1 / sdd1. It would not make sense otherwise. I do not like that this does not appear explicitly in filtreg_hemisph. So replaced the two arguments sdd1 and sdd2 by a single argument sdd. The cost is computing a division instead of multplication for every 3-dimensional position. Besides, this changes the results.