EvaluationLaplacianI.hpp

Go to the documentation of this file.

//@HEADER
/*
*******************************************************************************

   Copyright (C) 2004, 2005, 2007 EPFL, Politecnico di Milano, INRIA
   Copyright (C) 2010 EPFL, Politecnico di Milano, Emory UNiversity

   This file is part of the LifeV library

   LifeV is free software; you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published by
   the Free Software Foundation; either version 3 of the License, or
   (at your option) any later version.

   LifeV 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
   Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public
   License along with this library; if not, see <http://www.gnu.org/licenses/>


*******************************************************************************
*/
//@HEADER

/*!
 *   @file
     @brief This file contains the definition of the EvaluationLaplacianI class.

     @date 11/2014
     @author Davide Forti <davide.forti@epfl.ch>
 */

#ifndef EVALUATION_LAPLACIAN_I_HPP
#define EVALUATION_LAPLACIAN_I_HPP

#include <lifev/core/LifeV.hpp>

#include <lifev/eta/fem/ETCurrentFE.hpp>
#include <lifev/eta/fem/ETCurrentFlag.hpp>
#include <lifev/core/fem/QuadratureRule.hpp>

#include <lifev/core/array/VectorSmall.hpp>

#include <lifev/eta/expression/ExpressionLaplacianI.hpp>


namespace LifeV
{

namespace ExpressionAssembly
{

//! Evaluation of the laplacian of the basis function phi_i in the case of a vectorial FE.
/*!
  @author Davide Forti <davide.forti@epfl.ch>

  This class aims at representing the laplacian of a test function in the assembly.

  This class is an Evaluation class, and therefore, has all the methods
  required to work within the Evaluation trees.
 */
template <UInt fieldDim, UInt spaceDim>
class EvaluationLaplacianI
{
public:

    //! @name Public Types
    //@{

    //! Type of the values returned by this class
    typedef VectorSmall<fieldDim> return_Type;

    //@}


    //! @name Static constants
    //@{

    //! Flag for the global current FE
    const static flag_Type S_globalUpdateFlag;

    //! Flag for the test current FE
    const static flag_Type S_testUpdateFlag;

    //! Flag for the solution current FE
    const static flag_Type S_solutionUpdateFlag;

    //@}


    //! @name Constructors, destructor
    //@{

    //! Empty constructor
    EvaluationLaplacianI() {}

    //! Copy constructor
    EvaluationLaplacianI (const EvaluationLaplacianI& provider) : M_valuesPtr (provider.M_valuesPtr) {}

    //! Expression-based constructor
    explicit EvaluationLaplacianI (const ExpressionLaplacianI& /*expression*/) {}

    //! Destructor
    ~EvaluationLaplacianI() {}

    //@}


    //! @name Methods
    //@{

    //! Do nothing internal update
    void update (const UInt& /*iElement*/) {}

    //! Display method
    static void display (std::ostream& out = std::cout)
    {
        out << "laplacian(phi_i[3])[" << spaceDim << "]";
    }

    //@}


    //! @name Set Methods
    //@{

    //! Do nothing setter for the global current FE
    template< typename CFEType >
    void setGlobalCFE (const CFEType* /*globalCFE*/) {}

    //! Setter for the test current FE
    template< typename CFEType >
    void setTestCFE (const CFEType* testCFE)
    {
        ASSERT (testCFE != 0, "Nul pointer to the testCFE cannot be set");
        M_valuesPtr = & (testCFE->M_laplacian);
    }

    //! Do nothing setter for the solution current FE
    template< typename CFEType >
    void setSolutionCFE (const CFEType* /*solutionCFE*/) {}

    //! Do nothing setter for the quadrature rule
    void setQuadrature (const QuadratureRule&) {}

    //@}


    //! @name Get Methods
    //@{

    //! Getter for the value for a vector
    const return_Type& value_qi (const UInt& q, const UInt& i) const
    {
        ASSERT ( q < M_valuesPtr->size(), "Quadrature point index invalid");
        ASSERT ( i < (*M_valuesPtr) [q].size(), "Dof index invalid");
        return (*M_valuesPtr) [q][i];
    }

    //! Getter for the value for a matrix
    const return_Type& value_qij (const UInt& q, const UInt& i, const UInt& /*j*/) const
    {
        ASSERT ( q < M_valuesPtr->size(), "Quadrature point index invalid");
        ASSERT ( i < (*M_valuesPtr) [q].size(), "Dof index invalid");
        return (*M_valuesPtr) [q][i];
    }

    //@}

private:

