C++ et éléments finis Note de cours de DEA (version provisoire)

28 CHAPITRE 4. MÉTHODES D’ÉLÉMENTS FINIS P1
4.2.6 Le programme principale
Listing 14 (sfemGC.cpp)
#include
#include
#include
#include
#include
#include "sfem.hpp"
#include "RNM.hpp"
#include "GC.hpp"
R f(const R2 & ){return 1 ;} // right hand side
R g(const R2 & ){return 0 ;} // boundary condition
R u0(const R2 & ){return 0 ;} // initialization
class Laplacien2d: VirtualMatrice { public:
const Mesh & Th ;
typedef VirtualMatrice::plusAx plusAx ;
Laplacien2d(const Mesh & T) : Th(T) {} ;
void addMatMul(const KN_ & x, KN_ & Ax) const {
Axi+=
(∇w
K∈Th K
i |K, ∇w j |K) xj si i n’est pas sur le bord
for (int k=0 ;k