zig/test/behavior/floatop.zig
Andrew Kelley 4307436b99 move behavior tests from test/stage1/ to test/
And fix test cases to make them pass. This is in preparation for
starting to pass behavior tests with self-hosted.
2021-04-29 15:54:04 -07:00

466 lines
13 KiB
Zig

const std = @import("std");
const expect = std.testing.expect;
const math = std.math;
const pi = std.math.pi;
const e = std.math.e;
const Vector = std.meta.Vector;
const epsilon = 0.000001;
test "@sqrt" {
comptime testSqrt();
testSqrt();
}
fn testSqrt() void {
{
var a: f16 = 4;
expect(@sqrt(a) == 2);
}
{
var a: f32 = 9;
expect(@sqrt(a) == 3);
var b: f32 = 1.1;
expect(math.approxEqAbs(f32, @sqrt(b), 1.0488088481701516, epsilon));
}
{
var a: f64 = 25;
expect(@sqrt(a) == 5);
}
{
const a: comptime_float = 25.0;
expect(@sqrt(a) == 5.0);
}
// TODO https://github.com/ziglang/zig/issues/4026
//{
// var a: f128 = 49;
// expect(@sqrt(a) == 7);
//}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 3.3, 4.4 };
var result = @sqrt(v);
expect(math.approxEqAbs(f32, @sqrt(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @sqrt(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @sqrt(@as(f32, 3.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @sqrt(@as(f32, 4.4)), result[3], epsilon));
}
}
test "more @sqrt f16 tests" {
// TODO these are not all passing at comptime
expect(@sqrt(@as(f16, 0.0)) == 0.0);
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 2.0)), 1.414214, epsilon));
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 3.6)), 1.897367, epsilon));
expect(@sqrt(@as(f16, 4.0)) == 2.0);
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 7.539840)), 2.745877, epsilon));
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 19.230934)), 4.385309, epsilon));
expect(@sqrt(@as(f16, 64.0)) == 8.0);
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 64.1)), 8.006248, epsilon));
expect(math.approxEqAbs(f16, @sqrt(@as(f16, 8942.230469)), 94.563370, epsilon));
// special cases
expect(math.isPositiveInf(@sqrt(@as(f16, math.inf(f16)))));
expect(@sqrt(@as(f16, 0.0)) == 0.0);
expect(@sqrt(@as(f16, -0.0)) == -0.0);
expect(math.isNan(@sqrt(@as(f16, -1.0))));
expect(math.isNan(@sqrt(@as(f16, math.nan(f16)))));
}
test "@sin" {
comptime testSin();
testSin();
}
fn testSin() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 0;
expect(@sin(a) == 0);
}
{
var a: f32 = 0;
expect(@sin(a) == 0);
}
{
var a: f64 = 0;
expect(@sin(a) == 0);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 3.3, 4.4 };
var result = @sin(v);
expect(math.approxEqAbs(f32, @sin(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @sin(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @sin(@as(f32, 3.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @sin(@as(f32, 4.4)), result[3], epsilon));
}
}
test "@cos" {
comptime testCos();
testCos();
}
fn testCos() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 0;
expect(@cos(a) == 1);
}
{
var a: f32 = 0;
expect(@cos(a) == 1);
}
{
var a: f64 = 0;
expect(@cos(a) == 1);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 3.3, 4.4 };
var result = @cos(v);
expect(math.approxEqAbs(f32, @cos(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @cos(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @cos(@as(f32, 3.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @cos(@as(f32, 4.4)), result[3], epsilon));
}
}
test "@exp" {
comptime testExp();
testExp();
}
fn testExp() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 0;
expect(@exp(a) == 1);
}
{
var a: f32 = 0;
expect(@exp(a) == 1);
}
{
var a: f64 = 0;
expect(@exp(a) == 1);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 0.3, 0.4 };
var result = @exp(v);
expect(math.approxEqAbs(f32, @exp(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @exp(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @exp(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @exp(@as(f32, 0.4)), result[3], epsilon));
}
}
test "@exp2" {
comptime testExp2();
testExp2();
}
fn testExp2() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 2;
expect(@exp2(a) == 4);
}
{
var a: f32 = 2;
expect(@exp2(a) == 4);
}
{
var a: f64 = 2;
expect(@exp2(a) == 4);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 0.3, 0.4 };
var result = @exp2(v);
expect(math.approxEqAbs(f32, @exp2(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @exp2(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @exp2(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @exp2(@as(f32, 0.4)), result[3], epsilon));
}
}
test "@log" {
// Old musl (and glibc?), and our current math.ln implementation do not return 1
// so also accept those values.
comptime testLog();
testLog();
}
fn testLog() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = e;
expect(math.approxEqAbs(f16, @log(a), 1, epsilon));
}
{
var a: f32 = e;
expect(@log(a) == 1 or @log(a) == @bitCast(f32, @as(u32, 0x3f7fffff)));
}
{
var a: f64 = e;
expect(@log(a) == 1 or @log(a) == @bitCast(f64, @as(u64, 0x3ff0000000000000)));
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 0.3, 0.4 };
var result = @log(v);
expect(math.approxEqAbs(f32, @log(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @log(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @log(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @log(@as(f32, 0.4)), result[3], epsilon));
}
}
test "@log2" {
comptime testLog2();
testLog2();
}
fn testLog2() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 4;
expect(@log2(a) == 2);
}
{
var a: f32 = 4;
expect(@log2(a) == 2);
}
{
var a: f64 = 4;
expect(@log2(a) == 2);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 0.3, 0.4 };
var result = @log2(v);
expect(math.approxEqAbs(f32, @log2(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @log2(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @log2(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @log2(@as(f32, 0.4)), result[3], epsilon));
}
}
test "@log10" {
comptime testLog10();
testLog10();
}
fn testLog10() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 100;
expect(@log10(a) == 2);
}
{
var a: f32 = 100;
expect(@log10(a) == 2);
}
{
var a: f64 = 1000;
expect(@log10(a) == 3);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, 2.2, 0.3, 0.4 };
var result = @log10(v);
expect(math.approxEqAbs(f32, @log10(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @log10(@as(f32, 2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @log10(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @log10(@as(f32, 0.4)), result[3], epsilon));
}
}
test "@fabs" {
comptime testFabs();
testFabs();
}
fn testFabs() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = -2.5;
var b: f16 = 2.5;
expect(@fabs(a) == 2.5);
expect(@fabs(b) == 2.5);
}
{
var a: f32 = -2.5;
var b: f32 = 2.5;
expect(@fabs(a) == 2.5);
expect(@fabs(b) == 2.5);
}
{
var a: f64 = -2.5;
var b: f64 = 2.5;
expect(@fabs(a) == 2.5);
expect(@fabs(b) == 2.5);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, -2.2, 0.3, -0.4 };
var result = @fabs(v);
expect(math.approxEqAbs(f32, @fabs(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @fabs(@as(f32, -2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @fabs(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @fabs(@as(f32, -0.4)), result[3], epsilon));
}
}
test "@floor" {
comptime testFloor();
testFloor();
}
fn testFloor() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 2.1;
expect(@floor(a) == 2);
}
{
var a: f32 = 2.1;
expect(@floor(a) == 2);
}
{
var a: f64 = 3.5;
expect(@floor(a) == 3);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, -2.2, 0.3, -0.4 };
var result = @floor(v);
expect(math.approxEqAbs(f32, @floor(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @floor(@as(f32, -2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @floor(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @floor(@as(f32, -0.4)), result[3], epsilon));
}
}
test "@ceil" {
comptime testCeil();
testCeil();
}
fn testCeil() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 2.1;
expect(@ceil(a) == 3);
}
{
var a: f32 = 2.1;
expect(@ceil(a) == 3);
}
{
var a: f64 = 3.5;
expect(@ceil(a) == 4);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, -2.2, 0.3, -0.4 };
var result = @ceil(v);
expect(math.approxEqAbs(f32, @ceil(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @ceil(@as(f32, -2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @ceil(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @ceil(@as(f32, -0.4)), result[3], epsilon));
}
}
test "@trunc" {
comptime testTrunc();
testTrunc();
}
fn testTrunc() void {
// TODO test f128, and c_longdouble
// https://github.com/ziglang/zig/issues/4026
{
var a: f16 = 2.1;
expect(@trunc(a) == 2);
}
{
var a: f32 = 2.1;
expect(@trunc(a) == 2);
}
{
var a: f64 = -3.5;
expect(@trunc(a) == -3);
}
{
var v: Vector(4, f32) = [_]f32{ 1.1, -2.2, 0.3, -0.4 };
var result = @trunc(v);
expect(math.approxEqAbs(f32, @trunc(@as(f32, 1.1)), result[0], epsilon));
expect(math.approxEqAbs(f32, @trunc(@as(f32, -2.2)), result[1], epsilon));
expect(math.approxEqAbs(f32, @trunc(@as(f32, 0.3)), result[2], epsilon));
expect(math.approxEqAbs(f32, @trunc(@as(f32, -0.4)), result[3], epsilon));
}
}
test "floating point comparisons" {
testFloatComparisons();
comptime testFloatComparisons();
}
fn testFloatComparisons() void {
inline for ([_]type{ f16, f32, f64, f128 }) |ty| {
// No decimal part
{
const x: ty = 1.0;
expect(x == 1);
expect(x != 0);
expect(x > 0);
expect(x < 2);
expect(x >= 1);
expect(x <= 1);
}
// Non-zero decimal part
{
const x: ty = 1.5;
expect(x != 1);
expect(x != 2);
expect(x > 1);
expect(x < 2);
expect(x >= 1);
expect(x <= 2);
}
}
}
test "different sized float comparisons" {
testDifferentSizedFloatComparisons();
comptime testDifferentSizedFloatComparisons();
}
fn testDifferentSizedFloatComparisons() void {
var a: f16 = 1;
var b: f64 = 2;
expect(a < b);
}
// TODO This is waiting on library support for the Windows build (not sure why the other's don't need it)
//test "@nearbyint" {
// comptime testNearbyInt();
// testNearbyInt();
//}
//fn testNearbyInt() void {
// // TODO test f16, f128, and c_longdouble
// // https://github.com/ziglang/zig/issues/4026
// {
// var a: f32 = 2.1;
// expect(@nearbyint(a) == 2);
// }
// {
// var a: f64 = -3.75;
// expect(@nearbyint(a) == -4);
// }
//}