Post date: 2018/02/14 6:08:44, Last update: 2024/09/08

adaptmesh() を実行し，メッシュが変更された場合，本来は明示的な再代入が必要

• varf+matrix の場合はadaptmesh() 実行後に，明示的に再代入が必要

fespace Vh(Th,P1);

Vh u, uu; 

varf Mat(u,uu) = ...  

adaptmesh(Th); 

u=u; // !! 

matrix A=Mat(Vh,Vh); 

u=A^-1*u; 

• • ベクトル変数 [u1,u2] の場合は、ポインタで指定 u1[]=u1[];

• problem の場合は再代入しなくても大丈夫らしい

Vh u, uu; 

problem Mat(u,uu) = ...  

adaptmesh(Th); 

Mat;