EqSolver Class Reference

#include <EqSolver.h>

Inheritance diagram for EqSolver:

[legend]

Collaboration diagram for EqSolver:

[legend]

List of all members.

---

Detailed Description

The simplest concrete implementation of IEqSolver interface.

See also:

GaudiGSL/IEqSolver.h

Author:

Kirill Miklyaev kirillm@iris1.itep.ru

Date:

2003-07-07

Definition at line 28 of file EqSolver.h.

Public Types
typedef std::vector< Equations >	Jacobi
Public Member Functions
virtual StatusCode	solver (const Equations &funcs, Arg &arg) const
	Solving nonlinear systemwith N equations in N unknowns of the function "GenFunc".
virtual StatusCode	initialize ()
	Overriding initialize.
virtual StatusCode	finalize ()
	standard finalization method
virtual	~EqSolver ()
	Destructor.
	EqSolver (const std::string &type, const std::string &name, const IInterface *parent)
	Standard constructor.
Private Member Functions
	EqSolver ()
	default constructor is private
	EqSolver (const EqSolver &)
	copy constructor is private
EqSolver &	operator= (const EqSolver &)
	assignement operator is
Private Attributes
std::string	m_algType
double	m_max_iter
double	m_norm_residual
const gsl_multiroot_fdfsolver_type *	m_type
	AlgTool type (concrete class name).
Friends
class	ToolFactory< EqSolver >
Classes
class	EqSolverMisc

---

Member Typedef Documentation

typedef std::vector<Equations> EqSolver::Jacobi

Definition at line 35 of file EqSolver.h.

---

Constructor & Destructor Documentation

EqSolver::~EqSolver

(

)

[virtual]

Destructor.

Definition at line 344 of file EqSolver.cpp.

{}

EqSolver::EqSolver	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Standard constructor.

See also:

GaudiTool

Parameters:

	type	tool type
	name	tool name
	parent	parent of the tool

declare my special interface

declare type of the algorithm for root finding

declare maximum of iteration

declare the absolute error bound for the residual value

Definition at line 75 of file EqSolver.cpp.

  : GaudiTool       ( type, name , parent )
  , m_algType       ( "fdfsolver_hybridsj" ) 
  , m_max_iter      ( 1000                 ) 
  , m_norm_residual ( 1.0e-7               )
  , m_type          ( 0                    )
  
{
  declareInterface <IEqSolver>(this);
  declareProperty ( "Algorithm", m_algType       );
  declareProperty ( "Iteration", m_max_iter      );
  declareProperty ( "Residual" , m_norm_residual );
}

EqSolver::EqSolver

(

)

[private]

default constructor is private

EqSolver::EqSolver

(

const EqSolver &

)

[private]

copy constructor is private

---

Member Function Documentation

StatusCode EqSolver::solver	(	const Equations &	funcs,
		Arg &	arg
	)			const [virtual]

Solving nonlinear systemwith N equations in N unknowns of the function "GenFunc".

Solving nonlinear system of the function "GenFunc".

See also:

IEqSolver.h

Returns:

StatusCode

Definition at line 193 of file EqSolver.cpp.

{
  using namespace Genfun;
    
  gsl_vector_view vect = gsl_vector_view_array ( &arg[0] ,
                                                 arg.dimension() );
  MsgStream log( msgSvc(), name() );
  
  EqSolverMisc local (funcs, arg);

  const gsl_multiroot_fdfsolver_type *T = m_type;
  gsl_multiroot_fdfsolver *s;
  int    status = 0 ;
  size_t iter   = 0 ;

  gsl_multiroot_function_fdf function;
  
  function.f = &fun_gsl;
  function.df = &dfun_gsl;
  function.fdf = &fdfun_gsl;
  function.n = vect.vector.size;
  function.params = (void*) &local;

  s = gsl_multiroot_fdfsolver_alloc(T, vect.vector.size);
  gsl_multiroot_fdfsolver_set (s, &function, &vect.vector);
   
  for (iter = 0; iter < m_max_iter; ++iter)
    {
      status = gsl_multiroot_fdfsolver_iterate (s);
      if (status)
        {
          return Error 
            ("Error from gsl_gsl_multiroot_fdfsolver_iterate '"
             + std::string(gsl_strerror(status)) + "'") ;
          break;
        }
  
      status = gsl_multiroot_test_residual (s->f, 
                                            m_norm_residual);
      
      if ( status != GSL_CONTINUE ) { break; }
    }
  
  for (unsigned int i = 0; i < vect.vector.size; ++i)
    {
      gsl_vector_set (&vect.vector, i, gsl_vector_get (s->x, i));
    }
  
  if (status == GSL_SUCCESS)
    {
      log << MSG::DEBUG 
          << "We stopped in the method on the " << iter 
          << " iteration (we have maximum "     << m_max_iter 
          << " iterations)"                     << endreq;
    }
  else if (status == GSL_CONTINUE && iter <= m_max_iter )
    {
      return Error ( "Method finished with '" 
                     + std::string(gsl_strerror(status)) 
                     +  "' error" );
    }
  else
    {
      return Error ( "Method finished with '" +
                     std::string(gsl_strerror(status)) 
                     +  "' error" );
    }
  
  gsl_multiroot_fdfsolver_free (s);
  
  if (status) 
    { 
      return Error ( "Method finished with '" 
                     + std::string(gsl_strerror(status)) 
                     +  "' error" );
    }
  
  return GSL_SUCCESS;
}

