codes-workload-mpi-replay.c 15.1 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
/*
 * Copyright (C) 2013 University of Chicago.
 * See COPYRIGHT notice in top-level directory.
 *
 */

/* SUMMARY:
 *
 *  MPI replay tool for replaying workloads from the codes workload API. 
 *
 */

#include <stdio.h>
#include <unistd.h>
#include <assert.h>
#include <inttypes.h>
#include <math.h>
#include <getopt.h>
#include <fcntl.h>
20
#include <sys/stat.h>
21 22 23 24
#include <mpi.h>

#include "codes/codes-workload.h"
#include "codes/quickhash.h"
25
#include "codes/configuration.h"
26

27
/* hash table entry for looking up file descriptor of a workload file id */
28
struct file_info
29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44
{
    struct qlist_head hash_link;
    uint64_t file_hash;
    int file_descriptor;
};

int replay_workload_op(struct codes_workload_op replay_op, int rank, long long int op_number);
int hash_file_compare(void *key, struct qlist_head *link);

/* command line options */
static int opt_verbose = 0;
static int opt_noop = 0;

/* hash table for storing file descriptors of opened files */
static struct qhash_table *fd_table = NULL;

45 46 47
/* file stream to log rank events to, if verbose turned on */
static FILE *log_stream = NULL;

48 49
void usage(char *exename)
{
50 51
    fprintf(stderr, "Usage: %s [OPTIONS] --conf <conf_file_path>\n       "
            "--test-dir <workload_test_dir>\n\n", exename);
52
    fprintf(stderr, "\t<conf_file_path> : (absolute) path to a valid workload configuration file\n");
53 54
    fprintf(stderr, "\t<workload_test_dir> : the directory to replay the workload I/O in\n");
    fprintf(stderr, "\n\t[OPTIONS] includes:\n");
55 56 57 58 59 60
    fprintf(stderr, "\t\t--noop : do not perform i/o\n");
    fprintf(stderr, "\t\t    -v : verbose (output i/o details)\n");   

    exit(1);
}

61
void parse_args(int argc, char **argv, char **conf_path, char **test_dir)
62 63 64 65
{
    int index;
    static struct option long_opts[] =
    {
66
        {"conf", 1, NULL, 'c'},
67
        {"test-dir", 1, NULL, 'd'},
68
        {"noop", 0, NULL, 'n'},
69 70 71 72
        {"help", 0, NULL, 0},
        {0, 0, 0, 0}
    };

73
    *conf_path = NULL;
74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89
    *test_dir = NULL;
    while (1)
    {
        int c = getopt_long(argc, argv, "v", long_opts, &index);

        if (c == -1)
            break;

        switch (c)
        {
            case 'v':
                opt_verbose = 1;
                break;
            case 'n':
                opt_noop = 1;
                break;
90 91
            case 'c':
                *conf_path = optarg;
92 93 94 95 96 97 98 99 100 101 102 103
                break;
            case 'd':
                *test_dir = optarg;
                break;
            case 0:
            case '?':
            default:
                usage(argv[0]);
                break;
        }
    }

104
    if (optind < argc || !(*conf_path) || !(*test_dir))
105 106 107 108 109 110 111
    {
        usage(argv[0]);
    }

    return;
}

