def METHOD_NAME(self, params, residual_global):
    """Compute the cost of GeoIMC optimization problem
    Args:
        params (Iterator): An iterator containing the manifold point at which
        the cost needs to be evaluated.
        residual_global (csr_matrix): Residual matrix.
    """
    U = params[0]
    B = params[1]
    V = params[2]
    regularizer = 0.5 * self.lambda1 * np.sum(B**2)
    IMCProblem._computeLoss_csrmatrix(
        self.X.dot(U.dot(B)),
        V.T.dot(self.Z.T),
        self.Y.data,
        self.Y.indices,
        self.Y.indptr,
        residual_global,
    )
    cost = 0.5 * np.sum((residual_global) ** 2) / self.nSamples + regularizer
    return cost