user.cpp

OFELI is an object oriented library of C++ classes for development of finite element codes. Its main

/*==============================================================================

                                    O  F  E  L  I

                            Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2004 Rachid Touzani

   This program 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; Version 2 of the License.

   This program 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 this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================
                        User Defined functions to prescribe data
  ==============================================================================*/


#include "User.h"


complex<real_t> User::SurfaceForce(const Point<real_t> &x, int code, real_t time, size_t dof)
{
   real_t z = -2*OFELI_PI*cos(OFELI_PI*x.y);
   if (code==1)
     return complex<real_t>(0,z);
   else if (code==2)
     return complex<real_t>(z,0);
   else
     return complex<real_t>(0,0);
}


complex<real_t> User::ExactSolution(const Point<real_t> &x)
{
   real_t z = cos(OFELI_PI*x.y);
   return complex<real_t>(cos(2*OFELI_PI*x.x)*z,sin(2*OFELI_PI*x.x)*z);
}