CshapeFunction_FEM_Q4.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 CSHAPEFUNCTION_FEM_Q4
#define CSHAPEFUNCTION_FEM_Q4

#include "CshapeFunction_FEM.h"

template<typename T,typename Timg> 
class CshapeFunction_FEM_Q4 : public CshapeFunction_FEM<T,Timg>
{
  
  public:

    CshapeFunction_FEM_Q4() : CshapeFunction_FEM<T,Timg>()
    {//constructor

      this->class_name = "Shape function : Finite Element Method 4 nodes";
      
    }
    
    virtual void init(CParameterNetCDF &fp)
    {
      this->m_dof = 8;
    }
 
    virtual void exec(const CImg<Timg> &box, const std::vector<T> &x, CImgList<T> &N)
    {// assign CImgList<T> N[dimensions][modes] for required Q4P1 modes at coordonate (x,y,z)
      
      float ratio = box.height()*box.width();

      cimglist_for(N, dim)
      {// trigonometric orientation for nodes
        N(dim,0) = (box.width()-x[0])*(box.height()-x[1])/ratio;//node 0
        N(dim,1) = x[0]*(box.height()-x[1])/ratio;//node 1
        N(dim,3) = x[0]*x[1]/ratio;//node 2
        N(dim,2) = (box.width()-x[0])*x[1]/ratio;//node 3
      }
      
    }
       
};

#endif