Cmesh_global.h

Go to the documentation of this file.

/**********************************************************************
 * Copyright (C) 2012, The YaDICs Project Developers. 
 * See the COPYRIGHT file at the top-level directory of this distribution ./COPYRIGHT.
 * See ./COPYING file for copying and redistribution conditions.
 * 
 * This file is part of YaDICs.
 * YaDICs is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * YaDICs is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with YaDICs.  If not, see <http://www.gnu.org/licenses/>.
 * 
 * Information about how to use the software are provided at http://yadic.univ-lille1.fr/
 **********************************************************************/


#ifndef CMESH_GLOBAL
#define CMESH_GLOBAL

#include "Cmesh_generate.h"


template<typename T,typename Timg> 
class Cmesh_global : public Cmesh_generate<T,Timg>
{
 
  public:
    
    
  
  Cmesh_global() : Cmesh_generate<T,Timg>()
  {// constructor
    
    this->class_name = "Cmesh_global";

  }

  virtual int assign(const Cimage<Timg> &oImage)
  {// initialize attribute dimensions
    
    for (int i=0; i<this->num_var(); i++){this->m_grid.push_back(oImage.cur_dim(i));}
    
    this->m_node.assign(this->num_var(), 1,1,1,1, 0);
    
    this->m_connect.assign(1, this->m_dof);
    if (this->m_verbose){printf("\tConnectivity table : [%i](%i)\n",this->m_connect.size(),this->m_connect[0].width());}

    this->m_win.assign(3,0);
    this->m_pitch.assign(3,0);
    for (int i=0; i<3; i++){this->m_win[i] = oImage.cur_dim(i); this->m_pitch[i] = oImage.cur_dim(i);}
    
    this->m_numPix = oImage.m_curImg[0].size();
    
    return 0;

  }

  virtual void set_nodes(const Cimage<Timg> &oImage, CImgList<T> &nodes)
  {// node table definition
  
    cimglist_for(nodes, dim)
    {
      nodes[dim][0] = ((float)oImage.cur_dim(dim))/2.;
    }
    
  }
  
  virtual void set_connectivity(const Cimage<Timg> &oImage, const CImgList<T> &nodes, CImgList<int> &elements)
  {// element and connectivity table definition
    
    cimg_forX(elements[0],nd)
    {
      elements[0][nd] = nd;
    }
    
  }
  
  virtual int grid_dims(const int &i) const
  {
  
    switch(i)
    {
      case 0: return this->m_grid[0]; break;
      case 1: return this->m_grid[1]; break;
      case 2: 
        if (this->_3D)
        {return this->m_grid[2];}
        else
        {1;}
        break;
      case 3: return 1; break;
      default:
      {
        std::cerr<<this->class_name<<"::"<<__func__<<": error: id="<<i<<" which is over [0-3]\n"<<std::flush;
        return 0;
      }
    }
  
  }
 
 
};

#endif