moments.h

// Copyright (c) 2020 Chris Richardson & Matthew Scroggs
// FEniCS Project
// SPDX-License-Identifier:    MIT
 
#pragma once
 
#include "cell.h"
#include <utility>
#include <vector>
#include <xtensor/xtensor.hpp>
 
namespace basix
{
 
class FiniteElement;
 
namespace moments
{
 
std::pair<std::vector<xt::xtensor<double, 2>>,
          std::vector<xt::xtensor<double, 3>>>
make_integral_moments(const FiniteElement& moment_space, cell::type celltype,
                      std::size_t value_size, int q_deg);
 
std::pair<std::vector<xt::xtensor<double, 2>>,
          std::vector<xt::xtensor<double, 3>>>
make_dot_integral_moments(const FiniteElement& V, cell::type celltype,
                          std::size_t value_size, int q_deg);
 
std::pair<std::vector<xt::xtensor<double, 2>>,
          std::vector<xt::xtensor<double, 3>>>
make_tangent_integral_moments(const FiniteElement& V, cell::type celltype,
                              std::size_t value_size, int q_deg);
 
std::pair<std::vector<xt::xtensor<double, 2>>,
          std::vector<xt::xtensor<double, 3>>>
make_normal_integral_moments(const FiniteElement& V, cell::type celltype,
                             std::size_t value_size, int q_deg);
 
} // namespace moments
} // namespace basix