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>(gel_rand())/GEL_RAND_MAX - 0.5f;}

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

int 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;

    gel_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;

	return 0;

}