雨谷の日和

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

Javaでの行列式計算(ループ、動的)

今度は、Javaで実測を行う。
ループ展開したものを元に、実際に動くソースコードを書くと、以下のようになる。

public class Test {

  public static long expand(long[ ] arraybase, int dim){
    long[ ][ ] array = new long[dim][ ];
    long[ ] return_value = new long[dim];
    int[ ] i = new int[dim];
    int d = dim;
    int flg = 0;

    array[dim - 1] = arraybase;
    for(int index = 0; index < dim; index++)return_value[index] = 0;
    for(int index = 0; index < dim; index++)i[index] = 0;

    do {
      if(d < 2){
        return_value[d - 1] = array[d - 1][0];
        d++;
        flg = 1;
      }else{
        if(flg == 0){
          if(i[d - 1] == 0){
            array[d - 2] = new long[(d - 1) * (d - 1)];
            return_value[d - 1] = 0;
          }
          if(i[d - 1] < d){
            int l = 0;
            for(int j = 0; j < d; ++j){
              if(j != i[d - 1]){
              for(int k = 1; k < d; ++k)array[d - 2][(l * (d - 1)) + k - 1] = array[d - 1][(j * d) + k];
                ++l;
              }
            }
            d--;
            i[d - 1] = 0;
          }
        }else{
          if(i[d - 1] < d){
            if(i[d - 1] % 2 == 0){
              return_value[d - 1] += (array[d - 1][i[d - 1] * d] * return_value[d - 2]);
            }else{
              return_value[d - 1] -= (array[d - 1][i[d - 1] * d] * return_value[d - 2]);
            }
            i[d - 1]++;
            flg = 0;
          }
          if(i[d - 1] == d){
            d++;
            flg = 1;
          }
        }
      }
    }while(d <= dim);
    return return_value[dim - 1];
  }

  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,
    };
    System.out.println("result: " + expand(array, 11));
  }

}

次回以降、これを最適化無し・最適化有り・AOTの順に実測を行い、以前の結果と比較しよう。