从用numpy.save（...）保存的文件将numpy数组加载到C中

import numpy as np

arr = np.array([1, 2, 3, 4, 5])
np.save('my_array.npy', arr)

#include<stdio.h>
#include <stdlib.h>
#include <stdint.h>

int main() {
    FILE *file;
    char header[6];
    int32_t ndims, dims[32], type_num;
    size_t size;
    void *data;

    file = fopen("my_array.npy", "rb");
    if (file == NULL) {
        printf("Error: could not open file.\n");
        return 1;
    }

    // Read header
    fread(header, 1, 6, file);
    if (strncmp(header, "\x93NUMPY", 6) != 0) {
        printf("Error: not a valid NumPy file.\n");
        fclose(file);
        return 1;
    }

    // Read metadata
    fread(&ndims, 1, 1, file);
    fread(dims, 4, ndims, file);
    fread(&type_num, 1, 1, file);
    fseek(file, 1, SEEK_CUR);

    // Calculate size of data
    size = 1;
    for (int i = 0; i < ndims; i++) {
        size *= dims[i];
    }

    // Allocate memory for data
    data = malloc(size * sizeof(int32_t));
    if (data == NULL) {
        printf("Error: could not allocate memory.\n");
        fclose(file);
        return 1;
    }

    // Read data
    fread(data, size, sizeof(int32_t), file);

    // Print data
    for (int i = 0; i< size; i++) {
        printf("%d ", ((int32_t *)data)[i]);
    }
    printf("\n");

    // Clean up
    free(data);
    fclose(file);

    return 0;
}