An Object Finite Element Library

Laplace1DL2.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 Laplace1DL2
              for 1-D Laplace equation using 2-node line element (P1)

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


#ifndef __LAPLACE_1DL2_H
#define __LAPLACE_1DL2_H

#include <float.h>
#include <stdlib.h>
#include <math.h>

#include "OFELI_Config.h"
#include "Equation.h"
#include "Element.h"
#include "Line2.h"

namespace OFELI {

class Laplace1DL2
{

public :

    Laplace1DL2(Element *el);

    ~Laplace1DL2()
    {
       delete _ln;
#ifdef _OFELI_DEBUG
       std::clog << "An instance of class Laplace1DL2 is destructed.\n";
       std::clog << "File : " << __FILE__ << ", Line : " << __LINE__ << endl;
#endif
    }

    void Matrix(double coef=1.);

    void BodyRHS(const Vect<double> &f);

    void BoundaryRHS(int n, double p);

    double *A() { return _a.getArray(); }

    double *b() { return _b.getArray(); }

private :

   Line2                  *_ln;
   Element                *_el;
   LocalMatrix<double,2,2> _a;
   LocalVect<double,2>     _b;
};

} /* namespace OFELI */

#endif

Copyright © 1998-2007 Rachid Touzani