Quelques EDP simples r´esolues avec FreeFem++, Astuces et Trucs

More documents

Recommendations

Info

eal tgv = 1e30 ; // a hude value of exact penalisation of boundary condition // the variatonnal form of the problem: varf a(uh,vh) = // definion of the problem int2d(Th)( dx(uh)*dx(vh) + dy(uh)*dy(vh) ) // bilinear form - int2d(Th)( f*vh ) // linear form + on(1,2,3,4,uh=fd) ; // boundary condition form matrix A=a(Vh,Vh,tgv=tgv,solver=CG) ; matrix AA=a(Vh,Vh) ; // two version of the problem // the mass Matrix construction: varf vM(uh,vh) = int2d(Th)(uh*vh) ; matrix M=vM(Vh,Vh) ; // to do a fast computing of L2 u’*(w=M*u)) norm : sqrt(
Aii=A.diag ; // g<strong>et</strong> the diagonal of the matrix rhs = a(0,Vh,tgv=tgv) ; Ik =0 ; uhp=0 ; Vh lambda=0 ; for(int iter=0 ;iter
Page 1 and 2: Quelques EDP simples résolues avec
Page 3 and 4: HISTORIQUE Idée suivre le cours de
Page 5 and 6: The main characteristics of FreeFem
Page 7 and 8: The main characteristics of FreeFem
Page 9 and 10: trouver u solution de Equation de l
Page 11 and 12: Equation de la chaleur/ table de co
Page 13 and 14: Equation de la chaleur avec de matr
Page 15 and 16: Un problème de nappe phréatique T
Page 17 and 18: The new domain is : Ω k+1 = F(Ω
Page 19: Inequation Variationnelle mesh Th=s
Page 23 and 24: Methode de Joint Methode de joint ,
Page 25 and 26: Methode de Joint macro defspace(i)
Page 27 and 28: Methode de Joint lh[]=0 ; varf vDD(

Quelques EDP simples r´esolues avec FreeFem++, Astuces et Trucs

Create successful ePaper yourself

Delete template?

Save as template?