core/testsuite/fe_function/user_fun.cpp

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 LifeV.

    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 LifeV.  If not, see <http://www.gnu.org/licenses/>.

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

/**
   @file
   @author A. Fumagalli <alessio.fumagalli@mail.polimi.it>
   @date 2010-07-29
*/

#include "user_fun.hpp"

namespace dataProblem
{

Real MyFun::eval ( const UInt& iElem, const point_Type& P, const Real& time ) const
{
    const Real scalar1 = scalarField (0).eval ( iElem, P, time );
    const Real scalar2 = scalarField (1).eval ( iElem, P, time );

    const Vector vector = vectorField (0).eval ( iElem, P, time );

    return ( std::sin (scalar1) + scalar2 * scalar2 ) /  vector (0);
}

} // Namespace DataProblem