Skeleton for complex elementary functions #146
Closed
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
- Renamed scaled/unscaled to multiplied/divided - Adopted new Real module name - Adopted expMinusOne, log(onePlus:) names
Rebased on master, then cleaned up exp and expm1. New algorithm for expm1 that gives good componentwise accuracy as well as normwise. I _think_ that this is fundamentally new; I can't recall seeing it anywhere before, though it's simple enough that I expect someone has thought of it before. Also beefed up testing for exp and expm1. Left sketches of log/pow/sqrt/root in place, though I'm still working on filling out the testing for those.
NevinBR
reviewed
Aug 27, 2020
| // exp((log(gfm) + log(9/8), π/4) = exp((log(gfm*9/8), π/4)) | ||
| // = gfm*9/8 * (1/sqrt(2), 1/(sqrt(2)) | ||
| let x = T.log(.greatestFiniteMagnitude) + T.log(9/8) | ||
| let huge = Complex<T>.expMinusOne(Complex(x, .pi/4)) |
There was a problem hiding this comment.
This is in the testExp method, so it should be testing exp, not expMinusOne, right?
There was a problem hiding this comment.
Ah, yeah. Cleaned it up for expMinusOne and copied it back for improvements =)
|
I've created a branch in the main Swift Numerics repo for further work on this, complex-elfns. If you'd like to contribute, please raise PRs against that branch. |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR sets up some basic testing, and adds implementations for the following ElementaryFunctions on Complex:
expexpMinusOnelogpowrootsqrtexpandexpMinusOneare reasonably well tested. Additional testing is needed for the latter four. I also added stubs for all the other functions (which simply fatalError), except for the functions that can be defined by trivial symmetries. I expect to work through the remaining operations fairly quickly, but contributions are welcome and encouraged.If you want to make contributions, I would prefer PRs that add only a single function (or its obvious counterparts, e.g. sin and cos could come together), plus testing, as I expect to subject all of these implementations to fairly careful scrutiny.