Ну в gabx это реализовано поинтереснее, производительность заметно выше, чем в Вашем примере... Тем более, что там использованы инструкции SSE/SSE2 (как алтернатива).