start a bit C
alloc matrices
matrix indexing by iy, ix
correct Sx, Sy init