An Object Finite Element Library

EC2D2T3.h

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

                                    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 EC2D2T3
     for Eddy Current Problems in Two-Dimensions with a vector magnetic field
                           using the 3-Node triangle

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


#ifndef __EC2D2T3_H
#define __EC2D2T3_H


#include "Equa_Electromagnetics.h"
#include "Triang3.h"
#include "SpMatrix.h"
#include "SkMatrix.h"

#define GAUSS2 0.577350269189625764509148780501957

namespace OFELI {


class EC2D2T3 : virtual public Equa_Electromagnetics<double,3,6,2,4>
{

 public :

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

    EC2D2T3(const Element *el);

    EC2D2T3(const Side *sd);

    EC2D2T3(const Side *sd1, const Side *sd2);

    ~EC2D2T3()
    {
        if (_tr) {
           delete _tr;
           _tr = NULL;
        }
    }

    void RHS(double coef=1.);

    void FEBlock(double omega);

    void BEBlocks(size_t n1, size_t n2, SpMatrix<double> &L, SpMatrix<double> &U, SpMatrix<double> &D);

    std::complex<double> Const(double omega, const Vect<double>  &u, std::complex<double>  &I);

    double MagneticPressure(const Vect<double>  &u);

 private :

   double        _a3, _ll1, _ll2;
   Triang3       *_tr;
   Point<double> _N1, _N2, _M1, _M2;
   size_t        _ns, _nt, _i1, _j1, _i2, _j2;
   size_t _log_det(const Point<double> &ck, const Point<double> &cl);
   complex<double> _ablog(size_t det, complex<double> a, complex<double> b, double t);
   void _Lkl(const Point<double> &uk, const Point<double> &vk,
             const Point<double> &ul, const Point<double> &vl,
             double &d1, double &d2);
   void _Dkl(const Point<double> &uk, const Point<double> &vk,
             const Point<double> &ul, const Point<double> &vl, double &d);
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani