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.

Sign up for GitHub

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

Jump to bottom

Cayley/bessel5 #39

Closed
wants to merge 10 commits into from
Closed

Cayley/bessel5 #39

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
23f7e34
cyl_seq
jtlap Nov 4, 2023
b0ff4ca
y0 without allocation
jtlap Nov 5, 2023
07bb1ca
besselyn
jtlap Nov 5, 2023
123e25b
bessel_i
jtlap Nov 5, 2023
600ca66
bessel_k
jtlap Nov 5, 2023
793cc22
bessely
jtlap Nov 5, 2023
5a83221
bessel_h
jtlap Nov 5, 2023
4fa737d
besse_y with_alloca
jtlap Nov 5, 2023
0e8f684
missing include
jtlap Nov 5, 2023
1adb49f
suppressing cyl_seq temporarily
jtlap Nov 6, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
1 change: 1 addition & 0 deletions include/kyosu/details/cayleyify.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23,6 +23,7 @@ auto cayley_extend_rev(auto f, auto z1, auto z2)
return kyosu::real(c) + kyosu::ipart(c)*kyosu::sign(p);
}


auto cayley_extend2(auto f, auto z, auto... rs)
{
auto p = kyosu::pure(z);
Expand Down
13 changes: 6 additions & 7 deletions include/kyosu/types/impl/detail/besseli0.hpp → include/kyosu/details/with_alloca.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,17 +6,16 @@
*/
//======================================================================================================================
#pragma once
#include <kyosu/types/impl/detail/besselj0.hpp>
#include <concepts>

namespace kyosu::_
{

//===-------------------------------------------------------------------------------------------
// cyl_bessel_i0
//===-------------------------------------------------------------------------------------------
template<typename Z>
auto dispatch(eve::tag_of<kyosu::cyl_bessel_i0>, Z z) noexcept
template<typename T, std::invocable<T*> F>
decltype(auto) with_alloca(auto size, F f)
{
return cyl_bessel_j0(complex(-ipart(z), real(z)));
T* p = (T*)(__builtin_alloca_with_align(size*sizeof(T), 8*alignof(T)));
return f(p);
}

}
17 changes: 6 additions & 11 deletions include/kyosu/types/impl/bessel.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,14 +13,9 @@
#include <complex>

#include <kyosu/types/impl/detail/bessel_utils.hpp>
#include <kyosu/types/impl/detail/besselj0.hpp>
#include <kyosu/types/impl/detail/besselj1.hpp>
#include <kyosu/types/impl/detail/besseljn.hpp>
#include <kyosu/types/impl/detail/besseli0.hpp>
#include <kyosu/types/impl/detail/besseli1.hpp>
#include <kyosu/types/impl/detail/besselin.hpp>
#include <kyosu/types/impl/detail/bessely0.hpp>
#include <kyosu/types/impl/detail/bessely1.hpp>
#include <kyosu/types/impl/detail/besselyn.hpp>
#include <kyosu/types/impl/detail/besselhn.hpp>
#include <kyosu/types/impl/detail/besselkn.hpp>

#include <kyosu/types/impl/detail/bessel_j.hpp>
#include <kyosu/types/impl/detail/bessel_i.hpp>
#include <kyosu/types/impl/detail/bessel_y.hpp>
#include <kyosu/types/impl/detail/bessel_h.hpp>
#include <kyosu/types/impl/detail/bessel_k.hpp>
5 changes: 2 additions & 3 deletions include/kyosu/types/impl/detail/besselhn.hpp → include/kyosu/types/impl/detail/bessel_h.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,9 +6,8 @@
*/
//======================================================================================================================
#pragma once
#include <kyosu/types/impl/detail/besselj0.hpp>
#include <kyosu/types/impl/detail/bessely0.hpp>
#include <kyosu/types/impl/detail/bessely1.hpp>
#include <kyosu/types/impl/detail/bessel_j.hpp>
#include <kyosu/types/impl/detail/bessel_y.hpp>

namespace kyosu::_
{
Expand Down
29 changes: 27 additions & 2 deletions include/kyosu/types/impl/detail/besselin.hpp → include/kyosu/types/impl/detail/bessel_i.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,11 +6,36 @@
*/
//======================================================================================================================
#pragma once
#include <kyosu/types/impl/detail/besseljn.hpp>

namespace kyosu::_
{

//===-------------------------------------------------------------------------------------------
// cyl_bessel_i0
//===-------------------------------------------------------------------------------------------
template<typename Z>
auto dispatch(eve::tag_of<kyosu::cyl_bessel_i0>, Z z) noexcept
{
return cyl_bessel_j0(complex(-ipart(z), real(z)));
}

//===-------------------------------------------------------------------------------------------
// cyl_bessel_i1
//===-------------------------------------------------------------------------------------------
template<typename Z>
auto dispatch(eve::tag_of<kyosu::cyl_bessel_i1>, Z z) noexcept
{
if constexpr(concepts::complex<Z> )
{
auto [r, i] = cyl_bessel_j1(complex(-ipart(z), real(z)));
return complex(i, -r);
}
else
{
return cayley_extend(cyl_bessel_i1, z);
}
}

//===-------------------------------------------------------------------------------------------
// cyl_bessel_in
//===-------------------------------------------------------------------------------------------
Expand All @@ -33,5 +58,5 @@ namespace kyosu::_
{
return cayley_extend_rev(cyl_bessel_in, n, z);
}
}
}
}