C++ et éléments finis Note de cours de DEA (version provisoire)
C++ et éléments finis Note de cours de DEA (version provisoire)
C++ et éléments finis Note de cours de DEA (version provisoire)
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
24 CHAPITRE 4. MÉTHODES D’ÉLÉMENTS FINIS P1<br />
– H(i) qui calcule directement ∇w i par la formule (??).<br />
Listing 7 (sfem.hpp - la classe Triangle)<br />
class Triangle: public Label {<br />
Vertex *vertices[3] ; // tableau <strong>de</strong> trois pointeurs <strong>de</strong> type Vertex<br />
public:<br />
R area ;<br />
Triangle(){} ; // constructeur par défaut vi<strong>de</strong><br />
Vertex & operator[](int i) const {<br />
ASSERTION(i>=0 && i r<strong>et</strong>ourne un somm<strong>et</strong><br />
void s<strong>et</strong>(Vertex * v0,int i0,int i1,int i2,int r) {<br />
vertices[0]=v0+i0 ; vertices[1]=v0+i1; vertices[2]=v0+i2;<br />
R2 AB(*vertices[0],*vertices[1]);<br />
R2 AC(*vertices[0],*vertices[2]);<br />
area = (ABˆAC)*0.5 ;<br />
lab=r ;<br />
ASSERTION(area>=0) ; }<br />
R2 Edge(int i) const {<br />
ASSERTION(i>=0 && i =0 && i =0 && i