EqSolver.cpp File Reference

---

Detailed Description

Implementation file for class : EqSolver.

See also:

EqSolver.h

Author:

Kirill Miklyaev kirillm@iris1.itep.ru

Date:

2002-09-14

Definition in file EqSolver.cpp.

#include <stdlib.h>
#include <stdio.h>
#include "GaudiKernel/ToolFactory.h"
#include "GaudiKernel/MsgStream.h"
#include "gsl/gsl_vector.h"
#include "gsl/gsl_multimin.h"
#include "gsl/gsl_math.h"
#include "gsl/gsl_matrix.h"
#include "gsl/gsl_linalg.h"
#include "gsl/gsl_blas.h"
#include "gsl/gsl_errno.h"
#include "EqSolver.h"

Include dependency graph for EqSolver.cpp:

Go to the source code of this file.

Functions
int	fun_gsl (const gsl_vector *v, void *params, gsl_vector *f)
	The function which we solver.
int	dfun_gsl (const gsl_vector *v, void *params, gsl_matrix *df)
	The gradient of the function "fun_gsl".
int	fdfun_gsl (const gsl_vector *v, void *params, gsl_vector *f, gsl_matrix *df)
	Compute both the function "fun_gsl" and gradient of this function "dfun_gsl" together.

---

Function Documentation

int @120::dfun_gsl	(	const gsl_vector *	v,
		void *	params,
		gsl_matrix *	df
	)			[static]

The gradient of the function "fun_gsl".

Definition at line 125 of file EqSolver.cpp.

  {
    EqSolver::EqSolverMisc* local =
      static_cast <EqSolver::EqSolverMisc*> (params);
    const EqSolver::Jacobi&     jac = local -> jacobi()    ;
    Argument&                   arg = local -> argument()  ;

    {for ( unsigned int i = 0; i < v->size; ++i)
      {
        arg[i] = gsl_vector_get (v, i);
      }
    }

    {for( unsigned int i = 0 ; i < v->size ; ++i )
      {
        for (unsigned int j = 0; j < v->size; ++j)
          {
            Genfun::GENFUNCTION  f = *(jac[i][j]) ;
            gsl_matrix_set      ( df , i , j , f ( arg ) ) ;
          }
      }
    }
    return GSL_SUCCESS;
  }

int @120::fdfun_gsl	(	const gsl_vector *	v,
		void *	params,
		gsl_vector *	f,
		gsl_matrix *	df
	)			[static]

Compute both the function "fun_gsl" and gradient of this function "dfun_gsl" together.

Definition at line 155 of file EqSolver.cpp.

  {
    EqSolver::EqSolverMisc* local =
      static_cast <EqSolver::EqSolverMisc*> (params);
    const EqSolver::Equations&  eqs = *(local->equations()) ;
    const EqSolver::Jacobi&     jac = local->jacobi()    ;
    Argument&                   arg = local->argument()  ;

    {for ( unsigned int i = 0; i < v->size; ++i)
      {
        arg[i] = gsl_vector_get (v, i);
      }
    }

    {for( unsigned int i = 0 ; i < v->size ; ++i )
      {
        gsl_vector_set(f, i, (*eqs[i])(arg ) ) ;
        for (unsigned int j = 0; j < v->size; ++j)
          {
            Genfun::GENFUNCTION  f = *(jac[i][j]) ;
            gsl_matrix_set      ( df , i , j , f(arg) ) ;
          }
      }
    }
    return GSL_SUCCESS;
  }

int @120::fun_gsl	(	const gsl_vector *	v,
		void *	params,
		gsl_vector *	f
	)			[static]

The function which we solver.

Definition at line 99 of file EqSolver.cpp.

  {
    EqSolver::EqSolverMisc* local =
      static_cast <EqSolver::EqSolverMisc*> (params);
    const EqSolver::Equations& eqs = *(local->equations()) ;
    Argument&                  arg = local->argument()  ;

    {for ( unsigned int i = 0; i < v->size; ++i)
      {
        arg[i] = gsl_vector_get (v, i);
      }
    }

    {for ( unsigned int i = 0; i < v->size; ++i)
      {
        gsl_vector_set(f, i, (*eqs[i])(arg));
      }
    }
    return GSL_SUCCESS;
  }