Subversion Repositories gelsvn

#include <iostream>

#include "CGLA/Vec3f.h"

#include "CGLA/Mat3x3f.h"

#include "CGLA/statistics.h"

#include "CGLA/eigensolution.h"

using namespace std;

using namespace CGLA;

float frand() { return static_cast<float>(rand())/RAND_MAX - 0.5f;}

/* This is a non-exhaustive test program for CGLA */

main()

{

	Vec3f norm(1,0,0);

	norm.normalize();

	Vec3f a,b;

	orthogonal(norm,a,b);

	cout << "The frame " << endl;

	cout << norm << a << b << endl;

	vector<Vec3f> vec;

	srand(0);

	for(int i=0;i<100000;++i)

	  {

	    Vec3f p = 1.0f*frand()*norm + 2.02f*frand()*a + 3.05f*frand()*b;

	    vec.push_back(p);

	  }

	Mat3x3f A(0);

	Vec3f m = covariance(vec, A);

	cout << "Mean and covariance " << endl;

	cout << m << "\n" << A << endl;

	Mat3x3f Q,L;

	int n = power_eigensolution(A, Q, L);

	cout << "The " << n << " eigensolutions are ";

	cout << Q << L << endl;

	cout << "Dot products " 

			 << dot(Q[0], Q[1]) << " "

			 << dot(Q[0], Q[2]) << " " 

			 << dot(Q[1], Q[2]) << endl;

}