   if (solver_handle) cusolverDnDestroy(solver_handle);
 
   cuda::device_reset();
-
-  return 0;
 }

    -1,  0,  5, 2,
     0,  1,  2, 5,
  };

  std::cout << "Input matrix:\n";
  print_matrix(matrix, n);

#include <cusolverDn.h>

#include <cudalicious/cudalicious.h>

constexpr auto K_MAX = 50;

template<typename T>
void print_matrix(const std::vector<T>& matrix, const int n)
