inference.h

Go to the documentation of this file.

/* =====================================================================================
 
                      The UQ Toolkit (UQTk) version 3.1.5
                          Copyright (2024) NTESS
                        https://www.sandia.gov/UQToolkit/
                        https://github.com/sandialabs/UQTk
 
     Copyright 2024 National Technology & Engineering Solutions of Sandia, LLC (NTESS).
     Under the terms of Contract DE-NA0003525 with NTESS, the U.S. Government
     retains certain rights in this software.
 
     This file is part of The UQ Toolkit (UQTk)
 
     UQTk is open source software: you can redistribute it and/or modify
     it under the terms of BSD 3-Clause License
 
     UQTk 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
     BSD 3 Clause License for more details.
 
     You should have received a copy of the BSD 3 Clause License
     along with UQTk. If not, see https://choosealicense.com/licenses/bsd-3-clause/.
 
     Questions? Contact the UQTk Developers at https://github.com/sandialabs/UQTk/discussions
     Sandia National Laboratories, Livermore, CA, USA
===================================================================================== */
 
#include <string>
 
#include "Array1D.h"
#include "Array2D.h"
 
 
void infer_model(Array1D< Array2D<double> (*)(Array2D<double>&, Array2D<double>&, Array2D<double>&, void *) > forwardFuncs, void* funcInfo,
        string likType,
        string priorType, double priora, double priorb,
        Array2D<double>& xdata,Array1D<Array1D<double>>& ydata, Array2D<double>& xgrid,
        int dataNoiseInference, Array1D<double>& datanoise_array,
        int pdim,int order,Array1D<int>& rndInd,Array2D<double>& fixIndNom,string pdfType,string pcType,
        int seed, int nmcmc, double mcmcgamma, bool optimflag, Array1D<double>& chstart, Array1D<double>& chsig,
        double likParam, int likParam_int,
        Array2D<double>& pgrid,Array2D<double>& pchain, int nburn, int nstep,
        Array1D<double>& mapparam, Array1D<double>& datavar_map,
        Array1D<double>& pmean_map, Array1D<double>& pvar_map,
        Array1D<double>& fmean_map, Array1D<double>& fvar_map,
        Array1D<double>& postave_datavar,
        Array1D<double>& p_postave_mean, Array1D<double>& p_postave_var, Array1D<double>& p_postvar_mean,
        Array2D<double>& f_postsam_mean, Array1D<double>& f_postave_mean, Array1D<double>& f_postave_var, Array1D<double>& f_postvar_mean,
        Array2D<double>& paramPCcfs);
 
double LogPosterior(Array1D<double>& m, void* mypost_void);