OpenMP "parallel for"在特定程序中很奇怪

我大约一个月前开始编程，最近我一直在尝试用c++程序学习使用OpenMP进行多核开发。在我编写的大型程序中，我无法使OpenMP正确地为我工作。我想并行执行下面的循环。

for(int iAngle=0; iAngle<int(nAngles); iAngle++){
    for(int k=0; k<int(nRadSamples); k++){
        int x=coordX[k][iAngle];
        int y=coordY[k][iAngle];
        for(int i=yMin; i<int(yMax); i++){
            int iy=i+y;
            for(int j=xMin; j<int(xMax); j++){
                response[iAngle][i][j] += inputSlice[iy][j+x];
            }
        }
    }
}

coordX、coordY、response和inputSlice

vector<vector<int> >
vector<vector<int> >
vector<vector<vector<float> > > 
vector<vector<float> >

。添加

#pragma omp parallel for

作为循环上方的行。

我不认为我的问题来自访问响应和inputSlice共享变量，因为即使是基本的OpenMP代码在这个特定的程序中执行奇怪。例如基本程序

#include <stdio.h>
#include <omp.h>
int main() {
//////////////////////
#pragma omp parallel for
for(int i = 0; i<2; i++){
    int thread = omp_get_thread_num() + 1 ;
    int numThreads = omp_get_num_threads();
    printf("Thread %d of %d printing %dn", thread, numThreads,i);
}
 ////////////////////

return 0;
}

Thread 2 of 8 printing 1
Thread 1 of 8 printing 0

但是当我复制并粘贴///////////////边框内的代码到

int main(int argc, char* argv[])

函数，它输出

Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1
Thread 1 of 1 printing 0
Thread 1 of 1 printing 1   

就好像每个线程单独执行整个for循环，彼此看不到，这可以理解地减慢了程序的运行时间。'i'没有在我的程序的main函数的作用域中声明。

当编译我的程序时，我正在链接到其他库，包括添加-pthreads， -lguide和-ltiff，除了-fopenmp，当编译gcc4.4时。

任何帮助与这个特定的问题，或我的编码风格一般将非常感激!我的头一直在敲击键盘。