MatOpC

cimport cython


cdef class MatOpC:
    #no need to declare the matrices anymore, they're declared inside the pxd file
    #cdef public double[:, ::1] matA, matB
    
    def __init__(self, A, B):
        self.matA = A
        self.matB = B
        
    @cython.boundscheck(False)
    @cython.wraparound(False)
    cpdef cumsumC(self):
        cdef int i, j, lA, lB
        cdef float result = 0
        
        lA = len(self.matA)
        lB = len(self.matB)
        for i in range(0,lA):
            for j in range(0,lB):
                result = result + self.matA[i,j] * self.matB[i,j]
        return(result)