linear-algebra

1熱

0答えて

Qはスパークで配布行マトリックスとする、私はその転置Q'でQの外積を計算します。ただし、行行列にはmultiply()メソッドがありますが、引数としてローカル行列のみを受け入れることができます。コードイラスト（スカラ）：私はRowMatrixローカル一つと分散ではない分散のドット積を実行する val phi = new RowMatrix(phiRDD) // phiRDD is an

-6熱

1答えて

Scenekitを使用して、アイデンティティ四元数を取得するにはどうすればよいですか？

ああ、続けてください。私は常に世界のx、に向いているSCNNodeようSCNTransformConstraint姿勢制約を作成したいのy &のz軸に関係なく、ノードの親はについて/移動を向いていますか。私はこのような世界空間での作業姿勢の制約を作成することができます let orientationConstraint = SCNTransformConstraint.orientation

3熱

2答えて

Goで並列saxpyを実装してもコア全体でうまくスケーリングできない

saxpyの実装を実装しようとしていますが、両方ともブロックされていて、自分のマシンで使用できる8コアを使って並列計算できます。私は、自分のマシンのL1キャッシュに収まる小さいサイズのベクトルxとy（分割256kB - 128kBデータ、128kBコード）を連続して計算できるという前提から始めました。この仮定をテストするために、saxpy（BSS）のブロックされたシリアルバージョンとsaxpy（B

1熱

1答えて

Tridiagonal Gauss-JordanエリミネーションMatlabを試そう

タイトルから推測したように、私は三重対角ガウスジョーダン消去を試みています。私はデフォルトのソルバーなしでそれをしようとしています。私の答えが正しく出てこないので、私のコードにどこにエラーがあるのか援助が必要です。私が持っているコードを使ってA/bとxの値が違ってきています。 n = 4; #Range for diagonals ranged = [15 20]; rangesd

3熱

1答えて

このアルゴリズムを正しく実装していますか？

私は非可逆残差行列からの崩壊の影響を受けないブロック共役勾配アルゴリズムを実装しようとしています。しかし、私は無意味な結果を得ています（各反復において、Rcurrentのランクは、小さくならず、増加しなくてはなりません）。これは、Hao JiおよびYaohang Liによる論文「A breakdown-free block conjugate gradient法」に示されています。これはジュリア

-1熱

1答えて

np.eigとnp.svdの異なる固有値

np.eig（np.dot（A.T、A））** 2とnp.eig（A）で計算された固有値を比較しています。いくつかの値は同じですが、いくつかは同じではありません。誰が私になぜこれが起こるのか教えてもらえますか？

3熱

1答えて

netlib-javaのソースコードはどこにありますか？

スパークのMLLibにはhttps://github.com/fommil/netlib-javaがdependencyであると記載されています。ただし、このリポジトリはgithubから削除されました（上記のリンクは無効です）。 Googleのキャッシュされたバージョンのページはまだ存在し、資金不足のために図書館が廃止されるという注意があります。ライブラリはもはやメンテナンスされていませんか？

-2熱

1答えて

数式をPythonに変換するには？

私はPaperの2つの公式をPythonで使いたいと思います。この式：私は、このコードにそれを翻訳： P = A #(size 24x24) Q = B #(size 24x24) sum_of_all = 0 for row in range(0,P.shape[0]): for column in range (0,P.shape[1]):