dot.h		dot.h
			// -*- C++ -*-
	/************************************************************************** *		/************************************************************************** *
	* blitz/meta/dot.h Tiny vector dot product metaprogram		* blitz/meta/dot.h Tiny vector dot product metaprogram
	*		*
	* $Id: dot.h,v 1.2 2001/01/24 20:22:51 tveldhui Exp $		* $Id: dot.h,v 1.5 2005/05/07 04:17:57 julianc Exp $
	*		*
	* Copyright (C) 1997-2001 Todd Veldhuizen <tveldhui@oonumerics.org>		* Copyright (C) 1997-2001 Todd Veldhuizen <tveldhui@oonumerics.org>
	*		*
	* This program is free software; you can redistribute it and/or		* This program is free software; you can redistribute it and/or
	* modify it under the terms of the GNU General Public License		* modify it under the terms of the GNU General Public License
	* as published by the Free Software Foundation; either version 2		* as published by the Free Software Foundation; either version 2
	* of the License, or (at your option) any later version.		* of the License, or (at your option) any later version.
	*		*
	* This program is distributed in the hope that it will be useful,		* This program is distributed in the hope that it will be useful,
	* but WITHOUT ANY WARRANTY; without even the implied warranty of		* but WITHOUT ANY WARRANTY; without even the implied warranty of
	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the		* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	* GNU General Public License for more details.		* GNU General Public License for more details.
	*		*
	* Suggestions: blitz-dev@oonumerics.org		* Suggestions: blitz-dev@oonumerics.org
	* Bugs: blitz-bugs@oonumerics.org		* Bugs: blitz-bugs@oonumerics.org
	*		*
	* For more information, please see the Blitz++ Home Page:		* For more information, please see the Blitz++ Home Page:
	* http://oonumerics.org/blitz/		* http://oonumerics.org/blitz/
	*		*
	**************************************************************************		**************************************************************************
	*		*/
	* $Log: dot.h,v $
	* Revision 1.2 2001/01/24 20:22:51 tveldhui
	* Updated copyright date in headers.
	*
	* Revision 1.1.1.1 2000/06/19 12:26:13 tveldhui
	* Imported sources
	*
	* Revision 1.2 1998/03/14 00:08:44 tveldhui
	* 0.2-alpha-05
	*
	* Revision 1.1 1997/07/16 14:51:20 tveldhui
	* Update: Alpha release 0.2 (Arrays)
	*
	*/
	#ifndef BZ_META_DOT_H		#ifndef BZ_META_DOT_H
	#define BZ_META_DOT_H		#define BZ_META_DOT_H
	#ifndef BZ_PROMOTE_H		#ifndef BZ_PROMOTE_H
	#include <blitz/promote.h>		#include <blitz/promote.h>
	#endif		#endif
	#ifndef BZ_METAPROG_H		#ifndef BZ_METAPROG_H
	#include <blitz/meta/metaprog.h>		#include <blitz/meta/metaprog.h>
	#endif		#endif
	BZ_NAMESPACE(blitz)		BZ_NAMESPACE(blitz)
	template<int N, int I>		template<int N, int I>
	class _bz_meta_vectorDot {		class _bz_meta_vectorDot {
	public:		public:
	enum { loopFlag = (I < N-1) ? 1 : 0 };		static const int loopFlag = (I < N-1) ? 1 : 0;
	template<class T_expr1, class T_expr2>		template<typename T_expr1, typename T_expr2>
	static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,		static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,
	_bz_typename T_expr2::T_numtype)		_bz_typename T_expr2::T_numtype)
	f(const T_expr1& a, const T_expr2& b)		f(const T_expr1& a, const T_expr2& b)
	{		{
	return a[I] * b[I]		return a[I] * b[I]
	+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);		+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);
	}		}
	template<class T_expr1, class T_expr2>		template<typename T_expr1, typename T_expr2>
	static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,		static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,
	_bz_typename T_expr2::T_numtype)		_bz_typename T_expr2::T_numtype)
	f_value_ref(T_expr1 a, const T_expr2& b)		f_value_ref(T_expr1 a, const T_expr2& b)
	{		{
	return a[I] * b[I]		return a[I] * b[I]
	+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);		+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);
	}		}
	template<class T_expr1, class T_expr2>		template<typename T_expr1, typename T_expr2>
	static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,		static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,
	_bz_typename T_expr2::T_numtype)		_bz_typename T_expr2::T_numtype)
	f_ref_value(const T_expr1& a, T_expr2 b)		f_ref_value(const T_expr1& a, T_expr2 b)
	{		{
	return a[I] * b[I]		return a[I] * b[I]
	+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);		+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::f(a,b);
	}		}
	template<class T_expr1, class P_numtype2>		template<typename T_expr1, typename P_numtype2>
	static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,		static inline BZ_PROMOTE(_bz_typename T_expr1::T_numtype,
	P_numtype2)		P_numtype2)
	dotWithArgs(const T_expr1& a, P_numtype2 i1, P_numtype2 i2=0,		dotWithArgs(const T_expr1& a, P_numtype2 i1, P_numtype2 i2=0,
	P_numtype2 i3=0, P_numtype2 i4=0, P_numtype2 i5=0, P_numtype2 i6=0,		P_numtype2 i3=0, P_numtype2 i4=0, P_numtype2 i5=0, P_numtype2 i6=0,
	P_numtype2 i7=0, P_numtype2 i8=0, P_numtype2 i9=0, P_numtype2 i10=0 )		P_numtype2 i7=0, P_numtype2 i8=0, P_numtype2 i9=0, P_numtype2 i10=0 )
	{		{
	return a[I] * i1		return a[I] * i1
	+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::dotWithAr gs		+ _bz_meta_vectorDot<loopFlag * N, loopFlag * (I+1)>::dotWithAr gs
	(a, i2, i3, i4, i5, i6, i7, i8, i9);		(a, i2, i3, i4, i5, i6, i7, i8, i9);
	}		}
	};		};
	template<>		template<>
	class _bz_meta_vectorDot<0,0> {		class _bz_meta_vectorDot<0,0> {
	public:		public:
	template<class T_expr1, class T_expr2>		template<typename T_expr1, typename T_expr2>
	static inline _bz_meta_nullOperand f(const T_expr1&, const T_expr2&)		static inline _bz_meta_nullOperand f(const T_expr1&, const T_expr2&)
	{ return _bz_meta_nullOperand(); }		{ return _bz_meta_nullOperand(); }
	template<class T_expr1, class P_numtype2>		template<typename T_expr1, typename P_numtype2>
	static inline _bz_meta_nullOperand		static inline _bz_meta_nullOperand
	dotWithArgs(const T_expr1& a, P_numtype2 i1, P_numtype2 i2=0,		dotWithArgs(const T_expr1& a, P_numtype2 i1, P_numtype2 i2=0,
	P_numtype2 i3=0, P_numtype2 i4=0, P_numtype2 i5=0, P_numtype2 i6=0,		P_numtype2 i3=0, P_numtype2 i4=0, P_numtype2 i5=0, P_numtype2 i6=0,
	P_numtype2 i7=0, P_numtype2 i8=0, P_numtype2 i9=0, P_numtype2 i10=0 )		P_numtype2 i7=0, P_numtype2 i8=0, P_numtype2 i9=0, P_numtype2 i10=0 )
	{		{
	return _bz_meta_nullOperand();		return _bz_meta_nullOperand();
	}		}
	};		};
End of changes. 10 change blocks.
	24 lines changed or deleted		11 lines changed or added
This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/