multiindex.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 MULTIINDEX_H
#define MULTIINDEX_H
 
#include "Array1D.h"
#include "Array2D.h"
 
 
extern "C" void heap_ext_(const int *,const int *, int *, int *, int *);
 
int computeNPCTerms(int ndim,int norder);
 
int computeMultiIndex(int ndim,int norder, Array2D<int> &mi);
 
int computeMultiIndexT(int ndim,int norder, int *mi);
 
int computeMultiIndex(int ndim,int norder, Array2D<int> &mi, string ordtype);
 
 
int computeMultiIndexTP(Array1D<int>& maxorders, Array2D<int>& mindex);
 
int computeNPCTermsHDMR(int ndim, Array1D<int>& maxorders);
 
int computeMultiIndexHDMR(int ndim, Array1D<int>& maxorders,Array2D<int>& mindex);
 
void decodeMindex(Array1D< Array2D<int> >& sp_mindex, int ndim, Array2D<int>& mindex);
 
 
void upOrder(Array2D<int>& mindex,Array2D<int>& new_mindex);
 
bool is_admis(Array1D<int>& mindex_try,Array2D<int>& mindex);
 
void getOrders(Array2D<int>& mindex,Array1D<int>& orders);
 
int get_invmindex(Array1D<int> mi);
 
int get_invmindex_ord(Array1D<int> mi);
 
 
#endif // MULTIINDEX_H