スピードアップforループwith .push_back

Question

OpenMPまたはCUDAを使用してこのループをスピードアップできるかどうか疑問に思っています。現在のところ、それは逐次処理と細かい動作しますが、私のコーディングを最適化しようとしている私は「M：私は効果なしで、次の試してみた処理を並列化するための努力において

for (int curCol = 0; curCol < numRows; ++curCol){  //Long Loop 
     int lb = csc_colIndices[curCol]; 
     int ub = csc_colIndices[curCol + 1]; 

     // push back the diagonal value to L matrix 
     vec_L_val[curCol].push_back(1.0f); 
     vec_L_indices[curCol].push_back(curCol); 

     for (int curIndex = lb; curIndex < ub; ++curIndex){ 
      int curRow = csc_indices[curIndex]; 
      float curVal = csc_val[curIndex]; 

      if (!Equal(curVal, 0) && curRow <= curCol){// U entry 
       vec_U_val[curCol].push_back(curVal); 
       vec_U_indices[curCol].push_back(curRow); 
      } 
      else if (!Equal(curVal, 0) && curRow > curCol){// L entry 
       vec_L_val[curCol].push_back(curVal); 
       vec_L_indices[curCol].push_back(curRow); 
      } 
     } 
    } 

：

#pragma omp parallel for private(curCol) shared(curIndex) 

マイ？

どのように私はこのコードを向上させることができます...疑いが.push_backの使用であるが、私は間違っている可能性が

Answer 1

最初の事前割当て全部：

for (int curCol = 0; curCol < numRows; ++curCol) 
{ 
    vec_L_val[curCol].resize(SIZE_OF_THE_INNER_VECTOR); 
    vec_L_indices[curCol].resize(SIZE_OF_THE_INNER_VECTOR); 
    vec_U_val[curCol].resize(SIZE_OF_THE_INNER_VECTOR) 
    vec_U_indices[curCol].resize(SIZE_OF_INNER_VECTOR) 
} 

内部ループ構造の内部に再アロケーションが必要ないため、内側ループがおそらくより速く動作するでしょう。

for (int curCol = 0; curCol < numRows; ++curCol){  //Long Loop 
    int lb = csc_colIndices[curCol]; 
    int ub = csc_colIndices[curCol + 1]; 

    // push back the diagonal value to L matrix 
    vec_L_val[curCol].push_back(1.0f); 
    vec_L_indices[curCol].push_back(curCol); 

    for (int curIndex = lb; curIndex < ub; ++curIndex){ 
     int curRow = csc_indices[curIndex]; 
     float curVal = csc_val[curIndex]; 

     if (!Equal(curVal, 0) && curRow <= curCol){// U entry 
      vec_U_val[curCol].push_back(curVal); 
      vec_U_indices[curCol].push_back(curRow); 
     } 
     else if (!Equal(curVal, 0) && curRow > curCol){// L entry 
      vec_L_val[curCol].push_back(curVal); 
      vec_L_indices[curCol].push_back(curRow); 
     } 
    } 
}