An Object Finite Element Library

HelmholtzBT3.h

Go to the documentation of this file.

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

                                    O  F  E  L  I

                           Object  Finite  Element  Library

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

   Copyright (C) 1998 - 2008 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

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

                         Definition of class HelmholtzBT3
               for Helmholtz Equation in a bounded domain using
                         3-node triangular finite element

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


#ifndef __HelmholtzBT3_H
#define __HelmholtzBT3_H


#include "Equa_Electromagnetics.h"
#include "Triang3.h"
#include "Line2.h"
#include "UserData.h"

using std::complex;

namespace OFELI {

class HelmholtzBT3 : virtual public Equa_Electromagnetics<complex<double> ,3,3,2,2>
{

 public :

    HelmholtzBT3()
    {
       _tr = NULL; _ln = NULL;
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class HelmholtzBT3 is constructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    HelmholtzBT3(Element *el);

    HelmholtzBT3(Side *sd);

    ~HelmholtzBT3();

    void LHS(const double &wave_nb);

    void BoundaryRHS(UserData<complex<double> > &ud);

 private :

   double        _rho, _mu;
   Triang3      *_tr;
   Line2        *_ln;
   double        _lx[3], _ly[3];
   Point<double> _x[3];
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani