xref: /freebsd/contrib/llvm-project/libcxx/src/valarray.cpp (revision 0b57cec536236d46e3dba9bd041533462f33dbb7)
1*0b57cec5SDimitry Andric //===------------------------ valarray.cpp --------------------------------===//
2*0b57cec5SDimitry Andric //
3*0b57cec5SDimitry Andric // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4*0b57cec5SDimitry Andric // See https://llvm.org/LICENSE.txt for license information.
5*0b57cec5SDimitry Andric // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6*0b57cec5SDimitry Andric //
7*0b57cec5SDimitry Andric //===----------------------------------------------------------------------===//
8*0b57cec5SDimitry Andric 
9*0b57cec5SDimitry Andric #include "valarray"
10*0b57cec5SDimitry Andric 
11*0b57cec5SDimitry Andric _LIBCPP_BEGIN_NAMESPACE_STD
12*0b57cec5SDimitry Andric 
13*0b57cec5SDimitry Andric // These two symbols are part of the v1 ABI but not part of the >=v2 ABI.
14*0b57cec5SDimitry Andric #if _LIBCPP_ABI_VERSION == 1
15*0b57cec5SDimitry Andric template valarray<size_t>::valarray(size_t);
16*0b57cec5SDimitry Andric template valarray<size_t>::~valarray();
17*0b57cec5SDimitry Andric #endif
18*0b57cec5SDimitry Andric 
19*0b57cec5SDimitry Andric template void valarray<size_t>::resize(size_t, size_t);
20*0b57cec5SDimitry Andric 
21*0b57cec5SDimitry Andric void
22*0b57cec5SDimitry Andric gslice::__init(size_t __start)
23*0b57cec5SDimitry Andric {
24*0b57cec5SDimitry Andric     valarray<size_t> __indices(__size_.size());
25*0b57cec5SDimitry Andric     size_t __k = __size_.size() != 0;
26*0b57cec5SDimitry Andric     for (size_t __i = 0; __i < __size_.size(); ++__i)
27*0b57cec5SDimitry Andric         __k *= __size_[__i];
28*0b57cec5SDimitry Andric     __1d_.resize(__k);
29*0b57cec5SDimitry Andric     if (__1d_.size())
30*0b57cec5SDimitry Andric     {
31*0b57cec5SDimitry Andric         __k = 0;
32*0b57cec5SDimitry Andric         __1d_[__k] = __start;
33*0b57cec5SDimitry Andric         while (true)
34*0b57cec5SDimitry Andric         {
35*0b57cec5SDimitry Andric             size_t __i = __indices.size() - 1;
36*0b57cec5SDimitry Andric             while (true)
37*0b57cec5SDimitry Andric             {
38*0b57cec5SDimitry Andric                 if (++__indices[__i] < __size_[__i])
39*0b57cec5SDimitry Andric                 {
40*0b57cec5SDimitry Andric                     ++__k;
41*0b57cec5SDimitry Andric                     __1d_[__k] = __1d_[__k-1] + __stride_[__i];
42*0b57cec5SDimitry Andric                     for (size_t __j = __i + 1; __j != __indices.size(); ++__j)
43*0b57cec5SDimitry Andric                         __1d_[__k] -= __stride_[__j] * (__size_[__j] - 1);
44*0b57cec5SDimitry Andric                     break;
45*0b57cec5SDimitry Andric                 }
46*0b57cec5SDimitry Andric                 else
47*0b57cec5SDimitry Andric                 {
48*0b57cec5SDimitry Andric                     if (__i == 0)
49*0b57cec5SDimitry Andric                         return;
50*0b57cec5SDimitry Andric                     __indices[__i--] = 0;
51*0b57cec5SDimitry Andric                 }
52*0b57cec5SDimitry Andric             }
53*0b57cec5SDimitry Andric         }
54*0b57cec5SDimitry Andric     }
55*0b57cec5SDimitry Andric }
56*0b57cec5SDimitry Andric 
57*0b57cec5SDimitry Andric _LIBCPP_END_NAMESPACE_STD
58