Как работает алгоритм Ричардсона-Люси? Пример кода?

Я пытаюсь понять, как работает деконволюция. Я понимаю идею, стоящую за этим, но я хочу понять некоторые из реальных алгоритмов, которые его реализуют - алгоритмы, которые принимают в качестве входных данных размытое изображение с его функцией точечной выборки (ядро размытия) и производят на выходе скрытое изображение.

На данный момент я нашел алгоритм Ричардсона–Люси, в котором математика не кажется такой уж сложной, однако я не могу понять, как работает настоящий алгоритм. В Википедии сказано:

Это приводит к уравнению, которое может быть решено итеративносогласно...

однако он не показывает фактическую петлю. Может ли кто-нибудь указать мне ресурс, где объясняется фактический алгоритм. В Google мне удалось найти только методы, которые используют Ричардсона-Люси в качестве одного из шагов, но не настоящий алгоритм Ричардсона-Люси.

Алгоритм на любом языке или псевдокод был бы хорош, однако, если бы он был доступен на Python, это было бы потрясающе.

Заранее спасибо.

Редактировать

По сути, то, что я хочу выяснить, это размытое изображение (nxm):

x00 x01 x02 x03 .. x0n
x10 x11 x12 x13 .. x1n
...
xm0 xm1 xm2 xm3 .. xmn

и ядро ​​(ixj), которое использовалось для получения размытого изображения:

p00 p01 p02 .. p0i
p10 p11 p12 .. p1i
...
pj0 pj1 pj2 .. pji

Каковы точные шаги в алгоритме Ричардсона-Люси, чтобы выяснить исходное изображение.