RuntimeError: β = 1050 > d = 1049 for n=1030, q=2^64

[ElenaKirshanova]

Executing

params=LWE.Parameters(n=1030, q=2^64, Xs=ND.UniformMod(2), Xe=ND.DiscreteGaussian(stddev=3.19)) primal_usvp_cost = LWE.primal_usvp(params, red_cost_model=RC.BDGL16)

does not pass sanity_check():

RuntimeError: β = 1050 > d = 1049

Oddly, for the default red_cost_model (red_cost_model = RC.MATZOV) the issue does not appear.

Potential fix: checking inside _solve_for_d

lattice-estimator/estimator/lwe_primal.py

Line 40 in eae2fdb

def _solve_for_d(params, m, beta, tau, xi):

that the returned d is at least beta.

[ElenaKirshanova]

Just noticed that it can be related to #87 (comment)