Context Navigation

source: XMLIO_V2/external/include/blitz/matuops.h @ 80

Last change on this file since 80 was 80, checked in by ymipsl, 14 years ago
ajout lib externe
Property svn:eol-style set to `native`
File size: 147.0 KB

Line
1	// Generated source file. Do not edit.
2	// Created by: genmatuops.cpp Dec 10 2003 17:58:05
3
4	#ifndef BZ_MATUOPS_H
5	#define BZ_MATUOPS_H
6
7	BZ_NAMESPACE(blitz)
8
9	#ifndef BZ_MATEXPR_H
10	#error <blitz/matuops.h> must be included via <blitz/matexpr.h>
11	#endif
12
13	/****************************************************************************
14	* abs
15	****************************************************************************/
16
17	template<class P_numtype1, class P_struct1>
18	inline
19	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
20	_bz_abs<P_numtype1> > >
21	abs(const Matrix<P_numtype1, P_struct1>& d1)
22	{
23	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
24	_bz_abs<P_numtype1> > T_expr;
25
26	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
27	}
28
29	template<class P_expr1>
30	inline
31	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
32	_bz_abs<typename P_expr1::T_numtype> > >
33	abs(_bz_MatExpr<P_expr1> d1)
34	{
35	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
36	_bz_abs<typename P_expr1::T_numtype> > T_expr;
37
38	return _bz_MatExpr<T_expr>(T_expr(d1));
39	}
40
41
42	/****************************************************************************
43	* acos
44	****************************************************************************/
45
46	template<class P_numtype1, class P_struct1>
47	inline
48	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
49	_bz_acos<P_numtype1> > >
50	acos(const Matrix<P_numtype1, P_struct1>& d1)
51	{
52	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
53	_bz_acos<P_numtype1> > T_expr;
54
55	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
56	}
57
58	template<class P_expr1>
59	inline
60	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
61	_bz_acos<typename P_expr1::T_numtype> > >
62	acos(_bz_MatExpr<P_expr1> d1)
63	{
64	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
65	_bz_acos<typename P_expr1::T_numtype> > T_expr;
66
67	return _bz_MatExpr<T_expr>(T_expr(d1));
68	}
69
70
71	/****************************************************************************
72	* acosh
73	****************************************************************************/
74
75	#ifdef BZ_HAVE_IEEE_MATH
76	template<class P_numtype1, class P_struct1>
77	inline
78	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
79	_bz_acosh<P_numtype1> > >
80	acosh(const Matrix<P_numtype1, P_struct1>& d1)
81	{
82	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
83	_bz_acosh<P_numtype1> > T_expr;
84
85	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
86	}
87
88	template<class P_expr1>
89	inline
90	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
91	_bz_acosh<typename P_expr1::T_numtype> > >
92	acosh(_bz_MatExpr<P_expr1> d1)
93	{
94	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
95	_bz_acosh<typename P_expr1::T_numtype> > T_expr;
96
97	return _bz_MatExpr<T_expr>(T_expr(d1));
98	}
99
100	#endif
101
102	/****************************************************************************
103	* asin
104	****************************************************************************/
105
106	template<class P_numtype1, class P_struct1>
107	inline
108	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
109	_bz_asin<P_numtype1> > >
110	asin(const Matrix<P_numtype1, P_struct1>& d1)
111	{
112	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
113	_bz_asin<P_numtype1> > T_expr;
114
115	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
116	}
117
118	template<class P_expr1>
119	inline
120	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
121	_bz_asin<typename P_expr1::T_numtype> > >
122	asin(_bz_MatExpr<P_expr1> d1)
123	{
124	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
125	_bz_asin<typename P_expr1::T_numtype> > T_expr;
126
127	return _bz_MatExpr<T_expr>(T_expr(d1));
128	}
129
130
131	/****************************************************************************
132	* asinh
133	****************************************************************************/
134
135	#ifdef BZ_HAVE_IEEE_MATH
136	template<class P_numtype1, class P_struct1>
137	inline
138	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
139	_bz_asinh<P_numtype1> > >
140	asinh(const Matrix<P_numtype1, P_struct1>& d1)
141	{
142	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
143	_bz_asinh<P_numtype1> > T_expr;
144
145	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
146	}
147
148	template<class P_expr1>
149	inline
150	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
151	_bz_asinh<typename P_expr1::T_numtype> > >
152	asinh(_bz_MatExpr<P_expr1> d1)
153	{
154	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
155	_bz_asinh<typename P_expr1::T_numtype> > T_expr;
156
157	return _bz_MatExpr<T_expr>(T_expr(d1));
158	}
159
160	#endif
161
162	/****************************************************************************
163	* atan
164	****************************************************************************/
165
166	template<class P_numtype1, class P_struct1>
167	inline
168	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
169	_bz_atan<P_numtype1> > >
170	atan(const Matrix<P_numtype1, P_struct1>& d1)
171	{
172	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
173	_bz_atan<P_numtype1> > T_expr;
174
175	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
176	}
177
178	template<class P_expr1>
179	inline
180	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
181	_bz_atan<typename P_expr1::T_numtype> > >
182	atan(_bz_MatExpr<P_expr1> d1)
183	{
184	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
185	_bz_atan<typename P_expr1::T_numtype> > T_expr;
186
187	return _bz_MatExpr<T_expr>(T_expr(d1));
188	}
189
190
191	/****************************************************************************
192	* atan2
193	****************************************************************************/
194
195	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
196	inline
197	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
198	_bz_atan2<P_numtype1,P_numtype2> > >
199	atan2(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
200	{
201	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
202	_bz_atan2<P_numtype1,P_numtype2> > T_expr;
203
204	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
205	}
206
207	template<class P_numtype1, class P_struct1, class P_expr2>
208	inline
209	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
210	_bz_atan2<P_numtype1,typename P_expr2::T_numtype> > >
211	atan2(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
212	{
213	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
214	_bz_atan2<P_numtype1,typename P_expr2::T_numtype> > T_expr;
215
216	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
217	}
218
219	template<class P_numtype1, class P_struct1>
220	inline
221	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
222	_bz_atan2<P_numtype1,int> > >
223	atan2(const Matrix<P_numtype1, P_struct1>& d1, int d2)
224	{
225	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
226	_bz_atan2<P_numtype1,int> > T_expr;
227
228	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
229	}
230
231	template<class P_numtype1, class P_struct1>
232	inline
233	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
234	_bz_atan2<P_numtype1,float> > >
235	atan2(const Matrix<P_numtype1, P_struct1>& d1, float d2)
236	{
237	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
238	_bz_atan2<P_numtype1,float> > T_expr;
239
240	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
241	}
242
243	template<class P_numtype1, class P_struct1>
244	inline
245	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
246	_bz_atan2<P_numtype1,double> > >
247	atan2(const Matrix<P_numtype1, P_struct1>& d1, double d2)
248	{
249	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
250	_bz_atan2<P_numtype1,double> > T_expr;
251
252	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
253	}
254
255	template<class P_numtype1, class P_struct1>
256	inline
257	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
258	_bz_atan2<P_numtype1,long double> > >
259	atan2(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
260	{
261	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
262	_bz_atan2<P_numtype1,long double> > T_expr;
263
264	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
265	}
266
267	template<class P_numtype1, class P_struct1, class T2>
268	inline
269	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
270	_bz_atan2<P_numtype1,complex<T2> > > >
271	atan2(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
272	{
273	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
274	_bz_atan2<P_numtype1,complex<T2> > > T_expr;
275
276	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
277	}
278
279	template<class P_expr1, class P_numtype2, class P_struct2>
280	inline
281	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
282	_bz_atan2<typename P_expr1::T_numtype,P_numtype2> > >
283	atan2(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
284	{
285	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
286	_bz_atan2<typename P_expr1::T_numtype,P_numtype2> > T_expr;
287
288	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
289	}
290
291	template<class P_expr1, class P_expr2>
292	inline
293	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
294	_bz_atan2<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
295	atan2(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
296	{
297	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
298	_bz_atan2<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
299
300	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
301	}
302
303	template<class P_expr1>
304	inline
305	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
306	_bz_atan2<typename P_expr1::T_numtype,int> > >
307	atan2(_bz_MatExpr<P_expr1> d1, int d2)
308	{
309	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
310	_bz_atan2<typename P_expr1::T_numtype,int> > T_expr;
311
312	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
313	}
314
315	template<class P_expr1>
316	inline
317	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
318	_bz_atan2<typename P_expr1::T_numtype,float> > >
319	atan2(_bz_MatExpr<P_expr1> d1, float d2)
320	{
321	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
322	_bz_atan2<typename P_expr1::T_numtype,float> > T_expr;
323
324	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
325	}
326
327	template<class P_expr1>
328	inline
329	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
330	_bz_atan2<typename P_expr1::T_numtype,double> > >
331	atan2(_bz_MatExpr<P_expr1> d1, double d2)
332	{
333	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
334	_bz_atan2<typename P_expr1::T_numtype,double> > T_expr;
335
336	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
337	}
338
339	template<class P_expr1>
340	inline
341	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
342	_bz_atan2<typename P_expr1::T_numtype,long double> > >
343	atan2(_bz_MatExpr<P_expr1> d1, long double d2)
344	{
345	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
346	_bz_atan2<typename P_expr1::T_numtype,long double> > T_expr;
347
348	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
349	}
350
351	template<class P_expr1, class T2>
352	inline
353	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
354	_bz_atan2<typename P_expr1::T_numtype,complex<T2> > > >
355	atan2(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
356	{
357	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
358	_bz_atan2<typename P_expr1::T_numtype,complex<T2> > > T_expr;
359
360	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
361	}
362
363	template<class P_numtype2, class P_struct2>
364	inline
365	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
366	_bz_atan2<int,P_numtype2> > >
367	atan2(int d1, const Matrix<P_numtype2, P_struct2>& d2)
368	{
369	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
370	_bz_atan2<int,P_numtype2> > T_expr;
371
372	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
373	}
374
375	template<class P_expr2>
376	inline
377	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
378	_bz_atan2<int,typename P_expr2::T_numtype> > >
379	atan2(int d1, _bz_MatExpr<P_expr2> d2)
380	{
381	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
382	_bz_atan2<int,typename P_expr2::T_numtype> > T_expr;
383
384	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
385	}
386
387	template<class P_numtype2, class P_struct2>
388	inline
389	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
390	_bz_atan2<float,P_numtype2> > >
391	atan2(float d1, const Matrix<P_numtype2, P_struct2>& d2)
392	{
393	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
394	_bz_atan2<float,P_numtype2> > T_expr;
395
396	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
397	}
398
399	template<class P_expr2>
400	inline
401	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
402	_bz_atan2<float,typename P_expr2::T_numtype> > >
403	atan2(float d1, _bz_MatExpr<P_expr2> d2)
404	{
405	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
406	_bz_atan2<float,typename P_expr2::T_numtype> > T_expr;
407
408	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
409	}
410
411	template<class P_numtype2, class P_struct2>
412	inline
413	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
414	_bz_atan2<double,P_numtype2> > >
415	atan2(double d1, const Matrix<P_numtype2, P_struct2>& d2)
416	{
417	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
418	_bz_atan2<double,P_numtype2> > T_expr;
419
420	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
421	}
422
423	template<class P_expr2>
424	inline
425	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
426	_bz_atan2<double,typename P_expr2::T_numtype> > >
427	atan2(double d1, _bz_MatExpr<P_expr2> d2)
428	{
429	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
430	_bz_atan2<double,typename P_expr2::T_numtype> > T_expr;
431
432	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
433	}
434
435	template<class P_numtype2, class P_struct2>
436	inline
437	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
438	_bz_atan2<long double,P_numtype2> > >
439	atan2(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
440	{
441	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
442	_bz_atan2<long double,P_numtype2> > T_expr;
443
444	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
445	}
446
447	template<class P_expr2>
448	inline
449	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
450	_bz_atan2<long double,typename P_expr2::T_numtype> > >
451	atan2(long double d1, _bz_MatExpr<P_expr2> d2)
452	{
453	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
454	_bz_atan2<long double,typename P_expr2::T_numtype> > T_expr;
455
456	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
457	}
458
459	template<class T1, class P_numtype2, class P_struct2>
460	inline
461	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
462	_bz_atan2<complex<T1> ,P_numtype2> > >
463	atan2(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
464	{
465	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
466	_bz_atan2<complex<T1> ,P_numtype2> > T_expr;
467
468	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
469	}
470
471	template<class T1, class P_expr2>
472	inline
473	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
474	_bz_atan2<complex<T1> ,typename P_expr2::T_numtype> > >
475	atan2(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
476	{
477	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
478	_bz_atan2<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
479
480	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
481	}
482
483
484	/****************************************************************************
485	* atanh
486	****************************************************************************/
487
488	#ifdef BZ_HAVE_IEEE_MATH
489	template<class P_numtype1, class P_struct1>
490	inline
491	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
492	_bz_atanh<P_numtype1> > >
493	atanh(const Matrix<P_numtype1, P_struct1>& d1)
494	{
495	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
496	_bz_atanh<P_numtype1> > T_expr;
497
498	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
499	}
500
501	template<class P_expr1>
502	inline
503	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
504	_bz_atanh<typename P_expr1::T_numtype> > >
505	atanh(_bz_MatExpr<P_expr1> d1)
506	{
507	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
508	_bz_atanh<typename P_expr1::T_numtype> > T_expr;
509
510	return _bz_MatExpr<T_expr>(T_expr(d1));
511	}
512
513	#endif
514
515	/****************************************************************************
516	* _class
517	****************************************************************************/
518
519	#ifdef BZ_HAVE_SYSTEM_V_MATH
520	template<class P_numtype1, class P_struct1>
521	inline
522	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
523	_bz__class<P_numtype1> > >
524	_class(const Matrix<P_numtype1, P_struct1>& d1)
525	{
526	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
527	_bz__class<P_numtype1> > T_expr;
528
529	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
530	}
531
532	template<class P_expr1>
533	inline
534	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
535	_bz__class<typename P_expr1::T_numtype> > >
536	_class(_bz_MatExpr<P_expr1> d1)
537	{
538	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
539	_bz__class<typename P_expr1::T_numtype> > T_expr;
540
541	return _bz_MatExpr<T_expr>(T_expr(d1));
542	}
543
544	#endif
545
546	/****************************************************************************
547	* cbrt
548	****************************************************************************/
549
550	#ifdef BZ_HAVE_IEEE_MATH
551	template<class P_numtype1, class P_struct1>
552	inline
553	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
554	_bz_cbrt<P_numtype1> > >
555	cbrt(const Matrix<P_numtype1, P_struct1>& d1)
556	{
557	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
558	_bz_cbrt<P_numtype1> > T_expr;
559
560	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
561	}
562
563	template<class P_expr1>
564	inline
565	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
566	_bz_cbrt<typename P_expr1::T_numtype> > >
567	cbrt(_bz_MatExpr<P_expr1> d1)
568	{
569	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
570	_bz_cbrt<typename P_expr1::T_numtype> > T_expr;
571
572	return _bz_MatExpr<T_expr>(T_expr(d1));
573	}
574
575	#endif
576
577	/****************************************************************************
578	* ceil
579	****************************************************************************/
580
581	template<class P_numtype1, class P_struct1>
582	inline
583	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
584	_bz_ceil<P_numtype1> > >
585	ceil(const Matrix<P_numtype1, P_struct1>& d1)
586	{
587	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
588	_bz_ceil<P_numtype1> > T_expr;
589
590	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
591	}
592
593	template<class P_expr1>
594	inline
595	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
596	_bz_ceil<typename P_expr1::T_numtype> > >
597	ceil(_bz_MatExpr<P_expr1> d1)
598	{
599	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
600	_bz_ceil<typename P_expr1::T_numtype> > T_expr;
601
602	return _bz_MatExpr<T_expr>(T_expr(d1));
603	}
604
605
606	/****************************************************************************
607	* copysign
608	****************************************************************************/
609
610	#ifdef BZ_HAVE_SYSTEM_V_MATH
611	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
612	inline
613	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
614	_bz_copysign<P_numtype1,P_numtype2> > >
615	copysign(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
616	{
617	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
618	_bz_copysign<P_numtype1,P_numtype2> > T_expr;
619
620	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
621	}
622
623	template<class P_numtype1, class P_struct1, class P_expr2>
624	inline
625	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
626	_bz_copysign<P_numtype1,typename P_expr2::T_numtype> > >
627	copysign(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
628	{
629	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
630	_bz_copysign<P_numtype1,typename P_expr2::T_numtype> > T_expr;
631
632	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
633	}
634
635	template<class P_numtype1, class P_struct1>
636	inline
637	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
638	_bz_copysign<P_numtype1,int> > >
639	copysign(const Matrix<P_numtype1, P_struct1>& d1, int d2)
640	{
641	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
642	_bz_copysign<P_numtype1,int> > T_expr;
643
644	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
645	}
646
647	template<class P_numtype1, class P_struct1>
648	inline
649	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
650	_bz_copysign<P_numtype1,float> > >
651	copysign(const Matrix<P_numtype1, P_struct1>& d1, float d2)
652	{
653	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
654	_bz_copysign<P_numtype1,float> > T_expr;
655
656	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
657	}
658
659	template<class P_numtype1, class P_struct1>
660	inline
661	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
662	_bz_copysign<P_numtype1,double> > >
663	copysign(const Matrix<P_numtype1, P_struct1>& d1, double d2)
664	{
665	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
666	_bz_copysign<P_numtype1,double> > T_expr;
667
668	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
669	}
670
671	template<class P_numtype1, class P_struct1>
672	inline
673	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
674	_bz_copysign<P_numtype1,long double> > >
675	copysign(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
676	{
677	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
678	_bz_copysign<P_numtype1,long double> > T_expr;
679
680	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
681	}
682
683	template<class P_numtype1, class P_struct1, class T2>
684	inline
685	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
686	_bz_copysign<P_numtype1,complex<T2> > > >
687	copysign(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
688	{
689	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
690	_bz_copysign<P_numtype1,complex<T2> > > T_expr;
691
692	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
693	}
694
695	template<class P_expr1, class P_numtype2, class P_struct2>
696	inline
697	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
698	_bz_copysign<typename P_expr1::T_numtype,P_numtype2> > >
699	copysign(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
700	{
701	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
702	_bz_copysign<typename P_expr1::T_numtype,P_numtype2> > T_expr;
703
704	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
705	}
706
707	template<class P_expr1, class P_expr2>
708	inline
709	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
710	_bz_copysign<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
711	copysign(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
712	{
713	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
714	_bz_copysign<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
715
716	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
717	}
718
719	template<class P_expr1>
720	inline
721	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
722	_bz_copysign<typename P_expr1::T_numtype,int> > >
723	copysign(_bz_MatExpr<P_expr1> d1, int d2)
724	{
725	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
726	_bz_copysign<typename P_expr1::T_numtype,int> > T_expr;
727
728	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
729	}
730
731	template<class P_expr1>
732	inline
733	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
734	_bz_copysign<typename P_expr1::T_numtype,float> > >
735	copysign(_bz_MatExpr<P_expr1> d1, float d2)
736	{
737	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
738	_bz_copysign<typename P_expr1::T_numtype,float> > T_expr;
739
740	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
741	}
742
743	template<class P_expr1>
744	inline
745	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
746	_bz_copysign<typename P_expr1::T_numtype,double> > >
747	copysign(_bz_MatExpr<P_expr1> d1, double d2)
748	{
749	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
750	_bz_copysign<typename P_expr1::T_numtype,double> > T_expr;
751
752	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
753	}
754
755	template<class P_expr1>
756	inline
757	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
758	_bz_copysign<typename P_expr1::T_numtype,long double> > >
759	copysign(_bz_MatExpr<P_expr1> d1, long double d2)
760	{
761	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
762	_bz_copysign<typename P_expr1::T_numtype,long double> > T_expr;
763
764	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
765	}
766
767	template<class P_expr1, class T2>
768	inline
769	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
770	_bz_copysign<typename P_expr1::T_numtype,complex<T2> > > >
771	copysign(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
772	{
773	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
774	_bz_copysign<typename P_expr1::T_numtype,complex<T2> > > T_expr;
775
776	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
777	}
778
779	template<class P_numtype2, class P_struct2>
780	inline
781	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
782	_bz_copysign<int,P_numtype2> > >
783	copysign(int d1, const Matrix<P_numtype2, P_struct2>& d2)
784	{
785	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
786	_bz_copysign<int,P_numtype2> > T_expr;
787
788	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
789	}
790
791	template<class P_expr2>
792	inline
793	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
794	_bz_copysign<int,typename P_expr2::T_numtype> > >
795	copysign(int d1, _bz_MatExpr<P_expr2> d2)
796	{
797	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
798	_bz_copysign<int,typename P_expr2::T_numtype> > T_expr;
799
800	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
801	}
802
803	template<class P_numtype2, class P_struct2>
804	inline
805	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
806	_bz_copysign<float,P_numtype2> > >
807	copysign(float d1, const Matrix<P_numtype2, P_struct2>& d2)
808	{
809	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
810	_bz_copysign<float,P_numtype2> > T_expr;
811
812	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
813	}
814
815	template<class P_expr2>
816	inline
817	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
818	_bz_copysign<float,typename P_expr2::T_numtype> > >
819	copysign(float d1, _bz_MatExpr<P_expr2> d2)
820	{
821	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
822	_bz_copysign<float,typename P_expr2::T_numtype> > T_expr;
823
824	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
825	}
826
827	template<class P_numtype2, class P_struct2>
828	inline
829	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
830	_bz_copysign<double,P_numtype2> > >
831	copysign(double d1, const Matrix<P_numtype2, P_struct2>& d2)
832	{
833	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
834	_bz_copysign<double,P_numtype2> > T_expr;
835
836	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
837	}
838
839	template<class P_expr2>
840	inline
841	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
842	_bz_copysign<double,typename P_expr2::T_numtype> > >
843	copysign(double d1, _bz_MatExpr<P_expr2> d2)
844	{
845	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
846	_bz_copysign<double,typename P_expr2::T_numtype> > T_expr;
847
848	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
849	}
850
851	template<class P_numtype2, class P_struct2>
852	inline
853	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
854	_bz_copysign<long double,P_numtype2> > >
855	copysign(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
856	{
857	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
858	_bz_copysign<long double,P_numtype2> > T_expr;
859
860	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
861	}
862
863	template<class P_expr2>
864	inline
865	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
866	_bz_copysign<long double,typename P_expr2::T_numtype> > >
867	copysign(long double d1, _bz_MatExpr<P_expr2> d2)
868	{
869	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
870	_bz_copysign<long double,typename P_expr2::T_numtype> > T_expr;
871
872	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
873	}
874
875	template<class T1, class P_numtype2, class P_struct2>
876	inline
877	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
878	_bz_copysign<complex<T1> ,P_numtype2> > >
879	copysign(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
880	{
881	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
882	_bz_copysign<complex<T1> ,P_numtype2> > T_expr;
883
884	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
885	}
886
887	template<class T1, class P_expr2>
888	inline
889	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
890	_bz_copysign<complex<T1> ,typename P_expr2::T_numtype> > >
891	copysign(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
892	{
893	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
894	_bz_copysign<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
895
896	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
897	}
898
899	#endif
900
901	/****************************************************************************
902	* cos
903	****************************************************************************/
904
905	template<class P_numtype1, class P_struct1>
906	inline
907	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
908	_bz_cos<P_numtype1> > >
909	cos(const Matrix<P_numtype1, P_struct1>& d1)
910	{
911	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
912	_bz_cos<P_numtype1> > T_expr;
913
914	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
915	}
916
917	template<class P_expr1>
918	inline
919	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
920	_bz_cos<typename P_expr1::T_numtype> > >
921	cos(_bz_MatExpr<P_expr1> d1)
922	{
923	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
924	_bz_cos<typename P_expr1::T_numtype> > T_expr;
925
926	return _bz_MatExpr<T_expr>(T_expr(d1));
927	}
928
929
930	/****************************************************************************
931	* cosh
932	****************************************************************************/
933
934	template<class P_numtype1, class P_struct1>
935	inline
936	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
937	_bz_cosh<P_numtype1> > >
938	cosh(const Matrix<P_numtype1, P_struct1>& d1)
939	{
940	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
941	_bz_cosh<P_numtype1> > T_expr;
942
943	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
944	}
945
946	template<class P_expr1>
947	inline
948	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
949	_bz_cosh<typename P_expr1::T_numtype> > >
950	cosh(_bz_MatExpr<P_expr1> d1)
951	{
952	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
953	_bz_cosh<typename P_expr1::T_numtype> > T_expr;
954
955	return _bz_MatExpr<T_expr>(T_expr(d1));
956	}
957
958
959	/****************************************************************************
960	* drem
961	****************************************************************************/
962
963	#ifdef BZ_HAVE_SYSTEM_V_MATH
964	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
965	inline
966	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
967	_bz_drem<P_numtype1,P_numtype2> > >
968	drem(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
969	{
970	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
971	_bz_drem<P_numtype1,P_numtype2> > T_expr;
972
973	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
974	}
975
976	template<class P_numtype1, class P_struct1, class P_expr2>
977	inline
978	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
979	_bz_drem<P_numtype1,typename P_expr2::T_numtype> > >
980	drem(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
981	{
982	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
983	_bz_drem<P_numtype1,typename P_expr2::T_numtype> > T_expr;
984
985	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
986	}
987
988	template<class P_numtype1, class P_struct1>
989	inline
990	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
991	_bz_drem<P_numtype1,int> > >
992	drem(const Matrix<P_numtype1, P_struct1>& d1, int d2)
993	{
994	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
995	_bz_drem<P_numtype1,int> > T_expr;
996
997	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
998	}
999
1000	template<class P_numtype1, class P_struct1>
1001	inline
1002	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1003	_bz_drem<P_numtype1,float> > >
1004	drem(const Matrix<P_numtype1, P_struct1>& d1, float d2)
1005	{
1006	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1007	_bz_drem<P_numtype1,float> > T_expr;
1008
1009	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
1010	}
1011
1012	template<class P_numtype1, class P_struct1>
1013	inline
1014	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1015	_bz_drem<P_numtype1,double> > >
1016	drem(const Matrix<P_numtype1, P_struct1>& d1, double d2)
1017	{
1018	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1019	_bz_drem<P_numtype1,double> > T_expr;
1020
1021	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
1022	}
1023
1024	template<class P_numtype1, class P_struct1>
1025	inline
1026	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1027	_bz_drem<P_numtype1,long double> > >
1028	drem(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
1029	{
1030	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1031	_bz_drem<P_numtype1,long double> > T_expr;
1032
1033	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
1034	}
1035
1036	template<class P_numtype1, class P_struct1, class T2>
1037	inline
1038	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1039	_bz_drem<P_numtype1,complex<T2> > > >
1040	drem(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
1041	{
1042	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1043	_bz_drem<P_numtype1,complex<T2> > > T_expr;
1044
1045	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
1046	}
1047
1048	template<class P_expr1, class P_numtype2, class P_struct2>
1049	inline
1050	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1051	_bz_drem<typename P_expr1::T_numtype,P_numtype2> > >
1052	drem(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
1053	{
1054	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1055	_bz_drem<typename P_expr1::T_numtype,P_numtype2> > T_expr;
1056
1057	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
1058	}
1059
1060	template<class P_expr1, class P_expr2>
1061	inline
1062	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1063	_bz_drem<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
1064	drem(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
1065	{
1066	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1067	_bz_drem<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
1068
1069	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
1070	}
1071
1072	template<class P_expr1>
1073	inline
1074	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1075	_bz_drem<typename P_expr1::T_numtype,int> > >
1076	drem(_bz_MatExpr<P_expr1> d1, int d2)
1077	{
1078	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1079	_bz_drem<typename P_expr1::T_numtype,int> > T_expr;
1080
1081	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
1082	}
1083
1084	template<class P_expr1>
1085	inline
1086	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1087	_bz_drem<typename P_expr1::T_numtype,float> > >
1088	drem(_bz_MatExpr<P_expr1> d1, float d2)
1089	{
1090	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1091	_bz_drem<typename P_expr1::T_numtype,float> > T_expr;
1092
1093	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
1094	}
1095
1096	template<class P_expr1>
1097	inline
1098	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1099	_bz_drem<typename P_expr1::T_numtype,double> > >
1100	drem(_bz_MatExpr<P_expr1> d1, double d2)
1101	{
1102	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1103	_bz_drem<typename P_expr1::T_numtype,double> > T_expr;
1104
1105	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
1106	}
1107
1108	template<class P_expr1>
1109	inline
1110	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1111	_bz_drem<typename P_expr1::T_numtype,long double> > >
1112	drem(_bz_MatExpr<P_expr1> d1, long double d2)
1113	{
1114	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1115	_bz_drem<typename P_expr1::T_numtype,long double> > T_expr;
1116
1117	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
1118	}
1119
1120	template<class P_expr1, class T2>
1121	inline
1122	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1123	_bz_drem<typename P_expr1::T_numtype,complex<T2> > > >
1124	drem(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
1125	{
1126	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1127	_bz_drem<typename P_expr1::T_numtype,complex<T2> > > T_expr;
1128
1129	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
1130	}
1131
1132	template<class P_numtype2, class P_struct2>
1133	inline
1134	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1135	_bz_drem<int,P_numtype2> > >
1136	drem(int d1, const Matrix<P_numtype2, P_struct2>& d2)
1137	{
1138	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1139	_bz_drem<int,P_numtype2> > T_expr;
1140
1141	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
1142	}
1143
1144	template<class P_expr2>
1145	inline
1146	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1147	_bz_drem<int,typename P_expr2::T_numtype> > >
1148	drem(int d1, _bz_MatExpr<P_expr2> d2)
1149	{
1150	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1151	_bz_drem<int,typename P_expr2::T_numtype> > T_expr;
1152
1153	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
1154	}
1155
1156	template<class P_numtype2, class P_struct2>
1157	inline
1158	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1159	_bz_drem<float,P_numtype2> > >
1160	drem(float d1, const Matrix<P_numtype2, P_struct2>& d2)
1161	{
1162	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1163	_bz_drem<float,P_numtype2> > T_expr;
1164
1165	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
1166	}
1167
1168	template<class P_expr2>
1169	inline
1170	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1171	_bz_drem<float,typename P_expr2::T_numtype> > >
1172	drem(float d1, _bz_MatExpr<P_expr2> d2)
1173	{
1174	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1175	_bz_drem<float,typename P_expr2::T_numtype> > T_expr;
1176
1177	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
1178	}
1179
1180	template<class P_numtype2, class P_struct2>
1181	inline
1182	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1183	_bz_drem<double,P_numtype2> > >
1184	drem(double d1, const Matrix<P_numtype2, P_struct2>& d2)
1185	{
1186	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1187	_bz_drem<double,P_numtype2> > T_expr;
1188
1189	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
1190	}
1191
1192	template<class P_expr2>
1193	inline
1194	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1195	_bz_drem<double,typename P_expr2::T_numtype> > >
1196	drem(double d1, _bz_MatExpr<P_expr2> d2)
1197	{
1198	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1199	_bz_drem<double,typename P_expr2::T_numtype> > T_expr;
1200
1201	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
1202	}
1203
1204	template<class P_numtype2, class P_struct2>
1205	inline
1206	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1207	_bz_drem<long double,P_numtype2> > >
1208	drem(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
1209	{
1210	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1211	_bz_drem<long double,P_numtype2> > T_expr;
1212
1213	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
1214	}
1215
1216	template<class P_expr2>
1217	inline
1218	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1219	_bz_drem<long double,typename P_expr2::T_numtype> > >
1220	drem(long double d1, _bz_MatExpr<P_expr2> d2)
1221	{
1222	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1223	_bz_drem<long double,typename P_expr2::T_numtype> > T_expr;
1224
1225	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
1226	}
1227
1228	template<class T1, class P_numtype2, class P_struct2>
1229	inline
1230	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
1231	_bz_drem<complex<T1> ,P_numtype2> > >
1232	drem(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
1233	{
1234	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
1235	_bz_drem<complex<T1> ,P_numtype2> > T_expr;
1236
1237	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
1238	}
1239
1240	template<class T1, class P_expr2>
1241	inline
1242	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
1243	_bz_drem<complex<T1> ,typename P_expr2::T_numtype> > >
1244	drem(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
1245	{
1246	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
1247	_bz_drem<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
1248
1249	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
1250	}
1251
1252	#endif
1253
1254	/****************************************************************************
1255	* exp
1256	****************************************************************************/
1257
1258	template<class P_numtype1, class P_struct1>
1259	inline
1260	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1261	_bz_exp<P_numtype1> > >
1262	exp(const Matrix<P_numtype1, P_struct1>& d1)
1263	{
1264	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1265	_bz_exp<P_numtype1> > T_expr;
1266
1267	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1268	}
1269
1270	template<class P_expr1>
1271	inline
1272	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1273	_bz_exp<typename P_expr1::T_numtype> > >
1274	exp(_bz_MatExpr<P_expr1> d1)
1275	{
1276	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1277	_bz_exp<typename P_expr1::T_numtype> > T_expr;
1278
1279	return _bz_MatExpr<T_expr>(T_expr(d1));
1280	}
1281
1282
1283	/****************************************************************************
1284	* expm1
1285	****************************************************************************/
1286
1287	#ifdef BZ_HAVE_IEEE_MATH
1288	template<class P_numtype1, class P_struct1>
1289	inline
1290	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1291	_bz_expm1<P_numtype1> > >
1292	expm1(const Matrix<P_numtype1, P_struct1>& d1)
1293	{
1294	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1295	_bz_expm1<P_numtype1> > T_expr;
1296
1297	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1298	}
1299
1300	template<class P_expr1>
1301	inline
1302	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1303	_bz_expm1<typename P_expr1::T_numtype> > >
1304	expm1(_bz_MatExpr<P_expr1> d1)
1305	{
1306	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1307	_bz_expm1<typename P_expr1::T_numtype> > T_expr;
1308
1309	return _bz_MatExpr<T_expr>(T_expr(d1));
1310	}
1311
1312	#endif
1313
1314	/****************************************************************************
1315	* erf
1316	****************************************************************************/
1317
1318	#ifdef BZ_HAVE_IEEE_MATH
1319	template<class P_numtype1, class P_struct1>
1320	inline
1321	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1322	_bz_erf<P_numtype1> > >
1323	erf(const Matrix<P_numtype1, P_struct1>& d1)
1324	{
1325	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1326	_bz_erf<P_numtype1> > T_expr;
1327
1328	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1329	}
1330
1331	template<class P_expr1>
1332	inline
1333	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1334	_bz_erf<typename P_expr1::T_numtype> > >
1335	erf(_bz_MatExpr<P_expr1> d1)
1336	{
1337	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1338	_bz_erf<typename P_expr1::T_numtype> > T_expr;
1339
1340	return _bz_MatExpr<T_expr>(T_expr(d1));
1341	}
1342
1343	#endif
1344
1345	/****************************************************************************
1346	* erfc
1347	****************************************************************************/
1348
1349	#ifdef BZ_HAVE_IEEE_MATH
1350	template<class P_numtype1, class P_struct1>
1351	inline
1352	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1353	_bz_erfc<P_numtype1> > >
1354	erfc(const Matrix<P_numtype1, P_struct1>& d1)
1355	{
1356	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1357	_bz_erfc<P_numtype1> > T_expr;
1358
1359	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1360	}
1361
1362	template<class P_expr1>
1363	inline
1364	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1365	_bz_erfc<typename P_expr1::T_numtype> > >
1366	erfc(_bz_MatExpr<P_expr1> d1)
1367	{
1368	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1369	_bz_erfc<typename P_expr1::T_numtype> > T_expr;
1370
1371	return _bz_MatExpr<T_expr>(T_expr(d1));
1372	}
1373
1374	#endif
1375
1376	/****************************************************************************
1377	* fabs
1378	****************************************************************************/
1379
1380	template<class P_numtype1, class P_struct1>
1381	inline
1382	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1383	_bz_abs<P_numtype1> > >
1384	fabs(const Matrix<P_numtype1, P_struct1>& d1)
1385	{
1386	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1387	_bz_abs<P_numtype1> > T_expr;
1388
1389	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1390	}
1391
1392	template<class P_expr1>
1393	inline
1394	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1395	_bz_abs<typename P_expr1::T_numtype> > >
1396	fabs(_bz_MatExpr<P_expr1> d1)
1397	{
1398	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1399	_bz_abs<typename P_expr1::T_numtype> > T_expr;
1400
1401	return _bz_MatExpr<T_expr>(T_expr(d1));
1402	}
1403
1404
1405	/****************************************************************************
1406	* floor
1407	****************************************************************************/
1408
1409	template<class P_numtype1, class P_struct1>
1410	inline
1411	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1412	_bz_floor<P_numtype1> > >
1413	floor(const Matrix<P_numtype1, P_struct1>& d1)
1414	{
1415	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
1416	_bz_floor<P_numtype1> > T_expr;
1417
1418	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
1419	}
1420
1421	template<class P_expr1>
1422	inline
1423	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1424	_bz_floor<typename P_expr1::T_numtype> > >
1425	floor(_bz_MatExpr<P_expr1> d1)
1426	{
1427	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
1428	_bz_floor<typename P_expr1::T_numtype> > T_expr;
1429
1430	return _bz_MatExpr<T_expr>(T_expr(d1));
1431	}
1432
1433
1434	/****************************************************************************
1435	* fmod
1436	****************************************************************************/
1437
1438	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
1439	inline
1440	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1441	_bz_fmod<P_numtype1,P_numtype2> > >
1442	fmod(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
1443	{
1444	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1445	_bz_fmod<P_numtype1,P_numtype2> > T_expr;
1446
1447	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
1448	}
1449
1450	template<class P_numtype1, class P_struct1, class P_expr2>
1451	inline
1452	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
1453	_bz_fmod<P_numtype1,typename P_expr2::T_numtype> > >
1454	fmod(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
1455	{
1456	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
1457	_bz_fmod<P_numtype1,typename P_expr2::T_numtype> > T_expr;
1458
1459	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
1460	}
1461
1462	template<class P_numtype1, class P_struct1>
1463	inline
1464	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
1465	_bz_fmod<P_numtype1,int> > >
1466	fmod(const Matrix<P_numtype1, P_struct1>& d1, int d2)
1467	{
1468	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
1469	_bz_fmod<P_numtype1,int> > T_expr;
1470
1471	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
1472	}
1473
1474	template<class P_numtype1, class P_struct1>
1475	inline
1476	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1477	_bz_fmod<P_numtype1,float> > >
1478	fmod(const Matrix<P_numtype1, P_struct1>& d1, float d2)
1479	{
1480	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1481	_bz_fmod<P_numtype1,float> > T_expr;
1482
1483	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
1484	}
1485
1486	template<class P_numtype1, class P_struct1>
1487	inline
1488	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1489	_bz_fmod<P_numtype1,double> > >
1490	fmod(const Matrix<P_numtype1, P_struct1>& d1, double d2)
1491	{
1492	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1493	_bz_fmod<P_numtype1,double> > T_expr;
1494
1495	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
1496	}
1497
1498	template<class P_numtype1, class P_struct1>
1499	inline
1500	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1501	_bz_fmod<P_numtype1,long double> > >
1502	fmod(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
1503	{
1504	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1505	_bz_fmod<P_numtype1,long double> > T_expr;
1506
1507	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
1508	}
1509
1510	template<class P_numtype1, class P_struct1, class T2>
1511	inline
1512	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1513	_bz_fmod<P_numtype1,complex<T2> > > >
1514	fmod(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
1515	{
1516	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1517	_bz_fmod<P_numtype1,complex<T2> > > T_expr;
1518
1519	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
1520	}
1521
1522	template<class P_expr1, class P_numtype2, class P_struct2>
1523	inline
1524	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1525	_bz_fmod<typename P_expr1::T_numtype,P_numtype2> > >
1526	fmod(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
1527	{
1528	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1529	_bz_fmod<typename P_expr1::T_numtype,P_numtype2> > T_expr;
1530
1531	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
1532	}
1533
1534	template<class P_expr1, class P_expr2>
1535	inline
1536	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1537	_bz_fmod<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
1538	fmod(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
1539	{
1540	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1541	_bz_fmod<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
1542
1543	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
1544	}
1545
1546	template<class P_expr1>
1547	inline
1548	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1549	_bz_fmod<typename P_expr1::T_numtype,int> > >
1550	fmod(_bz_MatExpr<P_expr1> d1, int d2)
1551	{
1552	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1553	_bz_fmod<typename P_expr1::T_numtype,int> > T_expr;
1554
1555	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
1556	}
1557
1558	template<class P_expr1>
1559	inline
1560	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1561	_bz_fmod<typename P_expr1::T_numtype,float> > >
1562	fmod(_bz_MatExpr<P_expr1> d1, float d2)
1563	{
1564	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1565	_bz_fmod<typename P_expr1::T_numtype,float> > T_expr;
1566
1567	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
1568	}
1569
1570	template<class P_expr1>
1571	inline
1572	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1573	_bz_fmod<typename P_expr1::T_numtype,double> > >
1574	fmod(_bz_MatExpr<P_expr1> d1, double d2)
1575	{
1576	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1577	_bz_fmod<typename P_expr1::T_numtype,double> > T_expr;
1578
1579	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
1580	}
1581
1582	template<class P_expr1>
1583	inline
1584	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1585	_bz_fmod<typename P_expr1::T_numtype,long double> > >
1586	fmod(_bz_MatExpr<P_expr1> d1, long double d2)
1587	{
1588	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1589	_bz_fmod<typename P_expr1::T_numtype,long double> > T_expr;
1590
1591	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
1592	}
1593
1594	template<class P_expr1, class T2>
1595	inline
1596	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1597	_bz_fmod<typename P_expr1::T_numtype,complex<T2> > > >
1598	fmod(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
1599	{
1600	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1601	_bz_fmod<typename P_expr1::T_numtype,complex<T2> > > T_expr;
1602
1603	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
1604	}
1605
1606	template<class P_numtype2, class P_struct2>
1607	inline
1608	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1609	_bz_fmod<int,P_numtype2> > >
1610	fmod(int d1, const Matrix<P_numtype2, P_struct2>& d2)
1611	{
1612	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1613	_bz_fmod<int,P_numtype2> > T_expr;
1614
1615	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
1616	}
1617
1618	template<class P_expr2>
1619	inline
1620	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1621	_bz_fmod<int,typename P_expr2::T_numtype> > >
1622	fmod(int d1, _bz_MatExpr<P_expr2> d2)
1623	{
1624	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1625	_bz_fmod<int,typename P_expr2::T_numtype> > T_expr;
1626
1627	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
1628	}
1629
1630	template<class P_numtype2, class P_struct2>
1631	inline
1632	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1633	_bz_fmod<float,P_numtype2> > >
1634	fmod(float d1, const Matrix<P_numtype2, P_struct2>& d2)
1635	{
1636	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1637	_bz_fmod<float,P_numtype2> > T_expr;
1638
1639	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
1640	}
1641
1642	template<class P_expr2>
1643	inline
1644	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1645	_bz_fmod<float,typename P_expr2::T_numtype> > >
1646	fmod(float d1, _bz_MatExpr<P_expr2> d2)
1647	{
1648	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1649	_bz_fmod<float,typename P_expr2::T_numtype> > T_expr;
1650
1651	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
1652	}
1653
1654	template<class P_numtype2, class P_struct2>
1655	inline
1656	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1657	_bz_fmod<double,P_numtype2> > >
1658	fmod(double d1, const Matrix<P_numtype2, P_struct2>& d2)
1659	{
1660	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1661	_bz_fmod<double,P_numtype2> > T_expr;
1662
1663	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
1664	}
1665
1666	template<class P_expr2>
1667	inline
1668	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1669	_bz_fmod<double,typename P_expr2::T_numtype> > >
1670	fmod(double d1, _bz_MatExpr<P_expr2> d2)
1671	{
1672	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1673	_bz_fmod<double,typename P_expr2::T_numtype> > T_expr;
1674
1675	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
1676	}
1677
1678	template<class P_numtype2, class P_struct2>
1679	inline
1680	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1681	_bz_fmod<long double,P_numtype2> > >
1682	fmod(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
1683	{
1684	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1685	_bz_fmod<long double,P_numtype2> > T_expr;
1686
1687	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
1688	}
1689
1690	template<class P_expr2>
1691	inline
1692	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1693	_bz_fmod<long double,typename P_expr2::T_numtype> > >
1694	fmod(long double d1, _bz_MatExpr<P_expr2> d2)
1695	{
1696	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1697	_bz_fmod<long double,typename P_expr2::T_numtype> > T_expr;
1698
1699	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
1700	}
1701
1702	template<class T1, class P_numtype2, class P_struct2>
1703	inline
1704	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
1705	_bz_fmod<complex<T1> ,P_numtype2> > >
1706	fmod(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
1707	{
1708	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
1709	_bz_fmod<complex<T1> ,P_numtype2> > T_expr;
1710
1711	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
1712	}
1713
1714	template<class T1, class P_expr2>
1715	inline
1716	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
1717	_bz_fmod<complex<T1> ,typename P_expr2::T_numtype> > >
1718	fmod(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
1719	{
1720	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
1721	_bz_fmod<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
1722
1723	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
1724	}
1725
1726
1727	/****************************************************************************
1728	* hypot
1729	****************************************************************************/
1730
1731	#ifdef BZ_HAVE_SYSTEM_V_MATH
1732	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
1733	inline
1734	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1735	_bz_hypot<P_numtype1,P_numtype2> > >
1736	hypot(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
1737	{
1738	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1739	_bz_hypot<P_numtype1,P_numtype2> > T_expr;
1740
1741	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
1742	}
1743
1744	template<class P_numtype1, class P_struct1, class P_expr2>
1745	inline
1746	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
1747	_bz_hypot<P_numtype1,typename P_expr2::T_numtype> > >
1748	hypot(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
1749	{
1750	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
1751	_bz_hypot<P_numtype1,typename P_expr2::T_numtype> > T_expr;
1752
1753	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
1754	}
1755
1756	template<class P_numtype1, class P_struct1>
1757	inline
1758	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
1759	_bz_hypot<P_numtype1,int> > >
1760	hypot(const Matrix<P_numtype1, P_struct1>& d1, int d2)
1761	{
1762	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
1763	_bz_hypot<P_numtype1,int> > T_expr;
1764
1765	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
1766	}
1767
1768	template<class P_numtype1, class P_struct1>
1769	inline
1770	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1771	_bz_hypot<P_numtype1,float> > >
1772	hypot(const Matrix<P_numtype1, P_struct1>& d1, float d2)
1773	{
1774	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
1775	_bz_hypot<P_numtype1,float> > T_expr;
1776
1777	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
1778	}
1779
1780	template<class P_numtype1, class P_struct1>
1781	inline
1782	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1783	_bz_hypot<P_numtype1,double> > >
1784	hypot(const Matrix<P_numtype1, P_struct1>& d1, double d2)
1785	{
1786	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
1787	_bz_hypot<P_numtype1,double> > T_expr;
1788
1789	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
1790	}
1791
1792	template<class P_numtype1, class P_struct1>
1793	inline
1794	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1795	_bz_hypot<P_numtype1,long double> > >
1796	hypot(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
1797	{
1798	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
1799	_bz_hypot<P_numtype1,long double> > T_expr;
1800
1801	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
1802	}
1803
1804	template<class P_numtype1, class P_struct1, class T2>
1805	inline
1806	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1807	_bz_hypot<P_numtype1,complex<T2> > > >
1808	hypot(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
1809	{
1810	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
1811	_bz_hypot<P_numtype1,complex<T2> > > T_expr;
1812
1813	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
1814	}
1815
1816	template<class P_expr1, class P_numtype2, class P_struct2>
1817	inline
1818	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1819	_bz_hypot<typename P_expr1::T_numtype,P_numtype2> > >
1820	hypot(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
1821	{
1822	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
1823	_bz_hypot<typename P_expr1::T_numtype,P_numtype2> > T_expr;
1824
1825	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
1826	}
1827
1828	template<class P_expr1, class P_expr2>
1829	inline
1830	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1831	_bz_hypot<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
1832	hypot(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
1833	{
1834	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
1835	_bz_hypot<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
1836
1837	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
1838	}
1839
1840	template<class P_expr1>
1841	inline
1842	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1843	_bz_hypot<typename P_expr1::T_numtype,int> > >
1844	hypot(_bz_MatExpr<P_expr1> d1, int d2)
1845	{
1846	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
1847	_bz_hypot<typename P_expr1::T_numtype,int> > T_expr;
1848
1849	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
1850	}
1851
1852	template<class P_expr1>
1853	inline
1854	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1855	_bz_hypot<typename P_expr1::T_numtype,float> > >
1856	hypot(_bz_MatExpr<P_expr1> d1, float d2)
1857	{
1858	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
1859	_bz_hypot<typename P_expr1::T_numtype,float> > T_expr;
1860
1861	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
1862	}
1863
1864	template<class P_expr1>
1865	inline
1866	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1867	_bz_hypot<typename P_expr1::T_numtype,double> > >
1868	hypot(_bz_MatExpr<P_expr1> d1, double d2)
1869	{
1870	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
1871	_bz_hypot<typename P_expr1::T_numtype,double> > T_expr;
1872
1873	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
1874	}
1875
1876	template<class P_expr1>
1877	inline
1878	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1879	_bz_hypot<typename P_expr1::T_numtype,long double> > >
1880	hypot(_bz_MatExpr<P_expr1> d1, long double d2)
1881	{
1882	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
1883	_bz_hypot<typename P_expr1::T_numtype,long double> > T_expr;
1884
1885	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
1886	}
1887
1888	template<class P_expr1, class T2>
1889	inline
1890	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1891	_bz_hypot<typename P_expr1::T_numtype,complex<T2> > > >
1892	hypot(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
1893	{
1894	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
1895	_bz_hypot<typename P_expr1::T_numtype,complex<T2> > > T_expr;
1896
1897	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
1898	}
1899
1900	template<class P_numtype2, class P_struct2>
1901	inline
1902	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1903	_bz_hypot<int,P_numtype2> > >
1904	hypot(int d1, const Matrix<P_numtype2, P_struct2>& d2)
1905	{
1906	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
1907	_bz_hypot<int,P_numtype2> > T_expr;
1908
1909	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
1910	}
1911
1912	template<class P_expr2>
1913	inline
1914	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1915	_bz_hypot<int,typename P_expr2::T_numtype> > >
1916	hypot(int d1, _bz_MatExpr<P_expr2> d2)
1917	{
1918	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
1919	_bz_hypot<int,typename P_expr2::T_numtype> > T_expr;
1920
1921	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
1922	}
1923
1924	template<class P_numtype2, class P_struct2>
1925	inline
1926	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1927	_bz_hypot<float,P_numtype2> > >
1928	hypot(float d1, const Matrix<P_numtype2, P_struct2>& d2)
1929	{
1930	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
1931	_bz_hypot<float,P_numtype2> > T_expr;
1932
1933	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
1934	}
1935
1936	template<class P_expr2>
1937	inline
1938	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1939	_bz_hypot<float,typename P_expr2::T_numtype> > >
1940	hypot(float d1, _bz_MatExpr<P_expr2> d2)
1941	{
1942	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
1943	_bz_hypot<float,typename P_expr2::T_numtype> > T_expr;
1944
1945	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
1946	}
1947
1948	template<class P_numtype2, class P_struct2>
1949	inline
1950	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1951	_bz_hypot<double,P_numtype2> > >
1952	hypot(double d1, const Matrix<P_numtype2, P_struct2>& d2)
1953	{
1954	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1955	_bz_hypot<double,P_numtype2> > T_expr;
1956
1957	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
1958	}
1959
1960	template<class P_expr2>
1961	inline
1962	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1963	_bz_hypot<double,typename P_expr2::T_numtype> > >
1964	hypot(double d1, _bz_MatExpr<P_expr2> d2)
1965	{
1966	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
1967	_bz_hypot<double,typename P_expr2::T_numtype> > T_expr;
1968
1969	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
1970	}
1971
1972	template<class P_numtype2, class P_struct2>
1973	inline
1974	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1975	_bz_hypot<long double,P_numtype2> > >
1976	hypot(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
1977	{
1978	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
1979	_bz_hypot<long double,P_numtype2> > T_expr;
1980
1981	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
1982	}
1983
1984	template<class P_expr2>
1985	inline
1986	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1987	_bz_hypot<long double,typename P_expr2::T_numtype> > >
1988	hypot(long double d1, _bz_MatExpr<P_expr2> d2)
1989	{
1990	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
1991	_bz_hypot<long double,typename P_expr2::T_numtype> > T_expr;
1992
1993	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
1994	}
1995
1996	template<class T1, class P_numtype2, class P_struct2>
1997	inline
1998	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
1999	_bz_hypot<complex<T1> ,P_numtype2> > >
2000	hypot(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
2001	{
2002	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
2003	_bz_hypot<complex<T1> ,P_numtype2> > T_expr;
2004
2005	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
2006	}
2007
2008	template<class T1, class P_expr2>
2009	inline
2010	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2011	_bz_hypot<complex<T1> ,typename P_expr2::T_numtype> > >
2012	hypot(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
2013	{
2014	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2015	_bz_hypot<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
2016
2017	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
2018	}
2019
2020	#endif
2021
2022	/****************************************************************************
2023	* ilogb
2024	****************************************************************************/
2025
2026	#ifdef BZ_HAVE_IEEE_MATH
2027	template<class P_numtype1, class P_struct1>
2028	inline
2029	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2030	_bz_ilogb<P_numtype1> > >
2031	ilogb(const Matrix<P_numtype1, P_struct1>& d1)
2032	{
2033	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2034	_bz_ilogb<P_numtype1> > T_expr;
2035
2036	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2037	}
2038
2039	template<class P_expr1>
2040	inline
2041	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2042	_bz_ilogb<typename P_expr1::T_numtype> > >
2043	ilogb(_bz_MatExpr<P_expr1> d1)
2044	{
2045	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2046	_bz_ilogb<typename P_expr1::T_numtype> > T_expr;
2047
2048	return _bz_MatExpr<T_expr>(T_expr(d1));
2049	}
2050
2051	#endif
2052
2053	/****************************************************************************
2054	* blitz_isnan
2055	****************************************************************************/
2056
2057	#ifdef BZ_HAVE_IEEE_MATH
2058	template<class P_numtype1, class P_struct1>
2059	inline
2060	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2061	_bz_blitz_isnan<P_numtype1> > >
2062	blitz_isnan(const Matrix<P_numtype1, P_struct1>& d1)
2063	{
2064	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2065	_bz_blitz_isnan<P_numtype1> > T_expr;
2066
2067	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2068	}
2069
2070	template<class P_expr1>
2071	inline
2072	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2073	_bz_blitz_isnan<typename P_expr1::T_numtype> > >
2074	blitz_isnan(_bz_MatExpr<P_expr1> d1)
2075	{
2076	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2077	_bz_blitz_isnan<typename P_expr1::T_numtype> > T_expr;
2078
2079	return _bz_MatExpr<T_expr>(T_expr(d1));
2080	}
2081
2082	#endif
2083
2084	/****************************************************************************
2085	* itrunc
2086	****************************************************************************/
2087
2088	#ifdef BZ_HAVE_SYSTEM_V_MATH
2089	template<class P_numtype1, class P_struct1>
2090	inline
2091	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2092	_bz_itrunc<P_numtype1> > >
2093	itrunc(const Matrix<P_numtype1, P_struct1>& d1)
2094	{
2095	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2096	_bz_itrunc<P_numtype1> > T_expr;
2097
2098	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2099	}
2100
2101	template<class P_expr1>
2102	inline
2103	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2104	_bz_itrunc<typename P_expr1::T_numtype> > >
2105	itrunc(_bz_MatExpr<P_expr1> d1)
2106	{
2107	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2108	_bz_itrunc<typename P_expr1::T_numtype> > T_expr;
2109
2110	return _bz_MatExpr<T_expr>(T_expr(d1));
2111	}
2112
2113	#endif
2114
2115	/****************************************************************************
2116	* j0
2117	****************************************************************************/
2118
2119	#ifdef BZ_HAVE_IEEE_MATH
2120	template<class P_numtype1, class P_struct1>
2121	inline
2122	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2123	_bz_j0<P_numtype1> > >
2124	j0(const Matrix<P_numtype1, P_struct1>& d1)
2125	{
2126	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2127	_bz_j0<P_numtype1> > T_expr;
2128
2129	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2130	}
2131
2132	template<class P_expr1>
2133	inline
2134	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2135	_bz_j0<typename P_expr1::T_numtype> > >
2136	j0(_bz_MatExpr<P_expr1> d1)
2137	{
2138	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2139	_bz_j0<typename P_expr1::T_numtype> > T_expr;
2140
2141	return _bz_MatExpr<T_expr>(T_expr(d1));
2142	}
2143
2144	#endif
2145
2146	/****************************************************************************
2147	* j1
2148	****************************************************************************/
2149
2150	#ifdef BZ_HAVE_IEEE_MATH
2151	template<class P_numtype1, class P_struct1>
2152	inline
2153	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2154	_bz_j1<P_numtype1> > >
2155	j1(const Matrix<P_numtype1, P_struct1>& d1)
2156	{
2157	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2158	_bz_j1<P_numtype1> > T_expr;
2159
2160	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2161	}
2162
2163	template<class P_expr1>
2164	inline
2165	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2166	_bz_j1<typename P_expr1::T_numtype> > >
2167	j1(_bz_MatExpr<P_expr1> d1)
2168	{
2169	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2170	_bz_j1<typename P_expr1::T_numtype> > T_expr;
2171
2172	return _bz_MatExpr<T_expr>(T_expr(d1));
2173	}
2174
2175	#endif
2176
2177	/****************************************************************************
2178	* lgamma
2179	****************************************************************************/
2180
2181	#ifdef BZ_HAVE_IEEE_MATH
2182	template<class P_numtype1, class P_struct1>
2183	inline
2184	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2185	_bz_lgamma<P_numtype1> > >
2186	lgamma(const Matrix<P_numtype1, P_struct1>& d1)
2187	{
2188	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2189	_bz_lgamma<P_numtype1> > T_expr;
2190
2191	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2192	}
2193
2194	template<class P_expr1>
2195	inline
2196	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2197	_bz_lgamma<typename P_expr1::T_numtype> > >
2198	lgamma(_bz_MatExpr<P_expr1> d1)
2199	{
2200	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2201	_bz_lgamma<typename P_expr1::T_numtype> > T_expr;
2202
2203	return _bz_MatExpr<T_expr>(T_expr(d1));
2204	}
2205
2206	#endif
2207
2208	/****************************************************************************
2209	* log
2210	****************************************************************************/
2211
2212	template<class P_numtype1, class P_struct1>
2213	inline
2214	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2215	_bz_log<P_numtype1> > >
2216	log(const Matrix<P_numtype1, P_struct1>& d1)
2217	{
2218	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2219	_bz_log<P_numtype1> > T_expr;
2220
2221	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2222	}
2223
2224	template<class P_expr1>
2225	inline
2226	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2227	_bz_log<typename P_expr1::T_numtype> > >
2228	log(_bz_MatExpr<P_expr1> d1)
2229	{
2230	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2231	_bz_log<typename P_expr1::T_numtype> > T_expr;
2232
2233	return _bz_MatExpr<T_expr>(T_expr(d1));
2234	}
2235
2236
2237	/****************************************************************************
2238	* logb
2239	****************************************************************************/
2240
2241	#ifdef BZ_HAVE_IEEE_MATH
2242	template<class P_numtype1, class P_struct1>
2243	inline
2244	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2245	_bz_logb<P_numtype1> > >
2246	logb(const Matrix<P_numtype1, P_struct1>& d1)
2247	{
2248	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2249	_bz_logb<P_numtype1> > T_expr;
2250
2251	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2252	}
2253
2254	template<class P_expr1>
2255	inline
2256	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2257	_bz_logb<typename P_expr1::T_numtype> > >
2258	logb(_bz_MatExpr<P_expr1> d1)
2259	{
2260	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2261	_bz_logb<typename P_expr1::T_numtype> > T_expr;
2262
2263	return _bz_MatExpr<T_expr>(T_expr(d1));
2264	}
2265
2266	#endif
2267
2268	/****************************************************************************
2269	* log1p
2270	****************************************************************************/
2271
2272	#ifdef BZ_HAVE_IEEE_MATH
2273	template<class P_numtype1, class P_struct1>
2274	inline
2275	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2276	_bz_log1p<P_numtype1> > >
2277	log1p(const Matrix<P_numtype1, P_struct1>& d1)
2278	{
2279	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2280	_bz_log1p<P_numtype1> > T_expr;
2281
2282	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2283	}
2284
2285	template<class P_expr1>
2286	inline
2287	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2288	_bz_log1p<typename P_expr1::T_numtype> > >
2289	log1p(_bz_MatExpr<P_expr1> d1)
2290	{
2291	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2292	_bz_log1p<typename P_expr1::T_numtype> > T_expr;
2293
2294	return _bz_MatExpr<T_expr>(T_expr(d1));
2295	}
2296
2297	#endif
2298
2299	/****************************************************************************
2300	* log10
2301	****************************************************************************/
2302
2303	template<class P_numtype1, class P_struct1>
2304	inline
2305	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2306	_bz_log10<P_numtype1> > >
2307	log10(const Matrix<P_numtype1, P_struct1>& d1)
2308	{
2309	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2310	_bz_log10<P_numtype1> > T_expr;
2311
2312	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2313	}
2314
2315	template<class P_expr1>
2316	inline
2317	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2318	_bz_log10<typename P_expr1::T_numtype> > >
2319	log10(_bz_MatExpr<P_expr1> d1)
2320	{
2321	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2322	_bz_log10<typename P_expr1::T_numtype> > T_expr;
2323
2324	return _bz_MatExpr<T_expr>(T_expr(d1));
2325	}
2326
2327
2328	/****************************************************************************
2329	* nearest
2330	****************************************************************************/
2331
2332	#ifdef BZ_HAVE_SYSTEM_V_MATH
2333	template<class P_numtype1, class P_struct1>
2334	inline
2335	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2336	_bz_nearest<P_numtype1> > >
2337	nearest(const Matrix<P_numtype1, P_struct1>& d1)
2338	{
2339	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
2340	_bz_nearest<P_numtype1> > T_expr;
2341
2342	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
2343	}
2344
2345	template<class P_expr1>
2346	inline
2347	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2348	_bz_nearest<typename P_expr1::T_numtype> > >
2349	nearest(_bz_MatExpr<P_expr1> d1)
2350	{
2351	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
2352	_bz_nearest<typename P_expr1::T_numtype> > T_expr;
2353
2354	return _bz_MatExpr<T_expr>(T_expr(d1));
2355	}
2356
2357	#endif
2358
2359	/****************************************************************************
2360	* nextafter
2361	****************************************************************************/
2362
2363	#ifdef BZ_HAVE_SYSTEM_V_MATH
2364	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
2365	inline
2366	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2367	_bz_nextafter<P_numtype1,P_numtype2> > >
2368	nextafter(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
2369	{
2370	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2371	_bz_nextafter<P_numtype1,P_numtype2> > T_expr;
2372
2373	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
2374	}
2375
2376	template<class P_numtype1, class P_struct1, class P_expr2>
2377	inline
2378	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2379	_bz_nextafter<P_numtype1,typename P_expr2::T_numtype> > >
2380	nextafter(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
2381	{
2382	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2383	_bz_nextafter<P_numtype1,typename P_expr2::T_numtype> > T_expr;
2384
2385	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
2386	}
2387
2388	template<class P_numtype1, class P_struct1>
2389	inline
2390	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2391	_bz_nextafter<P_numtype1,int> > >
2392	nextafter(const Matrix<P_numtype1, P_struct1>& d1, int d2)
2393	{
2394	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2395	_bz_nextafter<P_numtype1,int> > T_expr;
2396
2397	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
2398	}
2399
2400	template<class P_numtype1, class P_struct1>
2401	inline
2402	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2403	_bz_nextafter<P_numtype1,float> > >
2404	nextafter(const Matrix<P_numtype1, P_struct1>& d1, float d2)
2405	{
2406	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2407	_bz_nextafter<P_numtype1,float> > T_expr;
2408
2409	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
2410	}
2411
2412	template<class P_numtype1, class P_struct1>
2413	inline
2414	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
2415	_bz_nextafter<P_numtype1,double> > >
2416	nextafter(const Matrix<P_numtype1, P_struct1>& d1, double d2)
2417	{
2418	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
2419	_bz_nextafter<P_numtype1,double> > T_expr;
2420
2421	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
2422	}
2423
2424	template<class P_numtype1, class P_struct1>
2425	inline
2426	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
2427	_bz_nextafter<P_numtype1,long double> > >
2428	nextafter(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
2429	{
2430	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
2431	_bz_nextafter<P_numtype1,long double> > T_expr;
2432
2433	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
2434	}
2435
2436	template<class P_numtype1, class P_struct1, class T2>
2437	inline
2438	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
2439	_bz_nextafter<P_numtype1,complex<T2> > > >
2440	nextafter(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
2441	{
2442	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
2443	_bz_nextafter<P_numtype1,complex<T2> > > T_expr;
2444
2445	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
2446	}
2447
2448	template<class P_expr1, class P_numtype2, class P_struct2>
2449	inline
2450	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2451	_bz_nextafter<typename P_expr1::T_numtype,P_numtype2> > >
2452	nextafter(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
2453	{
2454	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2455	_bz_nextafter<typename P_expr1::T_numtype,P_numtype2> > T_expr;
2456
2457	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
2458	}
2459
2460	template<class P_expr1, class P_expr2>
2461	inline
2462	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
2463	_bz_nextafter<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
2464	nextafter(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
2465	{
2466	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
2467	_bz_nextafter<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
2468
2469	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
2470	}
2471
2472	template<class P_expr1>
2473	inline
2474	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
2475	_bz_nextafter<typename P_expr1::T_numtype,int> > >
2476	nextafter(_bz_MatExpr<P_expr1> d1, int d2)
2477	{
2478	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
2479	_bz_nextafter<typename P_expr1::T_numtype,int> > T_expr;
2480
2481	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
2482	}
2483
2484	template<class P_expr1>
2485	inline
2486	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
2487	_bz_nextafter<typename P_expr1::T_numtype,float> > >
2488	nextafter(_bz_MatExpr<P_expr1> d1, float d2)
2489	{
2490	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
2491	_bz_nextafter<typename P_expr1::T_numtype,float> > T_expr;
2492
2493	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
2494	}
2495
2496	template<class P_expr1>
2497	inline
2498	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
2499	_bz_nextafter<typename P_expr1::T_numtype,double> > >
2500	nextafter(_bz_MatExpr<P_expr1> d1, double d2)
2501	{
2502	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
2503	_bz_nextafter<typename P_expr1::T_numtype,double> > T_expr;
2504
2505	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
2506	}
2507
2508	template<class P_expr1>
2509	inline
2510	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
2511	_bz_nextafter<typename P_expr1::T_numtype,long double> > >
2512	nextafter(_bz_MatExpr<P_expr1> d1, long double d2)
2513	{
2514	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
2515	_bz_nextafter<typename P_expr1::T_numtype,long double> > T_expr;
2516
2517	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
2518	}
2519
2520	template<class P_expr1, class T2>
2521	inline
2522	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
2523	_bz_nextafter<typename P_expr1::T_numtype,complex<T2> > > >
2524	nextafter(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
2525	{
2526	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
2527	_bz_nextafter<typename P_expr1::T_numtype,complex<T2> > > T_expr;
2528
2529	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
2530	}
2531
2532	template<class P_numtype2, class P_struct2>
2533	inline
2534	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
2535	_bz_nextafter<int,P_numtype2> > >
2536	nextafter(int d1, const Matrix<P_numtype2, P_struct2>& d2)
2537	{
2538	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
2539	_bz_nextafter<int,P_numtype2> > T_expr;
2540
2541	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
2542	}
2543
2544	template<class P_expr2>
2545	inline
2546	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
2547	_bz_nextafter<int,typename P_expr2::T_numtype> > >
2548	nextafter(int d1, _bz_MatExpr<P_expr2> d2)
2549	{
2550	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
2551	_bz_nextafter<int,typename P_expr2::T_numtype> > T_expr;
2552
2553	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
2554	}
2555
2556	template<class P_numtype2, class P_struct2>
2557	inline
2558	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
2559	_bz_nextafter<float,P_numtype2> > >
2560	nextafter(float d1, const Matrix<P_numtype2, P_struct2>& d2)
2561	{
2562	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
2563	_bz_nextafter<float,P_numtype2> > T_expr;
2564
2565	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
2566	}
2567
2568	template<class P_expr2>
2569	inline
2570	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
2571	_bz_nextafter<float,typename P_expr2::T_numtype> > >
2572	nextafter(float d1, _bz_MatExpr<P_expr2> d2)
2573	{
2574	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
2575	_bz_nextafter<float,typename P_expr2::T_numtype> > T_expr;
2576
2577	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
2578	}
2579
2580	template<class P_numtype2, class P_struct2>
2581	inline
2582	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2583	_bz_nextafter<double,P_numtype2> > >
2584	nextafter(double d1, const Matrix<P_numtype2, P_struct2>& d2)
2585	{
2586	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2587	_bz_nextafter<double,P_numtype2> > T_expr;
2588
2589	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
2590	}
2591
2592	template<class P_expr2>
2593	inline
2594	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
2595	_bz_nextafter<double,typename P_expr2::T_numtype> > >
2596	nextafter(double d1, _bz_MatExpr<P_expr2> d2)
2597	{
2598	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
2599	_bz_nextafter<double,typename P_expr2::T_numtype> > T_expr;
2600
2601	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
2602	}
2603
2604	template<class P_numtype2, class P_struct2>
2605	inline
2606	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2607	_bz_nextafter<long double,P_numtype2> > >
2608	nextafter(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
2609	{
2610	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2611	_bz_nextafter<long double,P_numtype2> > T_expr;
2612
2613	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
2614	}
2615
2616	template<class P_expr2>
2617	inline
2618	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
2619	_bz_nextafter<long double,typename P_expr2::T_numtype> > >
2620	nextafter(long double d1, _bz_MatExpr<P_expr2> d2)
2621	{
2622	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
2623	_bz_nextafter<long double,typename P_expr2::T_numtype> > T_expr;
2624
2625	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
2626	}
2627
2628	template<class T1, class P_numtype2, class P_struct2>
2629	inline
2630	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
2631	_bz_nextafter<complex<T1> ,P_numtype2> > >
2632	nextafter(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
2633	{
2634	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
2635	_bz_nextafter<complex<T1> ,P_numtype2> > T_expr;
2636
2637	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
2638	}
2639
2640	template<class T1, class P_expr2>
2641	inline
2642	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2643	_bz_nextafter<complex<T1> ,typename P_expr2::T_numtype> > >
2644	nextafter(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
2645	{
2646	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2647	_bz_nextafter<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
2648
2649	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
2650	}
2651
2652	#endif
2653
2654	/****************************************************************************
2655	* pow
2656	****************************************************************************/
2657
2658	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
2659	inline
2660	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2661	_bz_pow<P_numtype1,P_numtype2> > >
2662	pow(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
2663	{
2664	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2665	_bz_pow<P_numtype1,P_numtype2> > T_expr;
2666
2667	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
2668	}
2669
2670	template<class P_numtype1, class P_struct1, class P_expr2>
2671	inline
2672	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2673	_bz_pow<P_numtype1,typename P_expr2::T_numtype> > >
2674	pow(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
2675	{
2676	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2677	_bz_pow<P_numtype1,typename P_expr2::T_numtype> > T_expr;
2678
2679	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
2680	}
2681
2682	template<class P_numtype1, class P_struct1>
2683	inline
2684	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2685	_bz_pow<P_numtype1,int> > >
2686	pow(const Matrix<P_numtype1, P_struct1>& d1, int d2)
2687	{
2688	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2689	_bz_pow<P_numtype1,int> > T_expr;
2690
2691	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
2692	}
2693
2694	template<class P_numtype1, class P_struct1>
2695	inline
2696	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2697	_bz_pow<P_numtype1,float> > >
2698	pow(const Matrix<P_numtype1, P_struct1>& d1, float d2)
2699	{
2700	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2701	_bz_pow<P_numtype1,float> > T_expr;
2702
2703	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
2704	}
2705
2706	template<class P_numtype1, class P_struct1>
2707	inline
2708	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
2709	_bz_pow<P_numtype1,double> > >
2710	pow(const Matrix<P_numtype1, P_struct1>& d1, double d2)
2711	{
2712	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
2713	_bz_pow<P_numtype1,double> > T_expr;
2714
2715	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
2716	}
2717
2718	template<class P_numtype1, class P_struct1>
2719	inline
2720	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
2721	_bz_pow<P_numtype1,long double> > >
2722	pow(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
2723	{
2724	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
2725	_bz_pow<P_numtype1,long double> > T_expr;
2726
2727	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
2728	}
2729
2730	template<class P_numtype1, class P_struct1, class T2>
2731	inline
2732	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
2733	_bz_pow<P_numtype1,complex<T2> > > >
2734	pow(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
2735	{
2736	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
2737	_bz_pow<P_numtype1,complex<T2> > > T_expr;
2738
2739	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
2740	}
2741
2742	template<class P_expr1, class P_numtype2, class P_struct2>
2743	inline
2744	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2745	_bz_pow<typename P_expr1::T_numtype,P_numtype2> > >
2746	pow(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
2747	{
2748	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2749	_bz_pow<typename P_expr1::T_numtype,P_numtype2> > T_expr;
2750
2751	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
2752	}
2753
2754	template<class P_expr1, class P_expr2>
2755	inline
2756	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
2757	_bz_pow<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
2758	pow(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
2759	{
2760	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
2761	_bz_pow<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
2762
2763	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
2764	}
2765
2766	template<class P_expr1>
2767	inline
2768	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
2769	_bz_pow<typename P_expr1::T_numtype,int> > >
2770	pow(_bz_MatExpr<P_expr1> d1, int d2)
2771	{
2772	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
2773	_bz_pow<typename P_expr1::T_numtype,int> > T_expr;
2774
2775	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
2776	}
2777
2778	template<class P_expr1>
2779	inline
2780	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
2781	_bz_pow<typename P_expr1::T_numtype,float> > >
2782	pow(_bz_MatExpr<P_expr1> d1, float d2)
2783	{
2784	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
2785	_bz_pow<typename P_expr1::T_numtype,float> > T_expr;
2786
2787	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
2788	}
2789
2790	template<class P_expr1>
2791	inline
2792	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
2793	_bz_pow<typename P_expr1::T_numtype,double> > >
2794	pow(_bz_MatExpr<P_expr1> d1, double d2)
2795	{
2796	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
2797	_bz_pow<typename P_expr1::T_numtype,double> > T_expr;
2798
2799	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
2800	}
2801
2802	template<class P_expr1>
2803	inline
2804	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
2805	_bz_pow<typename P_expr1::T_numtype,long double> > >
2806	pow(_bz_MatExpr<P_expr1> d1, long double d2)
2807	{
2808	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
2809	_bz_pow<typename P_expr1::T_numtype,long double> > T_expr;
2810
2811	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
2812	}
2813
2814	template<class P_expr1, class T2>
2815	inline
2816	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
2817	_bz_pow<typename P_expr1::T_numtype,complex<T2> > > >
2818	pow(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
2819	{
2820	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
2821	_bz_pow<typename P_expr1::T_numtype,complex<T2> > > T_expr;
2822
2823	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
2824	}
2825
2826	template<class P_numtype2, class P_struct2>
2827	inline
2828	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
2829	_bz_pow<int,P_numtype2> > >
2830	pow(int d1, const Matrix<P_numtype2, P_struct2>& d2)
2831	{
2832	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
2833	_bz_pow<int,P_numtype2> > T_expr;
2834
2835	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
2836	}
2837
2838	template<class P_expr2>
2839	inline
2840	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
2841	_bz_pow<int,typename P_expr2::T_numtype> > >
2842	pow(int d1, _bz_MatExpr<P_expr2> d2)
2843	{
2844	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
2845	_bz_pow<int,typename P_expr2::T_numtype> > T_expr;
2846
2847	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
2848	}
2849
2850	template<class P_numtype2, class P_struct2>
2851	inline
2852	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
2853	_bz_pow<float,P_numtype2> > >
2854	pow(float d1, const Matrix<P_numtype2, P_struct2>& d2)
2855	{
2856	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
2857	_bz_pow<float,P_numtype2> > T_expr;
2858
2859	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
2860	}
2861
2862	template<class P_expr2>
2863	inline
2864	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
2865	_bz_pow<float,typename P_expr2::T_numtype> > >
2866	pow(float d1, _bz_MatExpr<P_expr2> d2)
2867	{
2868	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
2869	_bz_pow<float,typename P_expr2::T_numtype> > T_expr;
2870
2871	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
2872	}
2873
2874	template<class P_numtype2, class P_struct2>
2875	inline
2876	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2877	_bz_pow<double,P_numtype2> > >
2878	pow(double d1, const Matrix<P_numtype2, P_struct2>& d2)
2879	{
2880	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2881	_bz_pow<double,P_numtype2> > T_expr;
2882
2883	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
2884	}
2885
2886	template<class P_expr2>
2887	inline
2888	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
2889	_bz_pow<double,typename P_expr2::T_numtype> > >
2890	pow(double d1, _bz_MatExpr<P_expr2> d2)
2891	{
2892	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
2893	_bz_pow<double,typename P_expr2::T_numtype> > T_expr;
2894
2895	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
2896	}
2897
2898	template<class P_numtype2, class P_struct2>
2899	inline
2900	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2901	_bz_pow<long double,P_numtype2> > >
2902	pow(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
2903	{
2904	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
2905	_bz_pow<long double,P_numtype2> > T_expr;
2906
2907	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
2908	}
2909
2910	template<class P_expr2>
2911	inline
2912	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
2913	_bz_pow<long double,typename P_expr2::T_numtype> > >
2914	pow(long double d1, _bz_MatExpr<P_expr2> d2)
2915	{
2916	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
2917	_bz_pow<long double,typename P_expr2::T_numtype> > T_expr;
2918
2919	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
2920	}
2921
2922	template<class T1, class P_numtype2, class P_struct2>
2923	inline
2924	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
2925	_bz_pow<complex<T1> ,P_numtype2> > >
2926	pow(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
2927	{
2928	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
2929	_bz_pow<complex<T1> ,P_numtype2> > T_expr;
2930
2931	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
2932	}
2933
2934	template<class T1, class P_expr2>
2935	inline
2936	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2937	_bz_pow<complex<T1> ,typename P_expr2::T_numtype> > >
2938	pow(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
2939	{
2940	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
2941	_bz_pow<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
2942
2943	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
2944	}
2945
2946
2947	/****************************************************************************
2948	* remainder
2949	****************************************************************************/
2950
2951	#ifdef BZ_HAVE_SYSTEM_V_MATH
2952	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
2953	inline
2954	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2955	_bz_remainder<P_numtype1,P_numtype2> > >
2956	remainder(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
2957	{
2958	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
2959	_bz_remainder<P_numtype1,P_numtype2> > T_expr;
2960
2961	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
2962	}
2963
2964	template<class P_numtype1, class P_struct1, class P_expr2>
2965	inline
2966	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2967	_bz_remainder<P_numtype1,typename P_expr2::T_numtype> > >
2968	remainder(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
2969	{
2970	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
2971	_bz_remainder<P_numtype1,typename P_expr2::T_numtype> > T_expr;
2972
2973	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
2974	}
2975
2976	template<class P_numtype1, class P_struct1>
2977	inline
2978	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2979	_bz_remainder<P_numtype1,int> > >
2980	remainder(const Matrix<P_numtype1, P_struct1>& d1, int d2)
2981	{
2982	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
2983	_bz_remainder<P_numtype1,int> > T_expr;
2984
2985	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
2986	}
2987
2988	template<class P_numtype1, class P_struct1>
2989	inline
2990	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2991	_bz_remainder<P_numtype1,float> > >
2992	remainder(const Matrix<P_numtype1, P_struct1>& d1, float d2)
2993	{
2994	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
2995	_bz_remainder<P_numtype1,float> > T_expr;
2996
2997	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
2998	}
2999
3000	template<class P_numtype1, class P_struct1>
3001	inline
3002	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3003	_bz_remainder<P_numtype1,double> > >
3004	remainder(const Matrix<P_numtype1, P_struct1>& d1, double d2)
3005	{
3006	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3007	_bz_remainder<P_numtype1,double> > T_expr;
3008
3009	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
3010	}
3011
3012	template<class P_numtype1, class P_struct1>
3013	inline
3014	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3015	_bz_remainder<P_numtype1,long double> > >
3016	remainder(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
3017	{
3018	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3019	_bz_remainder<P_numtype1,long double> > T_expr;
3020
3021	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
3022	}
3023
3024	template<class P_numtype1, class P_struct1, class T2>
3025	inline
3026	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3027	_bz_remainder<P_numtype1,complex<T2> > > >
3028	remainder(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
3029	{
3030	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3031	_bz_remainder<P_numtype1,complex<T2> > > T_expr;
3032
3033	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
3034	}
3035
3036	template<class P_expr1, class P_numtype2, class P_struct2>
3037	inline
3038	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3039	_bz_remainder<typename P_expr1::T_numtype,P_numtype2> > >
3040	remainder(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
3041	{
3042	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3043	_bz_remainder<typename P_expr1::T_numtype,P_numtype2> > T_expr;
3044
3045	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
3046	}
3047
3048	template<class P_expr1, class P_expr2>
3049	inline
3050	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3051	_bz_remainder<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
3052	remainder(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
3053	{
3054	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3055	_bz_remainder<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
3056
3057	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
3058	}
3059
3060	template<class P_expr1>
3061	inline
3062	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3063	_bz_remainder<typename P_expr1::T_numtype,int> > >
3064	remainder(_bz_MatExpr<P_expr1> d1, int d2)
3065	{
3066	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3067	_bz_remainder<typename P_expr1::T_numtype,int> > T_expr;
3068
3069	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
3070	}
3071
3072	template<class P_expr1>
3073	inline
3074	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3075	_bz_remainder<typename P_expr1::T_numtype,float> > >
3076	remainder(_bz_MatExpr<P_expr1> d1, float d2)
3077	{
3078	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3079	_bz_remainder<typename P_expr1::T_numtype,float> > T_expr;
3080
3081	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
3082	}
3083
3084	template<class P_expr1>
3085	inline
3086	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3087	_bz_remainder<typename P_expr1::T_numtype,double> > >
3088	remainder(_bz_MatExpr<P_expr1> d1, double d2)
3089	{
3090	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3091	_bz_remainder<typename P_expr1::T_numtype,double> > T_expr;
3092
3093	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
3094	}
3095
3096	template<class P_expr1>
3097	inline
3098	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3099	_bz_remainder<typename P_expr1::T_numtype,long double> > >
3100	remainder(_bz_MatExpr<P_expr1> d1, long double d2)
3101	{
3102	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3103	_bz_remainder<typename P_expr1::T_numtype,long double> > T_expr;
3104
3105	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
3106	}
3107
3108	template<class P_expr1, class T2>
3109	inline
3110	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3111	_bz_remainder<typename P_expr1::T_numtype,complex<T2> > > >
3112	remainder(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
3113	{
3114	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3115	_bz_remainder<typename P_expr1::T_numtype,complex<T2> > > T_expr;
3116
3117	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
3118	}
3119
3120	template<class P_numtype2, class P_struct2>
3121	inline
3122	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3123	_bz_remainder<int,P_numtype2> > >
3124	remainder(int d1, const Matrix<P_numtype2, P_struct2>& d2)
3125	{
3126	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3127	_bz_remainder<int,P_numtype2> > T_expr;
3128
3129	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
3130	}
3131
3132	template<class P_expr2>
3133	inline
3134	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3135	_bz_remainder<int,typename P_expr2::T_numtype> > >
3136	remainder(int d1, _bz_MatExpr<P_expr2> d2)
3137	{
3138	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3139	_bz_remainder<int,typename P_expr2::T_numtype> > T_expr;
3140
3141	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
3142	}
3143
3144	template<class P_numtype2, class P_struct2>
3145	inline
3146	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
3147	_bz_remainder<float,P_numtype2> > >
3148	remainder(float d1, const Matrix<P_numtype2, P_struct2>& d2)
3149	{
3150	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
3151	_bz_remainder<float,P_numtype2> > T_expr;
3152
3153	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
3154	}
3155
3156	template<class P_expr2>
3157	inline
3158	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
3159	_bz_remainder<float,typename P_expr2::T_numtype> > >
3160	remainder(float d1, _bz_MatExpr<P_expr2> d2)
3161	{
3162	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
3163	_bz_remainder<float,typename P_expr2::T_numtype> > T_expr;
3164
3165	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
3166	}
3167
3168	template<class P_numtype2, class P_struct2>
3169	inline
3170	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3171	_bz_remainder<double,P_numtype2> > >
3172	remainder(double d1, const Matrix<P_numtype2, P_struct2>& d2)
3173	{
3174	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3175	_bz_remainder<double,P_numtype2> > T_expr;
3176
3177	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
3178	}
3179
3180	template<class P_expr2>
3181	inline
3182	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
3183	_bz_remainder<double,typename P_expr2::T_numtype> > >
3184	remainder(double d1, _bz_MatExpr<P_expr2> d2)
3185	{
3186	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
3187	_bz_remainder<double,typename P_expr2::T_numtype> > T_expr;
3188
3189	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
3190	}
3191
3192	template<class P_numtype2, class P_struct2>
3193	inline
3194	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3195	_bz_remainder<long double,P_numtype2> > >
3196	remainder(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
3197	{
3198	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3199	_bz_remainder<long double,P_numtype2> > T_expr;
3200
3201	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
3202	}
3203
3204	template<class P_expr2>
3205	inline
3206	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
3207	_bz_remainder<long double,typename P_expr2::T_numtype> > >
3208	remainder(long double d1, _bz_MatExpr<P_expr2> d2)
3209	{
3210	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
3211	_bz_remainder<long double,typename P_expr2::T_numtype> > T_expr;
3212
3213	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
3214	}
3215
3216	template<class T1, class P_numtype2, class P_struct2>
3217	inline
3218	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
3219	_bz_remainder<complex<T1> ,P_numtype2> > >
3220	remainder(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
3221	{
3222	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
3223	_bz_remainder<complex<T1> ,P_numtype2> > T_expr;
3224
3225	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
3226	}
3227
3228	template<class T1, class P_expr2>
3229	inline
3230	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
3231	_bz_remainder<complex<T1> ,typename P_expr2::T_numtype> > >
3232	remainder(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
3233	{
3234	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
3235	_bz_remainder<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
3236
3237	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
3238	}
3239
3240	#endif
3241
3242	/****************************************************************************
3243	* rint
3244	****************************************************************************/
3245
3246	#ifdef BZ_HAVE_IEEE_MATH
3247	template<class P_numtype1, class P_struct1>
3248	inline
3249	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3250	_bz_rint<P_numtype1> > >
3251	rint(const Matrix<P_numtype1, P_struct1>& d1)
3252	{
3253	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3254	_bz_rint<P_numtype1> > T_expr;
3255
3256	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3257	}
3258
3259	template<class P_expr1>
3260	inline
3261	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3262	_bz_rint<typename P_expr1::T_numtype> > >
3263	rint(_bz_MatExpr<P_expr1> d1)
3264	{
3265	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3266	_bz_rint<typename P_expr1::T_numtype> > T_expr;
3267
3268	return _bz_MatExpr<T_expr>(T_expr(d1));
3269	}
3270
3271	#endif
3272
3273	/****************************************************************************
3274	* rsqrt
3275	****************************************************************************/
3276
3277	#ifdef BZ_HAVE_SYSTEM_V_MATH
3278	template<class P_numtype1, class P_struct1>
3279	inline
3280	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3281	_bz_rsqrt<P_numtype1> > >
3282	rsqrt(const Matrix<P_numtype1, P_struct1>& d1)
3283	{
3284	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3285	_bz_rsqrt<P_numtype1> > T_expr;
3286
3287	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3288	}
3289
3290	template<class P_expr1>
3291	inline
3292	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3293	_bz_rsqrt<typename P_expr1::T_numtype> > >
3294	rsqrt(_bz_MatExpr<P_expr1> d1)
3295	{
3296	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3297	_bz_rsqrt<typename P_expr1::T_numtype> > T_expr;
3298
3299	return _bz_MatExpr<T_expr>(T_expr(d1));
3300	}
3301
3302	#endif
3303
3304	/****************************************************************************
3305	* scalb
3306	****************************************************************************/
3307
3308	#ifdef BZ_HAVE_SYSTEM_V_MATH
3309	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
3310	inline
3311	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3312	_bz_scalb<P_numtype1,P_numtype2> > >
3313	scalb(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
3314	{
3315	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3316	_bz_scalb<P_numtype1,P_numtype2> > T_expr;
3317
3318	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
3319	}
3320
3321	template<class P_numtype1, class P_struct1, class P_expr2>
3322	inline
3323	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
3324	_bz_scalb<P_numtype1,typename P_expr2::T_numtype> > >
3325	scalb(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
3326	{
3327	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
3328	_bz_scalb<P_numtype1,typename P_expr2::T_numtype> > T_expr;
3329
3330	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
3331	}
3332
3333	template<class P_numtype1, class P_struct1>
3334	inline
3335	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
3336	_bz_scalb<P_numtype1,int> > >
3337	scalb(const Matrix<P_numtype1, P_struct1>& d1, int d2)
3338	{
3339	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
3340	_bz_scalb<P_numtype1,int> > T_expr;
3341
3342	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
3343	}
3344
3345	template<class P_numtype1, class P_struct1>
3346	inline
3347	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
3348	_bz_scalb<P_numtype1,float> > >
3349	scalb(const Matrix<P_numtype1, P_struct1>& d1, float d2)
3350	{
3351	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
3352	_bz_scalb<P_numtype1,float> > T_expr;
3353
3354	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
3355	}
3356
3357	template<class P_numtype1, class P_struct1>
3358	inline
3359	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3360	_bz_scalb<P_numtype1,double> > >
3361	scalb(const Matrix<P_numtype1, P_struct1>& d1, double d2)
3362	{
3363	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3364	_bz_scalb<P_numtype1,double> > T_expr;
3365
3366	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
3367	}
3368
3369	template<class P_numtype1, class P_struct1>
3370	inline
3371	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3372	_bz_scalb<P_numtype1,long double> > >
3373	scalb(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
3374	{
3375	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3376	_bz_scalb<P_numtype1,long double> > T_expr;
3377
3378	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
3379	}
3380
3381	template<class P_numtype1, class P_struct1, class T2>
3382	inline
3383	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3384	_bz_scalb<P_numtype1,complex<T2> > > >
3385	scalb(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
3386	{
3387	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3388	_bz_scalb<P_numtype1,complex<T2> > > T_expr;
3389
3390	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
3391	}
3392
3393	template<class P_expr1, class P_numtype2, class P_struct2>
3394	inline
3395	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3396	_bz_scalb<typename P_expr1::T_numtype,P_numtype2> > >
3397	scalb(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
3398	{
3399	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3400	_bz_scalb<typename P_expr1::T_numtype,P_numtype2> > T_expr;
3401
3402	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
3403	}
3404
3405	template<class P_expr1, class P_expr2>
3406	inline
3407	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3408	_bz_scalb<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
3409	scalb(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
3410	{
3411	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3412	_bz_scalb<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
3413
3414	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
3415	}
3416
3417	template<class P_expr1>
3418	inline
3419	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3420	_bz_scalb<typename P_expr1::T_numtype,int> > >
3421	scalb(_bz_MatExpr<P_expr1> d1, int d2)
3422	{
3423	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3424	_bz_scalb<typename P_expr1::T_numtype,int> > T_expr;
3425
3426	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
3427	}
3428
3429	template<class P_expr1>
3430	inline
3431	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3432	_bz_scalb<typename P_expr1::T_numtype,float> > >
3433	scalb(_bz_MatExpr<P_expr1> d1, float d2)
3434	{
3435	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3436	_bz_scalb<typename P_expr1::T_numtype,float> > T_expr;
3437
3438	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
3439	}
3440
3441	template<class P_expr1>
3442	inline
3443	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3444	_bz_scalb<typename P_expr1::T_numtype,double> > >
3445	scalb(_bz_MatExpr<P_expr1> d1, double d2)
3446	{
3447	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3448	_bz_scalb<typename P_expr1::T_numtype,double> > T_expr;
3449
3450	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
3451	}
3452
3453	template<class P_expr1>
3454	inline
3455	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3456	_bz_scalb<typename P_expr1::T_numtype,long double> > >
3457	scalb(_bz_MatExpr<P_expr1> d1, long double d2)
3458	{
3459	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3460	_bz_scalb<typename P_expr1::T_numtype,long double> > T_expr;
3461
3462	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
3463	}
3464
3465	template<class P_expr1, class T2>
3466	inline
3467	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3468	_bz_scalb<typename P_expr1::T_numtype,complex<T2> > > >
3469	scalb(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
3470	{
3471	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3472	_bz_scalb<typename P_expr1::T_numtype,complex<T2> > > T_expr;
3473
3474	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
3475	}
3476
3477	template<class P_numtype2, class P_struct2>
3478	inline
3479	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3480	_bz_scalb<int,P_numtype2> > >
3481	scalb(int d1, const Matrix<P_numtype2, P_struct2>& d2)
3482	{
3483	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3484	_bz_scalb<int,P_numtype2> > T_expr;
3485
3486	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
3487	}
3488
3489	template<class P_expr2>
3490	inline
3491	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3492	_bz_scalb<int,typename P_expr2::T_numtype> > >
3493	scalb(int d1, _bz_MatExpr<P_expr2> d2)
3494	{
3495	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3496	_bz_scalb<int,typename P_expr2::T_numtype> > T_expr;
3497
3498	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
3499	}
3500
3501	template<class P_numtype2, class P_struct2>
3502	inline
3503	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
3504	_bz_scalb<float,P_numtype2> > >
3505	scalb(float d1, const Matrix<P_numtype2, P_struct2>& d2)
3506	{
3507	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
3508	_bz_scalb<float,P_numtype2> > T_expr;
3509
3510	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
3511	}
3512
3513	template<class P_expr2>
3514	inline
3515	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
3516	_bz_scalb<float,typename P_expr2::T_numtype> > >
3517	scalb(float d1, _bz_MatExpr<P_expr2> d2)
3518	{
3519	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
3520	_bz_scalb<float,typename P_expr2::T_numtype> > T_expr;
3521
3522	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
3523	}
3524
3525	template<class P_numtype2, class P_struct2>
3526	inline
3527	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3528	_bz_scalb<double,P_numtype2> > >
3529	scalb(double d1, const Matrix<P_numtype2, P_struct2>& d2)
3530	{
3531	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3532	_bz_scalb<double,P_numtype2> > T_expr;
3533
3534	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
3535	}
3536
3537	template<class P_expr2>
3538	inline
3539	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
3540	_bz_scalb<double,typename P_expr2::T_numtype> > >
3541	scalb(double d1, _bz_MatExpr<P_expr2> d2)
3542	{
3543	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
3544	_bz_scalb<double,typename P_expr2::T_numtype> > T_expr;
3545
3546	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
3547	}
3548
3549	template<class P_numtype2, class P_struct2>
3550	inline
3551	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3552	_bz_scalb<long double,P_numtype2> > >
3553	scalb(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
3554	{
3555	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
3556	_bz_scalb<long double,P_numtype2> > T_expr;
3557
3558	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
3559	}
3560
3561	template<class P_expr2>
3562	inline
3563	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
3564	_bz_scalb<long double,typename P_expr2::T_numtype> > >
3565	scalb(long double d1, _bz_MatExpr<P_expr2> d2)
3566	{
3567	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
3568	_bz_scalb<long double,typename P_expr2::T_numtype> > T_expr;
3569
3570	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
3571	}
3572
3573	template<class T1, class P_numtype2, class P_struct2>
3574	inline
3575	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
3576	_bz_scalb<complex<T1> ,P_numtype2> > >
3577	scalb(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
3578	{
3579	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
3580	_bz_scalb<complex<T1> ,P_numtype2> > T_expr;
3581
3582	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
3583	}
3584
3585	template<class T1, class P_expr2>
3586	inline
3587	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
3588	_bz_scalb<complex<T1> ,typename P_expr2::T_numtype> > >
3589	scalb(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
3590	{
3591	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
3592	_bz_scalb<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
3593
3594	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
3595	}
3596
3597	#endif
3598
3599	/****************************************************************************
3600	* sin
3601	****************************************************************************/
3602
3603	template<class P_numtype1, class P_struct1>
3604	inline
3605	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3606	_bz_sin<P_numtype1> > >
3607	sin(const Matrix<P_numtype1, P_struct1>& d1)
3608	{
3609	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3610	_bz_sin<P_numtype1> > T_expr;
3611
3612	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3613	}
3614
3615	template<class P_expr1>
3616	inline
3617	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3618	_bz_sin<typename P_expr1::T_numtype> > >
3619	sin(_bz_MatExpr<P_expr1> d1)
3620	{
3621	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3622	_bz_sin<typename P_expr1::T_numtype> > T_expr;
3623
3624	return _bz_MatExpr<T_expr>(T_expr(d1));
3625	}
3626
3627
3628	/****************************************************************************
3629	* sinh
3630	****************************************************************************/
3631
3632	template<class P_numtype1, class P_struct1>
3633	inline
3634	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3635	_bz_sinh<P_numtype1> > >
3636	sinh(const Matrix<P_numtype1, P_struct1>& d1)
3637	{
3638	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3639	_bz_sinh<P_numtype1> > T_expr;
3640
3641	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3642	}
3643
3644	template<class P_expr1>
3645	inline
3646	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3647	_bz_sinh<typename P_expr1::T_numtype> > >
3648	sinh(_bz_MatExpr<P_expr1> d1)
3649	{
3650	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3651	_bz_sinh<typename P_expr1::T_numtype> > T_expr;
3652
3653	return _bz_MatExpr<T_expr>(T_expr(d1));
3654	}
3655
3656
3657	/****************************************************************************
3658	* sqr
3659	****************************************************************************/
3660
3661	template<class P_numtype1, class P_struct1>
3662	inline
3663	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3664	_bz_sqr<P_numtype1> > >
3665	sqr(const Matrix<P_numtype1, P_struct1>& d1)
3666	{
3667	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3668	_bz_sqr<P_numtype1> > T_expr;
3669
3670	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3671	}
3672
3673	template<class P_expr1>
3674	inline
3675	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3676	_bz_sqr<typename P_expr1::T_numtype> > >
3677	sqr(_bz_MatExpr<P_expr1> d1)
3678	{
3679	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3680	_bz_sqr<typename P_expr1::T_numtype> > T_expr;
3681
3682	return _bz_MatExpr<T_expr>(T_expr(d1));
3683	}
3684
3685
3686	/****************************************************************************
3687	* sqrt
3688	****************************************************************************/
3689
3690	template<class P_numtype1, class P_struct1>
3691	inline
3692	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3693	_bz_sqrt<P_numtype1> > >
3694	sqrt(const Matrix<P_numtype1, P_struct1>& d1)
3695	{
3696	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3697	_bz_sqrt<P_numtype1> > T_expr;
3698
3699	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3700	}
3701
3702	template<class P_expr1>
3703	inline
3704	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3705	_bz_sqrt<typename P_expr1::T_numtype> > >
3706	sqrt(_bz_MatExpr<P_expr1> d1)
3707	{
3708	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3709	_bz_sqrt<typename P_expr1::T_numtype> > T_expr;
3710
3711	return _bz_MatExpr<T_expr>(T_expr(d1));
3712	}
3713
3714
3715	/****************************************************************************
3716	* tan
3717	****************************************************************************/
3718
3719	template<class P_numtype1, class P_struct1>
3720	inline
3721	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3722	_bz_tan<P_numtype1> > >
3723	tan(const Matrix<P_numtype1, P_struct1>& d1)
3724	{
3725	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3726	_bz_tan<P_numtype1> > T_expr;
3727
3728	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3729	}
3730
3731	template<class P_expr1>
3732	inline
3733	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3734	_bz_tan<typename P_expr1::T_numtype> > >
3735	tan(_bz_MatExpr<P_expr1> d1)
3736	{
3737	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3738	_bz_tan<typename P_expr1::T_numtype> > T_expr;
3739
3740	return _bz_MatExpr<T_expr>(T_expr(d1));
3741	}
3742
3743
3744	/****************************************************************************
3745	* tanh
3746	****************************************************************************/
3747
3748	template<class P_numtype1, class P_struct1>
3749	inline
3750	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3751	_bz_tanh<P_numtype1> > >
3752	tanh(const Matrix<P_numtype1, P_struct1>& d1)
3753	{
3754	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3755	_bz_tanh<P_numtype1> > T_expr;
3756
3757	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3758	}
3759
3760	template<class P_expr1>
3761	inline
3762	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3763	_bz_tanh<typename P_expr1::T_numtype> > >
3764	tanh(_bz_MatExpr<P_expr1> d1)
3765	{
3766	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3767	_bz_tanh<typename P_expr1::T_numtype> > T_expr;
3768
3769	return _bz_MatExpr<T_expr>(T_expr(d1));
3770	}
3771
3772
3773	/****************************************************************************
3774	* uitrunc
3775	****************************************************************************/
3776
3777	#ifdef BZ_HAVE_SYSTEM_V_MATH
3778	template<class P_numtype1, class P_struct1>
3779	inline
3780	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3781	_bz_uitrunc<P_numtype1> > >
3782	uitrunc(const Matrix<P_numtype1, P_struct1>& d1)
3783	{
3784	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
3785	_bz_uitrunc<P_numtype1> > T_expr;
3786
3787	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
3788	}
3789
3790	template<class P_expr1>
3791	inline
3792	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3793	_bz_uitrunc<typename P_expr1::T_numtype> > >
3794	uitrunc(_bz_MatExpr<P_expr1> d1)
3795	{
3796	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
3797	_bz_uitrunc<typename P_expr1::T_numtype> > T_expr;
3798
3799	return _bz_MatExpr<T_expr>(T_expr(d1));
3800	}
3801
3802	#endif
3803
3804	/****************************************************************************
3805	* unordered
3806	****************************************************************************/
3807
3808	#ifdef BZ_HAVE_SYSTEM_V_MATH
3809	template<class P_numtype1, class P_struct1, class P_numtype2, class P_struct2>
3810	inline
3811	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3812	_bz_unordered<P_numtype1,P_numtype2> > >
3813	unordered(const Matrix<P_numtype1, P_struct1>& d1, const Matrix<P_numtype2, P_struct2>& d2)
3814	{
3815	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3816	_bz_unordered<P_numtype1,P_numtype2> > T_expr;
3817
3818	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2._bz_getRef()));
3819	}
3820
3821	template<class P_numtype1, class P_struct1, class P_expr2>
3822	inline
3823	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
3824	_bz_unordered<P_numtype1,typename P_expr2::T_numtype> > >
3825	unordered(const Matrix<P_numtype1, P_struct1>& d1, _bz_MatExpr<P_expr2> d2)
3826	{
3827	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExpr<P_expr2>,
3828	_bz_unordered<P_numtype1,typename P_expr2::T_numtype> > T_expr;
3829
3830	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), d2));
3831	}
3832
3833	template<class P_numtype1, class P_struct1>
3834	inline
3835	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
3836	_bz_unordered<P_numtype1,int> > >
3837	unordered(const Matrix<P_numtype1, P_struct1>& d1, int d2)
3838	{
3839	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<int>,
3840	_bz_unordered<P_numtype1,int> > T_expr;
3841
3842	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<int>(d2)));
3843	}
3844
3845	template<class P_numtype1, class P_struct1>
3846	inline
3847	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
3848	_bz_unordered<P_numtype1,float> > >
3849	unordered(const Matrix<P_numtype1, P_struct1>& d1, float d2)
3850	{
3851	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<float>,
3852	_bz_unordered<P_numtype1,float> > T_expr;
3853
3854	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<float>(d2)));
3855	}
3856
3857	template<class P_numtype1, class P_struct1>
3858	inline
3859	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3860	_bz_unordered<P_numtype1,double> > >
3861	unordered(const Matrix<P_numtype1, P_struct1>& d1, double d2)
3862	{
3863	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<double>,
3864	_bz_unordered<P_numtype1,double> > T_expr;
3865
3866	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<double>(d2)));
3867	}
3868
3869	template<class P_numtype1, class P_struct1>
3870	inline
3871	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3872	_bz_unordered<P_numtype1,long double> > >
3873	unordered(const Matrix<P_numtype1, P_struct1>& d1, long double d2)
3874	{
3875	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<long double>,
3876	_bz_unordered<P_numtype1,long double> > T_expr;
3877
3878	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<long double>(d2)));
3879	}
3880
3881	template<class P_numtype1, class P_struct1, class T2>
3882	inline
3883	_bz_MatExpr<_bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3884	_bz_unordered<P_numtype1,complex<T2> > > >
3885	unordered(const Matrix<P_numtype1, P_struct1>& d1, complex<T2> d2)
3886	{
3887	typedef _bz_MatExprOp<_bz_MatrixRef<P_numtype1, P_struct1>, _bz_MatExprConstant<complex<T2> > ,
3888	_bz_unordered<P_numtype1,complex<T2> > > T_expr;
3889
3890	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef(), _bz_MatExprConstant<complex<T2> > (d2)));
3891	}
3892
3893	template<class P_expr1, class P_numtype2, class P_struct2>
3894	inline
3895	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3896	_bz_unordered<typename P_expr1::T_numtype,P_numtype2> > >
3897	unordered(_bz_MatExpr<P_expr1> d1, const Matrix<P_numtype2, P_struct2>& d2)
3898	{
3899	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatrixRef<P_numtype2, P_struct2>,
3900	_bz_unordered<typename P_expr1::T_numtype,P_numtype2> > T_expr;
3901
3902	return _bz_MatExpr<T_expr>(T_expr(d1, d2._bz_getRef()));
3903	}
3904
3905	template<class P_expr1, class P_expr2>
3906	inline
3907	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3908	_bz_unordered<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > >
3909	unordered(_bz_MatExpr<P_expr1> d1, _bz_MatExpr<P_expr2> d2)
3910	{
3911	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExpr<P_expr2>,
3912	_bz_unordered<typename P_expr1::T_numtype,typename P_expr2::T_numtype> > T_expr;
3913
3914	return _bz_MatExpr<T_expr>(T_expr(d1, d2));
3915	}
3916
3917	template<class P_expr1>
3918	inline
3919	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3920	_bz_unordered<typename P_expr1::T_numtype,int> > >
3921	unordered(_bz_MatExpr<P_expr1> d1, int d2)
3922	{
3923	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<int>,
3924	_bz_unordered<typename P_expr1::T_numtype,int> > T_expr;
3925
3926	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<int>(d2)));
3927	}
3928
3929	template<class P_expr1>
3930	inline
3931	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3932	_bz_unordered<typename P_expr1::T_numtype,float> > >
3933	unordered(_bz_MatExpr<P_expr1> d1, float d2)
3934	{
3935	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<float>,
3936	_bz_unordered<typename P_expr1::T_numtype,float> > T_expr;
3937
3938	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<float>(d2)));
3939	}
3940
3941	template<class P_expr1>
3942	inline
3943	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3944	_bz_unordered<typename P_expr1::T_numtype,double> > >
3945	unordered(_bz_MatExpr<P_expr1> d1, double d2)
3946	{
3947	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<double>,
3948	_bz_unordered<typename P_expr1::T_numtype,double> > T_expr;
3949
3950	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<double>(d2)));
3951	}
3952
3953	template<class P_expr1>
3954	inline
3955	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3956	_bz_unordered<typename P_expr1::T_numtype,long double> > >
3957	unordered(_bz_MatExpr<P_expr1> d1, long double d2)
3958	{
3959	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<long double>,
3960	_bz_unordered<typename P_expr1::T_numtype,long double> > T_expr;
3961
3962	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<long double>(d2)));
3963	}
3964
3965	template<class P_expr1, class T2>
3966	inline
3967	_bz_MatExpr<_bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3968	_bz_unordered<typename P_expr1::T_numtype,complex<T2> > > >
3969	unordered(_bz_MatExpr<P_expr1> d1, complex<T2> d2)
3970	{
3971	typedef _bz_MatExprOp<_bz_MatExpr<P_expr1>, _bz_MatExprConstant<complex<T2> > ,
3972	_bz_unordered<typename P_expr1::T_numtype,complex<T2> > > T_expr;
3973
3974	return _bz_MatExpr<T_expr>(T_expr(d1, _bz_MatExprConstant<complex<T2> > (d2)));
3975	}
3976
3977	template<class P_numtype2, class P_struct2>
3978	inline
3979	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3980	_bz_unordered<int,P_numtype2> > >
3981	unordered(int d1, const Matrix<P_numtype2, P_struct2>& d2)
3982	{
3983	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatrixRef<P_numtype2, P_struct2>,
3984	_bz_unordered<int,P_numtype2> > T_expr;
3985
3986	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2._bz_getRef()));
3987	}
3988
3989	template<class P_expr2>
3990	inline
3991	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3992	_bz_unordered<int,typename P_expr2::T_numtype> > >
3993	unordered(int d1, _bz_MatExpr<P_expr2> d2)
3994	{
3995	typedef _bz_MatExprOp<_bz_MatExprConstant<int>, _bz_MatExpr<P_expr2>,
3996	_bz_unordered<int,typename P_expr2::T_numtype> > T_expr;
3997
3998	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<int>(d1), d2));
3999	}
4000
4001	template<class P_numtype2, class P_struct2>
4002	inline
4003	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
4004	_bz_unordered<float,P_numtype2> > >
4005	unordered(float d1, const Matrix<P_numtype2, P_struct2>& d2)
4006	{
4007	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatrixRef<P_numtype2, P_struct2>,
4008	_bz_unordered<float,P_numtype2> > T_expr;
4009
4010	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2._bz_getRef()));
4011	}
4012
4013	template<class P_expr2>
4014	inline
4015	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
4016	_bz_unordered<float,typename P_expr2::T_numtype> > >
4017	unordered(float d1, _bz_MatExpr<P_expr2> d2)
4018	{
4019	typedef _bz_MatExprOp<_bz_MatExprConstant<float>, _bz_MatExpr<P_expr2>,
4020	_bz_unordered<float,typename P_expr2::T_numtype> > T_expr;
4021
4022	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<float>(d1), d2));
4023	}
4024
4025	template<class P_numtype2, class P_struct2>
4026	inline
4027	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
4028	_bz_unordered<double,P_numtype2> > >
4029	unordered(double d1, const Matrix<P_numtype2, P_struct2>& d2)
4030	{
4031	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatrixRef<P_numtype2, P_struct2>,
4032	_bz_unordered<double,P_numtype2> > T_expr;
4033
4034	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2._bz_getRef()));
4035	}
4036
4037	template<class P_expr2>
4038	inline
4039	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
4040	_bz_unordered<double,typename P_expr2::T_numtype> > >
4041	unordered(double d1, _bz_MatExpr<P_expr2> d2)
4042	{
4043	typedef _bz_MatExprOp<_bz_MatExprConstant<double>, _bz_MatExpr<P_expr2>,
4044	_bz_unordered<double,typename P_expr2::T_numtype> > T_expr;
4045
4046	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<double>(d1), d2));
4047	}
4048
4049	template<class P_numtype2, class P_struct2>
4050	inline
4051	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
4052	_bz_unordered<long double,P_numtype2> > >
4053	unordered(long double d1, const Matrix<P_numtype2, P_struct2>& d2)
4054	{
4055	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatrixRef<P_numtype2, P_struct2>,
4056	_bz_unordered<long double,P_numtype2> > T_expr;
4057
4058	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2._bz_getRef()));
4059	}
4060
4061	template<class P_expr2>
4062	inline
4063	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
4064	_bz_unordered<long double,typename P_expr2::T_numtype> > >
4065	unordered(long double d1, _bz_MatExpr<P_expr2> d2)
4066	{
4067	typedef _bz_MatExprOp<_bz_MatExprConstant<long double>, _bz_MatExpr<P_expr2>,
4068	_bz_unordered<long double,typename P_expr2::T_numtype> > T_expr;
4069
4070	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<long double>(d1), d2));
4071	}
4072
4073	template<class T1, class P_numtype2, class P_struct2>
4074	inline
4075	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
4076	_bz_unordered<complex<T1> ,P_numtype2> > >
4077	unordered(complex<T1> d1, const Matrix<P_numtype2, P_struct2>& d2)
4078	{
4079	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatrixRef<P_numtype2, P_struct2>,
4080	_bz_unordered<complex<T1> ,P_numtype2> > T_expr;
4081
4082	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2._bz_getRef()));
4083	}
4084
4085	template<class T1, class P_expr2>
4086	inline
4087	_bz_MatExpr<_bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
4088	_bz_unordered<complex<T1> ,typename P_expr2::T_numtype> > >
4089	unordered(complex<T1> d1, _bz_MatExpr<P_expr2> d2)
4090	{
4091	typedef _bz_MatExprOp<_bz_MatExprConstant<complex<T1> > , _bz_MatExpr<P_expr2>,
4092	_bz_unordered<complex<T1> ,typename P_expr2::T_numtype> > T_expr;
4093
4094	return _bz_MatExpr<T_expr>(T_expr(_bz_MatExprConstant<complex<T1> > (d1), d2));
4095	}
4096
4097	#endif
4098
4099	/****************************************************************************
4100	* y0
4101	****************************************************************************/
4102
4103	#ifdef BZ_HAVE_IEEE_MATH
4104	template<class P_numtype1, class P_struct1>
4105	inline
4106	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
4107	_bz_y0<P_numtype1> > >
4108	y0(const Matrix<P_numtype1, P_struct1>& d1)
4109	{
4110	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
4111	_bz_y0<P_numtype1> > T_expr;
4112
4113	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
4114	}
4115
4116	template<class P_expr1>
4117	inline
4118	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
4119	_bz_y0<typename P_expr1::T_numtype> > >
4120	y0(_bz_MatExpr<P_expr1> d1)
4121	{
4122	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
4123	_bz_y0<typename P_expr1::T_numtype> > T_expr;
4124
4125	return _bz_MatExpr<T_expr>(T_expr(d1));
4126	}
4127
4128	#endif
4129
4130	/****************************************************************************
4131	* y1
4132	****************************************************************************/
4133
4134	#ifdef BZ_HAVE_IEEE_MATH
4135	template<class P_numtype1, class P_struct1>
4136	inline
4137	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
4138	_bz_y1<P_numtype1> > >
4139	y1(const Matrix<P_numtype1, P_struct1>& d1)
4140	{
4141	typedef _bz_MatExprUnaryOp<_bz_MatrixRef<P_numtype1, P_struct1>,
4142	_bz_y1<P_numtype1> > T_expr;
4143
4144	return _bz_MatExpr<T_expr>(T_expr(d1._bz_getRef()));
4145	}
4146
4147	template<class P_expr1>
4148	inline
4149	_bz_MatExpr<_bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
4150	_bz_y1<typename P_expr1::T_numtype> > >
4151	y1(_bz_MatExpr<P_expr1> d1)
4152	{
4153	typedef _bz_MatExprUnaryOp<_bz_MatExpr<P_expr1>,
4154	_bz_y1<typename P_expr1::T_numtype> > T_expr;
4155
4156	return _bz_MatExpr<T_expr>(T_expr(d1));
4157	}
4158
4159	#endif
4160
4161
4162	BZ_NAMESPACE_END
4163
4164	#endif // BZ_MATUOPS_H

Note: See TracBrowser for help on using the repository browser.

Download in other formats: