Asdfgの日記

freefemを使いこなそう

freefem

2次元TMモード 電磁波動方程式 弱形式 その1(とfreefemをはじめた経緯など)

今日は2次元の場合のTMモード電磁波方程式の弱形式について記録しておきます 以下 /* から */まではコメントなのでお急ぎの方は下へ/************************************************************** 私がfreefemを取り組み始めたのは 2014年の春から…

【注意】複素関数空間の場合の転置コマンド「’」はエルミート転置!

私みたいにマニュアルに目を通さないで freefemを扱ってる初心者の方に向けての 注意喚起です周波数空間でマクスウェル方程式を解きたい場合 転置記号「’」を使いながら内積を使用することになる 例えば、 という数式をコード化すると Vh<complex> [Ex,Ey,Ez],[Wx,Wy,</complex>…

サーモグラフィックのような色諧調

初めての計算はコードを間違える可能性が大なので かならずplotで結果のグラフを確認しながらプログラムを行うことが多いです 大体グラフの変になり方からバグの位置が想定されることが多数ですただfreefemのデフォルトのplotの色諧調はオシャレとはいいがた…

freefem並列計算ことはじめ

これまで結局2次元計算ばかりでかなり良い結果を得てきていたので なかなか3Dを使ってこなかったのですが とうとう波動光学も本質的に2Dで近似しにくい系を 扱う羽目になってしまいました さてとりあえずexampleを動かすことから http://www.freefem.org/…

【エラー報告】3D 周期 missing face

3次元モデルでフォトニックバンドを計算しようとしている際のエラー (*下記はお試しの例) missing face 340 146 147 150 missing face 341 147 150 151 missing face 342 145 146 149 missing face 343 146 149 150 missing face 344 144 145 148 missing…