AnalyticalExprPython.hpp

Go to the documentation of this file.

/**
 * \file AnalyticalExprPython.hpp
 * \brief Analytical expression Python bridge
 */
 
#ifndef __ANALYTICAL_EXPR_PYTHON_HPP__
#define __ANALYTICAL_EXPR_PYTHON_HPP__
 
#ifdef ENABLE_PYTHON_BINDING
  #include <boost/python.hpp>
  #include <boost/python/def.hpp>
  #include <boost/python/numpy.hpp>
 
namespace bp = boost::python;
namespace np = boost::python::numpy;
 
namespace EshelbianPlasticity {
 
struct AnalyticalExprPython {
  AnalyticalExprPython() = default;
  virtual ~AnalyticalExprPython() = default;
 
  MoFEMErrorCode analyticalExprInit(const std::string py_file);
  MoFEMErrorCode evalAnalyticalDisp(double delta_t, double t, np::ndarray x,
                                    np::ndarray y, np::ndarray z,
                                    np::ndarray nx, np::ndarray ny,
                                    np::ndarray nz,
                                    const std::string &block_name,
                                    np::ndarray &analytical_expr);
  MoFEMErrorCode evalAnalyticalTraction(double delta_t, double t, np::ndarray x,
                                        np::ndarray y, np::ndarray z,
                                        np::ndarray nx, np::ndarray ny,
                                        np::ndarray nz,
                                        const std::string &block_name,
                                        np::ndarray &analytical_expr);
 
  MoFEMErrorCode evalAnalyticalExternalStrain(double delta_t, double t,
                                              np::ndarray x, np::ndarray y,
                                              np::ndarray z,
                                              const std::string &block_name,
                                              np::ndarray &analytical_expr);
 
  MoFEMErrorCode evalAnalyticalElastic(double delta_t, double t,
                                       np::ndarray x, np::ndarray y,
                                       np::ndarray z,
                                       const std::string &block_name,
                                       np::ndarray &analytical_expr);
 
  template <typename T>
  inline std::vector<T>
  py_list_to_std_vector(const boost::python::object &iterable) {
    return std::vector<T>(boost::python::stl_input_iterator<T>(iterable),
                          boost::python::stl_input_iterator<T>());
  }
 
private:
  bp::object mainNamespace;
  bp::object analyticalDispFun;
  bp::object analyticalTractionFun;
  bp::object analyticalExternalStrainFun;
  bp::object analyticalElasticFun;
};
 
extern boost::weak_ptr<AnalyticalExprPython> AnalyticalExprPythonWeakPtr;
 
} // namespace EshelbianPlasticity
#endif // ENABLE_PYTHON_BINDING
 
#endif // __ANALYTICAL_EXPR_PYTHON_HPP__