#include <stdio.h>
#include <stdlib.h>
#include <gsl/gsl_spmatrix.h>
#include <gsl/gsl_version.h>

int main(void) {
 size_t v = 1;
 v <<= 32;
 v += 1;
 printf("gsl header version : %s\n", GSL_VERSION);
 printf("gsl library version: %s\n", gsl_version);
 printf("size of _size_t is %d bytes\n", sizeof(size_t));
 if(v == 1) printf("WTF\n");
 gsl_spmatrix *m = gsl_spmatrix_alloc(1,1);
 gsl_spmatrix_set( m, 1, 0, 1.5 );
 printf("1: got: %f\n", gsl_spmatrix_get( m, 1, 0));
 printf("1: nnz: %d\n", gsl_spmatrix_nnz( m ));
 gsl_spmatrix_set( m, v, 0, 1.75 );
 if (v >= m->size1) printf("will be reported as 'out of range'\n");
 printf("2: got: %f\n", gsl_spmatrix_get( m, v, 0 ));
 printf("2: nnz: %d\n", gsl_spmatrix_nnz( m ));
 return 0;
}