TFiniteDifferences.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 TFiniteDifferences_H
#define TFiniteDifferences_H

// Uranie
#include "TStochAttribut.h"
#include "TSensitivity.h"

using std::cout;
using std::endl;
using namespace URANIE::DataServer;

namespace URANIE
{
namespace Sensitivity
{

//void SensitivityInteractiveMethod(Double_t *param, Double_t *out);      

class TFiniteDifferences: public TSensitivity
{

protected:

    TMatrixD _sMatrix;                               
    TString _sSensitiveAtt;                
    TString _sSamplingOption;              
    Int_t _nbIn;                             
    Int_t _nbAux;                        
protected:

public:
    //---------------------------------------------
    TFiniteDifferences();


#ifdef ROOT5
    TFiniteDifferences(TDataServer *tdsNominal, void *fcn, TString outputAtt,
            TString sensitiveAtt = "", TString samplingOption = "steps=1%");
#else
    TFiniteDifferences(TDataServer *tdsNominal, void (*fcn)(Double_t*,Double_t*),
            TString sensitiveAtt, TString outputAtt,
            TString samplingOption = "steps=1%");
#endif


    TFiniteDifferences(TDataServer *tdsNominal, const char* fcn,
                       TString sensitiveAtt, TString outputAtt,
                       TString samplingOption );

 

    TFiniteDifferences(TDataServer *tdsNominal, URANIE::Launcher::TCode *fcode,
            TString sensitiveAtt = "", TString samplingOption = "steps=1%");


    TFiniteDifferences(TDataServer *tdsNominal, URANIE::Relauncher::TRun *run,
            TString sensitiveAtt = "", TString samplingOption = "steps=1%");


    TFiniteDifferences(TDataServer *tdsNominal, TString inputAtt,
            TString outputAtt, TString sensitiveAtt = TString(""));

    virtual ~TFiniteDifferences();

    //---------------------------------------------
    TMatrixD getSensitivityMatrix()
    {
        return _sMatrix;
    }

    TDataServer* getDataServer()
    {
        return _tds;
    }

    //---------------------------------------------

    virtual void generateSample(Option_t * option = "");


    //---------------------------------------------

//       void computeIndexes(Option_t * option = "");
    void evaluateIndexes(Option_t * option = "");


    void preTreatment();
    
    //---------------------------------------------
    virtual void printLog();
    ClassDef(URANIE::Sensitivity::TFiniteDifferences, ID_SENSITIVITY)
    //Definition of an variable

};

} // Fin du namespace Sensitivity
} // Fin du namespace URANIE

#endif