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