TOptimizerOpt.h

Go to the documentation of this file.

 
// Copyright (C) 2013-2024 CEA/DES
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published
// by the Free Software Foundation, either version 3 of the License, or any
// later version.
//
// This program 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 Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
// $Id$
// $Author$
// $Date$
// $Revision$
// $State$
 
#ifndef TOPTIMIZEROPT_H
#define TOPTIMIZEROPT_H
 
#include <iostream>
 
//ROOT
#include "Rtypes.h"
#include "TMethodCall.h"
 
// Uranie
#include "Optimizer.h"
#include "TDataServer.h"
// Newmat
 
#include "newmat.h"
#include "opt++.h"
//using namespace OPTPP;
#ifdef HAVE_NAMESPACES
using namespace NEWMAT;
#endif
 
// Opt++
//extern typedef void(* OPTPP::USERFCN2)(int, int, const NEWMAT::ColumnVector &, real &, NEWMAT::ColumnVector &, NEWMAT::SymmetricMatrix &, int &);
 
namespace URANIE
{
namespace Optimizer
{
class TOptimizerOpt: public TNamed
{
public:
    enum EOptimType
    {
        kFunction, kSumOfSquare
    };
private:
    Bool_t _blog;                          
    URANIE::DataServer::TDataServer *_tds;     
 
//fga 2015/01/12    void * _function;                       //! rem : Do not process by rootcint
//fga 2015/01/12    OPTPP::USERFCN2 _MyOptFunction;         //! rem : Do not process by rootcint
//fga 2015/01/12    OPTPP::OPTPPINITFCN _MyOptInitFunction; //! rem : Do not process by rootcint
 
    TString _sLibrary;
    TString _sFunctionName;
    TString _sInitFunctionName;
 
    Int_t _niter;
 
    //      OPTPP::USERFCN0 *_theFunction ;
    TMethodCall *_fMethodCall; 
    void (*_fMethod)(Double_t *gin, Double_t *f);
//       void         (*_fMethodOpt)(int mode, int n, const NEWMAT::ColumnVector& x, double& fx, NEWMAT::ColumnVector& g, NEWMAT::SymmetricMatrix& H, int &result);
    void (*_fMethodOpt)(int n, const ColumnVector& x, double& fx, int& result);
 
    TString _sname;
public:
    //---------------------------------------------
 
    TOptimizerOpt();
//       TOptimizerOpt (URANIE::DataServer::TDataServer *tds, void *function = NULL);
    TOptimizerOpt(URANIE::DataServer::TDataServer *tds, const char *functionname);
    TOptimizerOpt(URANIE::DataServer::TDataServer *tds, TString slibrary,
            TString sfunction, TString sinit = "");
    virtual ~TOptimizerOpt();
 
    //---------------------------------------------
 
    virtual void init(Option_t *option = "");
    //      virtual void optimize(void *fj, Option_t *option ="");
    virtual void optimize(Option_t *option = "");
    virtual void terminate(Option_t *option = "");
    virtual void clean(Option_t *option = "");
 
    void run(Option_t *option = "");
//      void run();
 
    TMethodCall *getMethodCall() const
    {
        return _fMethodCall;
    }
    void incrementIteration()
    {
        _niter++;
    }
    Int_t getIteration()
    {
        return _niter;
    }
 
    //---------------------------------------------
    void setLog()
    {
        _blog = kTRUE;
    }
    void unsetLog()
    {
        _blog = kFALSE;
    }
    void changeLog()
    {
        _blog = _blog ? kFALSE : kTRUE;
    }
    Bool_t getLog()
    {
        return _blog;
    }
    virtual void printLog(Option_t *option = "");
 
    ClassDef(URANIE::Optimizer::TOptimizerOpt, ID_OPTIMIZER)
    //Classe de
};
 
} // Fin du namespace Optimizer
} // Fin du namespace URANIE
 
R__EXTERN URANIE::Optimizer::TOptimizerOpt *gOptOptimizer;
#endif