quadData.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2024 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#ifndef FE_QUADDATA_H
#define FE_QUADDATA_H
 
#include <util/io.h>
#include <util/point.h> // definition of Point
#include <vector>
 
namespace fe {
 
struct QuadData {
 
  double d_w;
 
  util::Point d_p;
 
  std::vector<double> d_shapes;
 
  std::vector<std::vector<double>> d_derShapes;
 
  std::vector<std::vector<double>> d_J;
 
  double d_detJ;
 
  QuadData() : d_w(0.), d_p(util::Point()), d_detJ(0.){};
 
  std::string printStr(int nt = 0, int lvl = 0) const {
 
    auto tabS = util::io::getTabS(nt);
    std::ostringstream oss;
    oss << tabS << "------- QuadData --------" << std::endl << std::endl;
    oss << tabS << "Weight = " << d_w << std::endl;
    oss << tabS << "Point = " << d_p.printStr() << std::endl;
    oss << tabS << "Shapes = " << util::io::printStr(d_shapes, 0) << std::endl;
    oss << tabS << "Derivative = " << util::io::printStr(d_derShapes, 0)
        << std::endl;
    oss << tabS << "Jacobian = " << util::io::printStr(d_J, 0)
        << std::endl;
    oss << tabS << "Det(J) = " << d_detJ << std::endl;
    oss << std::endl;
 
    return oss.str();
  };
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); };
};
 
} // namespace fe
 
#endif // FE_QUADDATA_H