Subversion Repositories gelsvn

#include <algorithm>

#include "CGLA/Mat4x4d.h"

	cout << "Determinant of singular 4x4 matrix" << endl;

 	Mat4x4d m_odd(Vec4d(0, 373, 139129, 1),

 								Vec4d(-373, -125, 154754, 1),

 								Vec4d(-50.6168, 21.2469, 3013.5, 100),

 								Vec4d(-50.6168, 21.2469, 3013.5, 100));

	cout << m_odd << endl;

	cout << "Determinant: " << determinant(m_odd) << endl;

	Mat4x4d m_odd_t = transpose(m_odd);

	cout << m_odd_t << endl;

	cout << "Det of transpose:" << determinant(m_odd_t) << endl;

	cout << "Determinant of random matrix (MatLab says 0.0491)" << endl;

	Mat4x4d Mrand(Vec4d(0.9169,0.3529,0.2028,0.1988),

								Vec4d(0.4103,0.8132,0.1987,0.0153),

								Vec4d(0.8936,0.0099,0.6038,0.7468),

								Vec4d(0.0579,0.1389,0.2722,0.4451));

	cout << "det " << determinant(Mrand) << endl;

	cout << "det of transpose " << determinant(transpose(Mrand)) << endl;