LCPLaw.hpp

/*
 * This file is part of bogus, a C++ sparse block matrix library.
 *
 * Copyright 2013 Gilles Daviet <gdaviet@gmail.com>
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at http://mozilla.org/MPL/2.0/.
*/

#ifndef BOGUS_LCPLAW_HPP
#define BOGUS_LCPLAW_HPP

#include <cmath>

namespace bogus
{


template < typename Scalar >
class LCPLaw
{
public:
    enum{ dimension = 1 } ;

    typedef LocalProblemTraits< dimension, Scalar > Traits ;

    LCPLaw( ) {}

    Scalar eval( const unsigned problemIndex,
                 const typename Traits::Vector &x,
                 const typename Traits::Vector &y ) const ;


    bool solveLocal(
            const unsigned problemIndex,
            const typename Traits::Matrix &A,
            const typename Traits::Vector &b,
            typename Traits::Vector &x,
            const Scalar scaling
            ) const ;

    void projectOnConstraint( const unsigned problemIndex, typename Traits::Vector &x ) const ;

    template< typename Segment >
    void dualityCOV( const unsigned , const Segment& ,
                     typename Traits::Vector& s ) const
    { s->setZero() ; }
} ;


}

#endif