double gweights = 0.0f;          // 0.020s

for (int j = yT; j <= yB; ++j) {
      // 0.666s
  for (int k = xL; k <= xR; ++k) {  // 2.347s
    const double gwTt = gwT[k];
  // 1.048s

    const int di = s1[k] - s2[k];
  // 2.282s
    const int di2 = di * di;
       // 2.530s
    const double di3 = di2 * gwTt;
// 10.685s

    gweights += gwTt;
           // 1.146s
    diff += di3;
                // 7.087s
  }

}

ICCでも速度出ない…。

int * doubleって遅い?