source: XIOS/trunk/src/node/interpolate_domain.cpp @ 1014

#include "interpolate_domain.hpp"
#include "type.hpp"

namespace xios {

  /// ////////////////////// Définitions ////////////////////// ///

  CInterpolateDomain::CInterpolateDomain(void)
    : CObjectTemplate<CInterpolateDomain>(), CInterpolateDomainAttributes(), CTransformation<CDomain>()
  { /* Ne rien faire de plus */ }

  CInterpolateDomain::CInterpolateDomain(const StdString & id)
    : CObjectTemplate<CInterpolateDomain>(id), CInterpolateDomainAttributes(), CTransformation<CDomain>()
  { /* Ne rien faire de plus */ }

  CInterpolateDomain::~CInterpolateDomain(void)
  {}

  CTransformation<CDomain>* CInterpolateDomain::create(const StdString& id, xml::CXMLNode* node)
  {
    CInterpolateDomain* interpDomain = CInterpolateDomainGroup::get("interpolate_domain_definition")->createChild(id);
    if (node) interpDomain->parse(*node);
    return static_cast<CTransformation<CDomain>*>(interpDomain);
  }

  bool CInterpolateDomain::_dummyRegistered = CInterpolateDomain::registerTrans();
  bool CInterpolateDomain::registerTrans()
  {
    registerTransformation(TRANS_INTERPOLATE_DOMAIN, create);
  }

  //----------------------------------------------------------------

  StdString CInterpolateDomain::GetName(void)    { return StdString("interpolate_domain"); }
  StdString CInterpolateDomain::GetDefName(void) { return StdString("interpolate_domain"); }
  ENodeType CInterpolateDomain::GetType(void)    { return eInterpolateDomain; }

  void CInterpolateDomain::checkValid(CDomain* domainSrc)
  {
    int order = 2;
    if (!this->order.isEmpty()) order = this->order.getValue();
    else this->order.setValue(order);
    if (order < 1)
    {
       ERROR("void CInterpolateDomain::checkValid(CDomain* domainSrc)",
             << "Interpolation order is less than 1, it should be greater than 0."
             << "Please define a correct one") ;
    }

    if (this->mode.isEmpty()) this->mode.setValue(mode_attr::compute);
    if (this->write_weight.isEmpty()) this->write_weight.setValue(false);

    StdString weightFile;
    switch (this->mode)
    {
      case mode_attr::read:
        if (this->weight_filename.isEmpty())
        {
          if (!this->write_weight)
            ERROR("void CInterpolateDomain::checkValid(CDomain* domainSrc)",
                 << "Read mode is activated but there is no file specified." << std::endl
                 << "Please define a correct file containing interpolation weights with option 'file'. ");
        }
        else
        {
          weightFile = this->weight_filename;
          ifstream f(weightFile.c_str());
          if (!f.good())
            ERROR("void CInterpolateDomain::checkValid(CDomain* domainSrc)",
                  << "Read mode is activated but file "  << weightFile << " doesn't exist." << std::endl
                  << "Please check this file ");
        }
        break;
      case mode_attr::compute:
        break;
      case mode_attr::read_or_compute:
        if (!this->weight_filename.isEmpty() && !this->write_weight)
        {
          weightFile = this->weight_filename;
          ifstream f(weightFile.c_str());
          if (!f.good())
            ERROR("void CInterpolateDomain::checkValid(CDomain* domainSrc)",
                  << "read_or_compute mode is activated but file "  << weightFile << " doesn't exist." << std::endl
                  << "Please check this file ");
        }
        break;
      default:
        break;
    }

  }

}