雨谷の日和

過去19年で2,600を超えるアニメの第1話だけは見続けた僕のお勧めアニメがハズレなはずがない

Javaでの行列式計算(事前メモリ確保での実装)

それでは、Javaでも事前にメモリを確保してから処理を行うように書いてみよう。
以下のようになる。

class Test {

  private static long[ ][ ] array_buf;

  private static long expand(long[ ] array, int dim) {
    if(dim < 2)return array[0];
    long[ ] subarray = array_buf[dim - 2];
    long value = 0;
    for(int i = 0; i < dim; ++i) {
      int l = 0;
      for(int j = 0; j < dim; ++j){
        if(j != i){
          for(int k = 1; k < dim; ++k)subarray[(l * (dim - 1)) + k - 1] = array[(j * dim) + k];
          ++l;
        }
      }
      if(i % 2 == 0){
        value += (array[i * dim] * expand(subarray, (dim - 1)));
      }else{
        value -= (array[i * dim] * expand(subarray, (dim - 1)));
      }
    }
    return value;
  }

  public static void main(String[ ] args) throws Exception {
    long[ ] array = new long[ ]{
      1, 2, 3, 4, 6, 1, 7, 1, 9, 0, 0, 
      5, 1, 5, 7, 7, 3, 5, 3, 6, 1, 6, 
      7, 6, 4, 8, 4, 6, 1, 7, 1, 4, 8, 
      6, 5, 1, 0, 1, 5, 9, 1, 7, 5, 1, 
      2, 3, 1, 1, 5, 8, 4, 6, 3, 1, 5, 
      3, 1, 0, 3, 2, 9, 4, 8, 0, 3, 1, 
      3, 2, 8, 5, 0, 6, 3, 4, 4, 8, 5, 
      4, 9, 9, 4, 6, 1, 0, 1, 2, 1, 6, 
      9, 7, 6, 1, 6, 4, 5, 2, 9, 4, 4, 
      1, 2, 2, 2, 8, 3, 4, 9, 1, 2, 4, 
      6, 1, 6, 4, 4, 8, 6, 3, 4, 1, 4, 
    };

    array_buf = new long[10][11 * 11];
    System.out.println("result: " + expand(array, 11));
  }
}

次回と次々回をかけて、最適化無/有の実行時profileを見てみることにしたい。