Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Investigate the threshold for parallelism with the poles #22

Open
jmert opened this issue Jan 20, 2020 · 0 comments
Open

Investigate the threshold for parallelism with the poles #22

jmert opened this issue Jan 20, 2020 · 0 comments

Comments

@jmert
Copy link
Owner

jmert commented Jan 20, 2020

The parallelism check coded in

CMB.jl/src/sphere.jl

Lines 13 to 24 in e06c74b

rtepsone(::Type{T}) where T = sqrt(eps(one(T)))
"""
∥(u, v)
Test whether vector ``u`` is parallel to vector ``v``. Assumes that both are unit
normalized.
"""
function (u, v)
T = promote_type(eltype(u), eltype(v))
return (one(T) - abs(uv)) < rtepsone(T)
end
was chosen relatively arbitrarily. Whether the limit can be made more strict should be investigated.

@jmert jmert changed the title Investigate the threshold parallelism to the poles Investigate the threshold for parallelism with the poles Jan 20, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant