tests/geometry/test-linalg.cpp
author Aravindakshan Babu <akshan@stanford.edu>
Thu, 08 Jul 2010 23:50:39 -0700
branchdev
changeset 217 64a272a34cb2
parent 55 7e71f5984931
permissions -rw-r--r--
Added extra functionality to Point class( an iterator ) and PersistenceDiagram( dimension property and __len__ func ). persistence-diagram.h: Added a new read-only dimension member and member function to access it. With a new constructor that that takes in an int type to initialize dimension. persistence-diagram.cpp: Added new bp::init constructor. Takes in an integer type to initialize the dimension. Exposed the dimension property. Exposed the size property via a __len__ method. __init__.py: Added an iterator for Point objects. This iterates over the coords and then the data( if present ).

#include "geometry/linalg.h"
#include <iostream>

#include <synaps/upol.h>
#include <synaps/upol/gcd.h>
#include "geometry/rational-function.h"


typedef			UPolDse<double>					Polynomial;
typedef			RationalFunction<Polynomial>	RationalF;

//typedef 		LinearAlgebra<double>			LinAlg;
//typedef 		LinearAlgebra<Polynomial>		LinAlg;
typedef 		LinearAlgebra<RationalF>		LinAlg;

int main()
{
	LinAlg::MatrixType  a(2,2);
	a(0,0) = Polynomial("3*x"); a(1,0) = Polynomial(4); a(0,1) = Polynomial(0); a(1,1) = Polynomial(7);

	std::cout << a << std::endl;
	std::cout << LinAlg::determinant(a) << std::endl;
	std::cout << Polynomial("3*x^2 + 4*x") / Polynomial("3*x^2") << std::endl;

	LinAlg::VectorType b(2), x;
	b(0) = Polynomial(4); b(1) = Polynomial(3);
	LinAlg::solve(a, b, x);
	std::cout << x << std::endl;

	x(0).normalize(); x(1).normalize();
	std::cout << x << std::endl;
}