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