00001 #ifndef POINT_H
00002 #define POINT_H
00003
00004 #include <boost/numeric/ublas/matrix.hpp>
00005 #include <boost/numeric/ublas/vector.hpp>
00006 #include <boost/numeric/ublas/io.hpp>
00007 #include <boost/numeric/ublas/lu.hpp>
00008 #include <vector>
00009 #include <stdlib.h>
00010 #include <stdio.h>
00011
00012 #include "ads.h"
00013 #include "LogControl.h"
00014
00015 #define verbose 1
00016 #define TOL 1e-5
00017
00018 typedef double DataType;
00019 typedef boost::numeric::ublas::vector<DataType> BPType;
00020 typedef ads::FixedArray<DIM,DataType> PType;
00021 typedef ads::FixedArray<DIM,int> multi_index_type;
00022
00023 typedef boost::numeric::ublas::matrix<DataType> MType;
00024 typedef boost::numeric::ublas::identity_matrix<std::size_t> imatrix;
00025
00026 const DataType d2r = std::asin(DataType(1.))/DataType(90.);
00027 const DataType pi = DataType(2.)*std::asin(DataType(1.));
00028
00029 static std::ostream* MKlog;
00030 static inline std::ostream& mklog() { return *MKlog; }
00031
00032 #endif // POINT_H