112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138
int load_workload(char *conf_path, int rank)
{
    config_lpgroups_t paramconf;
    char workload_type[MAX_NAME_LENGTH_WKLD];

    /* load the config file across all ranks */
    configuration_load(conf_path, MPI_COMM_WORLD, &config);

    /* get the PARAMS group out of the config file */
    configuration_get_lpgroups(&config, "PARAMS", &paramconf);

    /* get the workload type out of PARAMS */
    configuration_get_value(&config, "PARAMS", "workload_type",
                            workload_type, MAX_NAME_LENGTH_WKLD);

    /* set up the workload parameters and load into the workload API */
    if (strcmp(workload_type, "darshan_io_workload") == 0)
    {
        struct darshan_params d_params;
        char aggregator_count[10];

        /* get the darshan params from the config file */
        configuration_get_value(&config, "PARAMS", "log_file_path",
                                d_params.log_file_path, MAX_NAME_LENGTH_WKLD);
        configuration_get_value(&config, "PARAMS", "aggregator_count", aggregator_count, 10);
        d_params.aggregator_cnt = atoi(aggregator_count);

139 140
        d_params.stream = NULL;
#if 0
141
        d_params.stream = log_stream;
142 143 144
        opt_verbose = 0;
#endif

145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170
        return codes_workload_load(workload_type, (char *)&d_params, rank);
    }
    else if (strcmp(workload_type, "bgp_io_workload") == 0)
    {
        struct bgp_params b_params;
        char rank_count[10];

        /* get the bgp i/o params from the config file */
        configuration_get_value(&config, "PARAMS", "io_kernel_meta_path",
                                b_params.io_kernel_meta_path, MAX_NAME_LENGTH_WKLD);
        configuration_get_value(&config, "PARAMS", "bgp_config_file", 
                                b_params.bgp_config_file, MAX_NAME_LENGTH_WKLD);
        configuration_get_value(&config, "PARAMS", "rank_count", rank_count, 10);
        strcpy(b_params.io_kernel_path, "");
        strcpy(b_params.io_kernel_def_path, "");
        b_params.num_cns = atoi(rank_count);

        return codes_workload_load(workload_type, (char *)&b_params, rank);
    }
    else
    {
        fprintf(stderr, "Error: Invalid workload type specified (%s)\n", workload_type);
        return -1;
    }
}

171 172
int main(int argc, char *argv[])
{
173
    char *conf_path;
174
    char *replay_test_path;
175 176
    char *log_dir = "log";
    char my_log_path[MAX_NAME_LENGTH_WKLD];
177 178 179 180 181 182 183 184
    int nprocs;
    int myrank;
    int workload_id;
    struct codes_workload_op next_op;
    long long int replay_op_number = 1;
    int ret = 0;

    /* parse command line args */
185
    parse_args(argc, argv, &conf_path, &replay_test_path);
186 187 188 189 190 191

    /* initialize MPI */
    MPI_Init(&argc, &argv);
    MPI_Comm_size(MPI_COMM_WORLD, &nprocs);
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank);

192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212
    /* change the working directory to be the test directory */
    ret = chdir(replay_test_path);
    if (ret < 0)
    {
        fprintf(stderr, "Unable to change to testing directory (%s)\n", strerror(errno));
        goto error_exit;
    }

    /* set the path for logging this rank's events, if verbose is turned on */
    if (opt_verbose)
    {
        mkdir(log_dir, 0755);
        snprintf(my_log_path, MAX_NAME_LENGTH_WKLD, "%s/rank-%d.log", log_dir, myrank);
        log_stream = fopen(my_log_path, "w");
        if (log_stream == NULL)
        {
            fprintf(stderr, "Unable to open log file %s\n", my_log_path);
            goto error_exit;
        }
    }

213 214
    /* initialize workload generator from config file */
    workload_id = load_workload(conf_path, myrank);
215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249
    if (workload_id < 0)
    {
        goto error_exit;
    }

    /* initialize hash table for storing file descriptors */
    fd_table = qhash_init(hash_file_compare, quickhash_64bit_hash, 29);
    if (!fd_table)
    {
        fprintf(stderr, "File descriptor hash table memory error\n");
        goto error_exit;
    }

    /* replay loop */
    while (1)
    {
        /* get the next replay operation from the workload generator */
        codes_workload_get_next(workload_id, myrank, &next_op);

        if (next_op.op_type != CODES_WK_END)
        {
            /* replay the next workload operation */
            ret = replay_workload_op(next_op, myrank, replay_op_number++);
            if (ret < 0)
            {
                break;
            }
        }
        else
        {
            /* workload replay for this rank is complete */
            break;
        }
    }

250
    if (log_stream)
251 252
        fclose(log_stream);

253
    /* destroy and finalize the file descriptor hash table */
