Morph MORPH
SIGN IN SIGN UP
ziglang / zig UNCLAIMED

Moved to Codeberg

0 0 0 Zig
Normal View History Raw
std.math: Add documentation for all functions and algorithm sources
2019-05-01 18:15:57 +12:00
// Ported from musl, which is licensed under the MIT license:
// https://git.musl-libc.org/cgit/musl/tree/COPYRIGHT
//
// https://git.musl-libc.org/cgit/musl/tree/src/complex/__cexpf.c
// https://git.musl-libc.org/cgit/musl/tree/src/complex/__cexp.c
rename std lib files to new convention
2019-03-02 16:46:04 -05:00
const std = @import("../../std.zig");
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const debug = std.debug;
const math = std.math;
std: Fix complex ldexp implementation Two bugs in the implementation ported from musl made all the complex functions relying on ldexp return incorrect results in some cases. Spotted in #9047
2021-06-12 15:35:52 +02:00
const testing = std.testing;
std, src, doc, test: remove unused variables
2021-06-09 21:35:42 -04:00
const cmath = math.complex;
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const Complex = cmath.Complex;
std.math: Add documentation for all functions and algorithm sources
2019-05-01 18:15:57 +12:00
/// Returns exp(z) scaled to avoid overflow.
use Peer Type Resolution for standalone complex fn use peer type resolution Update complex.zig Revert "use peer type resolution" This reverts commit 1bc681ca5b36d2b55b5efab5a5dbec7e8a17332e. Revert "Update pow.zig" This reverts commit 5487e8d3159f832b5a0bf29a06bd12575182464f. Update pow.zig Revert "Update pow.zig" This reverts commit 521153d1ef004d627c785f2d3fe5e6497dc15073. Update pow.zig
2023-11-08 04:10:43 +08:00
pub fn ldexp_cexp(z: anytype, expt: i32) Complex(@TypeOf(z.re, z.im)) {
const T = @TypeOf(z.re, z.im);
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
return switch (T) {
f32 => ldexp_cexp32(z, expt),
f64 => ldexp_cexp64(z, expt),
else => unreachable,
};
}
use * for pointer type instead of & See #770 To help automatically translate code, see the zig-fmt-pointer-reform-2 branch. This will convert all & into *. Due to the syntax ambiguity (which is why we are making this change), even address-of & will turn into *, so you'll have to manually fix thes instances. You will be guaranteed to get compile errors for them - expected 'type', found 'foo'
2018-05-31 10:56:59 -04:00
fn frexp_exp32(x: f32, expt: *i32) f32 {
more std lib to postfix deref with zig fmt
2018-05-17 00:56:14 -04:00
const k = 235; // reduction constant
const kln2 = 162.88958740; // k * ln2
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
compiler-rt: math functions reorg * unify the logic for exporting math functions from compiler-rt, with the appropriate suffixes and prefixes. - add all missing f128 and f80 exports. Functions with missing implementations call other functions and have TODO comments. - also add f16 functions * move math functions from freestanding libc to compiler-rt (#7265) * enable all the f128 and f80 code in the stage2 compiler and behavior tests (#11161). * update std lib to use builtins rather than `std.math`.
2022-04-26 10:13:55 -07:00
const exp_x = @exp(x - kln2);
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const hx = @as(u32, @bitCast(exp_x));
remove integer and float casting syntax * add `@intCast` * add `@floatCast` * add `@floatToInt` * add `@intToFloat` See #1061
2018-06-17 02:57:07 -04:00
// TODO zig should allow this cast implicitly because it should know the value is in range
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
expt.* = @as(i32, @intCast(hx >> 23)) - (0x7f + 127) + k;
return @as(f32, @bitCast((hx & 0x7fffff) | ((0x7f + 127) << 23)));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
}
remove integer and float casting syntax * add `@intCast` * add `@floatCast` * add `@floatToInt` * add `@intToFloat` See #1061
2018-06-17 02:57:07 -04:00
fn ldexp_cexp32(z: Complex(f32), expt: i32) Complex(f32) {
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
var ex_expt: i32 = undefined;
const exp_x = frexp_exp32(z.re, &ex_expt);
const exptf = expt + ex_expt;
const half_expt1 = @divTrunc(exptf, 2);
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const scale1 = @as(f32, @bitCast((0x7f + half_expt1) << 23));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const half_expt2 = exptf - half_expt1;
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const scale2 = @as(f32, @bitCast((0x7f + half_expt2) << 23));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
std: Fix complex ldexp implementation Two bugs in the implementation ported from musl made all the complex functions relying on ldexp return incorrect results in some cases. Spotted in #9047
2021-06-12 15:35:52 +02:00
return Complex(f32).init(
std: use float builtins instead of std.math
2022-04-27 19:35:28 -07:00
@cos(z.im) * exp_x * scale1 * scale2,
@sin(z.im) * exp_x * scale1 * scale2,
std: Fix complex ldexp implementation Two bugs in the implementation ported from musl made all the complex functions relying on ldexp return incorrect results in some cases. Spotted in #9047
2021-06-12 15:35:52 +02:00
);
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
}
use * for pointer type instead of & See #770 To help automatically translate code, see the zig-fmt-pointer-reform-2 branch. This will convert all & into *. Due to the syntax ambiguity (which is why we are making this change), even address-of & will turn into *, so you'll have to manually fix thes instances. You will be guaranteed to get compile errors for them - expected 'type', found 'foo'
2018-05-31 10:56:59 -04:00
fn frexp_exp64(x: f64, expt: *i32) f64 {
more std lib to postfix deref with zig fmt
2018-05-17 00:56:14 -04:00
const k = 1799; // reduction constant
const kln2 = 1246.97177782734161156; // k * ln2
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
compiler-rt: math functions reorg * unify the logic for exporting math functions from compiler-rt, with the appropriate suffixes and prefixes. - add all missing f128 and f80 exports. Functions with missing implementations call other functions and have TODO comments. - also add f16 functions * move math functions from freestanding libc to compiler-rt (#7265) * enable all the f128 and f80 code in the stage2 compiler and behavior tests (#11161). * update std lib to use builtins rather than `std.math`.
2022-04-26 10:13:55 -07:00
const exp_x = @exp(x - kln2);
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const fx = @as(u64, @bitCast(exp_x));
const hx = @as(u32, @intCast(fx >> 32));
const lx = @as(u32, @truncate(fx));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
expt.* = @as(i32, @intCast(hx >> 20)) - (0x3ff + 1023) + k;
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const high_word = (hx & 0xfffff) | ((0x3ff + 1023) << 20);
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
return @as(f64, @bitCast((@as(u64, high_word) << 32) | lx));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
}
remove integer and float casting syntax * add `@intCast` * add `@floatCast` * add `@floatToInt` * add `@intToFloat` See #1061
2018-06-17 02:57:07 -04:00
fn ldexp_cexp64(z: Complex(f64), expt: i32) Complex(f64) {
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
var ex_expt: i32 = undefined;
const exp_x = frexp_exp64(z.re, &ex_expt);
update more of the std lib to use `@as`
2019-11-07 18:52:09 -05:00
const exptf = @as(i64, expt + ex_expt);
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const half_expt1 = @divTrunc(exptf, 2);
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const scale1 = @as(f64, @bitCast((0x3ff + half_expt1) << (20 + 32)));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
const half_expt2 = exptf - half_expt1;
all: migrate code to new cast builtin syntax Most of this migration was performed automatically with `zig fmt`. There were a few exceptions which I had to manually fix: * `@alignCast` and `@addrSpaceCast` cannot be automatically rewritten * `@truncate`'s fixup is incorrect for vectors * Test cases are not formatted, and their error locations change
2023-06-22 18:46:56 +01:00
const scale2 = @as(f64, @bitCast((0x3ff + half_expt2) << (20 + 32)));
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
std.math.Complex: Change `new()` to `init()`
2021-05-18 02:57:51 +07:00
return Complex(f64).init(
std: use float builtins instead of std.math
2022-04-27 19:35:28 -07:00
@cos(z.im) * exp_x * scale1 * scale2,
@sin(z.im) * exp_x * scale1 * scale2,
Add initial complex-number support - Library type instead of builtin - All C complex functions implemented Partial WIP: Needs more tests for edge cases.
2018-04-24 19:18:31 +12:00
);
}
Reference in New Issue Copy Permalink