Cfield_intercor.h

/**********************************************************************
 * 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 CFIELD_INTERCOR
#define CFIELD_INTERCOR


//-----------------------------CIMG_NETCDF-------------------------------
#define cimg_plugin  "plugins/add_fileformat.h"
#define cimg_use_netcdf

#ifdef cimg_use_netcdf
 #include "../NetCDF.Tool/struct_parameter_NetCDF.h"
// #define     cimg_plugin1 "plugins/netcdf_file_format4CImg1.h"
// #define cimglist_plugin1 "plugins/netcdf_file_format4CImgList1.h"
 #define     cimg_plugin2 "plugins/netcdf_file_format4CImg2.h"
 #define cimglist_plugin2 "plugins/netcdf_file_format4CImgList2.h"
 #include "../CImg.Tool/CImg_NetCDF.h"
#endif
//-----------------------------CIMG_NETCDF-------------------------------

#include "Cfield_opticalFlow_fem.h"


template<typename T, typename Timg>
class Cfield_intercor: public Cfield_opticalFlow_fem<T,Timg>
{
  
  
  public:

    Cfield_intercor() : Cfield_opticalFlow_fem<T,Timg>()
    {// constructor
      
      this->class_name = "Field : Inter-correlation";
      
    }
    
    
    virtual void init(CParameterNetCDF &fp)
    {// init from parameter file

      std::string attribute_name = "store_correlogram";
      int error = fp.loadAttribute(attribute_name,this->m_storeCorrelo);
      if (this->m_verbose and this->m_storeCorrelo){printf("\t\t Correlograms are stored\n");}
      
    }
    
    
    virtual void exec(const Cimage<Timg> &oImage, const Cmesh<T,Timg> &oMesh, CImgList<T> &modes_prev, CshapeFunction<T,Timg> &oShape_prev, const std::vector<T> &ref_prev, const CImgList<T> &nodes_prev)
    {
      if (this->m_storeCorrelo){this->m_correlogram.assign(oImage.cur_dim(0), oImage.cur_dim(1), oImage.cur_dim(2), oImage.cur_dim(3), 0);}
      
      this->m_mode.assign(oMesh.num_var()+2, oMesh.grid_dims(0),oMesh.grid_dims(1),oMesh.grid_dims(2),1, 0);

      Cfield<T,Timg>::exec(oImage, oMesh, modes_prev, oShape_prev, ref_prev, nodes_prev); 
      
    }
    

      
    virtual void project(const Cmesh<T,Timg> &oMesh, const Cimage<Timg> &oImage, CImgList<T> &modes_prev, CshapeFunction<T,Timg> &oShape_prev, const std::vector<T> &ref_prev, const CImgList<T> &nodes_prev, CImgList<T> &modes_cur)
    {//project modes on new base
      
      
      if ((oShape_prev.m_name).compare(0,5,"OFFEM") == 0)
      {//basic linear interpolation from previous the current grid of modes [U, V, W]

        
        CImgList<T> U(oMesh.num_var(), modes_prev[0].width(), modes_prev[0].height(), modes_prev[0].depth(), 1, 0);

        this->Field2Field(modes_prev, U);
        
        cimglist_for(U, mode)
        {modes_cur[mode+2] = U[mode];}
        
        
      }
      else if((oShape_prev.m_name).compare(0,2,"IC") == 0)
      {//basic linear interpolation from previous the current grid of modes [U, V, W]
      
        this->Field2Field(modes_prev, modes_cur);
     
      }
      else if((oShape_prev.m_name).compare(0,4,"OFIB") == 0)
      {//basic linear interpolation from previous the current grid and building of U, V and W at previous element centroid
        
        CImgList<T> U0(oMesh.num_var(), modes_prev[0].width(), modes_prev[0].height(), modes_prev[0].depth(), 1, 0);
        CImgList<T> U1(oMesh.num_var(), modes_cur[0].width(), modes_cur[0].height(), modes_cur[0].depth(), 1, 0);
        
        this->OFIBU(oMesh, oImage, oShape_prev, ref_prev, nodes_prev, modes_prev, U0);
        
        this->Field2Field(U0, U1);
        
        cimglist_for(U1, mode)
        {modes_cur[mode+2] = U1[mode];}

        
      }
      else if((oShape_prev.m_name).compare(0,3,"OFI") == 0)
      {
        
        CImgList<T> U(oMesh.num_var(), modes_cur[0].width(), modes_cur[0].height(), modes_cur[0].depth(), 1, 0);
        
        this->OFIU(oMesh, oImage, oShape_prev, ref_prev, nodes_prev, modes_prev, U);
        
        cimglist_for(U, mode)
        {modes_cur[mode+2] = U[mode];}
        
      }
      else{printf("unknown name for previous correlation method");}
      
    }
    
  
    
};

#endif