src/MEDCalc/cmp/MEDPresentationContour.cxx

   1 // Copyright (C) 2016-2021  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, or (at your option) any later version.
   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
  20 #include "MEDPresentationContour.hxx"
  21
  22 #include <SALOME_KernelServices.hxx>
  23 #undef LOG  // should be fixed in KERNEL - double definition
  24 #include <Basics_Utils.hxx>
  25
  26 #include <sstream>
  27
  28 const std::string MEDPresentationContour::TYPE_NAME = "MEDPresentationContour";
  29 const std::string MEDPresentationContour::PROP_NB_CONTOUR = "nbContour";
  30
  31 MEDPresentationContour::MEDPresentationContour(const MEDCALC::ContourParameters& params,
  32                                                const MEDCALC::ViewModeType viewMode) :
  33         MEDPresentation(params.fieldHandlerId, TYPE_NAME, viewMode, params.colorMap, params.scalarBarRange),
  34         _params(params)
  35 {
  36   setIntProperty(MEDPresentationContour::PROP_NB_CONTOUR, params.nbContours);
  37 }
  38
  39 void
  40 MEDPresentationContour::initFieldMeshInfos()
  41 {
  42   MEDPresentation::initFieldMeshInfos();
  43   _colorByType = "POINTS";
  44 }
  45
  46 void
  47 MEDPresentationContour::setNumberContours()
  48 {
  49   std::ostringstream oss;
  50
  51   oss << "min_max = " << _srcObjVar << ".PointData.GetArray('" << _fieldName << "').GetRange();";
  52   pushAndExecPyLine(oss.str()); oss.str("");
  53   oss << "delta = (min_max[1]-min_max[0])/float(" << _params.nbContours << ");";
  54   pushAndExecPyLine(oss.str()); oss.str("");
  55   oss << _objVar << ".Isosurfaces = [min_max[0]+0.5*delta+i*delta for i in range(" << _params.nbContours << ")];";
  56   pushAndExecPyLine(oss.str()); oss.str("");
  57 }
  58
  59 void
  60 MEDPresentationContour::internalGeneratePipeline()
  61 {
  62   MEDPresentation::internalGeneratePipeline();
  63
  64   MEDPyLockWrapper lock;
  65
  66   createSource();
  67   setTimestamp();
  68
  69   // Populate internal array of available components:
  70   fillAvailableFieldComponents();
  71   if (getIntProperty(MEDPresentation::PROP_NB_COMPONENTS) > 1)
  72     {
  73       const char * msg = "Contour presentation only works for scalar field!"; // this message will appear in GUI too
  74       STDLOG(msg);
  75       throw KERNEL::createSalomeException(msg);
  76     }
  77   if (_params.nbContours < 1)
  78     {
  79       const char * mes = "Invalid number of contours!";
  80       STDLOG(mes);
  81       throw KERNEL::createSalomeException(mes);
  82     }
  83
  84   setOrCreateRenderView(); // instantiate __viewXXX, needs to be after the exception above otherwise previous elements in the view will be hidden.
  85
  86   // Contour needs point data:
  87   applyCellToPointIfNeeded();
  88
  89   std::ostringstream oss;
  90   oss << _objVar << " = pvs.Contour(Input=" << _srcObjVar << ");";
  91   pushAndExecPyLine(oss.str()); oss.str("");
  92
  93   showObject();
  94
  95   oss << _objVar << ".ContourBy = ['POINTS', '" << _fieldName << "'];";
  96   pushAndExecPyLine(oss.str()); oss.str("");
  97
  98   // Colorize contour
  99   oss << _objVar << ".ComputeScalars = 1;";
 100   pushAndExecPyLine(oss.str()); oss.str("");
 101
 102   // Set number of contours
 103   setNumberContours();
 104
 105   colorBy();    // see initFieldInfo() - necessarily POINTS because of the conversion above
 106   showScalarBar();
 107   selectColorMap();
 108   rescaleTransferFunction();
 109   resetCameraAndRender();
 110 }
 111
 112 void
 113 MEDPresentationContour::updatePipeline(const MEDCALC::ContourParameters& params)
 114 {
 115   if (params.fieldHandlerId != _params.fieldHandlerId)
 116     throw KERNEL::createSalomeException("Unexpected updatePipeline error! Mismatching fieldHandlerId!");
 117
 118   if (params.scalarBarRange != _params.scalarBarRange)
 119     updateScalarBarRange<MEDPresentationContour, MEDCALC::ContourParameters>(params.scalarBarRange);
 120   if (params.colorMap != _params.colorMap)
 121     updateColorMap<MEDPresentationContour, MEDCALC::ContourParameters>(params.colorMap);
 122
 123   if (params.nbContours != _params.nbContours)
 124     {
 125       if (params.nbContours < 1)
 126         {
 127           const char * mes = "Invalid number of contours!";
 128           STDLOG(mes);
 129           throw KERNEL::createSalomeException(mes);
 130         }
 131       updateNbContours(params.nbContours);
 132     }
 133 }
 134
 135 void
 136 MEDPresentationContour::updateNbContours(const int nbContours)
 137 {
 138   _params.nbContours = nbContours;
 139
 140   // GUI helper:
 141   setIntProperty(MEDPresentationContour::PROP_NB_CONTOUR, nbContours);
 142
 143   // Update the pipeline:
 144   {
 145     MEDPyLockWrapper lock;
 146     setNumberContours();
 147     pushAndExecPyLine("pvs.Render();");
 148   }
 149 }