00001 // Emacs will be in -*- Mode: c++ -*-
00002 //
00003 // ********** DO NOT REMOVE THIS BANNER **********
00004 //
00005 // SUMMARY: Language for a Finite Element Method
00006 //
00007 // AUTHORS: C. Prud'homme
00008 // ORG :
00009 // E-MAIL : prudhomm@users.sourceforge.net
00010 //
00011 // ORIG-DATE: June-94
00012 // LAST-MOD: 13-Aug-00 at 23:18:21 by Christophe Prud'homme
00013 //
00014 // DESCRIPTION:
00015 /*
00016 This program is free software; you can redistribute it and/or modify
00017 it under the terms of the GNU General Public License as published by
00018 the Free Software Foundation; either version 2 of the License, or
00019 (at your option) any later version.
00020
00021 This program is distributed in the hope that it will be useful,
00022 but WITHOUT ANY WARRANTY; without even the implied warranty of
00023 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
00024 GNU General Public License for more details.
00025
00026 You should have received a copy of the GNU General Public License
00027 along with this program; if not, write to the Free Software
00028 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
00029
00030 */
00031 // DESCRIP-END.
00032 //
00033
00034 #ifndef __FEM1_H
00035 #define __FEM1_H 1
00036
00037 namespace fem
00038 {
00039 float id(const float x);
00040 //float norm2( float a);
00041 float norm(const float x, const float y);
00042
00043 void rhsPDE(int quadra, float* fw, float* f, float* g);
00044
00045 float gaussband (float* a, float* x, long n, long bdthl, int first, float eps);
00046
00047 float pdeian(float* a, float* u, float* f, float* g, float* u0,
00048 float* alpha, float* rho11, float* rho12, float* rho21, float* rho22,
00049 float* u1, float* u2, float* beta, int quadra, int factorize);
00050 }
00051 #endif /* __FEM1_H */