Question

我正在尝试学习xeon phi编程。

我在cpu上运行此代码，我正在使用offload pragma来处理我想在协处理器上运行的部分。

由于我正在编译cpu并使用offloads，我正在使用：

export MIC_ENV_PREFIX=MIC
export MIC_OMP_NUM_THREADS=120

以指定线程数。

我的问题：

1）运行代码，显示总是使用40个线程。

2）一次又一次地运行代码而不编译，我得到了   不同的时间结果。

代码：

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <omp.h>
#include <sys/time.h>

#include <cilk/cilk.h>
#include <cilk/reducer_opadd.h>


typedef CILK_C_DECLARE_REDUCER(float) reducer;


double dtime()
{
    double tseconds = 0.0;
    struct timeval mytime;
    gettimeofday(&mytime,(struct timezone*)0);
    tseconds = (double)(mytime.tv_sec + mytime.tv_usec*1.0e-6);
    return( tseconds * 1000 );
}

float openMPIntegration(

    int N,
    float * const ioA )
{

    float res = 0;

#if DOFFLOAD
    #pragma offload target (mic) 
    {
#endif

    #pragma omp parallel for reduction(+:res)
    for ( int i = 0; i < N; i++ )
    {
        res += ioA[ i ]; 
    }

#if DOFFLOAD
}
#endif

    return res;

}

float CilkIntegration(

    int N , 
    float * const ioA )
{


float res = 0;
#if DOFFLOAD
    #pragma offload target (mic) 
    {
#endif

    CILK_C_REDUCER_OPADD( sum, float , 0);
    CILK_C_REGISTER_REDUCER(sum);

    cilk_for ( int i = 0; i < N; i++ )
    {
        REDUCER_VIEW(sum) += ioA[ i ];
    }

    res = sum.value;
    CILK_C_UNREGISTER_REDUCER(sum);

#if DOFFLOAD
}
#endif

    return res;
}    

int main()
{
    int NbOfThreads;
    double tstart, tstop, ttime;

    int N = 1000000;
    float * A = (float*) _mm_malloc( N * sizeof(*A) , 32 );

    //fill A
    for ( int i = 0; i < N; i++ )
        A[ i ] = i;

#if DOFFLOAD
    #pragma offload target (mic)
#endif

    #pragma omp parallel
    #pragma omp master
    NbOfThreads = omp_get_num_threads();

    printf("\nUsing %d threads\r\n",NbOfThreads);

    tstart = dtime();   

    float openMPRes = openMPIntegration( N , A );

    tstop = dtime();    
    ttime = tstop - tstart;
    printf("\nopenMP integration = %10.3lf msecs \t value = %10.3f", ttime ,openMPRes);


    tstart = dtime();   
    float CilkRes = CilkIntegration( N , A );

    tstop = dtime();    
    ttime = tstop - tstart;
    printf("\nCilk   integration = %10.3lf msecs \t value = %10.3f", ttime,CilkRes);

    printf("\n");
    _mm_free( A );

    return 0;

}

我正在编译：

icc -std=c99 -DOFFLOAD -openmp -qopt-report -O3 xeon.c -o xeon

Answer 1

这不是严格意义上的OpenMP问题，因为它涉及不同的并行运行时模型的推荐组合，我认为您没有使用openmp标准化的卸载语法。简短的回答是没有实现建议组合OpenMP和cilkplus并行运行时模型。除此之外的下一步是，默认情况下，典型的OpenMP模型在一段时间内将硬件线程上下文的可用性阻止到OpenMP外部的线程模型，通常默认为0.200秒。使用omp减少而不是cilkplus减速器在风格上似乎更加一致，但在目前的实现中可能不是停止显示的决定。我猜你可能正在使用英特尔卸载模型，以便同时提供openmp标准和非标准卸载语法。

线程数和时间结果的一致性

1 个答案: