quadData.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_QUADDATA_H
#define FE_QUADDATA_H
 
#include <util/point.h>           // definition of Point3
#include <util/utilIO.h>
#include <vector>
 
namespace fe {
 
struct QuadData {
 
  double d_w;
 
  util::Point3 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::Point3()), 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