c言語 配列の計算について

Question

bの配列に0以外の数字が入っていたらaの配列に加算するメソッドを、以下の繰り返し文で書きました。
もっと簡単な(計算量少ない)方法がある気がするのですがどなたか教えていただけないでしょうか？

for(i = 0; i < 10; i++){
for(j = 0; j < 10; j++){
if(b[i][j] != 0){
a[i][j] = a[i][j] + b[i][j];
}
}
}

Answer 1

加算なら判定 if(b[i][j] != 0) は必要無いでしょう。

Comments

[edamame11]

計算量としてもっと早くなる書き方はございますでしょうか?
例えば、繰り返しなど使わずに1回の処理で加算する方法などです。

[LouiS0616]

SIMD演算ですかね(多分)。
行列処理を提供しているライブラリのソースを読めば参考になるかもしれません。

計算量を減らしたい理由にも依るでしょう。
活用に支障が出るほどに処理が遅く、且つボトルネックが当該行列演算だと特定できているのでしょうか。

[edamame11]

単純にシンプルな書き方が他にありそうな気がして投稿させていただきました。
SIMD演算を調べてみます。
ありがとうございました。

[LouiS0616]

単に高速化したいだけなら、Eigenなど使ってみると良いかもしれません。
行列演算に係る様々な機能を提供しているようです。

Answer 2

a[10][10], b[10][10] 双方の各要素が連続領域(スキマなし)に並んでいるなら

  int a[10][10];
  int b[10][10];
  ...
  int* pa = &a[0][0];
  int* pb = &b[0][0];
  for ( i = 0; i < 100; ++i ) {
    if ( pb[i] != 0 ) pa[i] += pb[i];
  }

のように、二重loopなしに書けます。
※ 計算量は変わりませんけどね。

Comments

[edamame11]

計算量はあまり変わらないんですね。
ありがとうございます。

[episteme]

その計算量が 計算機科学で言うところの時間計算量 であるなら、
あまり変わらない どころではなく まったく変わらない です。
実際にふたつのコードを実行し、ナノ秒単位で比較すれば差がでる"かもしれません"。

Answer 3

ざんねんながら、C言語では要素を全部舐めていくしか方法はありません

＃ポインタを使うなどして記述量を削減する程度しかできませんなー

Comments

[edamame11]

そうなんですね。
ありがとうございます。