0

我是新来的，而且我在编程方面也相对较新。我已经用 C 语言编写了一个程序，我需要使用 pthread 对其进行加速。我曾尝试使用 OpenMP 这样做，但我不知道如何调试它。此外，我需要确定程序是否使用 pthreads 和时间更快，但我不知道如何在我的代码中编写它。这是我的代码

enter code here
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <pthread.h>
#define NTHREADS 2
#define FYLLO(komvos) ((komvos) * 2 + 1)

long factorial(long);
void heap_function (int [], int, int );
void make_heap(long [], int );
void pop_heap(long [], int );

struct thread_data
{
long int    n;
long int k;
long *b;
};

main()
{

  long int n,k,c,fact=1;
  long *a,*b,*d,p[k];
  int i,j,rc;
  int q[]={2,3,4,5,6,7,8,9,12,13,14,15,16};
  pthread_t thread[NTHREADS];
  struct thread_data threada;
  for(i=0;i<NTHREADS;i++)
  {
    threada.n=n;
    threada.k=k;
    threada.b=b;
    pthread_create (&thread[i], NULL, (void *)&threada);
  }
  for (i=0; i<NTHREADS; i++) 
    rc = pthread_join (thread[i], NULL);

  for(i=0;i<13;i++)
  {
    k=pow(2,q[i])-1;
    if(a=(long*)malloc(i*sizeof(long))==NULL);
    {
       printf("Den yparxei diathesimi mnimi gia desmeusi\n");
       exit(1);
    }
    a[i]=k;
    for(a[0];a[13];a[i]++)
    {
        n=(pow(2,q[i]))*k;
        if(d=(long*)malloc((i*i)*sizeof(long))==NULL);
        {
            printf("Den yparxei diathesimi mnimi gia desmeusi\n");
            exit(1);
        }
        d[i]=n;
    }  
  c=(factorial(n))/((factorial(k))*(factorial(n-k)));     
  } 
  if(b=(long*)malloc(((i*i)+i)*sizeof(long))==NULL)
    {
        printf("Den yparxei diathesimi mnimi gia desmeusi\n");
        exit(1);
    }
  for(i=0;i<13;i++)
  {
    b[i]=a[i];
  }
  for(i=13;i<182;i++)                     /* Gia i=13 exoume i^2=169 kai i^2+i=182*/
  {
    b[i]=d[i];
  }    
  long heap[sizeof(b)];
  make_heap( heap, sizeof(b) );
  printf("To heap einai:\n");        
    for ( i = sizeof(b); i >=0; i-- ) 
        {
            printf( "%d ", heap[0] );
            pop_heap( heap, i );
        }
    for(i=(n-k);i<=n;i++)
        for(j=0;j<k;j++)
        {
            p[j]=heap[i];
            printf("Ta %d mikrotera stoixeia eina ta %ld\n",k,p[j]);
        } 
  free((void*)b);             
  getch();
}



long factorial(long n)
{
   int a;
   long result=1;

   for( a=1;a<=n;a++ )
     result=result*a;

   return(result);
}



void heap_function( int a[], int i, int n )
{

  while ( FYLLO( i ) < n )                  /* Vazoume sto heap ta stoixeia san           ypodentra */
{

  int fyllo = FYLLO( i );


  if ( fyllo + 1 < n && a[fyllo] < a[fyllo + 1] )     /* Dialegoume to maegalytero apo ta dyo paidia komvous */
  ++fyllo;


  if ( a[i] < a[fyllo] )                      /* Metaferoume to megalytero komvo sti riza */
  {
    int k = a[i];
    a[i] = a[fyllo];
    a[fyllo] = k;
  }


  ++i;                                    /* Synexizoume ston epomeno komvo */
  }
}


void make_heap( long a[], int n )        /*Dhmioyrgoume ti sinartisi make_heap gia na  mporesoume na valoume ta
                                        stoixeia pou dwsame mesa sto heap kai na ta ta3inomisoume*/
{
  int i = n / 2;

  while ( i-- > 0 )
  heap_function( a, i, n );
}


void pop_heap( long heap[], int n )      /*Dhmiourgoume ti sinartisi pop_heap gia na  mporesoume na e3agoume 
                                        ta stoixeia apo to heap apo to megalytero sto mikrotero*/
{
  long k = heap[0];
  heap[0] = heap[n];
  heap[n] = k;
  heap_function( heap, 0, n );          /*Afou emfanistei to prwto stoixeio kaloume ti  sinartisi heap_function 
                                        gia na ta3inomisei ta stoixeia pou menoun sto heap*/
}

对不起，我的邮件搞砸了，但我是新来的，现在我要开始使用它了

4

1 回答 1

5

添加线程可能不会加速您的程序，它可以让您将工作组织成可以并行运行的执行单元（并且在多核系统上，通常可以并行运行）。如果您不在多核系统上，如果您的一个或多个线程必须阻塞等待慢速输入，您仍然可以获得优势，因为其他线程可以继续运行；这可能会也可能不会给您更快的运行时间，具体取决于您的实际程序。

调试线程通常比调试单个线程更困难，如何做到这一点取决于您可用的工具。如果您的调试器无法让您的工作更轻松，我建议您首先让您的程序串行运行 - 仍然使用线程模型将其分解，但让每个代码在主线程中运行并让它运行直到完成，如果您的模型允许这样做。许多线程应用程序不能这样编写，因为线程在运行时相互依赖，但这仅取决于您正在做什么。

现在谈谈你的具体情况——当你还不知道如何游泳时，你正在潜入深渊。我建议您首先学习使用线程，而不必考虑为什么需要它们的复杂性，否则您会使问题变得比需要的更复杂。http://cs.gmu.edu/~white/CS571/Examples/Pthread/create.c有一个简单的例子可以开始。特别注意pthread_create()调用的参数并与您所做的进行比较；您的代码缺少第三个参数——作为线程运行的函数。您似乎根本没有这样的功能，相反，您似乎相信调用后的代码pthread_create()是并行运行的。这就是fork()工作原理，但这是非常不同的。

这应该足以让你开始。http://cs.gmu.edu/~white/CS571/Examples/pthread_examples.html有其他示例，“pthread 教程”的谷歌可能会有所帮助。

于 2012-03-29T11:52:25.430 回答