tvmet - This is the Tiny Vector Matrix Expression Templat…

Users » madcowswe » Code » tvmet

This is the Tiny Vector Matrix Expression Templates library found at http://tvmet.sourceforge.net. It is the fastest and most compact matrix lib out there (for < 10x10 matricies). I have done some minor tweaks to make it compile for mbed. For examples and hints on how to use, see: http://tvmet.sourceforge.net/usage.html

Dependents: Eurobot_2012_Secondary

VectorImpl.h@0:feb4117d16d8, 2012-03-28 (annotated)

Committer:: madcowswe
Date:: Wed Mar 28 15:53:45 2012 +0000
Revision:: 0:feb4117d16d8

Who changed what in which revision?

User	Revision	Line number	New contents of line
madcowswe	0:feb4117d16d8	1	/*
madcowswe	0:feb4117d16d8	2	* Tiny Vector Matrix Library
madcowswe	0:feb4117d16d8	3	* Dense Vector Matrix Libary of Tiny size using Expression Templates
madcowswe	0:feb4117d16d8	4	*
madcowswe	0:feb4117d16d8	5	* Copyright (C) 2001 - 2007 Olaf Petzold <opetzold@users.sourceforge.net>
madcowswe	0:feb4117d16d8	6	*
madcowswe	0:feb4117d16d8	7	* This library is free software; you can redistribute it and/or
madcowswe	0:feb4117d16d8	8	* modify it under the terms of the GNU lesser General Public
madcowswe	0:feb4117d16d8	9	* License as published by the Free Software Foundation; either
madcowswe	0:feb4117d16d8	10	* version 2.1 of the License, or (at your option) any later version.
madcowswe	0:feb4117d16d8	11	*
madcowswe	0:feb4117d16d8	12	* This library is distributed in the hope that it will be useful,
madcowswe	0:feb4117d16d8	13	* but WITHOUT ANY WARRANTY; without even the implied warranty of
madcowswe	0:feb4117d16d8	14	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
madcowswe	0:feb4117d16d8	15	* lesser General Public License for more details.
madcowswe	0:feb4117d16d8	16	*
madcowswe	0:feb4117d16d8	17	* You should have received a copy of the GNU lesser General Public
madcowswe	0:feb4117d16d8	18	* License along with this library; if not, write to the Free Software
madcowswe	0:feb4117d16d8	19	* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
madcowswe	0:feb4117d16d8	20	*
madcowswe	0:feb4117d16d8	21	* $Id: VectorImpl.h,v 1.31 2007-06-23 15:58:58 opetzold Exp $
madcowswe	0:feb4117d16d8	22	*/
madcowswe	0:feb4117d16d8	23
madcowswe	0:feb4117d16d8	24	#ifndef TVMET_VECTOR_IMPL_H
madcowswe	0:feb4117d16d8	25	#define TVMET_VECTOR_IMPL_H
madcowswe	0:feb4117d16d8	26
madcowswe	0:feb4117d16d8	27	#include <iomanip> // setw
madcowswe	0:feb4117d16d8	28
madcowswe	0:feb4117d16d8	29	#include <tvmet/Functional.h>
madcowswe	0:feb4117d16d8	30	#include <tvmet/Io.h>
madcowswe	0:feb4117d16d8	31
madcowswe	0:feb4117d16d8	32
madcowswe	0:feb4117d16d8	33	namespace tvmet {
madcowswe	0:feb4117d16d8	34
madcowswe	0:feb4117d16d8	35
madcowswe	0:feb4117d16d8	36	/*
madcowswe	0:feb4117d16d8	37	* member operators for i/o
madcowswe	0:feb4117d16d8	38	*/
madcowswe	0:feb4117d16d8	39	template<class T, std::size_t Sz>
madcowswe	0:feb4117d16d8	40	std::ostream& Vector<T, Sz>::print_xpr(std::ostream& os, std::size_t l) const
madcowswe	0:feb4117d16d8	41	{
madcowswe	0:feb4117d16d8	42	os << IndentLevel(l++) << "Vector[" << ops << "]<"
madcowswe	0:feb4117d16d8	43	<< typeid(T).name() << ", " << Size << ">,"
madcowswe	0:feb4117d16d8	44	<< IndentLevel(--l)
madcowswe	0:feb4117d16d8	45	<< std::endl;
madcowswe	0:feb4117d16d8	46
madcowswe	0:feb4117d16d8	47	return os;
madcowswe	0:feb4117d16d8	48	}
madcowswe	0:feb4117d16d8	49
madcowswe	0:feb4117d16d8	50
madcowswe	0:feb4117d16d8	51	template<class T, std::size_t Sz>
madcowswe	0:feb4117d16d8	52	std::ostream& Vector<T, Sz>::print_on(std::ostream& os) const
madcowswe	0:feb4117d16d8	53	{
madcowswe	0:feb4117d16d8	54	enum {
madcowswe	0:feb4117d16d8	55	complex_type = NumericTraits<value_type>::is_complex
madcowswe	0:feb4117d16d8	56	};
madcowswe	0:feb4117d16d8	57
madcowswe	0:feb4117d16d8	58	std::streamsize w = IoPrintHelper<Vector>::width(dispatch<complex_type>(), *this);
madcowswe	0:feb4117d16d8	59
madcowswe	0:feb4117d16d8	60	os << std::setw(0) << "[\n ";
madcowswe	0:feb4117d16d8	61	for(std::size_t i = 0; i < (Size - 1); ++i) {
madcowswe	0:feb4117d16d8	62	os << std::setw(w) << m_data[i] << ", ";
madcowswe	0:feb4117d16d8	63	}
madcowswe	0:feb4117d16d8	64	os << std::setw(w) << m_data[Size - 1] << "\n]";
madcowswe	0:feb4117d16d8	65
madcowswe	0:feb4117d16d8	66	return os;
madcowswe	0:feb4117d16d8	67	}
madcowswe	0:feb4117d16d8	68
madcowswe	0:feb4117d16d8	69
madcowswe	0:feb4117d16d8	70	/*
madcowswe	0:feb4117d16d8	71	* member operators with scalars, per se element wise
madcowswe	0:feb4117d16d8	72	*/
madcowswe	0:feb4117d16d8	73	#define TVMET_IMPLEMENT_MACRO(NAME, OP) \
madcowswe	0:feb4117d16d8	74	template<class T, std::size_t Sz> \
madcowswe	0:feb4117d16d8	75	inline \
madcowswe	0:feb4117d16d8	76	Vector<T, Sz>& Vector<T, Sz>::operator OP (value_type rhs) { \
madcowswe	0:feb4117d16d8	77	typedef XprLiteral<value_type> expr_type; \
madcowswe	0:feb4117d16d8	78	this->M_##NAME(XprVector<expr_type, Size>(expr_type(rhs))); \
madcowswe	0:feb4117d16d8	79	return *this; \
madcowswe	0:feb4117d16d8	80	}
madcowswe	0:feb4117d16d8	81
madcowswe	0:feb4117d16d8	82	TVMET_IMPLEMENT_MACRO(add_eq, +=)
madcowswe	0:feb4117d16d8	83	TVMET_IMPLEMENT_MACRO(sub_eq, -=)
madcowswe	0:feb4117d16d8	84	TVMET_IMPLEMENT_MACRO(mul_eq, *=)
madcowswe	0:feb4117d16d8	85	TVMET_IMPLEMENT_MACRO(div_eq, /=)
madcowswe	0:feb4117d16d8	86	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	87
madcowswe	0:feb4117d16d8	88
madcowswe	0:feb4117d16d8	89	#define TVMET_IMPLEMENT_MACRO(NAME, OP) \
madcowswe	0:feb4117d16d8	90	template<class T, std::size_t Sz> \
madcowswe	0:feb4117d16d8	91	inline \
madcowswe	0:feb4117d16d8	92	Vector<T, Sz>& Vector<T, Sz>::operator OP (std::size_t rhs) { \
madcowswe	0:feb4117d16d8	93	typedef XprLiteral<value_type> expr_type; \
madcowswe	0:feb4117d16d8	94	this->M_##NAME(XprVector<expr_type, Size>(expr_type(rhs))); \
madcowswe	0:feb4117d16d8	95	return *this; \
madcowswe	0:feb4117d16d8	96	}
madcowswe	0:feb4117d16d8	97
madcowswe	0:feb4117d16d8	98	TVMET_IMPLEMENT_MACRO(mod_eq, %=)
madcowswe	0:feb4117d16d8	99	TVMET_IMPLEMENT_MACRO(xor_eq,^=)
madcowswe	0:feb4117d16d8	100	TVMET_IMPLEMENT_MACRO(and_eq, &=)
madcowswe	0:feb4117d16d8	101	TVMET_IMPLEMENT_MACRO(or_eq, \|=)
madcowswe	0:feb4117d16d8	102	TVMET_IMPLEMENT_MACRO(shl_eq, <<=)
madcowswe	0:feb4117d16d8	103	TVMET_IMPLEMENT_MACRO(shr_eq, >>=)
madcowswe	0:feb4117d16d8	104	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	105
madcowswe	0:feb4117d16d8	106
madcowswe	0:feb4117d16d8	107	/*
madcowswe	0:feb4117d16d8	108	* member functions (operators) with vectors, for use with +=,-= ... <<=
madcowswe	0:feb4117d16d8	109	*/
madcowswe	0:feb4117d16d8	110	#define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe	0:feb4117d16d8	111	template<class T1, std::size_t Sz> \
madcowswe	0:feb4117d16d8	112	template <class T2> \
madcowswe	0:feb4117d16d8	113	inline Vector<T1, Sz>& \
madcowswe	0:feb4117d16d8	114	Vector<T1, Sz>::M_##NAME (const Vector<T2, Size>& rhs) { \
madcowswe	0:feb4117d16d8	115	this->M_##NAME( XprVector<typename Vector<T2, Size>::ConstReference, Size>(rhs.const_ref()) ); \
madcowswe	0:feb4117d16d8	116	return *this; \
madcowswe	0:feb4117d16d8	117	}
madcowswe	0:feb4117d16d8	118
madcowswe	0:feb4117d16d8	119	TVMET_IMPLEMENT_MACRO(add_eq)
madcowswe	0:feb4117d16d8	120	TVMET_IMPLEMENT_MACRO(sub_eq)
madcowswe	0:feb4117d16d8	121	TVMET_IMPLEMENT_MACRO(mul_eq)
madcowswe	0:feb4117d16d8	122	TVMET_IMPLEMENT_MACRO(div_eq)
madcowswe	0:feb4117d16d8	123	TVMET_IMPLEMENT_MACRO(mod_eq)
madcowswe	0:feb4117d16d8	124	TVMET_IMPLEMENT_MACRO(xor_eq)
madcowswe	0:feb4117d16d8	125	TVMET_IMPLEMENT_MACRO(and_eq)
madcowswe	0:feb4117d16d8	126	TVMET_IMPLEMENT_MACRO(or_eq)
madcowswe	0:feb4117d16d8	127	TVMET_IMPLEMENT_MACRO(shl_eq)
madcowswe	0:feb4117d16d8	128	TVMET_IMPLEMENT_MACRO(shr_eq)
madcowswe	0:feb4117d16d8	129	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	130
madcowswe	0:feb4117d16d8	131
madcowswe	0:feb4117d16d8	132	/*
madcowswe	0:feb4117d16d8	133	* member functions (operators) with expressions, for use width +=,-= ... <<=
madcowswe	0:feb4117d16d8	134	*/
madcowswe	0:feb4117d16d8	135	#define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe	0:feb4117d16d8	136	template<class T, std::size_t Sz> \
madcowswe	0:feb4117d16d8	137	template <class E> \
madcowswe	0:feb4117d16d8	138	inline \
madcowswe	0:feb4117d16d8	139	Vector<T, Sz>& \
madcowswe	0:feb4117d16d8	140	Vector<T, Sz>::M_##NAME (const XprVector<E, Size>& rhs) { \
madcowswe	0:feb4117d16d8	141	rhs.assign_to(*this, Fcnl_##NAME<value_type, typename E::value_type>()); \
madcowswe	0:feb4117d16d8	142	return *this; \
madcowswe	0:feb4117d16d8	143	}
madcowswe	0:feb4117d16d8	144
madcowswe	0:feb4117d16d8	145	TVMET_IMPLEMENT_MACRO(add_eq)
madcowswe	0:feb4117d16d8	146	TVMET_IMPLEMENT_MACRO(sub_eq)
madcowswe	0:feb4117d16d8	147	TVMET_IMPLEMENT_MACRO(mul_eq)
madcowswe	0:feb4117d16d8	148	TVMET_IMPLEMENT_MACRO(div_eq)
madcowswe	0:feb4117d16d8	149	TVMET_IMPLEMENT_MACRO(mod_eq)
madcowswe	0:feb4117d16d8	150	TVMET_IMPLEMENT_MACRO(xor_eq)
madcowswe	0:feb4117d16d8	151	TVMET_IMPLEMENT_MACRO(and_eq)
madcowswe	0:feb4117d16d8	152	TVMET_IMPLEMENT_MACRO(or_eq)
madcowswe	0:feb4117d16d8	153	TVMET_IMPLEMENT_MACRO(shl_eq)
madcowswe	0:feb4117d16d8	154	TVMET_IMPLEMENT_MACRO(shr_eq)
madcowswe	0:feb4117d16d8	155	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	156
madcowswe	0:feb4117d16d8	157
madcowswe	0:feb4117d16d8	158	/*
madcowswe	0:feb4117d16d8	159	* aliased member functions (operators) with vectors,
madcowswe	0:feb4117d16d8	160	* for use with +=,-= ... <<=
madcowswe	0:feb4117d16d8	161	*/
madcowswe	0:feb4117d16d8	162	#define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe	0:feb4117d16d8	163	template<class T1, std::size_t Sz> \
madcowswe	0:feb4117d16d8	164	template <class T2> \
madcowswe	0:feb4117d16d8	165	inline \
madcowswe	0:feb4117d16d8	166	Vector<T1, Sz>& \
madcowswe	0:feb4117d16d8	167	Vector<T1, Sz>::alias_##NAME (const Vector<T2, Size>& rhs) { \
madcowswe	0:feb4117d16d8	168	this->alias_##NAME( XprVector<typename Vector<T2, Size>::ConstReference, Size>(rhs.const_ref()) ); \
madcowswe	0:feb4117d16d8	169	return *this; \
madcowswe	0:feb4117d16d8	170	}
madcowswe	0:feb4117d16d8	171
madcowswe	0:feb4117d16d8	172	TVMET_IMPLEMENT_MACRO(assign)
madcowswe	0:feb4117d16d8	173	TVMET_IMPLEMENT_MACRO(add_eq)
madcowswe	0:feb4117d16d8	174	TVMET_IMPLEMENT_MACRO(sub_eq)
madcowswe	0:feb4117d16d8	175	TVMET_IMPLEMENT_MACRO(mul_eq)
madcowswe	0:feb4117d16d8	176	TVMET_IMPLEMENT_MACRO(div_eq)
madcowswe	0:feb4117d16d8	177	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	178
madcowswe	0:feb4117d16d8	179
madcowswe	0:feb4117d16d8	180	/*
madcowswe	0:feb4117d16d8	181	* aliased member functions (operators) with expressions,
madcowswe	0:feb4117d16d8	182	* for use width +=,-= ... <<=
madcowswe	0:feb4117d16d8	183	*/
madcowswe	0:feb4117d16d8	184	#define TVMET_IMPLEMENT_MACRO(NAME) \
madcowswe	0:feb4117d16d8	185	template<class T, std::size_t Sz> \
madcowswe	0:feb4117d16d8	186	template <class E> \
madcowswe	0:feb4117d16d8	187	inline \
madcowswe	0:feb4117d16d8	188	Vector<T, Sz>& \
madcowswe	0:feb4117d16d8	189	Vector<T, Sz>::alias_##NAME (const XprVector<E, Size>& rhs) { \
madcowswe	0:feb4117d16d8	190	typedef Vector<T, Sz> temp_type; \
madcowswe	0:feb4117d16d8	191	temp_type(rhs).assign_to(*this, Fcnl_##NAME<value_type, typename E::value_type>()); \
madcowswe	0:feb4117d16d8	192	return *this; \
madcowswe	0:feb4117d16d8	193	}
madcowswe	0:feb4117d16d8	194
madcowswe	0:feb4117d16d8	195	TVMET_IMPLEMENT_MACRO(assign)
madcowswe	0:feb4117d16d8	196	TVMET_IMPLEMENT_MACRO(add_eq)
madcowswe	0:feb4117d16d8	197	TVMET_IMPLEMENT_MACRO(sub_eq)
madcowswe	0:feb4117d16d8	198	TVMET_IMPLEMENT_MACRO(mul_eq)
madcowswe	0:feb4117d16d8	199	TVMET_IMPLEMENT_MACRO(div_eq)
madcowswe	0:feb4117d16d8	200	#undef TVMET_IMPLEMENT_MACRO
madcowswe	0:feb4117d16d8	201
madcowswe	0:feb4117d16d8	202
madcowswe	0:feb4117d16d8	203	} // namespace tvmet
madcowswe	0:feb4117d16d8	204
madcowswe	0:feb4117d16d8	205	#endif // TVMET_VECTOR_IMPL_H
madcowswe	0:feb4117d16d8	206
madcowswe	0:feb4117d16d8	207	// Local Variables:
madcowswe	0:feb4117d16d8	208	// mode:C++
madcowswe	0:feb4117d16d8	209	// tab-width:8
madcowswe	0:feb4117d16d8	210	// End:

Repository toolbox

Export to desktop IDE

Repository details

Type:	Library
Created:	28 Mar 2012
Imports:	62
Forks:	0
Commits:	1
Dependents:	1
Dependencies:	0
Followers:	2

VectorImpl.h@0:feb4117d16d8, 2012-03-28 (annotated)

Who changed what in which revision?

Repository toolbox

Repository details

Important Information for this Arm website

Access Warning