gcc50: Disconnect from buildworld.
[dragonfly.git] / contrib / gcc-5.0 / libstdc++-v3 / include / bits / valarray_after.h
1 // The template and inlines for the -*- C++ -*- internal _Meta class.
2
3 // Copyright (C) 1997-2015 Free Software Foundation, Inc.
4 //
5 // This file is part of the GNU ISO C++ Library.  This library is free
6 // software; you can redistribute it and/or modify it under the
7 // terms of the GNU General Public License as published by the
8 // Free Software Foundation; either version 3, or (at your option)
9 // any later version.
10
11 // This library is distributed in the hope that it will be useful,
12 // but WITHOUT ANY WARRANTY; without even the implied warranty of
13 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 // GNU General Public License for more details.
15
16 // Under Section 7 of GPL version 3, you are granted additional
17 // permissions described in the GCC Runtime Library Exception, version
18 // 3.1, as published by the Free Software Foundation.
19
20 // You should have received a copy of the GNU General Public License and
21 // a copy of the GCC Runtime Library Exception along with this program;
22 // see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
23 // <http://www.gnu.org/licenses/>.
24
25 /** @file bits/valarray_after.h
26  *  This is an internal header file, included by other library headers.
27  *  Do not attempt to use it directly. @headername{valarray}
28  */
29
30 // Written by Gabriel Dos Reis <Gabriel.Dos-Reis@cmla.ens-cachan.fr>
31
32 #ifndef _VALARRAY_AFTER_H
33 #define _VALARRAY_AFTER_H 1
34
35 #pragma GCC system_header
36
37 namespace std _GLIBCXX_VISIBILITY(default)
38 {
39 _GLIBCXX_BEGIN_NAMESPACE_VERSION
40
41   //
42   // gslice_array closure.
43   //
44   template<class _Dom>
45     class _GBase
46     {
47     public:
48       typedef typename _Dom::value_type value_type;
49       
50       _GBase (const _Dom& __e, const valarray<size_t>& __i)
51       : _M_expr (__e), _M_index(__i) {}
52       
53       value_type
54       operator[] (size_t __i) const
55       { return _M_expr[_M_index[__i]]; }
56       
57       size_t
58       size () const
59       { return _M_index.size(); }
60
61     private:
62       const _Dom&             _M_expr;
63       const valarray<size_t>& _M_index;
64     };
65
66   template<typename _Tp>
67     class _GBase<_Array<_Tp> >
68     {
69     public:
70       typedef _Tp value_type;
71       
72       _GBase (_Array<_Tp> __a, const valarray<size_t>& __i)
73       : _M_array (__a), _M_index(__i) {}
74       
75       value_type
76       operator[] (size_t __i) const
77       { return _M_array._M_data[_M_index[__i]]; }
78       
79       size_t
80       size () const
81       { return _M_index.size(); }
82
83     private:
84       const _Array<_Tp>       _M_array;
85       const valarray<size_t>& _M_index;
86     };
87
88   template<class _Dom>
89     struct _GClos<_Expr, _Dom>
90     : _GBase<_Dom>
91     {
92       typedef _GBase<_Dom> _Base;
93       typedef typename _Base::value_type value_type;
94       
95       _GClos (const _Dom& __e, const valarray<size_t>& __i)
96       : _Base (__e, __i) {}
97     };
98
99   template<typename _Tp>
100     struct _GClos<_ValArray, _Tp>
101     : _GBase<_Array<_Tp> >
102     {
103       typedef _GBase<_Array<_Tp> > _Base;
104       typedef typename _Base::value_type value_type;
105       
106       _GClos (_Array<_Tp> __a, const valarray<size_t>& __i)
107       : _Base (__a, __i) {}
108     };
109
110   //
111   // indirect_array closure
112   //
113   template<class _Dom>
114     class _IBase
115     {
116     public:
117       typedef typename _Dom::value_type value_type;
118
119       _IBase (const _Dom& __e, const valarray<size_t>& __i)
120       : _M_expr (__e), _M_index (__i) {}
121       
122       value_type
123       operator[] (size_t __i) const
124       { return _M_expr[_M_index[__i]]; }
125       
126       size_t
127       size() const
128       { return _M_index.size(); }
129
130     private:
131       const _Dom&             _M_expr;
132       const valarray<size_t>& _M_index;
133     };
134
135   template<class _Dom>
136     struct _IClos<_Expr, _Dom>
137     : _IBase<_Dom>
138     {
139       typedef _IBase<_Dom> _Base;
140       typedef typename _Base::value_type value_type;
141       
142       _IClos (const _Dom& __e, const valarray<size_t>& __i)
143       : _Base (__e, __i) {}
144     };
145
146   template<typename _Tp>
147     struct _IClos<_ValArray, _Tp>
148     : _IBase<valarray<_Tp> >
149     {
150       typedef _IBase<valarray<_Tp> > _Base;
151       typedef _Tp value_type;
152       
153       _IClos (const valarray<_Tp>& __a, const valarray<size_t>& __i)
154       : _Base (__a, __i) {}
155     };
156   
157   //
158   // class _Expr
159   //
160   template<class _Clos, typename _Tp>
161     class _Expr
162     {
163     public:
164       typedef _Tp value_type;
165
166       _Expr(const _Clos&);
167
168       const _Clos& operator()() const;
169
170       value_type operator[](size_t) const;
171       valarray<value_type> operator[](slice) const;
172       valarray<value_type> operator[](const gslice&) const;
173       valarray<value_type> operator[](const valarray<bool>&) const;
174       valarray<value_type> operator[](const valarray<size_t>&) const;
175
176       _Expr<_UnClos<__unary_plus, std::_Expr, _Clos>, value_type>
177       operator+() const;
178
179       _Expr<_UnClos<__negate, std::_Expr, _Clos>, value_type>
180       operator-() const;
181
182       _Expr<_UnClos<__bitwise_not, std::_Expr, _Clos>, value_type>
183       operator~() const;
184
185       _Expr<_UnClos<__logical_not, std::_Expr, _Clos>, bool>
186       operator!() const;
187
188       size_t size() const;
189       value_type sum() const;
190
191       valarray<value_type> shift(int) const;
192       valarray<value_type> cshift(int) const;
193
194       value_type min() const;
195       value_type max() const;
196
197       valarray<value_type> apply(value_type (*)(const value_type&)) const;
198       valarray<value_type> apply(value_type (*)(value_type)) const;
199
200     private:
201       const _Clos _M_closure;
202     };
203
204   template<class _Clos, typename _Tp>
205     inline
206     _Expr<_Clos, _Tp>::_Expr(const _Clos& __c) : _M_closure(__c) {}
207
208   template<class _Clos, typename _Tp>
209     inline const _Clos&
210     _Expr<_Clos, _Tp>::operator()() const
211     { return _M_closure; }
212
213   template<class _Clos, typename _Tp>
214     inline _Tp
215     _Expr<_Clos, _Tp>::operator[](size_t __i) const
216     { return _M_closure[__i]; }
217
218   template<class _Clos, typename _Tp>
219     inline valarray<_Tp>
220     _Expr<_Clos, _Tp>::operator[](slice __s) const
221     {
222       valarray<_Tp> __v = valarray<_Tp>(*this)[__s];
223       return __v;
224     }
225
226   template<class _Clos, typename _Tp>
227     inline valarray<_Tp>
228     _Expr<_Clos, _Tp>::operator[](const gslice& __gs) const
229     {
230       valarray<_Tp> __v = valarray<_Tp>(*this)[__gs];
231       return __v;
232     }
233
234   template<class _Clos, typename _Tp>
235     inline valarray<_Tp>
236     _Expr<_Clos, _Tp>::operator[](const valarray<bool>& __m) const
237     {
238       valarray<_Tp> __v = valarray<_Tp>(*this)[__m];
239       return __v;
240     }
241
242   template<class _Clos, typename _Tp>
243     inline valarray<_Tp>
244     _Expr<_Clos, _Tp>::operator[](const valarray<size_t>& __i) const
245     {
246       valarray<_Tp> __v = valarray<_Tp>(*this)[__i];
247       return __v;
248     }
249
250   template<class _Clos, typename _Tp>
251     inline size_t
252     _Expr<_Clos, _Tp>::size() const
253     { return _M_closure.size(); }
254
255   template<class _Clos, typename _Tp>
256     inline valarray<_Tp>
257     _Expr<_Clos, _Tp>::shift(int __n) const
258     {
259       valarray<_Tp> __v = valarray<_Tp>(*this).shift(__n);
260       return __v;
261     }
262
263   template<class _Clos, typename _Tp>
264     inline valarray<_Tp>
265     _Expr<_Clos, _Tp>::cshift(int __n) const
266     {
267       valarray<_Tp> __v = valarray<_Tp>(*this).cshift(__n);
268       return __v;
269     }
270
271   template<class _Clos, typename _Tp>
272     inline valarray<_Tp>
273     _Expr<_Clos, _Tp>::apply(_Tp __f(const _Tp&)) const
274     {
275       valarray<_Tp> __v = valarray<_Tp>(*this).apply(__f);
276       return __v;
277     }
278
279   template<class _Clos, typename _Tp>
280     inline valarray<_Tp>
281     _Expr<_Clos, _Tp>::apply(_Tp __f(_Tp)) const
282     {
283       valarray<_Tp> __v = valarray<_Tp>(*this).apply(__f);
284       return __v;
285     }
286
287   // XXX: replace this with a more robust summation algorithm.
288   template<class _Clos, typename _Tp>
289     inline _Tp
290     _Expr<_Clos, _Tp>::sum() const
291     {
292       size_t __n = _M_closure.size();
293       if (__n == 0)
294         return _Tp();
295       else
296         {
297           _Tp __s = _M_closure[--__n];
298           while (__n != 0)
299             __s += _M_closure[--__n];
300           return __s;
301         }
302     }
303
304   template<class _Clos, typename _Tp>
305     inline _Tp
306     _Expr<_Clos, _Tp>::min() const
307     { return __valarray_min(_M_closure); }
308
309   template<class _Clos, typename _Tp>
310     inline _Tp
311     _Expr<_Clos, _Tp>::max() const
312     { return __valarray_max(_M_closure); }
313
314   template<class _Dom, typename _Tp>
315     inline _Expr<_UnClos<__logical_not, _Expr, _Dom>, bool>
316     _Expr<_Dom, _Tp>::operator!() const
317     {
318       typedef _UnClos<__logical_not, std::_Expr, _Dom> _Closure;
319       return _Expr<_Closure, bool>(_Closure(this->_M_closure));
320     }
321
322 #define _DEFINE_EXPR_UNARY_OPERATOR(_Op, _Name)                           \
323   template<class _Dom, typename _Tp>                                      \
324     inline _Expr<_UnClos<_Name, std::_Expr, _Dom>, _Tp>                   \
325     _Expr<_Dom, _Tp>::operator _Op() const                                \
326     {                                                                     \
327       typedef _UnClos<_Name, std::_Expr, _Dom> _Closure;                  \
328       return _Expr<_Closure, _Tp>(_Closure(this->_M_closure));            \
329     }
330
331     _DEFINE_EXPR_UNARY_OPERATOR(+, __unary_plus)
332     _DEFINE_EXPR_UNARY_OPERATOR(-, __negate)
333     _DEFINE_EXPR_UNARY_OPERATOR(~, __bitwise_not)
334
335 #undef _DEFINE_EXPR_UNARY_OPERATOR
336
337 #define _DEFINE_EXPR_BINARY_OPERATOR(_Op, _Name)                        \
338   template<class _Dom1, class _Dom2>                                    \
339     inline _Expr<_BinClos<_Name, _Expr, _Expr, _Dom1, _Dom2>,           \
340            typename __fun<_Name, typename _Dom1::value_type>::result_type> \
341     operator _Op(const _Expr<_Dom1, typename _Dom1::value_type>& __v,   \
342                  const _Expr<_Dom2, typename _Dom2::value_type>& __w)   \
343     {                                                                   \
344       typedef typename _Dom1::value_type _Arg;                          \
345       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
346       typedef _BinClos<_Name, _Expr, _Expr, _Dom1, _Dom2> _Closure;     \
347       return _Expr<_Closure, _Value>(_Closure(__v(), __w()));           \
348     }                                                                   \
349                                                                         \
350   template<class _Dom>                                                  \
351     inline _Expr<_BinClos<_Name, _Expr, _Constant, _Dom,                \
352                           typename _Dom::value_type>,                   \
353              typename __fun<_Name, typename _Dom::value_type>::result_type> \
354     operator _Op(const _Expr<_Dom, typename _Dom::value_type>& __v,     \
355                  const typename _Dom::value_type& __t)                  \
356     {                                                                   \
357       typedef typename _Dom::value_type _Arg;                           \
358       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
359       typedef _BinClos<_Name, _Expr, _Constant, _Dom, _Arg> _Closure;   \
360       return _Expr<_Closure, _Value>(_Closure(__v(), __t));             \
361     }                                                                   \
362                                                                         \
363   template<class _Dom>                                                  \
364     inline _Expr<_BinClos<_Name, _Constant, _Expr,                      \
365                           typename _Dom::value_type, _Dom>,             \
366              typename __fun<_Name, typename _Dom::value_type>::result_type> \
367     operator _Op(const typename _Dom::value_type& __t,                  \
368                  const _Expr<_Dom, typename _Dom::value_type>& __v)     \
369     {                                                                   \
370       typedef typename _Dom::value_type _Arg;                           \
371       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
372       typedef _BinClos<_Name, _Constant, _Expr, _Arg, _Dom> _Closure;   \
373       return _Expr<_Closure, _Value>(_Closure(__t, __v()));             \
374     }                                                                   \
375                                                                         \
376   template<class _Dom>                                                  \
377     inline _Expr<_BinClos<_Name, _Expr, _ValArray,                      \
378                           _Dom, typename _Dom::value_type>,             \
379              typename __fun<_Name, typename _Dom::value_type>::result_type> \
380     operator _Op(const _Expr<_Dom,typename _Dom::value_type>& __e,      \
381                  const valarray<typename _Dom::value_type>& __v)        \
382     {                                                                   \
383       typedef typename _Dom::value_type _Arg;                           \
384       typedef typename __fun<_Name, _Arg>::result_type _Value;          \
385       typedef _BinClos<_Name, _Expr, _ValArray, _Dom, _Arg> _Closure;   \
386       return _Expr<_Closure, _Value>(_Closure(__e(), __v));             \
387     }                                                                   \
388                                                                         \
389   template<class _Dom>                                                  \
390     inline _Expr<_BinClos<_Name, _ValArray, _Expr,                      \
391                  typename _Dom::value_type, _Dom>,                      \
392              typename __fun<_Name, typename _Dom::value_type>::result_type> \
393     operator _Op(const valarray<typename _Dom::value_type>& __v,        \
394                  const _Expr<_Dom, typename _Dom::value_type>& __e)     \
395     {                                                                   \
396       typedef typename _Dom::value_type _Tp;                            \
397       typedef typename __fun<_Name, _Tp>::result_type _Value;           \
398       typedef _BinClos<_Name, _ValArray, _Expr, _Tp, _Dom> _Closure;    \
399       return _Expr<_Closure, _Value>(_Closure(__v, __e ()));            \
400     }
401
402     _DEFINE_EXPR_BINARY_OPERATOR(+, __plus)
403     _DEFINE_EXPR_BINARY_OPERATOR(-, __minus)
404     _DEFINE_EXPR_BINARY_OPERATOR(*, __multiplies)
405     _DEFINE_EXPR_BINARY_OPERATOR(/, __divides)
406     _DEFINE_EXPR_BINARY_OPERATOR(%, __modulus)
407     _DEFINE_EXPR_BINARY_OPERATOR(^, __bitwise_xor)
408     _DEFINE_EXPR_BINARY_OPERATOR(&, __bitwise_and)
409     _DEFINE_EXPR_BINARY_OPERATOR(|, __bitwise_or)
410     _DEFINE_EXPR_BINARY_OPERATOR(<<, __shift_left)
411     _DEFINE_EXPR_BINARY_OPERATOR(>>, __shift_right)
412     _DEFINE_EXPR_BINARY_OPERATOR(&&, __logical_and)
413     _DEFINE_EXPR_BINARY_OPERATOR(||, __logical_or)
414     _DEFINE_EXPR_BINARY_OPERATOR(==, __equal_to)
415     _DEFINE_EXPR_BINARY_OPERATOR(!=, __not_equal_to)
416     _DEFINE_EXPR_BINARY_OPERATOR(<, __less)
417     _DEFINE_EXPR_BINARY_OPERATOR(>, __greater)
418     _DEFINE_EXPR_BINARY_OPERATOR(<=, __less_equal)
419     _DEFINE_EXPR_BINARY_OPERATOR(>=, __greater_equal)
420
421 #undef _DEFINE_EXPR_BINARY_OPERATOR
422
423 #define _DEFINE_EXPR_UNARY_FUNCTION(_Name, _UName)                       \
424   template<class _Dom>                                                   \
425     inline _Expr<_UnClos<_UName, _Expr, _Dom>,                           \
426                  typename _Dom::value_type>                              \
427     _Name(const _Expr<_Dom, typename _Dom::value_type>& __e)             \
428     {                                                                    \
429       typedef typename _Dom::value_type _Tp;                             \
430       typedef _UnClos<_UName, _Expr, _Dom> _Closure;                     \
431       return _Expr<_Closure, _Tp>(_Closure(__e()));                      \
432     }                                                                    \
433                                                                          \
434   template<typename _Tp>                                                 \
435     inline _Expr<_UnClos<_UName, _ValArray, _Tp>, _Tp>                   \
436     _Name(const valarray<_Tp>& __v)                                      \
437     {                                                                    \
438       typedef _UnClos<_UName, _ValArray, _Tp> _Closure;                  \
439       return _Expr<_Closure, _Tp>(_Closure(__v));                        \
440     }
441
442     _DEFINE_EXPR_UNARY_FUNCTION(abs, _Abs)
443     _DEFINE_EXPR_UNARY_FUNCTION(cos, _Cos)
444     _DEFINE_EXPR_UNARY_FUNCTION(acos, _Acos)
445     _DEFINE_EXPR_UNARY_FUNCTION(cosh, _Cosh)
446     _DEFINE_EXPR_UNARY_FUNCTION(sin, _Sin)
447     _DEFINE_EXPR_UNARY_FUNCTION(asin, _Asin)
448     _DEFINE_EXPR_UNARY_FUNCTION(sinh, _Sinh)
449     _DEFINE_EXPR_UNARY_FUNCTION(tan, _Tan)
450     _DEFINE_EXPR_UNARY_FUNCTION(tanh, _Tanh)
451     _DEFINE_EXPR_UNARY_FUNCTION(atan, _Atan)
452     _DEFINE_EXPR_UNARY_FUNCTION(exp, _Exp)
453     _DEFINE_EXPR_UNARY_FUNCTION(log, _Log)
454     _DEFINE_EXPR_UNARY_FUNCTION(log10, _Log10)
455     _DEFINE_EXPR_UNARY_FUNCTION(sqrt, _Sqrt)
456
457 #undef _DEFINE_EXPR_UNARY_FUNCTION
458
459 #define _DEFINE_EXPR_BINARY_FUNCTION(_Fun, _UFun)                      \
460   template<class _Dom1, class _Dom2>                                   \
461     inline _Expr<_BinClos<_UFun, _Expr, _Expr, _Dom1, _Dom2>,          \
462                  typename _Dom1::value_type>                           \
463     _Fun(const _Expr<_Dom1, typename _Dom1::value_type>& __e1,         \
464          const _Expr<_Dom2, typename _Dom2::value_type>& __e2)         \
465     {                                                                  \
466       typedef typename _Dom1::value_type _Tp;                          \
467       typedef _BinClos<_UFun, _Expr, _Expr, _Dom1, _Dom2> _Closure;    \
468       return _Expr<_Closure, _Tp>(_Closure(__e1(), __e2()));           \
469     }                                                                  \
470                                                                        \
471   template<class _Dom>                                                 \
472     inline _Expr<_BinClos<_UFun, _Expr, _ValArray, _Dom,               \
473                           typename _Dom::value_type>,                  \
474                  typename _Dom::value_type>                            \
475     _Fun(const _Expr<_Dom, typename _Dom::value_type>& __e,            \
476          const valarray<typename _Dom::value_type>& __v)               \
477     {                                                                  \
478       typedef typename _Dom::value_type _Tp;                           \
479       typedef _BinClos<_UFun, _Expr, _ValArray, _Dom, _Tp> _Closure;   \
480       return _Expr<_Closure, _Tp>(_Closure(__e(), __v));               \
481     }                                                                  \
482                                                                        \
483   template<class _Dom>                                                 \
484     inline _Expr<_BinClos<_UFun, _ValArray, _Expr,                     \
485                           typename _Dom::value_type, _Dom>,            \
486                  typename _Dom::value_type>                            \
487     _Fun(const valarray<typename _Dom::valarray>& __v,                 \
488          const _Expr<_Dom, typename _Dom::value_type>& __e)            \
489     {                                                                  \
490       typedef typename _Dom::value_type _Tp;                           \
491       typedef _BinClos<_UFun, _ValArray, _Expr, _Tp, _Dom> _Closure;   \
492       return _Expr<_Closure, _Tp>(_Closure(__v, __e()));               \
493     }                                                                  \
494                                                                        \
495   template<class _Dom>                                                 \
496     inline _Expr<_BinClos<_UFun, _Expr, _Constant, _Dom,               \
497                           typename _Dom::value_type>,                  \
498                  typename _Dom::value_type>                            \
499     _Fun(const _Expr<_Dom, typename _Dom::value_type>& __e,            \
500          const typename _Dom::value_type& __t)                         \
501     {                                                                  \
502       typedef typename _Dom::value_type _Tp;                           \
503       typedef _BinClos<_UFun, _Expr, _Constant, _Dom, _Tp> _Closure;   \
504       return _Expr<_Closure, _Tp>(_Closure(__e(), __t));               \
505     }                                                                  \
506                                                                        \
507   template<class _Dom>                                                 \
508     inline _Expr<_BinClos<_UFun, _Constant, _Expr,                     \
509                           typename _Dom::value_type, _Dom>,            \
510                  typename _Dom::value_type>                            \
511     _Fun(const typename _Dom::value_type& __t,                         \
512          const _Expr<_Dom, typename _Dom::value_type>& __e)            \
513     {                                                                  \
514       typedef typename _Dom::value_type _Tp;                           \
515       typedef _BinClos<_UFun, _Constant, _Expr, _Tp, _Dom> _Closure;   \
516       return _Expr<_Closure, _Tp>(_Closure(__t, __e()));               \
517     }                                                                  \
518                                                                        \
519   template<typename _Tp>                                               \
520     inline _Expr<_BinClos<_UFun, _ValArray, _ValArray, _Tp, _Tp>, _Tp> \
521     _Fun(const valarray<_Tp>& __v, const valarray<_Tp>& __w)           \
522     {                                                                  \
523       typedef _BinClos<_UFun, _ValArray, _ValArray, _Tp, _Tp> _Closure;\
524       return _Expr<_Closure, _Tp>(_Closure(__v, __w));                 \
525     }                                                                  \
526                                                                        \
527   template<typename _Tp>                                               \
528     inline _Expr<_BinClos<_UFun, _ValArray, _Constant, _Tp, _Tp>, _Tp> \
529     _Fun(const valarray<_Tp>& __v, const _Tp& __t)                     \
530     {                                                                  \
531       typedef _BinClos<_UFun, _ValArray, _Constant, _Tp, _Tp> _Closure;\
532       return _Expr<_Closure, _Tp>(_Closure(__v, __t));                 \
533     }                                                                  \
534                                                                        \
535   template<typename _Tp>                                               \
536     inline _Expr<_BinClos<_UFun, _Constant, _ValArray, _Tp, _Tp>, _Tp> \
537     _Fun(const _Tp& __t, const valarray<_Tp>& __v)                     \
538     {                                                                  \
539       typedef _BinClos<_UFun, _Constant, _ValArray, _Tp, _Tp> _Closure;\
540       return _Expr<_Closure, _Tp>(_Closure(__t, __v));                 \
541     }
542
543 _DEFINE_EXPR_BINARY_FUNCTION(atan2, _Atan2)
544 _DEFINE_EXPR_BINARY_FUNCTION(pow, _Pow)
545
546 #undef _DEFINE_EXPR_BINARY_FUNCTION
547
548 _GLIBCXX_END_NAMESPACE_VERSION
549 } // namespace
550
551 #endif /* _CPP_VALARRAY_AFTER_H */