NonSmoothNewton.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_NS_NETWON_HPP
#define BOGUS_NS_NETWON_HPP

namespace bogus {


template < typename NSFunction >
class NonSmoothNewton
{
public:
  typedef typename NSFunction::Traits Traits ;
  typedef typename Traits::Scalar Scalar ;
  typedef typename Traits::Vector Vector ;
  typedef typename Traits::Matrix Matrix ;

  NonSmoothNewton( const NSFunction &func, Scalar tol )
  : m_func( func ), m_tol( tol ), m_maxIters( 20 )
  {}


  Scalar solve ( Vector &x ) const ;

  void setMaxIters( unsigned maxIters ) { m_maxIters = maxIters ; }

  void setTol( Scalar tol ) { m_tol = tol ; }

private:
  const NSFunction& m_func ;
  Scalar m_tol ;
  unsigned m_maxIters ;

} ;


}

#endif