    //! Storage of the pointer to the data
    std::vector< std::vector < return_Type > > const* M_valuesPtr;

};


template<UInt fieldDim, UInt spaceDim>
const flag_Type EvaluationLaplacianI<fieldDim, spaceDim>::S_globalUpdateFlag = ET_UPDATE_NONE;

template<UInt fieldDim, UInt spaceDim>
const flag_Type EvaluationLaplacianI<fieldDim, spaceDim>::S_testUpdateFlag = ET_UPDATE_LAPLACIAN;

template<UInt fieldDim, UInt spaceDim>
const flag_Type EvaluationLaplacianI<fieldDim, spaceDim>::S_solutionUpdateFlag = ET_UPDATE_NONE;


//! Evaluation of the basis function dphi_i in the case of a scalar FE.
/*!
  @author Samuel Quinodoz <samuel.quinodoz@epfl.ch>

  This class aims at representing the gradient of a test function in the assembly.

  This class is an Evaluation class, and therefore, has all the methods
  required to work within the Evaluation trees.
 */
template <UInt spaceDim>
class EvaluationLaplacianI<1, spaceDim>
{
public:

    //! @name Public Types
    //@{

    //! Type of the values returned by this class
    typedef Real return_Type;

    //@}


    //! @name Static constants
    //@{

    //! Flag for the global current FE
    const static flag_Type S_globalUpdateFlag;

    //! Flag for the test current FE
    const static flag_Type S_testUpdateFlag;

    //! Flag for the solution current FE
    const static flag_Type S_solutionUpdateFlag;

    //@}


    //! @name Constructors, destructor
    //@{

    //! Empty constructor
    EvaluationLaplacianI() {}

    //! Copy constructor
    EvaluationLaplacianI (const EvaluationLaplacianI& provider) : M_valuesPtr (provider.M_valuesPtr) {}

    //! Expression-based constructor
    explicit EvaluationLaplacianI (const ExpressionLaplacianI& /*expression*/) {}

    //! Destructor
    ~EvaluationLaplacianI() {}

    //@}


    //! @name Methods
    //@{

    //! Do nothing internal update
    void update (const UInt& /*iElement*/) {}

    //! Display method
    static void display (std::ostream& out = std::cout)
    {
        out << "dphi_i[" << spaceDim << "]";
    }

    //@}


    //! @name Set Methods
    //@{

    //! Do nothing setter for the global current FE
    template< typename CFEType >
    void setGlobalCFE (const CFEType* /*globalCFE*/) {}

    //! Setter for the test current FE
    template< typename CFEType >
    void setTestCFE (const CFEType* testCFE)
    {
        ASSERT (testCFE != 0, "Nul pointer to the testCFE cannot be set");
        M_valuesPtr = & (testCFE->M_laplacian);
    }

    //! Do nothing setter for the solution current FE
    template< typename CFEType >
    void setSolutionCFE (const CFEType* /*solutionCFE*/) {}

    //! Do nothing setter for the quadrature rule
    void setQuadrature (const QuadratureRule&) {}

    //@}


    //! @name Get Methods
    //@{

    //! Getter for the value for a vector
    const return_Type& value_qi (const UInt& q, const UInt& i) const
    {
        ASSERT ( q < M_valuesPtr->size(), "Quadrature point index invalid");
        ASSERT ( i < (*M_valuesPtr) [q].size(), "Dof index invalid");
        return (*M_valuesPtr) [q][i];
    }

    //! Getter for the value for a matrix
    const return_Type& value_qij (const UInt& q, const UInt& i, const UInt& /*j*/) const
    {
        ASSERT ( q < M_valuesPtr->size(), "Quadrature point index invalid");
        ASSERT ( i < (*M_valuesPtr) [q].size(), "Dof index invalid");
        return (*M_valuesPtr) [q][i];
    }

    //@}

private:

    //! Storage of the pointer to the data
    std::vector< std::vector < return_Type > > const* M_valuesPtr;

};


template<UInt spaceDim>
const flag_Type EvaluationLaplacianI<1, spaceDim>::S_globalUpdateFlag = ET_UPDATE_NONE;

template<UInt spaceDim>
const flag_Type EvaluationLaplacianI<1, spaceDim>::S_testUpdateFlag = ET_UPDATE_LAPLACIAN;

template<UInt spaceDim>
const flag_Type EvaluationLaplacianI<1, spaceDim>::S_solutionUpdateFlag = ET_UPDATE_NONE;

} // Namespace ExpressionAssembly

} // Namespace LifeV
#endif