Моделирование световой волны клеточным автоматом

Глобальная мегазадача - моделирование распространения световой волны, для синтез вычислительной голограммы.
В принципе, таких программок валяется по углам 100500, но они все для голограммы Габора, а вот для голограммы Денисюка, по понятным причинам - увы.

Я хочу попробовать подойти со стороным моделирования клеточным автоматом(тогда вычислительная сложность становится зависима ТОЛЬКО от размеров объемной дифракционной решетки).

Вот простенький пример: https://www.youtube.com/watch?v=noUpBKY2rIg


В принципе, там все достаточно прозрачно. Каждая клетка имеет вектор сила-ускорение-скорость-амплитуду, которые на следующем шаге рассчитываются из предыдущего с 4-х связностью:
Аn = (dH0n-1 + dH1n-1 + dH2n-1 + dH4n-1) / m
Vn = Vn-1 + A * dT
Hn = Hn-1 + Vn * dT

Вполне моделируется двумерная дифракционная решетка.
Но сразу уперся в вопрос поглощения волны. Любая неоднородность становится отражателем
И если с отражением от краёв можно как-то справиться, то внутри модулируемой решетки волну нужно поглощать на одной клетке!

P.S. Я, увы, в математике компетентен ровно на столько чтобы не путать знаки препинания.