summaryrefslogtreecommitdiffstats
path: root/examples/h5dsm_dset_r1m.c
blob: ff92bcb2db97b82f23dafc7c2d98f438e6cd12d2 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
#include "h5dsm_example.h"
#include <time.h>

int main(int argc, char *argv[]) {
    uuid_t pool_uuid;
    char *pool_grp = NULL;
    hid_t file = -1, dset = -1, fapl = -1;
    hid_t nfile = -1, ndset = -1;
    hsize_t dims[1] = {256 * 1024};
    int *buf = NULL, *nbuf = NULL;
    H5VL_daosm_snap_id_t snap_id;
    int i;

    (void)MPI_Init(&argc, &argv);

    /* Seed random number generator */
    srand(time(NULL));

    if(argc < 4 || argc > 5)
        PRINTF_ERROR("argc must be 4 or 5\n");

    /* Parse UUID */
    if(0 != uuid_parse(argv[1], pool_uuid))
        ERROR;

    /* Initialize VOL */
    if(H5VLdaosm_init(MPI_COMM_WORLD, pool_uuid, pool_grp) < 0)
        ERROR;

    /* Set up FAPL */
    if((fapl = H5Pcreate(H5P_FILE_ACCESS)) < 0)
        ERROR;
    if(H5Pset_fapl_daosm(fapl, MPI_COMM_WORLD, MPI_INFO_NULL) < 0)
        ERROR;
    if(H5Pset_all_coll_metadata_ops(fapl, true) < 0)
        ERROR;

    /* Open snapshot if specified */
    if(argc == 5) {
        snap_id = (H5VL_daosm_snap_id_t)atoi(argv[4]);
        printf("Opening snapshot %llu\n", (long long unsigned)snap_id);
        if(H5Pset_daosm_snap_open(fapl, snap_id) < 0)
            ERROR;
    } /* end if */

    /* Open file */
    if((file = H5Fopen(argv[2], H5F_ACC_RDWR, fapl)) < 0)
        ERROR;

    /* Open native file */
    if((nfile = H5Fopen(argv[2], H5F_ACC_RDWR, H5P_DEFAULT)) < 0)
        ERROR;

    /* Open dataset */
    if((dset = H5Dopen2(file, argv[3], H5P_DEFAULT)) < 0)
        ERROR;

    /* Open native dataset */
    if((ndset = H5Dopen2(nfile, argv[3], H5P_DEFAULT)) < 0)
        ERROR;

    /* Allocate buffers */
    if(NULL == (buf = (int *)malloc(dims[0] * sizeof(int))))
        ERROR;
    if(NULL == (nbuf = (int *)malloc(dims[0] * sizeof(int))))
        ERROR;

    /* Read data */
    printf("Reading daos-m dataset\n");
    if(H5Dread(dset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, buf) < 0)
        ERROR;
    
    /* Read native data */
    printf("Reading native dataset\n");
    if(H5Dread(ndset, H5T_NATIVE_INT, H5S_ALL, H5S_ALL, H5P_DEFAULT, nbuf) < 0)
        ERROR;

    /* Compare data */
    if(memcmp(buf, nbuf, dims[0] * sizeof(int)))
        PRINTF_ERROR("Buffers differ\n");
    else
        printf("Buffers are equal\n");

    /* Close */
    if(H5Dclose(dset) < 0)
        ERROR;
    if(H5Dclose(ndset) < 0)
        ERROR;
    if(H5Fclose(file) < 0)
        ERROR;
    if(H5Fclose(nfile) < 0)
        ERROR;
    if(H5Pclose(fapl) < 0)
        ERROR;
    free(buf);
    buf = NULL;
    free(nbuf);
    nbuf = NULL;

    printf("Success\n");

    (void)MPI_Finalize();
    return 0;

error:
    H5E_BEGIN_TRY {
        H5Dclose(dset);
        H5Dclose(ndset);
        H5Fclose(file);
        H5Fclose(nfile);
        H5Pclose(fapl);
    } H5E_END_TRY;

    free(buf);
    free(nbuf);

    (void)MPI_Finalize();
    return 1;
}