funclib.h

Go to the documentation of this file.

// This program 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 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 General Public License for more details.

// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.

#ifndef _FUNCLIB_H
#define _FUNCLIB_H

// common definitions needed here:
#include "./../stf.h"

namespace stf {



    double fexp(double x, const Vector_double& p);
    

    Vector_double fexp_jac(double x, const Vector_double& p);


    void fexp_init(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt, Vector_double& pInit );


    void fexp_init2(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt, Vector_double& pInit );
    

    double fexpde(double x, const Vector_double& p);

#if 0


    Vector_double fexpde_jac(double x, const Vector_double& p);
#endif
    

    void fexpde_init(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt,  Vector_double& pInit );
    
    

    double fexpbde(double x, const Vector_double& p);

#if 0


    Vector_double fexpbde_jac(double x, const Vector_double& p);
#endif
    

    void fexpbde_init(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt,  Vector_double& pInit );


    double falpha(double x, const Vector_double& p);
    

    Vector_double falpha_jac(double x, const Vector_double& p);
    

    double fHH(double x, const Vector_double& p);


    double fgauss(double x, const Vector_double& p);

    Vector_double fgauss_jac(double x, const Vector_double& p);


    double fgnabiexp(double x, const Vector_double& p);


    Vector_double fgnabiexp_jac(double x, const Vector_double& p);


    void falpha_init(const Vector_double& data, double base, double peak, double RTLoHI, double HalfWidth, double dt, Vector_double& pInit );


    void fgauss_init(const Vector_double& data, double base, double peak, double RTLoHI, double HalfWidth, double dt, Vector_double& pInit );


    void fHH_init(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt, Vector_double& pInit );


    void fgnabiexp_init(const Vector_double& data, double base, double peak, double RTLoHi, double HalfWidth, double dt, Vector_double& pInit );


    double xscale(double param, double xscale, double xoff, double yscale, double yoff);


    double xunscale(double param, double xscale, double xoff, double yscale, double yoff);


    double yscale(double param, double xscale, double xoff, double yscale, double yoff);


    double yscaleoffset(double param, double xscale, double xoff, double yscale, double yoff);


    double yunscale(double param, double xscale, double xoff, double yscale, double yoff);


    double yunscaleoffset(double param, double xscale, double xoff, double yscale, double yoff);


    std::vector<stf::parInfo> getParInfoExp(int n_exp);
    

    stf::Table outputWTau(const Vector_double& p, const std::vector<stf::parInfo>& parsInfo, double chisqr);
    

    std::size_t whereis(const Vector_double& data, double value);


    std::vector<stf::storedFunc> GetFuncLib();

}

#endif