def evaluate_g1( l, s1 ):
    """
    Evaluate the first constraint equation and also return the jacobians

    :param float l: The value of the modulus lambda
    :param float s1: The constraint value
    """

    return l - s1**2, {'lambda':1., 's1':- 2 * s1 }