packages/vtk/src/bbvtkPiecewiseFunction.cxx

   1 /*=========================================================================
   2   Program:   bbtk
   3   Module:    $RCSfile: bbvtkPiecewiseFunction.cxx,v $
   4   Language:  C++
   5   Date:      $Date: 2009/05/15 14:58:01 $
   6   Version:   $Revision: 1.2 $
   7 =========================================================================*/
   8
   9 /* ---------------------------------------------------------------------
  10
  11 * Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale)
  12 * Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux
  13 *
  14 *  This software is governed by the CeCILL-B license under French law and
  15 *  abiding by the rules of distribution of free software. You can  use,
  16 *  modify and/ or redistribute the software under the terms of the CeCILL-B
  17 *  license as circulated by CEA, CNRS and INRIA at the following URL
  18 *  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  19 *  or in the file LICENSE.txt.
  20 *
  21 *  As a counterpart to the access to the source code and  rights to copy,
  22 *  modify and redistribute granted by the license, users are provided only
  23 *  with a limited warranty  and the software's author,  the holder of the
  24 *  economic rights,  and the successive licensors  have only  limited
  25 *  liability.
  26 *
  27 *  The fact that you are presently reading this means that you have had
  28 *  knowledge of the CeCILL-B license and that you accept its terms.
  29 * ------------------------------------------------------------------------ */
  30
  31 /**
  32  *  \file
  33  *  \brief
  34  */
  35
  36
  37
  38 #ifdef _USE_VTK_
  39
  40
  41 #include "bbvtkPiecewiseFunction.h"
  42 #include "bbvtkPackage.h"
  43
  44 namespace bbvtk
  45 {
  46
  47
  48    BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,PiecewiseFunction)
  49    BBTK_BLACK_BOX_IMPLEMENTATION(PiecewiseFunction,bbtk::AtomicBlackBox);
  50
  51
  52         //---------------------------------------------------------------------
  53    void PiecewiseFunction::bbUserSetDefaultValues()
  54    {
  55      // Identity function
  56      std::vector<float> x,y;
  57      x.push_back(0);
  58      x.push_back(1);
  59      y.push_back(0);
  60      y.push_back(1);
  61      bbSetInputX(x);
  62      bbSetInputY(y);
  63          bbSetOutputOut(0);
  64    }
  65
  66
  67         //---------------------------------------------------------------------
  68    void PiecewiseFunction::bbUserInitializeProcessing()
  69    {
  70      bbSetOutputOut(vtkPiecewiseFunction::New());
  71    }
  72
  73         //---------------------------------------------------------------------
  74    void PiecewiseFunction::bbUserFinalizeProcessing()
  75    {
  76            if (bbGetOutputOut()!=NULL)
  77            {
  78                    bbGetOutputOut()->Delete();
  79                    bbSetOutputOut(0);
  80            }
  81    }
  82 //---------------------------------------------------------------------
  83
  84    void PiecewiseFunction::Process()
  85    {
  86      if ( ( bbGetInputStatus("X") != bbtk::UPTODATE ) ||
  87           ( bbGetInputStatus("Y") != bbtk::UPTODATE ) )
  88        {
  89          if  ( bbGetInputX().size() != bbGetInputY().size() )
  90            bbtkError(bbGetFullName()
  91                      <<" : input vectors X and Y do not have the same size");
  92          if  ( bbGetInputX().size() < 2 )
  93            bbtkError(bbGetFullName()
  94                      <<" : input vectors X and Y have a size < 2");
  95          bbGetOutputOut()->RemoveAllPoints ();
  96          //      std::cout << bbGetInputX().size()<< ","<< bbGetInputY().size()<< std::endl;
  97          std::vector<float>::const_iterator x,y;
  98          std::vector<float> X = bbGetInputX();
  99          std::vector<float> Y = bbGetInputY();
 100          for (x = X.begin(), y = Y.begin();
 101               x != X.end(),  y != Y.end();
 102               ++x, ++y )
 103            {
 104              //              std::cout << *x << " -> " << *y << std::endl;
 105              bbGetOutputOut()->AddPoint(*x,*y);
 106            }
 107        }
 108    }
 109
 110 }//namespace bbtk
 111
 112 #endif // _USE_VTK_
 113