0

我正在从文件中读取浮点矩阵。矩阵尺寸为 4k * 4K。使用下面的程序，它只会导致 now() 函数奇怪地重置。如果我将矩阵大小减小到 1k * 1K，它不会重置。虽然它正确读取浮点数，但最后几个值是垃圾。我不知道这些垃圾值是从哪里来的。我选择了 BUFFSIZE 6 的大小，因为它在浮点数中的位数将在 5-6 左右。不确定是否正确。

#include <sys/time.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <sys/mman.h>
#include <fcntl.h>
#include <stdio.h>
#include <unistd.h>
#include <iostream>
#include <stdlib.h>
#include <sstream>


#define ROWS 4000
#define COLS 4000
#define BUFFSIZE 6

//#define USE_FREAD
#define USE_MMAP

double now()
{
    struct timeval tv;
    gettimeofday(&tv, NULL);
    return tv.tv_sec + tv.tv_usec / 1000000.;
}

int main()
{

    double end_time;
    double total_time;
    int i, x, y, k;
    for (k = 0; k < 1; k++)
    {
    double start_time = now();

    FILE* in = fopen("resistence_file", "rb");

    float arr[ROWS][COLS];

    char temp[BUFFSIZE];

    int val;
    std::stringstream ss;
    char* floats  = (char*)mmap(

            0,

            ROWS * COLS * sizeof(float),

            PROT_READ,

            MAP_FILE | MAP_PRIVATE,

            fileno(in),

            0

            );

    fclose(in);
    ss<<floats;

    for (int i =0; i < ROWS; i++)
    {
        for (int j = 0; j < COLS; j++)
        {
            if ((ss.getline(temp, BUFFSIZE, ' ')) )
            {
                arr[i][j] = atof((temp));
            }
        }
    }

    for (int i =0; i < ROWS; i++)
    {
        for (int j = 0; j < COLS; j++)
        {
            printf("%.1f ", arr[i][j]);
        }
        printf("\n");
    }
    munmap(floats, ROWS * COLS * sizeof(float));


    end_time = now();
    total_time = end_time - start_time;

    printf("It took %f seconds to read %d * %d matrix \n", total_time, ROWS, COLS);
    }

    return 0;
}
4

2 回答 2

6

float arr[4000][4000];将需要 56Mb（假设sizeof(float)=4）。这很可能比您的可用堆栈大。

您需要移动arr以获得静态持续时间

static float arr[ROWS][COLS];
...
int main()

或动态分配，free稍后记住

int main()
{
    float (*arr)[COLS] = malloc(sizeof(*arr) * ROWS);
    ....
    free(arr);
于 2013-06-18T15:29:23.500 回答
0
  1. 您正在映射一个sizeof(float)*K字节文件，就好像它包含二进制数据一样。如果它包含已知数量的二进制浮点数据，您为什么要稍后尝试使用它，就好像它包含文本一样？如果它包含文本，它与什么有关系sizeof(float)？您是否绝对确定您的浮点数平均占用text字节数，包括空格分隔符？sizeof(float)在大多数地方，每个数字大约有 3 个非空白字符。
  2. 您正在ss输入一个可能会或可能不会以 NUL 结尾的字符串。

与崩溃没有直接关系：这里的整个业务完全浪费了stringstream时空getline（好吧，至少浪费了 64 兆字节的空间和至少一些时间）。更不用说atof程序员不应该接触六英尺长的杆子了。我希望您不要认真计划在生产代码中这样做。只需std::strtod在原始数组上使用来代替所有这些。

于 2013-06-18T18:05:11.773 回答