forked from mirror/libcxx
You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
55 lines
1.4 KiB
C++
55 lines
1.4 KiB
C++
//===----------------------------------------------------------------------===//
|
|
//
|
|
// The LLVM Compiler Infrastructure
|
|
//
|
|
// This file is dual licensed under the MIT and the University of Illinois Open
|
|
// Source Licenses. See LICENSE.TXT for details.
|
|
//
|
|
//===----------------------------------------------------------------------===//
|
|
|
|
// <valarray>
|
|
|
|
// template<class T> class valarray;
|
|
|
|
// template<class T>
|
|
// valarray<T>
|
|
// atan2(const valarray<T>& x, const valarray<T>& y);
|
|
|
|
#include <valarray>
|
|
#include <cassert>
|
|
#include <sstream>
|
|
#include <cstddef>
|
|
|
|
bool is_about(double x, double y, int p)
|
|
{
|
|
std::ostringstream o;
|
|
o.precision(p);
|
|
scientific(o);
|
|
o << x;
|
|
std::string a = o.str();
|
|
o.str("");
|
|
o << y;
|
|
return a == o.str();
|
|
}
|
|
|
|
int main()
|
|
{
|
|
{
|
|
typedef double T;
|
|
T a1[] = {-.9, -.5, 0., .5, .75};
|
|
T a2[] = {-.8, .25, 0.375, -.5, .75};
|
|
T a3[] = {-2.2974386674766221e+00,
|
|
-1.1071487177940904e+00,
|
|
0.0000000000000000e+00,
|
|
2.3561944901923448e+00,
|
|
7.8539816339744828e-01};
|
|
const unsigned N = sizeof(a1)/sizeof(a1[0]);
|
|
std::valarray<T> v1(a1, N);
|
|
std::valarray<T> v2(a2, N);
|
|
std::valarray<T> v3 = atan2(v1, v2);
|
|
assert(v3.size() == v1.size());
|
|
for (std::size_t i = 0; i < v3.size(); ++i)
|
|
assert(is_about(v3[i], a3[i], 10));
|
|
}
|
|
}
|