OpenMP优化功能

时间:2012-11-23 05:37:45

标签: c optimization openmp

我们已经获得了使用openmp优化结构不良程序的任务。我是编程方面的新手,所以我希望任何人都可以了解如何优化这个功能(我需要优化的众多功能之一):

void
entry_type3(F2D *sData, F2D *ones, F2D *quat, F2D **pos, F2D **vel)
{
    //Observation

   F2D *t;

   t = fSetArray(1, 3, 0);
    asubsref(t,2) = -9.8;


    F2D *accl = fDeepCopyRange(sData, 0, 1, 0, 3);
    F2D *gtemp = fMtimes( ones, t);
    F2D *gravity = quatRot(gtemp, quat);



    fFreeHandle(gtemp);
    fFreeHandle(t);

   t = fSetArray(3,3,0);
    asubsref(t,0) = 1;
    asubsref(t,4) = 1;
    asubsref(t,8) = 1;

    int n = ones->height;
    int i;
    for(i=0; i<(t->height*t->width); i++)
            asubsref(t,i) = asubsref(t,i)/STDDEV_ACCL;

   F2D *w = mcl( gravity, accl, t);

    generateSample(w, quat, *vel, *pos);
    fFreeHandle(t);


    //Motion model
  t = fMtimes(ones, accl);
    fFreeHandle(accl);
    accl = fMinus(t, gravity);
    fFreeHandle(w);
    fFreeHandle(gravity);
    fFreeHandle(t);




    F2D *is;
    #pragma omp parallel sections
    {
        #pragma omp section
        {
            F2D *is = quatConj(quat);
            F2D *s = quatRot(*vel, is);
            fFreeHandle(is);
            for(i=0; i<(s->height*s->width); i++)
            {
                 asubsref(s,i) = asubsref(s,i)*acclTimeInterval;
            }
            is = fPlus(*pos, s);
            fFreeHandle(*pos);
            *pos = fDeepCopy(is);
            fFreeHandle(is);
            fFreeHandle(s);
         }


        /** pos_ above stores: pos+quatRot(vel,quatConj(quat))*acclTimeInterval **/

        #pragma omp section
        {
            F2D *is = quatConj(quat);
            F2D *s = quatRot(accl, is);
            F2D* t = fDeepCopy(s);

            for(i=0; i<(s->height*s->width); i++)
            {
                asubsref(t,i) = 1/2*asubsref(s,i)*acclTimeInterval*acclTimeInterval;
            }

            /** t_ above stores: 1/2*quatRot(accl,quatCong(quat))*acclTimeInterval^2 **/

            fFreeHandle(s);
            fFreeHandle(is);


 s = randnWrapper(n,3);

            for(i=0; i<(s->height*s->width); i++)
            {
                asubsref(s,i) = asubsref(s,i) * M_STDDEV_POS;
            }

            /** s_ above stores: randn(n,3)*M_STDDEV_POS **/

       is = fPlus(*pos, t);
           fFreeHandle(*pos);
         *pos = fPlus(is, s);

            fFreeHandle(s);
            fFreeHandle(t);
            fFreeHandle(is);
    } 

}        
        //vel=vel+accl*acclTimeInterval+randn(n,3)*M_STDDEV_VEL;
#pragma omp parallel sections
{
#pragma omp section
{
 F2D *t = fDeepCopy(accl);
#pragma omp parallel for
    for(i=0; i<(accl->height*accl->width); i++)
    {
            asubsref(t,i) = asubsref(accl,i) * acclTimeInterval;
    }

    is = fPlus(*vel, t);
    fFreeHandle(accl);
    fFreeHandle(t);
}
#pragma omp section
{

 F2D *s = randnWrapper(n,3);
#pragma omp parallel for
    for(i=0; i<(s->height*s->width); i++)
    {
            asubsref(s,i) = asubsref(s,i) * M_STDDEV_VEL;
    }

    fFreeHandle(*vel);
    *vel = fPlus(is, s);
    fFreeHandle(is);
    fFreeHandle(s);
}
}
}

我已经添加了一些openmp larallels但它仍然运行得很慢,所以我希望你是否经验丰富的老兵可以指出我应该关注哪些方面以提高性能。

2 个答案:

答案 0 :(得分:0)

现在,看起来只是将一些pragma粘贴在helter skelter中,而不关心代码的并行性质。现在,你并没有真正分开处理器上的工作,你实际上是通过让每个核心做同样的事情来重复工作(并且可能会产生一些严重的错误答案)。很高兴看到原版,看看你改变了什么。你实际上需要以编程方式分割任务(例如,如果你正在处理一个图像,你需要告诉,比方说,处理器1做奇数像素而处理器2做偶数像素,OpenMP不够聪明,无法解决这个问题)。我很抱歉,如果这似乎没有那么有用,但是因为这是为了一项任务,我认为你可能需要为自己确定一定的数额 - 但重要的是,OpenMP不能只是让代码并行,你需要告诉它如何。查看omp_get_num_threads()和omp_get_thread_num()。

答案 1 :(得分:0)

此代码开头: for(i=0; i<(t->height*t->width); i++) asubsref(t,i) = asubsref(t,i)/STDDEV_ACCL; 可以有效地并行化。每个i的计算完全独立。

寻找这些独立的代码部分。当然,您应该考虑要并行化的内容。