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

add Rémy Oudompheng's implementation of the BMSS algorithm #36285

Merged
merged 4 commits into from
Oct 8, 2023

Conversation

yyyyx4
Copy link
Member

@yyyyx4 yyyyx4 commented Sep 17, 2023

The BMSS algorithm computes (the kernel polynomial of) a normalized isogeny from its degree $\ell$ and its domain and codomain in time $\widetilde O(\ell)$. Currently Sage uses Stark's algorithm for the same task, which takes time in $\Omega(\ell^2)$. An implementation in Sage of the BMSS algorithm is available thanks to @remyoudompheng, and in this patch we add it to the Sage library (following Rémy's suggestion).

Benchmark results for isogenies over a ~280-bit prime field (red is Stark, blue is BMSS; axes are degree $\ell$ and seconds taken):

Benchmark results for {Stark, BMSS} algorithms

Co-authored-by: Rémy Oudompheng <remyoudompheng@gmail.com>
@kwankyu
Copy link
Collaborator

kwankyu commented Sep 17, 2023

The achievement is impressive.

I made one comment. Otherwise, lgtm.

Copy link
Collaborator

@kwankyu kwankyu left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks.

Copy link

@github-actions github-actions bot left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

kwankyu approved this PR

@remyoudompheng
Copy link
Contributor

Thanks for taking care of this. Please note that this algorithm requires that field characteristic is "large enough" (solve_linear_de needs to divide by integers up to requested precision). If I understand correctly, a similar limitation (p > 4l+6) exists for compute_isogeny_stark due to weierstrass_p.

@vbraun
Copy link
Member

vbraun commented Sep 21, 2023

merge conflict

SageMath version 10.2.beta4, Release Date: 2023-09-24
@yyyyx4
Copy link
Member Author

yyyyx4 commented Sep 26, 2023

Resolved merge conflict and added some new code to fail more gracefully when the limitation mentioned in #36285 (comment) is encountered.

@yyyyx4 yyyyx4 requested a review from kwankyu September 26, 2023 08:02
@github-actions
Copy link

Documentation preview for this PR (built with commit 81996bb; changes) is ready! 🎉

vbraun pushed a commit to vbraun/sage that referenced this pull request Sep 28, 2023
…gorithm

The BMSS algorithm computes (the kernel polynomial of) a *normalized*
isogeny from its degree $\ell$ and its domain and codomain in time
$\widetilde O(\ell)$. Currently Sage uses Stark's algorithm for the same
task, which takes time in $\Omega(\ell^2)$. An implementation in Sage of
the BMSS algorithm is available thanks to @remyoudompheng, and in this
patch we add it to the Sage library (following Rémy's suggestion).

Benchmark results for isogenies over a ~280-bit prime field (red is
Stark, blue is BMSS; axes are degree $\ell$ and seconds taken):

![Benchmark results for {Stark, BMSS} algorithms](https://github.com/sag
emath/sage/assets/84067835/2a9f47d3-f1f9-4bb1-b3f1-6539de2a7562)

URL: sagemath#36285
Reported by: Lorenz Panny
Reviewer(s): github-actions[bot], Kwankyu Lee, Lorenz Panny
vbraun pushed a commit to vbraun/sage that referenced this pull request Oct 1, 2023
…gorithm

    
The BMSS algorithm computes (the kernel polynomial of) a *normalized*
isogeny from its degree $\ell$ and its domain and codomain in time
$\widetilde O(\ell)$. Currently Sage uses Stark's algorithm for the same
task, which takes time in $\Omega(\ell^2)$. An implementation in Sage of
the BMSS algorithm is available thanks to @remyoudompheng, and in this
patch we add it to the Sage library (following Rémy's suggestion).

Benchmark results for isogenies over a ~280-bit prime field (red is
Stark, blue is BMSS; axes are degree $\ell$ and seconds taken):

![Benchmark results for {Stark, BMSS} algorithms](https://github.com/sag
emath/sage/assets/84067835/2a9f47d3-f1f9-4bb1-b3f1-6539de2a7562)
    
URL: sagemath#36285
Reported by: Lorenz Panny
Reviewer(s): github-actions[bot], Kwankyu Lee, Lorenz Panny
@vbraun vbraun merged commit 1194290 into sagemath:develop Oct 8, 2023
@yyyyx4 yyyyx4 deleted the public/bmss branch October 8, 2023 13:59
@mkoeppe mkoeppe added this to the sage-10.2 milestone Oct 8, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

5 participants