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;
}
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;
}