254
    qhash_destroy_and_finalize(fd_table, struct file_info, hash_link, free);
255 256 257

error_exit:
    MPI_Finalize();
258

259 260 261 262 263 264 265 266 267 268
    return ret;
}

int replay_workload_op(struct codes_workload_op replay_op, int rank, long long int op_number)
{
    unsigned int secs;
    unsigned int usecs;
    int open_flags = O_RDWR;
    char file_name[50];
    int fildes;
269
    struct file_info *tmp_list = NULL;
270 271 272 273 274 275 276 277
    struct qlist_head *hash_link = NULL;
    char *buf = NULL;
    int ret;

    switch (replay_op.op_type)
    {
        case CODES_WK_DELAY:
            if (opt_verbose)
278
                fprintf(log_stream, "[Rank %d] Operation %lld : DELAY %lf seconds\n",
279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306
                       rank, op_number, replay_op.u.delay.seconds);

            if (!opt_noop)
            {
                /* satisfy delay using second delay then microsecond delay */
                secs = floor(replay_op.u.delay.seconds);
                usecs = round((replay_op.u.delay.seconds - secs) * 1000 * 1000);
                ret = sleep(secs);
                if (ret)
                {
                    /* error in sleep */
                    errno = EINTR;
                    fprintf(stderr, "Rank %d failure on operation %lld [DELAY: %s]\n",
                            rank, op_number, strerror(errno));
                    return -1;
                }
                ret = usleep(usecs);
                if (ret < 0)
                {
                    /* error in usleep */
                    fprintf(stderr, "Rank %d failure on operation %lld [DELAY: %s]\n",
                            rank, op_number, strerror(errno));
                    return -1;
                }
            }
            return 0;
        case CODES_WK_BARRIER:
            if (opt_verbose)
307
                fprintf(log_stream, "[Rank %d] Operation %lld : BARRIER\n", rank, op_number);
308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323

            if (!opt_noop)
            {
                /* implement barrier using MPI global barrier on all ranks */
                ret = MPI_Barrier(MPI_COMM_WORLD);
                if (ret != MPI_SUCCESS)
                {
                    /* error in MPI_Barrier */
                    fprintf(stderr, "Rank %d failure on operation %lld [BARRIER: %s]\n",
                            rank, op_number, "Invalid communicator");
                    return -1;
                }
            }
            return 0;
        case CODES_WK_OPEN:
            if (opt_verbose)
324
                fprintf(log_stream, "[Rank %d] Operation %lld: %s file %"PRIu64"\n", rank, op_number,
325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346
                       (replay_op.u.open.create_flag) ? "CREATE" : "OPEN", replay_op.u.open.file_id);

            if (!opt_noop)
            {
                /* set the create flag, if necessary */
                if (replay_op.u.open.create_flag)
                    open_flags |= O_CREAT;

                /* write the file hash to string to be used as the actual file name */
                snprintf(file_name, sizeof(file_name), "%"PRIu64, replay_op.u.open.file_id);

                /* perform the open operation */
                fildes = open(file_name, open_flags, 0666);
                if (fildes < 0)
                {
                    fprintf(stderr, "Rank %d failure on operation %lld [%s: %s]\n",
                            rank, op_number, (replay_op.u.open.create_flag) ? "CREATE" : "OPEN",
                            strerror(errno));
                    return -1;
                }

                /* save the file descriptor for this file in a hash table to be retrieved later */
347
                tmp_list = malloc(sizeof(struct file_info));
348 349 350 351 352 353 354 355 356 357 358 359 360
                if (!tmp_list)
                {
                    fprintf(stderr, "No memory available for file hash entry\n");
                    return -1;
                }

                tmp_list->file_hash = replay_op.u.open.file_id;
                tmp_list->file_descriptor = fildes;
                qhash_add(fd_table, &(replay_op.u.open.file_id), &(tmp_list->hash_link));
            }
            return 0;
        case CODES_WK_CLOSE:
            if (opt_verbose)
361 362
                fprintf(log_stream, "[Rank %d] Operation %lld : CLOSE file %"PRIu64"\n",
                        rank, op_number, replay_op.u.close.file_id);
363 364 365 366 367 368

            if (!opt_noop)
            {
                /* search for the corresponding file descriptor in the hash table */
                hash_link = qhash_search_and_remove(fd_table, &(replay_op.u.close.file_id));
                assert(hash_link);
369
                tmp_list = qhash_entry(hash_link, struct file_info, hash_link);
370 371 372 373 374 375 376 377 378 379 380 381 382 383 384
                fildes = tmp_list->file_descriptor;
                free(tmp_list);
     
                /* perform the close operation */
                ret = close(fildes);
                if (ret < 0)
                {
                    fprintf(stderr, "Rank %d failure on operation %lld [CLOSE: %s]\n",
                            rank, op_number, strerror(errno));
                    return -1;
                }
            }
            return 0;
        case CODES_WK_WRITE:
            if (opt_verbose)
385
                fprintf(log_stream, "[Rank %d] Operation %lld : WRITE file %"PRIu64" (sz = %"PRId64
386 387 388 389 390 391 392 393 394
                       ", off = %"PRId64")\n",
                       rank, op_number, replay_op.u.write.file_id, replay_op.u.write.size,    
                       replay_op.u.write.offset);

            if (!opt_noop)
            {
                /* search for the corresponding file descriptor in the hash table */
                hash_link = qhash_search(fd_table, &(replay_op.u.write.file_id));
                assert(hash_link);
395
                tmp_list = qhash_entry(hash_link, struct file_info, hash_link);
396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416
                fildes = tmp_list->file_descriptor;

                /* perform the write operation */
                buf = malloc(replay_op.u.write.size);
                if (!buf)
                {
                    fprintf(stderr, "No memory available for write buffer\n");
                    return -1;
                }
                ret = pwrite(fildes, buf, replay_op.u.write.size, replay_op.u.write.offset);
                free(buf);
                if (ret < 0)
                {
                    fprintf(stderr, "Rank %d failure on operation %lld [WRITE: %s]\n",
                            rank, op_number, strerror(errno));
                    return -1;
                }
            }
            return 0;
        case CODES_WK_READ:
            if (opt_verbose)
417 418
                fprintf(log_stream, "[Rank %d] Operation %lld : READ file %"PRIu64" (sz = %"PRId64
                        ", off = %"PRId64")\n", rank, op_number, replay_op.u.read.file_id,
419 420 421 422 423 424 425
                       replay_op.u.read.size, replay_op.u.read.offset);

            if (!opt_noop)
            {
                /* search for the corresponding file descriptor in the hash table */
                hash_link = qhash_search(fd_table, &(replay_op.u.read.file_id));
                assert(hash_link);
426
                tmp_list = qhash_entry(hash_link, struct file_info, hash_link);
427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447
                fildes = tmp_list->file_descriptor;

                /* perform the write operation */
                buf = malloc(replay_op.u.read.size);
                if (!buf)
                {
                    fprintf(stderr, "No memory available for write buffer\n");
                    return -1;
                }
                ret = pread(fildes, buf, replay_op.u.read.size, replay_op.u.read.offset);
                free(buf);
                if (ret < 0)
                {
                    fprintf(stderr, "Rank %d failure on operation %lld [READ: %s]\n",
                            rank, op_number, strerror(errno));
                    return -1;
                }
            }
            return 0;
        default:
            fprintf(stderr, "** Rank %d: INVALID OPERATION (op count = %lld) **\n", rank, op_number);
448
            return 0;
449 450 451 452 453 454 455
    }

}

int hash_file_compare(void *key, struct qlist_head *link)
{
    uint64_t *in_file_hash = (uint64_t *)key;
456
    struct file_info *tmp_file;
457

458
    tmp_file = qlist_entry(link, struct file_info, hash_link);
459 460 461 462 463
    if (tmp_file->file_hash == *in_file_hash)
        return 1;

    return 0;
}