Цикл FOR наше всё.
Разумеется, тут может "напрягать", что ради вычисления каждой точки нужно целую матрицу обойти и сложить коэффициенты.
Но в случае небольших фильтров это не проблема.
Для ускорения используют преобразование Фурье и свёртку: https://habr.com/post/265781/