Wpis z mikrobloga

@Kernel_poza_kontrolo: Tak w ogóle to co to za typy użyte do utworzenia H i G ?
Array i Matrix różnią się bodaj interpretacją \*, array uznaje to za element-wise.

No i chyba G[np.arange(imin,imax,1)-i+w, np.arange(jmin,jmax,1)-i+w] , skoro masz indeksowanie od zera.
@Kernel_poza_kontrolo: Poczytałbyś dokumentację, a nie mnie męczył :).
np.arange(imin,imax) działa dość dziwnie:
np.arange(3,7)
array([3, 4, 5, 6])
Czyli pewnie musisz poprawić na coś jak:
G[np.arange(imin,imax+1)-i+w, np.arange(jmin,jmax+1)-i+w]

i to chyba rozwiąże problem z pierwszym wymiarem tzn imin,imax+1, natomiast w drugim jakby jmin==jmax.

Ale to mój ostatni strzał.