StatusCode EqSolver::initialize

(

)

[virtual]

Overriding initialize.

Reimplemented from GaudiTool.

Definition at line 276 of file EqSolver.cpp.

{
  StatusCode sc = GaudiTool::initialize() ;
  
  MsgStream log( msgSvc() , name() ) ;
  
  if( sc.isFailure() ) 
    {
      return Error ("Could not initiliaze base class GaudiTool", sc);
    }
  /* The algorithm for multiional root-finding 
     (solving nonlinear systems with n equations in n unknowns)*/
  
  if(       "fdfsolver_hybridsj" == m_algType )
    {
      m_type = gsl_multiroot_fdfsolver_hybridsj  ;
      log << MSG::DEBUG 
          << "Root findind algoritms to be used: " 
          << "'gsl_multiroot_fdfsolver_hybridsj'" 
          << endreq;
    }
  else if ( "fdfsolver_hybridj" == m_algType ) 
    {
      m_type = gsl_multiroot_fdfsolver_hybridj   ;
      log << MSG::DEBUG 
          << "Root findind algoritms to be used: " 
          << "'gsl_multiroot_fdfsolver_hybridj'"
          << endreq;
    }
  else if ( "fdfsolver_newton" == m_algType )
    {
      m_type = gsl_multiroot_fdfsolver_newton    ;
      log << MSG::DEBUG 
          << "Root findind algoritms to be used: " 
          << "'gsl_multiroot_fdfsolver_newton'"
          << endreq;
    }
  else if ( "fdfsolver_gnewton" == m_algType )
    {
      m_type = gsl_multiroot_fdfsolver_gnewton   ;
      log << MSG::DEBUG 
          << "Root findind algoritms to be used: " 
          << "'gsl_multiroot_fdfsolver_gnewton'"
          << endreq;
    }
  else 
    {
      return Error(" Unknown algorith type '" 
                   + std::string(m_algType) + "'");
     }
  
  return StatusCode::SUCCESS;
}  

StatusCode EqSolver::finalize

(

void

)

[virtual]

standard finalization method

See also:

AlgTool

IAlgTool

Returns:

status code

Reimplemented from GaudiTool.

Definition at line 331 of file EqSolver.cpp.

{
  StatusCode sc = GaudiTool::finalize() ;
  
  MsgStream log( msgSvc() , name() ) ;
  
  if ( sc.isFailure() )
    {
      return Error("Could not finaliaze base class GaudiTool", sc);
    }
  return StatusCode::SUCCESS;
};

EqSolver& EqSolver::operator=

(

const EqSolver &

)

[private]

assignement operator is

---

Friends And Related Function Documentation

friend class ToolFactory< EqSolver > [friend]

Definition at line 33 of file EqSolver.h.

---

Member Data Documentation

std::string EqSolver::m_algType [private]

Definition at line 102 of file EqSolver.h.

double EqSolver::m_max_iter [private]

Definition at line 103 of file EqSolver.h.

double EqSolver::m_norm_residual [private]

Definition at line 104 of file EqSolver.h.

const gsl_multiroot_fdfsolver_type* EqSolver::m_type [private]

AlgTool type (concrete class name).

Reimplemented from AlgTool.

Definition at line 105 of file EqSolver.h.

---

The documentation for this class was generated from the following files:

• /afs/cern.ch/sw/Gaudi/releases/GAUDI/GAUDI_v20r4/GaudiGSL/src/Components/EqSolver.h
• /afs/cern.ch/sw/Gaudi/releases/GAUDI/GAUDI_v20r4/GaudiGSL/src/Components/EqSolver.cpp