#include <algorithm>
#include <gsl/gsl_sort.h>
#include <gsl/gsl_vector.h>
#include <iostream>
#include <memory>
#include <numeric>
#include <vector>

int main()
{
	/* gsl vector */
	unsigned int vector_size = 3;

	std::unique_ptr<gsl_vector,void(*)(gsl_vector*)>
		v( gsl_vector_alloc( vector_size ), gsl_vector_free );

	for(unsigned int i = 0; i < vector_size; ++i)
		gsl_vector_set( v.get(), i, i );

	for(unsigned int i = 0; i < vector_size; ++i)
		std::cout << gsl_vector_get( v.get(), i ) << std::endl;

	std::cout << std::endl;

	/* c++11 vector */
	std::vector<double> vec = { 2, -1, 0 };

	gsl_sort( vec.data(), 1, vec.size() );

	for(double x : vec)
		std::cout << x << std::endl;

	std::cout << std::endl;

	/* using STL */
	std::cout << std::max( vec.front(), vec.back() ) << std::endl;

	return 0;
}