Fix strided FFT

This commit is contained in:
Ben Edgington 2021-02-10 11:15:35 +00:00
parent ea18b23bf5
commit 8986f1cde3
2 changed files with 2 additions and 2 deletions

View File

@ -61,7 +61,7 @@ C_KZG_RET fft_fr(blst_fr *out, const blst_fr *in, const FFTSettings *fs, bool in
fr_from_uint64(&inv_len, n);
blst_fr_eucl_inverse(&inv_len, &inv_len);
fft_fr_fast(out, in, 1, fs->reverse_roots_of_unity, stride, n);
for (uint64_t i = 0; i < fs->max_width; i++) {
for (uint64_t i = 0; i < n; i++) {
blst_fr_mul(&out[i], &out[i], &inv_len);
}
} else {

View File

@ -60,7 +60,7 @@ C_KZG_RET fft_g1(blst_p1 *out, blst_p1 *in, FFTSettings *fs, bool inv, uint64_t
fr_from_uint64(&inv_len, n);
blst_fr_eucl_inverse(&inv_len, &inv_len);
fft_g1_fast(out, in, 1, fs->reverse_roots_of_unity, stride, n);
for (uint64_t i = 0; i < fs->max_width; i++) {
for (uint64_t i = 0; i < n; i++) {
p1_mul(&out[i], &out[i], &inv_len);
}
} else {