int main(int argc, char **argv)
{
  program_name = basename(argv[0]);

  /* Initialize library */
  spx_init();
  spx_log_info_console();

  if (argc < 2) {
    fprintf(stderr, "%s: too few arguments\n", program_name);
    print_usage();
    exit(1);
  }

  /* Load matrix from binary file */
  spx_matrix_t *A = spx_mat_restore(argv[1]);

  /* Change some values in the matrix */
  spx_index_t row = 1, col = 1;
  spx_value_t value = 0.42;
  spx_mat_set_entry(A, row, col, value);
  /* ... */

  /* Create random x and y vectors */
  spx_partition_t *parts = spx_mat_get_partition(A);
  spx_vector_t *x = spx_vec_create_random(spx_mat_get_ncols(A), parts);
  spx_vector_t *y = spx_vec_create_random(spx_mat_get_nrows(A), parts);

  /* Run 128 loops of the SpMV kernel */
  spx_value_t alpha = 0.8, beta = 0.42;
  const size_t nr_loops = 128;
  spx_timer_t t;
  double elapsed_time, flops;
  size_t i;

  spx_timer_clear(&t);
  spx_timer_start(&t);
  for (i = 0; i < nr_loops; i++) {
    spx_matvec_kernel(alpha, A, x, beta, y);
  }
  spx_timer_pause(&t);

  elapsed_time = spx_timer_get_secs(&t);
  flops = (double) (2 * nr_loops * spx_mat_get_nnz(A)) /
    ((double) 1000 * 1000 * elapsed_time);
  printf("SPMV time: %lf secs\n", elapsed_time);
  printf("MFLOPS: %lf\n", flops);

  /* Cleanup */
  spx_mat_destroy(A);
  spx_partition_destroy(parts);
  spx_vec_destroy(x);
  spx_vec_destroy(y);

  /* Shutdown library */
  spx_finalize();

  return 0;
}
Esempio n. 2
0
int main(int argc, char **argv)
{
    /* Initialize library */
    spx_init();
    spx_log_info_console();

    /* Load matrix from MMF file */
    spx_input_t *input = spx_input_load_mmf(argv[1]);

    /* Set tuning options */
    spx_option_set("spx.rt.nr_threads", "2");
    spx_option_set("spx.rt.cpu_affinity", "0,1");
    spx_option_set("spx.preproc.xform", "all");
    spx_option_set("spx.preproc.sampling", "portion");
    spx_option_set("spx.preproc.sampling.nr_samples", "48");
    spx_option_set("spx.preproc.sampling.portion", "0.01");
    /* spx_option_set("spx.matrix.symmetric", "true"); */

    /* Transform to CSX */
    spx_matrix_t *A = spx_mat_tune(input);

    /* Create random x and y vectors */
    spx_partition_t *parts = spx_mat_get_partition(A);
    spx_vector_t *x = spx_vec_create_random(spx_mat_get_ncols(A), parts);
    spx_vector_t *y = spx_vec_create_random(spx_mat_get_nrows(A), parts);

    /* Run 128 loops of the SpMV kernel */
    spx_value_t alpha = 0.8, beta = 0.42;
    const size_t nr_loops = 128;
    spx_timer_t t;
    double elapsed_time, flops;
    size_t i;

    spx_timer_clear(&t);
    spx_timer_start(&t);
    for (i = 0; i < nr_loops; i++) {
        spx_matvec_kernel(alpha, A, x, beta, y);
    }
    spx_timer_pause(&t);

    elapsed_time = spx_timer_get_secs(&t);
    flops = (double) (2 * nr_loops * spx_mat_get_nnz(A)) /
        ((double) 1000 * 1000 * elapsed_time);
    printf("Elapsed time: %lf secs\n", elapsed_time);
    printf("FLOPS: %lf\n", flops);

    /* Cleanup */
    spx_input_destroy(input);
    spx_mat_destroy(A);
    spx_partition_destroy(parts);
    spx_vec_destroy(x);
    spx_vec_destroy(y);

    /* Shutdown library */
    spx_finalize();

    return 0;
}
int main(int argc, char **argv)
{
    program_name = basename(argv[0]);

    /* Initialize library */
    spx_init();
    spx_log_info_console();

    if (argc < 3) {
        fprintf(stderr, "%s: too few arguments\n", program_name);
        print_usage();
        exit(1);
    }

    /* Load matrix from MMF file */
    spx_input_t *input = spx_input_load_mmf(argv[1]);

    /* Set tuning options */
    spx_option_set("spx.rt.nr_threads", "2");
    spx_option_set("spx.rt.cpu_affinity", "0,1");

    /* Transform to CSX */
    spx_matrix_t *A = spx_mat_tune(input);

    /* Save matrix in a binary file for  future session */
    spx_mat_save(A, argv[2]);

    /* Create random x and y vectors */
    spx_partition_t *parts = spx_mat_get_partition(A);
    spx_vector_t *x = spx_vec_create_random(spx_mat_get_ncols(A), parts);
    spx_vector_t *y = spx_vec_create_random(spx_mat_get_nrows(A), parts);

    /* Run 128 loops of the SpMV kernel */
    spx_value_t alpha = 0.8, beta = 0.42;
    const size_t nr_loops = 128;
    spx_timer_t t;
    double elapsed_time, flops;
    size_t i;

    spx_timer_clear(&t);
    spx_timer_start(&t);
    for (i = 0; i < nr_loops; i++) {
        spx_matvec_kernel(alpha, A, x, beta, y);
    }
    spx_timer_pause(&t);

    elapsed_time = spx_timer_get_secs(&t);
    flops = (double) (2 * nr_loops * spx_mat_get_nnz(A)) /
        ((double) 1000 * 1000 * elapsed_time);
    printf("SPMV time: %lf secs\n", elapsed_time);
    printf("MFLOPS: %lf\n", flops);

    /* Cleanup */
    spx_input_destroy(input);
    spx_mat_destroy(A);
    spx_partition_destroy(parts);
    spx_vec_destroy(x);
    spx_vec_destroy(y);

    /* Shutdown library */
    spx_finalize();

    return 0;
}