src/MEDCoupling/MEDCouplingRemapper.hxx

   1 // Copyright (C) 2007-2012  CEA/DEN, EDF R&D
   2 //
   3 // This library is free software; you can redistribute it and/or
   4 // modify it under the terms of the GNU Lesser General Public
   5 // License as published by the Free Software Foundation; either
   6 // version 2.1 of the License.
   7 //
   8 // This library is distributed in the hope that it will be useful,
   9 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 // Lesser General Public License for more details.
  12 //
  13 // You should have received a copy of the GNU Lesser General Public
  14 // License along with this library; if not, write to the Free Software
  15 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 //
  17 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 //
  19 // Author : Anthony Geay (CEA/DEN)
  20
  21 #ifndef __PARAMEDMEM_MEDCOUPLINGREMAPPER_HXX__
  22 #define __PARAMEDMEM_MEDCOUPLINGREMAPPER_HXX__
  23
  24 #include "MEDCoupling.hxx"
  25 #include "MEDCouplingTimeLabel.hxx"
  26 #include "InterpolationOptions.hxx"
  27 #include "MEDCouplingNatureOfField.hxx"
  28 #include "InterpKernelException.hxx"
  29
  30 #include <map>
  31 #include <vector>
  32
  33 namespace ParaMEDMEM
  34 {
  35   class MEDCouplingMesh;
  36   class MEDCouplingUMesh;
  37   class MEDCouplingFieldDouble;
  38   class MEDCouplingFieldTemplate;
  39 }
  40
  41 namespace ParaMEDMEM
  42 {
  43   class MEDCouplingRemapper : public TimeLabel, public INTERP_KERNEL::InterpolationOptions
  44   {
  45   public:
  46     MEDCOUPLINGREMAPPER_EXPORT MEDCouplingRemapper();
  47     MEDCOUPLINGREMAPPER_EXPORT ~MEDCouplingRemapper();
  48     MEDCOUPLINGREMAPPER_EXPORT int prepare(const MEDCouplingMesh *srcMesh, const MEDCouplingMesh *targetMesh, const char *method) throw(INTERP_KERNEL::Exception);
  49     MEDCOUPLINGREMAPPER_EXPORT int prepareEx(const MEDCouplingFieldTemplate *src, const MEDCouplingFieldTemplate *target) throw(INTERP_KERNEL::Exception);
  50     MEDCOUPLINGREMAPPER_EXPORT void transfer(const MEDCouplingFieldDouble *srcField, MEDCouplingFieldDouble *targetField, double dftValue) throw(INTERP_KERNEL::Exception);
  51     MEDCOUPLINGREMAPPER_EXPORT void partialTransfer(const MEDCouplingFieldDouble *srcField, MEDCouplingFieldDouble *targetField) throw(INTERP_KERNEL::Exception);
  52     MEDCOUPLINGREMAPPER_EXPORT void reverseTransfer(MEDCouplingFieldDouble *srcField, const MEDCouplingFieldDouble *targetField, double dftValue) throw(INTERP_KERNEL::Exception);
  53     MEDCOUPLINGREMAPPER_EXPORT MEDCouplingFieldDouble *transferField(const MEDCouplingFieldDouble *srcField, double dftValue) throw(INTERP_KERNEL::Exception);
  54     MEDCOUPLINGREMAPPER_EXPORT MEDCouplingFieldDouble *reverseTransferField(const MEDCouplingFieldDouble *targetField, double dftValue) throw(INTERP_KERNEL::Exception);
  55     MEDCOUPLINGREMAPPER_EXPORT bool setOptionInt(const std::string& key, int value);
  56     MEDCOUPLINGREMAPPER_EXPORT bool setOptionDouble(const std::string& key, double value);
  57     MEDCOUPLINGREMAPPER_EXPORT bool setOptionString(const std::string& key, const std::string& value);
  58     //
  59     MEDCOUPLINGREMAPPER_EXPORT int nullifiedTinyCoeffInCrudeMatrixAbs(double maxValAbs) throw(INTERP_KERNEL::Exception);
  60     MEDCOUPLINGREMAPPER_EXPORT int nullifiedTinyCoeffInCrudeMatrix(double scaleFactor) throw(INTERP_KERNEL::Exception);
  61     MEDCOUPLINGREMAPPER_EXPORT double getMaxValueInCrudeMatrix() const throw(INTERP_KERNEL::Exception);
  62   public:
  63     MEDCOUPLINGREMAPPER_EXPORT const std::vector<std::map<int,double> >& getCrudeMatrix() const;
  64     MEDCOUPLINGREMAPPER_EXPORT static void PrintMatrix(const std::vector<std::map<int,double> >& m);
  65   private:
  66     int prepareUU(const char *method) throw(INTERP_KERNEL::Exception);
  67     int prepareEE(const char *method) throw(INTERP_KERNEL::Exception);
  68     int prepareUC(const char *method) throw(INTERP_KERNEL::Exception);
  69     int prepareCU(const char *method) throw(INTERP_KERNEL::Exception);
  70     int prepareCC(const char *method) throw(INTERP_KERNEL::Exception);
  71     void updateTime() const;
  72     void releaseData(bool matrixSuppression);
  73     void transferUnderground(const MEDCouplingFieldDouble *srcField, MEDCouplingFieldDouble *targetField, bool isDftVal, double dftValue) throw(INTERP_KERNEL::Exception);
  74     void computeDeno(NatureOfField nat, const MEDCouplingFieldDouble *srcField, const MEDCouplingFieldDouble *trgField);
  75     void computeDenoFromScratch(NatureOfField nat, const MEDCouplingFieldDouble *srcField, const MEDCouplingFieldDouble *trgField) throw(INTERP_KERNEL::Exception);
  76     void computeProduct(const double *inputPointer, int inputNbOfCompo, bool isDftVal, double dftValue, double *resPointer);
  77     void computeReverseProduct(const double *inputPointer, int inputNbOfCompo, double dftValue, double *resPointer);
  78     void buildFinalInterpolationMatrixByConvolution(const std::vector< std::map<int,double> >& m1D,
  79                                                     const std::vector< std::map<int,double> >& m2D,
  80                                                     const int *corrCellIdSrc, int nbOf2DCellsSrc, int nbOf1DCellsSrc,
  81                                                     const int *corrCellIdTrg);
  82     static void ReverseMatrix(const std::vector<std::map<int,double> >& matIn, int nbColsMatIn,
  83                               std::vector<std::map<int,double> >& matOut);
  84     static void ComputeRowSumAndColSum(const std::vector<std::map<int,double> >& matrixDeno,
  85                                        std::vector<std::map<int,double> >& deno, std::vector<std::map<int,double> >& denoReverse);
  86     static void ComputeColSumAndRowSum(const std::vector<std::map<int,double> >& matrixDeno,
  87                                        std::vector<std::map<int,double> >& deno, std::vector<std::map<int,double> >& denoReverse);
  88   private:
  89     MEDCouplingMesh *_src_mesh;
  90     MEDCouplingMesh *_target_mesh;
  91     std::string _src_method;
  92     std::string _target_method;
  93     NatureOfField _nature_of_deno;
  94     unsigned int _time_deno_update;
  95     std::vector<std::map<int,double> > _matrix;
  96     std::vector<std::map<int,double> > _deno_multiply;
  97     std::vector<std::map<int,double> > _deno_reverse_multiply;
  98   };
  99 }
 100
 101 #endif