gauss_hermite_qr.hpp

Go to the documentation of this file.

#pragma once

#include <tuple>
#include <vector>

#include "../types.hpp"

#include "tables_gausshermite.hpp"


namespace waveblocks {
    namespace innerproducts {
        template <dim_t ORDER>
        struct GaussHermiteQR
        {
            static const dim_t D = 1;
            static const dim_t order = ORDER;

            using NodeMatrix = Eigen::Matrix<real_t,1,Eigen::Dynamic>;
            using WeightVector = Eigen::Matrix<real_t,1,Eigen::Dynamic>;

            static dim_t number_nodes()
            {
                return order;
            }

            static NodeMatrix nodes()
            {
                const auto& nodes = gauss_hermite_rules[ORDER-1].nodes;
                return NodeMatrix::Map(nodes.data(), nodes.size());
            }

            static WeightVector weights()
            {
                const auto& weights = gauss_hermite_rules[ORDER-1].weights;
                return WeightVector::Map(weights.data(), weights.size());
            }

            static std::tuple<NodeMatrix,WeightVector> nodes_and_weights()
            {
                return std::make_tuple(nodes(), weights());
            }

            static void clear_cache() {}
        };
    }
}