#pragma omp parallel for collapse(2)
    for(int particleRow = 0; particleRow < numberOfDipolesI; particleRow++)
    {
        for(int particleColumn = 0; particleColumn < numberOfDipolesJ; particleColumn++)
        {
            int ijx , ijy , ijz;
            //
            // Find location
            //
            i_jlocation(objData.dimensions,dipolePositionsI,dipolePositionsJ,particleRow,particleColumn,&ijx,&ijy,&ijz);

            int truePos = ijx*2*2*objData.dimensions[1]*objData.dimensions[2]+ijy*2*objData.dimensions[2]+ijz;
            //#pragma omp critical (flipCrit)
            {
                temp[truePos] = true;
            }

        }
    }