// Emacs will be in -*- Mode: c++ -*-
//
// ********** DO NOT REMOVE THIS BANNER **********
//
// SUMMARY: Language for a Finite Element Method
// RELEASE: 2.0     
// USAGE  : You may copy freely these files and use it for    
//          teaching or research. These or part of these may   
//          not be sold or used for a commercial purpose with- 
//          out our consent : fax (33)1 44 27 44 11        
//
// AUTHORS:  D. Bernardi, Y. Darmaillac F. Hecht,    
//           P. Parole O. Pironneau C. Prud'homme
// ORG    :          
// E-MAIL :   pironneau@ann.jussieu.fr     
//
// ORIG-DATE:     June-94
// LAST-MOD:      3-Jan-96 at 13:40:55 by Prud'homme Christophe
//
// DESCRIPTION:  
// DESCRIP-END.
//

#ifndef __FEM1C_H
#define __FEM1C_H 1

void rhsPDE(int quadra, creal*  fw, creal*  f, creal*  g);
       
float gaussband (creal*  a, creal*  x, long n, long bdthl, int first, float eps);

float pdeian(creal*  a, creal*  u, creal*  f, creal*  g, creal*  u0, 
        creal*  alpha, creal*  rho11, creal*  rho12, creal*  rho21, creal*  rho22,
        creal*  u1, creal*  u2, creal*  beta, int quadra, int factorize);


#endif /* __FEM1C_H */
