initialCondition.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 LOADING_INITIALCONDITION_H
#define LOADING_INITIALCONDITION_H
 
#include "util/point.h"          // definition of Point3
#include <string>
#include <vector>
 
// forward declaration of initial condition deck
namespace inp {
struct InitialConditionDeck;
}
 
// forward declaration
namespace fe {
class Mesh;
}
 
namespace loading {
 
class InitialCondition {
 
public:
  explicit InitialCondition(inp::InitialConditionDeck *deck);
 
  void apply(std::vector<util::Point3> *u, std::vector<util::Point3> *v,
             fe::Mesh *mesh);
 
  std::string printStr(int nt = 0, int lvl = 0) const;
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); };
 
private:
  double getICFormula(const std::string &fn_type,
                      const std::vector<double> &params, const util::Point3 &x,
                      const size_t &dof, const size_t &dim);
 
  inp::InitialConditionDeck *d_deck_p;
};
 
} // namespace loading
 
#endif // LOADING_INITIALCONDITION_H