bcs.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
===================================================================================== */
 
#ifndef BCS_H
#define BCS_H
 
#include "Array1D.h"
#include "Array2D.h"
 
 
#define MAX_IT 1000
 
 
 
void WBCS(Array2D<double> &PHI, Array1D<double> &y, Array1D<double> &sigma2,
                 double eta, Array1D<double> &lambda_init,
                 int adaptive, int optimal, double scale, int verbose,
                 Array1D<double> &weights, Array1D<int> &used,
                 Array1D<double> &errbars, Array1D<double> &basis,
                 Array1D<double> &alpha, Array2D<double> &Sig);
 
 
 
void BCS(Array2D<double> &PHI, Array1D<double> &y, double &sigma2,
                 double eta, Array1D<double> &lambda_init,
                 int adaptive, int optimal, double scale, int verbose,
                 Array1D<double> &weights, Array1D<int> &used,
                 Array1D<double> &errbars, Array1D<double> &basis,
                 Array1D<double> &alpha, double &lambda) ;
 
// void BCS(Array2D<double> &PHI, Array1D<double> &y, Array1D<double> &sigma2,
//                 double eta, Array1D<double> &lambda_init,
//                 int adaptive, int optimal, double scale, int verbose,
//                 Array1D<double> &weights, Array1D<int> &used,
//                 Array1D<double> &errbars, Array1D<double> &basis,
//                 Array1D<double> &alpha, Array1D<double> &lambda) ;
 
 
 
#endif // BCS_H