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

Fix some errors with rbo

This commit is contained in:
Dankrad Feist 2022-10-03 19:35:34 +01:00
parent 4f7cf29e5f
commit b611c51874
No known key found for this signature in database
GPG Key ID: 6815E6A20BEBBABA
1 changed files with 10 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
min-src/c_kzg_4844.c
View File

@ -801,9 +801,16 @@ C_KZG_RET load_trusted_setup(KZGSettings *out, FILE *in) {
TRY(new_fft_settings((FFTSettings*)out->fs, max_scale)); TRY(new_fft_settings((FFTSettings*)out->fs, max_scale));
//g1_t *g1_origorder;
//TRY(new_g1_array(&g1_origorder, out->length));
TRY(fft_g1(out->g1_values, g1_projective, true, out->length, out->fs)); TRY(fft_g1(out->g1_values, g1_projective, true, out->length, out->fs));
TRY(reverse_bit_order(out->g1_values, sizeof(g1_t), out->length));
printf("Reverse bit ordered\n");
free(g1_projective); free(g1_projective);
//free(g1_origorder);
return C_KZG_OK; return C_KZG_OK;
@ -1014,19 +1021,19 @@ C_KZG_RET evaluate_polynomial_in_evaluation_form(BLSFieldElement *out, const Pol
TRY(new_fr_array(&inverses_in, p->length)); TRY(new_fr_array(&inverses_in, p->length));
TRY(new_fr_array(&inverses, p->length)); TRY(new_fr_array(&inverses, p->length));
for (i = 0; i < p->length; i++) { for (i = 0; i < p->length; i++) {
if (fr_equal(x, &roots_of_unity[i * stride])) { if (fr_equal(x, &roots_of_unity[i])) {
*out = p->values[i]; *out = p->values[i];
free(inverses_in); free(inverses_in);
free(inverses); free(inverses);
return C_KZG_OK; return C_KZG_OK;
} }
fr_sub(&inverses_in[i], x, &roots_of_unity[i * stride]); fr_sub(&inverses_in[i], x, &roots_of_unity[i]);
} }
TRY(fr_batch_inv(inverses, inverses_in, p->length)); TRY(fr_batch_inv(inverses, inverses_in, p->length));
*out = fr_zero; *out = fr_zero;
for (i = 0; i < p->length; i++) { for (i = 0; i < p->length; i++) {
fr_mul(&tmp, &inverses[i], &roots_of_unity[i * stride]); fr_mul(&tmp, &inverses[i], &roots_of_unity[i]);
fr_mul(&tmp, &tmp, &p->values[i]); fr_mul(&tmp, &tmp, &p->values[i]);
fr_add(out, out, &tmp); fr_add(out, out, &tmp);
} }