baseElem.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_BASEELEM_H
#define FE_BASEELEM_H
 
#include "quadData.h"   // definition of QuadData
#include "util/point.h" // definition of Point3
 
namespace fe {
 
class BaseElem {
 
public:
  BaseElem(size_t order, size_t element_type);
 
  size_t getElemType() { return d_elemType; }
 
  size_t getQuadOrder() { return d_quadOrder; }
 
  size_t getNumQuadPoints() { return d_numQuadPts; }
 
  virtual double elemSize(const std::vector<util::Point3> &nodes) = 0;
 
  virtual std::vector<double>
  getShapes(const util::Point3 &p, const std::vector<util::Point3> &nodes);
 
  virtual std::vector<std::vector<double>>
  getDerShapes(const util::Point3 &p,
               const std::vector<util::Point3> &nodes);
 
  virtual std::vector<fe::QuadData>
  getQuadDatas(const std::vector<util::Point3> &nodes) = 0;
 
  virtual std::vector<fe::QuadData>
  getQuadPoints(const std::vector<util::Point3> &nodes) = 0;
 
  std::string printStr(int nt = 0, int lvl = 0) const;
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); };
 
protected:
  virtual std::vector<double> getShapes(const util::Point3 &p) = 0;
 
  virtual std::vector<std::vector<double>>
  getDerShapes(const util::Point3 &p) = 0;
 
  virtual util::Point3
  mapPointToRefElem(const util::Point3 &p,
                    const std::vector<util::Point3> &nodes);
 
  virtual double getJacobian(const util::Point3 &p,
                             const std::vector<util::Point3> &nodes,
                             std::vector<std::vector<double>> *J) = 0;
 
  virtual void init() = 0;
 
  size_t d_quadOrder;
 
  size_t d_numQuadPts;
 
  size_t d_elemType;
 
  std::vector<fe::QuadData> d_quads;
};
 
} // namespace fe
 
#endif // FE_BASEELEM_H