Please forget the dividing by 8 because a) it's wrong, should be sqrt(8). 8 is for full 2D (so after 2 steps) b) it's not needed, it would only be useful if we did inverse, too c) one has to work with 16bit anyway, so there's no range problem unless: somebody tries a saturated 8bit for small values a-h Sorry again, gruel