問題タブ [dot-product]

私は SSE/SIMD に関する SO の質問をたくさん読みました (例: Getting started with SSE ) が、いまだに混乱しています。私が欲しいのは、C (C99 FWIW) の 2 つの倍精度浮動小数点ベクトル間のドット積だけです。私はGCCを使用しています。

double ベクトルを SSE 型に変換して元に戻す方法など、単純で完全な例を誰かが投稿できますか?

[2012 年 10 月 8 日編集]

これは、私がうまくまとめた SSE2 コードの一部です。

Pythonを使用して、スパース行列の行の内積をその行の転置で取得しようとしています。X2と呼ばれる巨大なスパース行列があります。そして、結果（単一の数値であると想定されています）をNjcというリストに保存しています。

ただし、プログラムを実行すると、結果は単一の数値ではありません。それらは次のようになります：（0、0）355

（0、0）295

（0、0）15

（0、0）204

（0、0）66

...。

残念ながら、私の疎行列は非常に大きいため、密行列にすることはできません（私の記憶は爆発します）。左側のカップルなしで右側の数字だけを取得する方法はありますか？

2つのオプションの上部にあるhttp://www.blackpawn.com/texts/pointinpoly/default.htmlの数式を使用しようとしています。（私はすでに下を使用しました）。

CrossProduct（ba、p1-a）などに頭を包むことができないようです。誰かが私のためにこれらを拡張してください。以下にあります。

私が理解しているように、彼らはこれに出てくるはずです。

しかし、Excelでこれを試してみると、再び間違った答えが返ってきます。

助けてください！8- |

一連の頂点 (ピンク) を回転させて、頂点のパターンの 1 つのエッジが三角形のエッジ (白) と一致するようにします。

これを行うには、最初にエッジを表す 2 つのベクトルを作成します。floretAB と triangleAB (緑)。次に、2 つのクロス積を見つけて、頂点を回転できる軸 (赤) を取得します。

次に、2 つのベクトルの間の角度を取得し、それを回転軸と共に使用して四元数を作成します。最後に、クォータニオンを中心にすべての頂点を回転させます。

回転前

_

どのような回転を生成する必要がありますか

_

ただし、頂点はクォータニオンを中心に正しく回転しますが、角度は次の図のように正しくありません。

これは、2 つのベクトル間の角度を取得するために使用しているコードです。私が間違っていることを理解していません：

これは本当に私を悩ませているので、これに関する助けは最も高く評価されます.

ありがとう、ジェームズ

編集：残りのコードは次のとおりです。

_

_

Heres 私の頂点シェーダー コード:

このコードでは、Learning Modern 3D Graphics Programming サイトのチュートリアルを使用して照明を計算しようとしています。ところで、これは素晴らしい本ですが、ここにあるコードは基本的にリストされたコードと同じですが、いくつかの名前が変更されていますが、いくつかの実験を通じて、非常に独特なものを発見しました. セグメンテーション違反は、

色 = vec3(強度);

行、その行を削除するとプログラムが正常に実行されるため、それはわかっています（もちろん、照明なしを除く）

ただし、行を置き換えると

float 強度 = dot(normal, toLight);

と：

フロート強度 = 3.14159;

それは 1.0f にクランプされ、すべて正常に動作します。立方体は白です。負の数を入力すると同じことが起こります。数値はゼロに固定され、立方体は黒です...

これは奇妙なことです。事実は次のとおりです。

• 出力ベクトルを要求された色に設定する行を削除すると、セグメンテーション違反はなくなります。
• 内積の計算を削除すると、セグメンテーション違反も消えます
• 強度が一定のフロートである場合、すべてが正常に機能します。
• 色を強度を含まないものに設定するとうまくいきます。
• 編集:強度変数を使用して色変数以外の変数を設定すると、たとえばvec3 temp = vec3(intensity);セグメンテーション違反が発生しないなど、興味深い情報が見つかりましたが、強度変数に関連しないものに色を設定しても問題なく動作しますcolor = vec3(1.0f);が、色を設定tempするとクラッシュしますが、これは非常に奇妙です。

どんな提案でも大歓迎です。ありがとう

また、opengl、c++ コードを含めますが、必要ではないと思います。間違っている場合は修正してください。

perlプログラムに次の2つの同じサイズの配列があるとします。

reduceコアモジュールで定義された関数を使用してドット積を取得しようとしていList::Utilますが、以下が機能しません。

予想される130ではなく50を出力として取得します。

以下を使用してRでマトリックスを生成しています。

この行列は、3D の点の座標を示します。Rで以下を計算する方法は?

たとえば、行列が次の場合、

次に出力=1*2 + 4*5

私は R を学ぼうとしています。

ありがとう

私は持っています

ここでdata_mat、は行列および列ベクトルと考えることができ、 withdata_vecのすべての列の内積を計算し、それを別のに格納する方法を探しています。data_matdata_vecvector < int > data_out (4)

とを使用した例http://liveworkspace.org/code/2bW3X5%241は、行列の列の合計を計算するために使用できます。for_eachtransform

同様の方法で（またはSTL関数を使用するわずかに異なる方法で）、ベクトルを使用して行列列の列内積を計算することは可能ですか？

問題は、'd2 = d1 + d2'トリックが、列の内積の場合には機能しないことです。これを解決するd3も含める方法がある場合（d3 = d3 + d1 * d2）、に三元関数が存在しないようですtransform。

パンダのデータフレームとシリーズの行列乗算をしたい

df は、

サーは、

私の望む結果は、行列の積です。

答えは、

numpy ドット演算子を使用してデータフレームを再構築することでこれを行うことができます

これを行うためのパンダのネイティブメソッドはありますか?

同じ次元の2つの行列の行方向の内積をできるだけ速く計算したいと思います。これは私がそれをしている方法です：

そしてもちろん、出力は次のとおりです。