#include <gsl/gsl_matrix.h>
#include <gsl/gsl_linalg.h>

int main(void)
{
  gsl_matrix *A,*V;
  gsl_vector *S,*work;

  A=gsl_matrix_calloc(3,3);
  V=gsl_matrix_calloc(3,3);
  S=gsl_vector_calloc(3);
  work=gsl_vector_calloc(3);

  gsl_matrix_set(A,0,0,0);
  gsl_matrix_set(A,0,1,1);
  gsl_matrix_set(A,0,2,0);

  gsl_matrix_set(A,1,0,0);
  gsl_matrix_set(A,1,1,2);
  gsl_matrix_set(A,1,2,1);

  gsl_matrix_set(A,2,0,0);
  gsl_matrix_set(A,2,1,-1);
  gsl_matrix_set(A,2,2,-1);

  printf("gets to here\n");
  /*gsl_linalg_SV_decomp(A,V,S,work);*/
  gsl_linalg_SV_decomp_jacobi(A,V,S);
  printf("doesn't get to here\n");

  return(0);
}