'Evaluate w = M^-H v'

def rsolve(self, v, tol=0):

    if (self.collapsed is not None):
        return solve(self.collapsed.T.conj(), v)
    return LowRankMatrix._solve(v, np.conj(self.alpha), self.ds, self.cs)
