def METHOD_NAME(c, d):
    """
    Least squares fit with ridge regularization.
    """
    alval = 1.e-6
    reg = Ridge(alpha = alval, fit_intercept = False)
    reg.fit(c, d)
    return reg.coef_.T # return transpose if using sklearn ridge