FreeFem 3.5.x
|
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 */