linear-algebra

2熱

1答えて

私はPythonコードで計算されたNxN対称および三重対角行列を持っています。対角化したいのです。具体的にはN = 6000を扱っていますが、マトリックスが大きくなることがあります。それは疎であるので、対角化する最良の方法は、アルゴリズムscipy.sparse.linalg.eigsh()を使用することであると仮定しました。これは、他の疎な対称行列（三重対角行列ではありませんが）でもうまくいき

1熱

1答えて

線形システムを解くsympy問題

私はJupyter QtConsole上で動作するPython v.3.6を使用しています。私は、予測を調査スコアと結びつける個人プロジェクトのために、Sympyを使ってデータセット上の線形代数を試みています。本質的には、N = 14の線形方程式とM = 5の未知数を持つ拡張行列を設定し、システムを解決しようとしています。私の問題は、私は私の拡大行列にsolve_linear_systemコマン

1熱

1答えて

テンソルフローの線形モデル

私はTensorflowで単純な線形モデルを生成しようとしていました。ここでコードが... N = 400 features = 100 nSteps = 1000 data = (np.random.randn(N, features), np.random.randint(0, 2, N)) W = tf.placeholder(tf.float32, shape=(featur

1熱

1答えて

ベクトルと顔のOpenGL交差

私は地形上を移動すべき地形とオブジェクトを持っているので、オブジェクトの起源を検出できる関数を作ったので、オブジェクトのYを設定する必要があります顔に垂直なオブジェクトの原点からのベクトルとの交点の正確な高さ、私は顔の3つの頂点を持つので、その法線とその起点と最大点と最小点（境界ボックス）を計算できます。 enter image description here 編集：別の形で質問：私はそ

0熱

1答えて

インテリジェントなmklライブラリを使ってダイレクトソルバーを使って線形システムを解くと、dss_createは例外をスローします

次は私のコードのFEM問題を解決するクラスです。式の線形システムの係数行列は疎なので、このシステムを解くにはintel mkl direct solver for sparse matricesを使用しています。クラスは線形システムを解くために使用されています。私は簡単な問題のためにこれを試しましたが、結果をmatlabと比較することでうまく動作します。私はすべてのリンクが正しいことを確信して

1熱

2答えて

JavaでNaNを生成するコレスキー分解

これはmaths.seかSOのどちらの質問か分かりませんが、私のJavaに関連していると思います私はGaussian Processes（R&W）のテキストブックに従い、Javaでいくつかの例を実装しています。いくつかの例の1つの共通ステップは、共分散行列のコレスキー分解を生成することです。私の試みでは、限られたサイズ（33x33）までのマトリックスの成功を得ることができます。しかし、より大きな

2熱

1答えて

座標がどこにあるかを見つける

私はPoints[]の配列を渡しているので、交点を求める必要があります。私のテストでは、常に「交差していない」ことを示す座標セットが見つかりました。これを考慮する方法が必要です。与えられた座標の下で、交点を計算するためにC＃の公式/方程式はどうなるでしょうか？ private void button1_Click(object sender, EventArgs e) { Point[

0熱

1答えて

Tensorflow：tf.gradients（）は正確に何を返しますか

私はここで混乱しています。のは、我々は簡単なグラフを持っているとしましょう： a = tf.Variable(tf.truncated_normal(shape=[200, 1], mean=0., stddev=.5)) b = tf.Variable(tf.truncated_normal(shape=[200, 100], mean=0., stddev=.5)) add = a+b

2熱

1答えて

OpenGLはオブジェクトの顔カメラを作る

OpenGL ES3、C++ 11、Android、私はオブジェクトとカメラを持っているので、私はオブジェクトをカメラに向かって見て、物体とカメラのビュー行列、カメラは私が次のことを試してみましたので、私は、画面をタッチすると、オブジェクトがカメラに向けて回転させたい、オブジェクトの周囲を回転： glm::vec3 up(0, 1, 0); glm::vec3 camFor(glm::cross

0熱

2答えて

特異な正方行列を解くためにlu分解を使用することについての知識はありますか？

私はランクが不足している（ランク= 9）特異行列A（10 * 10）を持っており、Aの範囲空間にあるベクトルbを持っています。具体性のためにここに私の array([[ 0. , 0. , 0. , 0.86826141, 0. , 0. , 0.88788426, 0. , 0.4089203 , 0.88134901], [ 0. , 0. , 0.