Hey there. I'm pursuing my Matrix class and would like to make it work together with a Vector class I'm developing as well. I have a design question however and a technical question as well.

First, here's the code for the Vector class right now:
Code:
#ifndef VECTOR_H_INCLUDED
#define VECTOR_H_INCLUDED

#include <vector>

template < class T >
class Vector
{
public:
	Vector(int n = 3);
	Vector(const Vector< T >&);

	// Not quite sure about the type it should return...
	float Angle( ) const;
	void Normalize( );
	float Length( ) const;

	Vector< T > VecProduct(const Vector< T >&) const;
	T DotProduct(const Vector< T >&) const;

	bool IsParallel(const Vector< T >&) const;
	bool IsOrthogonal(const Vector< T >&) const;

	bool operator == (const Vector< T >&) const;
	bool operator != (const Vector< T >&) const;
	Vector< T >& operator = (const Vector< T >&);
	Vector< T >& operator += (const Vector< T >&);
	Vector< T >& operator + (const Vector< T >&) const;
	Vector< T >& operator -= (const Vector< T >&);
	Vector< T >& operator - (const Vector< T >&) const;
	Vector< T >& operator *= (const T&);
	Vector< T >& operator * (const T&) const;

private:
	
};

#endif // VECTOR_H_INCLUDED
It isn't as complete as it should, but this is a good start for a beginner's Vector class. I want to be able to use that Vector class with matrices; therefore I need it to let the user use a 2, 3, 5 or even 10 dimension vector. I first thought about using std::auto_ptr< > and allocate enough memory for each dimension, but copying objects would be a mess because of pointer ownership. Then I thought of using a vector, but it seems like overhead to me since vectors don't usually get large enough for vectors to be efficient. I thought then to use a pointer and allocate memory dynamically but this seems a bit hack-ish to me. Do you guys know a good way to solve that problem ?

Also, the technical question is rather simple, in fact. What type should the Angle() function return ? I thought returning a T would be dangerous because if T happens to be 'int' then we lose precision (which is also a problem for Normalize() and Length()). Then I thought returning float would be correct because of its precision but I'm not 100&#37; sure.

Thanks a lot.