scientific-computing

0熱

1答えて

[0 0 1]に沿って、点[x y z]からオブジェクトの表面までの距離を計算する必要があります（この段階では、単純な直腸ですが、後で任意の形状になります）。単位ベクトルを使用してサーフェスを定義し、[0 0 1]に沿ってすべてのプレーンまでの距離を求める線形代数計算を行うことができますが、コーディングとJavaについてかなり新しい人として、ライブラリーがあるかどうか、長期的には複雑な凸オブジェ

8熱

4答えて

matplotlib/pyplot/numpy/pythonでPlot Mandelbrot

私はPythonには新しく、Python "Scientific lecture notes Release 2013.1 Release"のチュートリアルに続いて学習します。以下のsrc screenshot（Pg 71）でこのMandelbrotの問題を解決してください。可能であれば、プログラミングのコンセプトは私にとって新しいものなので、ステップごとの説明を説明に付けてください。 http:

0熱

1答えて

特定のスレッド数に対するOpenMPの最適化

私は科学的なコードを開発していますので、もちろん速度は本質です。今ではその移植性は本当に問題ではないので、プログラムをコンパイルするときにopenmpスレッドがいくつあるか知っています。この情報を使用して追加の最適化を実行することはできますか？はいの場合はどうすればよいですか？この質問は非常に広いと指摘されて以来、私はそれを自動的にコンパイルの最適化に少し減らしたいと思います。コンパイラのフラグ

6熱

1答えて

HDF5ストレージオーバーヘッド

私はHDF5ファイルに多数の小さなデータセットを書きます。その結果のファイルサイズは、入れているデータの素朴な表から予想されるものの約10倍です。私のデータは階層的に編成されています group 0 -> subgroup 0 -> dataset (dimensions: 100 x 4, datatype: float) -> dataset (dimensi

1熱

1答えて

Pypparse installation ubuntu

ubuntuにpysparseをインストールする際に問題があります。私のpython-スパースパッケージをインストールしたのですが、例のコードを実行するときに from pysparse.sparse import spmatrix from pysparse.direct import superlu import numpy n = 100 A = poisson2d_sym_blk(

0熱

1答えて

Ubuntuは私の計算を殺しています

私はubuntuマシンで重い科学計算を実行しようとしています。このコードはPythonで書かれており、実行可能ファイルとしてマークされています。しかし、約1時間後には、システムが "killall python"コマンドを実行したように、 "terminated"以外のエラーメッセージなしで計算が終了します。なぜubuntuで計算が終了するのですか？無限ループはなく、解はうまく収束します（ただし

1熱

3答えて

最適化でコンパイルしたときにのみNaNが表示されるときに、数値コードのデバッグを開始するにはどうすればよいですか？

数値コード、デバッグを開始するための一般的な戦略は何ですか：積極的な最適化フラグを指定してコンパイルコードは、実行する（-O0を意味する）-gしてコンパイル出力コードではNaNのとInfの者をoccational生成しますデバッガでNaNとInfをもう生成しないのですか？私はポートランドグループC++コンパイラ、PGCCで働いていた。この場合、および使用される最適化オプション、その後 -w

6熱

1答えて

NumPy'sのスピードアップ機能

いくつかのしきい値を超える数字の1D配列のすべての値のインデックスを抽出しようとしています。配列の長さは1e9です。私のアプローチはNumPyで次のようである： idxs = where(data>threshold) これは受け入れられない上向きに20分のものを取ります。この機能をどのようにスピードアップできますか？または、より速い選択肢がありますか？（具体的には、それは1.86 GHz

14熱

3答えて

scipyに多次元積分を組み込む

動機付け：私は多次元積分を持っています。これは完全性のために私が以下に再現したものです。大きな異方性がある場合には、第2ビリアル係数の計算から来ている：ここでWは、すべての変数の関数です。これは既知の関数で、Python関数を定義することができます。プログラミングの質問：scipyにこの式を統合するにはどうすればよいですか？私は2つのトリプルクワッド（scipy.integrate.tplqu

0熱

1答えて

scipy.interpolate.bisplevによって与えられるBスプラインの偏微分の意味はなんですか？

私は、z = f(x,y)と記述することができるので、scipy.interpolate.bisplrepとbisplevという良い結果が得られます。は、結果が評価点における次数の派生物であるように、パラメータdx=nおよび/またはdy=nと共に使用することができる。私は（MATLABでsurfatureと呼ばれる）の平均とガウス曲率を計算し、これを使用することを計画し、そしてそれが一度に部分的