下三角行列と上三角行列との分解 † 次のようにある行列Aを下三角行列Lと上三角行列Uに分解することを考えます.
下三角行列と上三角行列との分解 † 次のようにある行列Aを下三角行列Lと上三角行列Uに分解することを考えます.
load "msh3" // Parameters int nn = 20; // Mesh quality // Mesh int[int] labs = [1, 2, 2, 1, 1, 2]; // Label numbering mesh3 Th = cube(nn, nn, nn, label=labs); // Remove the ]0.5,1[^3 domain of the cube Th = trunc(Th, (x < 0.5) | (y < 0.5) | (z < 0.5), label=1); // Fespace fespace Vh(Th, P1); Vh u, v; // Macro macro Grad(u) [dx(u), dy(u), dz(u)] // // Define the weak form and solve solve Poisson(u,
FreeFEM 1950 年代に構造解析から誕生した「有限要素法」は偏微分方程式の標準的な数値解法の 一つとして定着しましたが、 複雑な境界形状に良く適合した要素分割ができる ポテンシャル分布に適応した自由な要素分割ができる 材料定数の取扱が容易 といった際だった利点の半面、構造データとプログラムの作成が面倒で、 電卓のように気軽に使えるものではありませんでした。 ところが、 FreeFEM という見事なアイデアで、この状況はずいぶん改善されたと思います。 このプログラムは、フランスの UPMC (当時 Pierre et Marie Curie 大学)の スタッフによって開発された もので、Gfemと呼ばれるPascal風の言語により、 1) 境界形状 2) 境界条件 3) 偏微分方程式 を記述するだけで、メッシュの自動生成から数値解を求める作業のすべてをやって くれますから、何とも気軽
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
処理を実行中です
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く