status-im
/
c-kzg-4844
mirror of https://github.com/status-im/c-kzg-4844.git
2
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Remove tmp2 

This is possible assuming the fr_ ops don't need to read their inputs
after writing to their output.
This commit is contained in:
Ramana Kumar 2022-09-16 14:28:10 +01:00
parent 8eda02ab72
commit 12c3cb7645
No known key found for this signature in database
GPG Key ID: ED471C788B900433
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/poly.c
View File

@ -104,17 +104,17 @@ void eval_poly(fr_t *out, const poly *p, const fr_t *x) {
// TODO: optimize via batch inversion
void eval_poly_l(fr_t *out, const poly_l *p, const fr_t *x, const FFTSettings *fs) {
fr_t tmp, tmp2;
fr_t tmp;
uint64_t i;
const uint64_t stride = fs->max_width / p->length;
*out = fr_zero;
for (i = 0; i < p->length; i++) {
fr_sub(&tmp, x, &fs->expanded_roots_of_unity[i * stride]);
fr_inv(&tmp2, &tmp);
fr_mul(&tmp, &tmp2, &fs->expanded_roots_of_unity[i * stride]);
fr_mul(&tmp2, &tmp, &p->values[i]);
fr_add(out, out, &tmp2);
fr_inv(&tmp, &tmp);
fr_mul(&tmp, &tmp, &fs->expanded_roots_of_unity[i * stride]);
fr_mul(&tmp, &tmp, &p->values[i]);
fr_add(out, out, &tmp);
}
fr_from_uint64(&tmp, p->length);
fr_div(out, out, &tmp);