Abastro/mathematical-algorithms-cpp
Archived
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Clearing corrects everything

Browse source
This commit is contained in:
Abastro 2025-03-25 22:31:07 +09:00
parent 7123942c9d
commit 8eff17c1be
1 changed files with 5 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

7
src/karatsuba.cpp
View file

@ -133,6 +133,10 @@ void poly_mult_Karatsuba_step(const size_t deg_bnd, span<R> &a, span<R> &b,
return; return;
} }
// Result may be reused, so this needs clearing
for (auto &entry : result)
entry = 0;
const auto next_bnd = deg_bnd >> 1; const auto next_bnd = deg_bnd >> 1;
auto a0 = a.subspan(0, next_bnd); auto a0 = a.subspan(0, next_bnd);
auto a1 = a.subspan(next_bnd, next_bnd); auto a1 = a.subspan(next_bnd, next_bnd);
@ -144,8 +148,7 @@ void poly_mult_Karatsuba_step(const size_t deg_bnd, span<R> &a, span<R> &b,
auto prod0 = result.subspan(0, deg_bnd); auto prod0 = result.subspan(0, deg_bnd);
auto prod1 = result.subspan(deg_bnd, deg_bnd); auto prod1 = result.subspan(deg_bnd, deg_bnd);
auto prod_add_vec = vector<R>(deg_bnd); auto prod_add = buffer.subspan(next_bnd * 2, next_bnd * 2);
auto prod_add = span(prod_add_vec);
auto next_buffer = buffer.subspan(next_bnd * 4, next_bnd * 4); auto next_buffer = buffer.subspan(next_bnd * 4, next_bnd * 4);