問題タブ [finite-element-analysis]

私はpython2.7プログラムを実行したい(これ)。Python 2.7 の依存関係がインストールされているため、私は多くの問題を抱えています (私はこれに午後全体を費やしています)。

設定

Debian ベースの Ubuntu 16.04 64 ビット ([Mint XFCE 18) を実行しています。私のコンピューターは、2011 年製の Dell Inspiron N5110 で、デュアル ブート W7/U16 を備えています。キーボードから画面へのインターフェイスは、この世界では本当に新しいものであり、これを単独で解決するには、おそらくそれについてもっと学ぶ必要があります.

収入

aptを使用してさまざまなプログラムをインストールすることから始めました。

sudo apt-get install -y git python-dev libpython-dev libevent-dev libsuperlu-dev libblas-dev liblapack-dev

必要なプログラムをgit cloneした後、依存関係をインストールしました。理由はわかりませんが、うまくいきませんsudo pip install pysparseでした。代わりに、sudo pip install csc-pysparseうまくいきました。

問題

プログラムを実行すると、次のように表示されfrom pysparse import superlu, itsolvers, preconますImportError: cannot import name superlu。

なんで ？pip は依存関係の問題を解決して superlu をインストールしていませんか? superlu を手動でインストールする必要がありますか、それとも csc-pysparse の代わりに pysparse をインストールする必要がありますか?

（そして、甘やかしてください。これは、githubでの私の最初の投稿だったので、stackoverflowでの私の最初の投稿です）

いくつかの機能を拡張する必要があるため、コメントが非常に少ないコードを扱っています。変更が影響するかどうかを知るために、コードのすべての部分が何をするかを理解する必要があります。

このコードは、2D メッシュ生成 (ライブラリ Triangle を使用) を処理し、その上で PDE を解決しています。

ここに私が理解できないコードがあります：

最後に、正確に phiMapMaster_ とは何か、それをどのように使用できるかを知りたいです!

これはテンプレート クラス FiniteElement 内のメソッドであり、ORDER=1 と仮定し、Integrator:

(Point は std::complex と同じです)。ここでは、入力として三角形 (単純に 3 つの点を反時計回りに並べたもの)、点 (三角形の内部)、およびフーリエの係数を受け取るメソッド evaluate_point を示します。三角形で定義された基底

私は、データを格納するために非常に大きな配列 ([1 000 000] 要素で構成される配列) を必要とする私の学士号のために、C で記述された FEM プログラムに取り組んでおり、それらを操作しています。同様の異常に巨大なサイズを持つ 2D 配列も使用します ([100 000][100 000] など)。

プログラムはエラーを表示せずに終了します。プログラムが (大きな 2D 配列を使用する) メッシュを生成するためのルーチンを実行しようとすると、プログラムがクラッシュすることがわかっています。

A、B、Cを頂点とする三角形要素があります。三角形の面の法線Nに圧力Pを加えました。ここで、三角形要素の節点に作用する力を計算する必要があります。そのFEMプロセスで行われます。数式でそれを行う方法を教えてください。三角形のバルクデータを次のようにします。

数学的な方法で教えてください。ご尽力いただきありがとうございます

メッシュ用の小さなアンラッパーを作成しようとしています。これは、有限要素法を使用して、平坦化された表面と未加工の表面の間の最小の線形応力を解決します。現時点では、結果を得るために固定された頂点がいくつかあります。これがないと、三角形はランダムに回転および移動されます...

しかし、この固定は問題には必要ないため、より良い解決策は特異行列を直接解くことです。Petsc は、ヌル空間に関する情報を提供することによって、特異なシステムを解くいくつかの方法を提供します。http://www.mcs.anl.gov/petsc/petsc-current/docs/manual.pdf#section.4.6 Eigen でこれに代わるものがあるのだろうか。そうでない場合、頂点を固定/固定せずにこの問題を解決する他の可能性はありますか?

ありがとう、よろしく

詳細については、このリンクも参照してください: 開発履歴

ここで三角形要素を使用した例を見つけました。次に、メッシュの生成方法を変更して、三角形の要素を長方形の要素に置き換えましたが、それらを統合する方法がわかりません。これが私のバージョンです：

三角形の要素を統合するために使用されたコードのセクションの編集を開始しましたが、どのように進めればよいのか、または長方形の要素についても同様の方法で行われるのかどうかはわかりません。

アップデート

だから私はドヒョンが提案したことを私の限られた理解で取り入れようとしました、そしてこれが私が今得たものです：

数値積分の代わりに定積分を使用できると思いますが、結果は元のプログラムのものと一致しません。

私のプログラムの結果

オリジナルプログラムの結果