Generates a unitary random vector in D dimensions

gistfile1.c

std::vector<double> unitaryVector(unsigned int dimension)
{
    std::vector<double> vector(dimension);
	double norm = 0;
	for (unsigned int d = 0; d < dimension; d++)
	{
		vector[d] = grandom();  // RNG call for a gaussian with mean=0 and sigma=1
		norm += vector[d]*vector[d];
	}
	norm = sqrt(norm);
	for (size_t d = 0; d < dimension; d++)
	{
		vector[d] /= norm;
	}
	return vector;
}