Рассматриваются объемные сингулярные интегральные уравнения, описывающие задачи рассеяния электромагнитных волн на ограниченных, неоднородных и анизотропных диэлектрических структурах. Для дискретизации уравнений используется метод коллокации на неравномерной сетке. Для быстрого умножения матрицы получающейся СЛАУ на вектор вводится равномерная сетка. С использованием алгоритмов интерполяции функций и быстрого дискретного преобразования Фурье построен алгоритм умножения матрицы СЛАУ на вектор с числом операций для его реализации практически пропорциональным размерности матрицы. Кроме того, объем памяти для реализации алгоритма также пропорционален размерности.