Map_headerTemplate.h File Reference

Variables
const char	cpp_file []
const char	hpp_header []
const char	hpp_endheader []
const char	header []
const char	evalOperatorBegin []
const char	evalOperatorEnd []
const char	evalOperatorDerBegin []
const char	evalOperatorDerEnd []
const char	evalODEBegin []
const char	evalODEEnd []
const char	evalODEDerBegin []
const char	evalODEDerEnd []

Variable Documentation

cpp_file

const char cpp_file[]

Initial value:

="\
#include \"%s.hpp\"\n\
\n\
template class %s::%s<double>;\n\
template class %s::%s<capd::interval>;\n\
\n\
\n"

evalODEBegin

const char evalODEBegin[]

Initial value:

="\
template<class Scalar>\n\
void CAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::computeODECoefficients(VectorType u[], unsigned order) const{\n\
  ScalarType* p = data;\n\
  int i;\n\
  for(typename VectorType::const_iterator b=u[0].begin();b!=u[0].end();++b,p+=timeJetSize)\n\
    *p=*b;\n\
\n\
  for(unsigned k=0;k<order;++k){\n\
"

evalODEDerBegin

const char evalODEDerBegin[]

Initial value:

="\
template<class Scalar>\n\
void CAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::computeODECoefficients(VectorType u[], MatrixType der[], unsigned order) const{\n\
  ScalarType* p = data;\n\
  typename VectorType::const_iterator b=u[0].begin();\n\
  typename MatrixType::const_iterator i=der[0].begin(), e = der[0].begin() + dimIn;\n\
  for(;b!=u[0].end();++b,p+=timeJetSize){\n\
    *p=*b;\n\
    ScalarType* q = p+maxOrder;\n\
    while(i!=e) { *q=*i; q+=maxOrder; ++i; }\n\
    e += dimIn;\n\
  }\n\
\n\
  for(unsigned k=0;k<order;++k){\n\
"

evalODEDerEnd

const char evalODEDerEnd[]

Initial value:

="\n\
    int j;\n\
    for(j=0, p=data+k+1;j<dimIn;++j, p+=timeJetSize){\n\
      u[k+1][j] = *p = data[pos[j]+k]/(k+1);\n\
      for(int c=0, r=maxOrder; c<dimIn; ++c, r+=maxOrder)\n\
        der[k+1](j+1,c+1) = *(p+r) = data[pos[j]+r+k]/(k+1);\n\
    }\n\
  }\n\
}\n\
"

evalODEEnd

const char evalODEEnd[]

Initial value:

="\n\
    for(i=0, p=data+k+1;i<dimIn;++i, p+=timeJetSize)\n\
      u[k+1][i] = *p = data[pos[i]+k]/(k+1);\n\
  }\n\
}\n\
"

evalOperatorBegin

const char evalOperatorBegin[]

Initial value:

="\
template<class Scalar>\n\
typename CAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::ImageVectorType\nCAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::operator()(const VectorType& u) const{\n\
  ScalarType* p = data;\n\
  for(typename VectorType::const_iterator b=u.begin();b!=u.end();++b,p+=timeJetSize)\n\
    *p=*b;\n\
  \n\
"

evalOperatorDerBegin

const char evalOperatorDerBegin[]

Initial value:

="\n\
template<class Scalar>\n\
typename CAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::ImageVectorType\nCAPD__TEMPLATE_MAP_CLASSNAME__ <Scalar>::operator()(const VectorType& u, MatrixType& der) const{\n\
  ScalarType* p = data;\n\
  int i,j;\n\
  for(i=0;i<u.dimension();++i,p+=timeJetSize){\n\
    *p=u[i];\n\
    ScalarType* q=p+maxOrder;\n\
    for(j=0;j<u.dimension();++j,q+=maxOrder)\n\
      *q = (i==j) ? capd::TypeTraits<ScalarType>::one() : capd::TypeTraits<ScalarType>::zero();\n\
  }\n\
  \n\
"

evalOperatorDerEnd

const char evalOperatorDerEnd[]

Initial value:

="\n\
  ImageVectorType result(this->imageDimension());\n\
  for(i=0;i<this->imageDimension();++i){\n\
    ScalarType* p = data+pos[i];\n\
    result[i] = *p;\n\
    p+=maxOrder;\n\
    for(j=0;j<this->dimension();++j,p+=maxOrder)\n\
      der(i+1,j+1) = *p;\n\
  }\n\
  return result;\n\
}\n\
\n\
"

evalOperatorEnd

const char evalOperatorEnd[]

Initial value:

="\n\
  ImageVectorType result;\n\
  for(int i=0;i<result.dimension();++i)\n\
    result[i] = data[pos[i]];\n\
  return result;\n\
}\n\
\n\
"

header

const char header[]

hpp_endheader

const char hpp_endheader[]

Initial value:

="\n\
} /* namespace */\n\
\n\
#undef CAPD__TEMPLATE_MAP_CLASSNAME__\n\
\n\
#endif\n\n\
"

hpp_header

const char hpp_header[]

Initial value:

="\
#ifndef __%s_%s_HPP__\n\
#define __%s_%s_HPP__\n\n\
#include \"%s.h\"\n\
#include \"capd/autodiff/eval.hpp\"\n\
#include \"capd/vectalg/Vector.hpp\"\n\
#include \"capd/vectalg/Matrix.hpp\"\n\
#include \"capd/vectalg/Norm.hpp\"\n\
#include \"capd/matrixAlgorithms/floatMatrixAlgorithms.hpp\"\n\
#include \"capd/map/Function.hpp\"\n\
#include \"capd/dynsys/TaylorHOE.hpp\"\n\
#include \"capd/poincare/TimeMap.hpp\"\n\
#include \"capd/poincare/PoincareMap.hpp\"\n\
\n\
#define CAPD__TEMPLATE_MAP_CLASSNAME__ %s\n\
\n\
namespace %s{\n\
\n"