FuncMinimum.cpp File Reference

---

Detailed Description

Implementation file for class : FuncMinimum.

See also:

FuncMinimum.h

Author:

Kirill Miklyaev kirillm@iris1.itep.ru

Date:

2002-09-14

Definition in file FuncMinimum.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 "FuncMinimum.h"

Include dependency graph for FuncMinimum.cpp:

Go to the source code of this file.

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

---

Function Documentation

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

The gradient of the function "fun_gsl".

Definition at line 110 of file FuncMinimum.cpp.

  {
    FuncMinimum::FuncMinimumMisc * local = 
      static_cast <FuncMinimum::FuncMinimumMisc*> (params);
    const FuncMinimum::Gradient& grad  = local->gradient();
    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 < df->size ; ++i )
      
      {
        Genfun::GENFUNCTION f = *(grad[i]);
        gsl_vector_set ( df, i, f(arg) );
      }
    }
    
  }

void @123::fdfun_gsl	(	const gsl_vector *	v,
		void *	params,
		double *	f,
		gsl_vector *	df
	)			[static]

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

Definition at line 136 of file FuncMinimum.cpp.

  {
    *f = fun_gsl( v , params );
    dfun_gsl ( v , params, df);
  }

double @123::fun_gsl	(	const gsl_vector *	v,
		void *	params
	)			[static]

The function which we minimize.

Definition at line 92 of file FuncMinimum.cpp.

  {
    FuncMinimum::FuncMinimumMisc * local = 
      static_cast <FuncMinimum::FuncMinimumMisc*> (params);
    const FuncMinimum::GenFunc& eq  = *(local->equation());
    Argument&                   arg = local->argument()   ;
         
    {for (unsigned int i = 0; i < v->size; ++i)
      {
        arg[i] = gsl_vector_get (v, i);
      }
    }
    
    return eq(arg);
  }