baseElem.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2026 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#ifndef FE_BASEELEM_H
#define FE_BASEELEM_H
 
#include "quadData.h"   // definition of QuadData
#include "util/point.h" // definition of Point
 
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_quads.size(); }
 
  virtual double elemSize(const std::vector<util::Point> &nodes) = 0;
 
  virtual std::vector<double>
  getShapes(const util::Point &p, const std::vector<util::Point> &nodes);
 
  virtual std::vector<std::vector<double>>
  getDerShapes(const util::Point &p,
               const std::vector<util::Point> &nodes);
 
  std::vector<fe::QuadData>
  getQuadDatas(const std::vector<util::Point> &nodes);
 
  std::vector<fe::QuadData>
  getQuadPoints(const std::vector<util::Point> &nodes);
 
protected:
  virtual std::vector<double> getShapes(const util::Point &p) = 0;
 
  virtual std::vector<std::vector<double>>
  getDerShapes(const util::Point &p) = 0;
 
  virtual util::Point mapPointToRefElem(const util::Point &p,
                    const std::vector<util::Point> &nodes);
 
  virtual double getJacobian(const util::Point &p,
                             const std::vector<util::Point> &nodes,
                             std::vector<std::vector<double>> *J) = 0;
 
  virtual void init() = 0;
 
  size_t d_quadOrder;
 
  size_t d_elemType;
 
  std::vector<fe::QuadData> d_quads;
};
 
} // namespace fe
 
#endif // FE_BASEELEM_H