QuadratureRule.h

Go to the documentation of this file.

/*------------------------------------------------------------------------*/
/*  Copyright 2014 Sandia Corporation.                                    */
/*  This software is released under the license detailed                  */
/*  in the file, LICENSE, which is located in the top-level Nalu          */
/*  directory structure                                                   */
/*------------------------------------------------------------------------*/
#ifndef QuadratureRule_h
#define QuadratureRule_h

#include <vector>
#include <Teuchos_SerialDenseVector.hpp>

namespace sierra{
namespace nalu{

  // <abscissae, weights>
  std::pair<std::vector<double>, std::vector<double>>
  gauss_legendre_rule(int order);

  // <abscissae, weights>
  std::pair<std::vector<double>, std::vector<double>>
  gauss_lobatto_legendre_rule(int order, double xleft = -1.0, double xright = +1.0);

  // <abscissae, weights>
  std::pair<std::vector<double>, std::vector<double>>
  SGL_quadrature_rule(int order, std::vector<double> scsEndLocations);

  // a vector with -1 added at the first entry and +1 added at the last entry
  std::vector<double> pad_end_points(std::vector<double> x, double xleft = -1.0, double xright = +1.0);

} // namespace nalu
} // namespace Sierra

#endif