UPC выделение динамического массива и сортировка - PullRequest
Новая
       31

UPC выделение динамического массива и сортировка

0 голосов
/ 27 января 2012

Я пытаюсь прочитать поплавки из файла, а затем отсортировать их. Сортировка должна быть параллельной UPC. Это текущий код: 

#include <upc_relaxed.h>
#include <upc_collective.h>
#include <stdio.h>
#include <stdlib.h>

int lt_int( shared void *x, shared void *y ) {
  int x_val = *(shared int *)x, 
    y_val = *(shared int *)y;
  return x_val > y_val ? -1 : x_val < y_val ? 1 : 0;
}

shared int size=0;

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


  FILE *f;
  int i=0;
  if (MYTHREAD == 0) {
    f = fopen ("dane.dat", "r");
    while (feof(f) == 0) {        
      fscanf (f, "%f\n");                
      ++size;    
    }
    fclose(f);
  }

  upc_barrier;

  /* allocation goes wrong! */
  shared [] float *array = upc_all_alloc(size, sizeof(float));
  /* printf("%d\n",sizeof(array)); // it returns 8! */ 

  upc_barrier;

  if (MYTHREAD == 0) {
    f = fopen ("dane.dat", "r");
    i=0;   
    while (feof(f) == 0) {
      printf("%d\n", i);
      /* segmentation fault! */
      fscanf (f, "%f\n", &array[i]);    
      printf("%f\n", array[i]);            
      i++;    
    }
    fclose(f);
  }

  upc_barrier;
  upc_all_sort( array, sizeof(float), size/THREADS, size, lt_int, UPC_IN_ALLSYNC);
  upc_barrier;  

  if (MYTHREAD == 0) {
    for (i = 0; i<=atoi(argv[1]) ; ++i) {
      printf("%f\n", array[atoi(argv[1]) + (size/atoi(argv[1]))]);
    }
  }

  return 0;
}

И я не знаю, что я делаю не так. Я получаю ошибку сегментации, потому что распределение памяти происходит неправильно. Вы можете мне помочь?

1 Ответ

1 голос
/ 27 января 2012

Звонок неправильный: 

fscanf (f, "%f\n");

Также ваш объект array является указателем на float. Обычно sizeof array возвращает размер типа указателя (8 в вашей реализации), а не размер выделенного вами объекта массива. Вам следует проверить возвращаемое значение upc_all_alloc, чтобы убедиться, что во время распределения не было ошибок (если возвращаемое значение == NULL, распределение не удалось).

...