/*! \file scalar_algorithm_extract_scalar.cpp \author Ha NGUYEN \since 23 June 2016 \date 23 June 2016 \brief Algorithm for extract an axis to a scalar */ #include "scalar_algorithm_extract_axis.hpp" #include "axis.hpp" #include "scalar.hpp" #include "extract_axis_to_scalar.hpp" #include "grid.hpp" #include "grid_transformation_factory_impl.hpp" #include "reduction.hpp" namespace xios { CGenericAlgorithmTransformation* CScalarAlgorithmExtractAxis::create(bool isSource, CGrid* gridDst, CGrid* gridSrc, CTransformation* transformation, int elementPositionInGrid, std::map& elementPositionInGridSrc2ScalarPosition, std::map& elementPositionInGridSrc2AxisPosition, std::map& elementPositionInGridSrc2DomainPosition, std::map& elementPositionInGridDst2ScalarPosition, std::map& elementPositionInGridDst2AxisPosition, std::map& elementPositionInGridDst2DomainPosition) TRY { std::vector scalarListDestP = gridDst->getScalars(); std::vector axisListSrcP = gridSrc->getAxis(); CExtractAxisToScalar* extractAxis = dynamic_cast (transformation); int scalarDstIndex = elementPositionInGridDst2ScalarPosition[elementPositionInGrid]; int axisSrcIndex = elementPositionInGridSrc2AxisPosition[elementPositionInGrid]; return (new CScalarAlgorithmExtractAxis(isSource, scalarListDestP[scalarDstIndex], axisListSrcP[axisSrcIndex], extractAxis)); } CATCH bool CScalarAlgorithmExtractAxis::dummyRegistered_ = CScalarAlgorithmExtractAxis::registerTrans(); bool CScalarAlgorithmExtractAxis::registerTrans() TRY { return CGridTransformationFactory::registerTransformation(TRANS_EXTRACT_AXIS_TO_SCALAR, create); } CATCH CScalarAlgorithmExtractAxis::CScalarAlgorithmExtractAxis(bool isSource, CScalar* scalarDestination, CAxis* axisSource, CExtractAxisToScalar* algo) : CAlgorithmTransformationTransfer(isSource) TRY { algo->checkValid(scalarDestination, axisSource); pos_ = algo->position; this->transformationMapping_[0]=pos_ ; scalarDestination->checkAttributes() ; this->computeAlgorithm(axisSource->getLocalView(CElementView::WORKFLOW), scalarDestination->getLocalView(CElementView::WORKFLOW)) ; } CATCH }