def METHOD_NAME(self, x: NDArray) -> NDArray:
    ncp = get_array_module(x)
    y = ncp.zeros(self.nops, dtype=self.dtype)
    for iop, oper in enumerate(self.ops):
        y += oper.matvec(x[self.mmops[iop] : self.mmops[iop + 1]]).squeeze()
    return y