AnaSol.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ANA_SOL_H
#define _ANA_SOL_H
 
#include <iostream>
#include <vector>
 
class RealFunc;
class VecFunc;
class TensorFunc;
class STensorFunc;
namespace tinyxml2 { class XMLElement; }
 
 
class AnaSol
{
public:
  explicit AnaSol(RealFunc* s1 = nullptr, VecFunc* s2 = nullptr,
                  VecFunc* v1 = nullptr, TensorFunc* v2 = nullptr,
                  STensorFunc* v3 = nullptr);
 
  AnaSol(RealFunc* s, STensorFunc* sigma);
 
  AnaSol(VecFunc* v, STensorFunc* sigma)
    : vecSol(v), vecSecSol(nullptr), stressSol(sigma) {}
 
  explicit AnaSol(STensorFunc* sigma)
    : vecSol(nullptr), vecSecSol(nullptr), stressSol(sigma) {}
 
  AnaSol(std::istream& is, const int nlines, bool scalarSol = true);
 
  explicit AnaSol(const tinyxml2::XMLElement* elem, bool scalarSol = true);
 
  AnaSol(const AnaSol&) = delete;
 
  virtual ~AnaSol();
 
  char hasScalarSol() const
  {
    if (stressSol && !vecSecSol && !vecSol)
      return 3;
    else if (!scalSecSol.empty())
      return 2;
    else if (!scalSol.empty())
      return 1;
    else
      return 0;
  }
 
  char hasVectorSol() const
  {
    if (stressSol)
      return 3;
    else if (vecSecSol)
      return 2;
    else if (vecSol)
      return 1;
    else
      return 0;
  }
 
  RealFunc* getScalarSol(size_t idx = 0) const
  { return scalSol.size() <= idx ? nullptr : scalSol[idx]; }
 
  VecFunc* getScalarSecSol(size_t idx = 0) const
  { return scalSecSol.size() <= idx ? nullptr : scalSecSol[idx]; }
 
  VecFunc* getVectorSol() const { return vecSol; }
 
  TensorFunc* getVectorSecSol() const { return vecSecSol; }
 
  STensorFunc* getStressSol() const { return stressSol; }
 
  void initPatch(size_t pIdx);
 
  virtual void setupSecondarySolutions();
 
private:
  template<class Scalar>
  void parseExpressionFunctions(const tinyxml2::XMLElement* elem, bool scalarSol);
 
  void parseFieldFunctions(const tinyxml2::XMLElement* elem, bool scalarSol);
 
protected:
  bool symmetric = false; 
 
  std::vector<RealFunc*>   scalSol; 
  std::vector<VecFunc*> scalSecSol; 
 
  VecFunc*     vecSol;    
  TensorFunc*  vecSecSol; 
  STensorFunc* stressSol; 
 
};
 
#endif