mu11 is not used, you could explicitly compute m00, m10 and m01:

cdef double m00, m01, m10

for
...
m00 += ...
m10 += ...
...