fe-vect.edp

来自「FreeFem++可以生成高质量的有限元网格。可以用于流体力学」· EDP 代码 · 共 15 行

mesh Th = square(10,10); 
fespace Wh(Th,[P2,P2]);	// R2->R2
Wh [F1,F2]; //def de l'inconnu et des fonct tests

F1 = 1;
// il y a un bug ici car dans un element fini vectoriel
// il est impossible defini juste une composante
//  il faut faire:

cout << F1(0.5, 0.6) << endl;
cout << F2(0.5, 0.6) << endl;
cout << F1(0.5, 0.6);
[F1,F2]=[1,2]; 
cout << " -- F1 " <<  F1(0.5, 0.6) << endl;
cout << " -- F2 " << F2(0.5, 0.6) << endl;