1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
|
// @(#)root/mathmore:$Id$
// Authors: L. Moneta Dec 2006
/**********************************************************************
* *
* Copyright (c) 2004 ROOT Foundation, CERN/PH-SFT *
* *
* This library is free software; you can redistribute it and/or *
* modify it under the terms of the GNU General Public License *
* as published by the Free Software Foundation; either version 2 *
* of the License, or (at your option) any later version. *
* *
* This library is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU *
* General Public License for more details. *
* *
* You should have received a copy of the GNU General Public License *
* along with this library (see file COPYING); if not, write *
* to the Free Software Foundation, Inc., 59 Temple Place, Suite *
* 330, Boston, MA 02111-1307 USA, or contact the author. *
* *
**********************************************************************/
// Header file for class GSLMultiMinFunctionWrapper
//
// Created by: moneta at Sat Nov 13 14:54:41 2004
//
// Last update: Sat Nov 13 14:54:41 2004
//
#ifndef ROOT_Math_GSLMultiFitFunctionWrapper
#define ROOT_Math_GSLMultiFitFunctionWrapper
#include "gsl/gsl_multifit.h"
#include "gsl/gsl_multifit_nlin.h" // patch JWu dec20
#include "GSLMultiFitFunctionAdapter.h"
#include <cassert>
namespace ROOT {
namespace Math {
typedef double ( * GSLMultiFitFPointer ) ( const gsl_vector *, void *, gsl_vector *);
typedef void ( * GSLMultiFitDfPointer ) ( const gsl_vector *, void *, gsl_matrix *);
typedef void ( * GSLMultiFitFdfPointer ) ( const gsl_vector *, void *, gsl_vector *, gsl_matrix *);
/**
wrapper to a multi-dim function withtout derivatives for multi-dimensional
minimization algorithm
@ingroup MultiMin
*/
class GSLMultiFitFunctionWrapper {
public:
GSLMultiFitFunctionWrapper()
{
fFunc.f = 0;
fFunc.df = 0;
fFunc.fdf = 0;
fFunc.n = 0;
fFunc.p = 0;
fFunc.params = 0;
#if GSL_MAJOR_VERSION > 1
fFunc.nevalf = 0;
fFunc.nevaldf = 0;
#endif
}
/// Fill gsl function structure from a C++ function iterator and size and number of residuals
template<class FuncVector>
void SetFunction(const FuncVector & f, unsigned int nres, unsigned int npar ) {
const void * p = &f;
assert (p != 0);
fFunc.f = &GSLMultiFitFunctionAdapter<FuncVector >::F;
fFunc.df = &GSLMultiFitFunctionAdapter<FuncVector >::Df;
fFunc.fdf = &GSLMultiFitFunctionAdapter<FuncVector >::FDf;
fFunc.n = nres;
fFunc.p = npar;
fFunc.params = const_cast<void *>(p);
}
gsl_multifit_function_fdf * GetFunc() { return &fFunc; }
private:
gsl_multifit_function_fdf fFunc;
};
} // namespace Math
} // namespace ROOT
#endif /* ROOT_Math_GSLMultiMinFunctionWrapper */
|
