tetElem.h

//  Copyright (c) 2019 Prashant K. Jha
//  Copyright (c) 2019 Patrick Diehl
//
//  Distributed under the Boost Software License, Version 1.0. (See accompanying
//  file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
 
#ifndef FE_TETELEM_H
#define FE_TETELEM_H
 
#include "baseElem.h" // base class BaseElem
 
namespace fe {
 
class TetElem : public BaseElem {
 
public:
  explicit TetElem(size_t order);
 
  double elemSize(const std::vector<util::Point3> &nodes) override;
 
  std::vector<double>
  getShapes(const util::Point3 &p,
            const std::vector<util::Point3> &nodes) override;
 
  std::vector<std::vector<double>>
  getDerShapes(const util::Point3 &p,
               const std::vector<util::Point3> &nodes) override;
 
  std::vector<fe::QuadData>
  getQuadDatas(const std::vector<util::Point3> &nodes) override;
 
  std::vector<fe::QuadData>
  getQuadPoints(const std::vector<util::Point3> &nodes) override;
 
  std::string printStr(int nt = 0, int lvl = 0) const;
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); };
 
private:
  std::vector<double> getShapes(const util::Point3 &p) override;
 
  std::vector<std::vector<double>> getDerShapes(const util::Point3 &p) override;
 
  util::Point3
  mapPointToRefElem(const util::Point3 &p,
                    const std::vector<util::Point3> &nodes) override;
 
  double getJacobian(const util::Point3 &p,
                     const std::vector<util::Point3> &nodes,
                     std::vector<std::vector<double>> *J) override;
 
  void init() override;
};
 
} // namespace fe
 
#endif // FE_TETELEM_H