1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
use super::{check_bit_consistency, extract_bit};
use crate::fft::{Evaluations, Polynomial};
use dusk_bls12_381::Scalar;
use dusk_jubjub::EDWARDS_D;
#[derive(Debug, Eq, PartialEq)]
pub struct ProverKey {
pub q_l: (Polynomial, Evaluations),
pub q_r: (Polynomial, Evaluations),
pub q_c: (Polynomial, Evaluations),
pub q_ecc: (Polynomial, Evaluations),
}
impl ProverKey {
#[allow(clippy::too_many_arguments)]
pub(crate) fn compute_quotient_i(
&self,
index: usize,
ecc_separation_challenge: &Scalar,
w_l_i: &Scalar,
w_l_i_next: &Scalar,
w_r_i: &Scalar,
w_r_i_next: &Scalar,
w_o_i: &Scalar,
w_4_i: &Scalar,
w_4_i_next: &Scalar,
) -> Scalar {
let q_ecc_i = &self.q_ecc.1[index];
let q_c_i = &self.q_c.1[index];
let kappa = ecc_separation_challenge.square();
let kappa_sq = kappa.square();
let kappa_cu = kappa_sq * kappa;
let x_beta = &self.q_l.1[index];
let y_beta = &self.q_r.1[index];
let acc_x = w_l_i;
let acc_x_next = w_l_i_next;
let acc_y = w_r_i;
let acc_y_next = w_r_i_next;
let xy_alpha = w_o_i;
let accumulated_bit = w_4_i;
let accumulated_bit_next = w_4_i_next;
let bit = extract_bit(accumulated_bit, accumulated_bit_next);
let bit_consistency = check_bit_consistency(bit);
let y_alpha = bit.square() * (y_beta - Scalar::one()) + Scalar::one();
let x_alpha = bit * x_beta;
let xy_consistency = ((bit * q_c_i) - xy_alpha) * kappa;
let x_3 = acc_x_next;
let lhs = x_3 + (x_3 * xy_alpha * acc_x * acc_y * EDWARDS_D);
let rhs = (acc_x * y_alpha) + (acc_y * x_alpha);
let x_acc_consistency = (lhs - rhs) * kappa_sq;
let y_3 = acc_y_next;
let lhs = y_3 - (y_3 * xy_alpha * acc_x * acc_y * EDWARDS_D);
let rhs = (acc_y * y_alpha) + (acc_x * x_alpha);
let y_acc_consistency = (lhs - rhs) * kappa_cu;
let identity = bit_consistency + x_acc_consistency + y_acc_consistency + xy_consistency;
identity * q_ecc_i * ecc_separation_challenge
}
#[allow(clippy::too_many_arguments)]
pub(crate) fn compute_linearisation(
&self,
ecc_separation_challenge: &Scalar,
a_eval: &Scalar,
a_next_eval: &Scalar,
b_eval: &Scalar,
b_next_eval: &Scalar,
c_eval: &Scalar,
d_eval: &Scalar,
d_next_eval: &Scalar,
q_l_eval: &Scalar,
q_r_eval: &Scalar,
q_c_eval: &Scalar,
) -> Polynomial {
let q_ecc_poly = &self.q_ecc.0;
let kappa = ecc_separation_challenge.square();
let kappa_sq = kappa.square();
let kappa_cu = kappa_sq * kappa;
let x_beta_eval = q_l_eval;
let y_beta_eval = q_r_eval;
let acc_x = a_eval;
let acc_x_next = a_next_eval;
let acc_y = b_eval;
let acc_y_next = b_next_eval;
let xy_alpha = c_eval;
let accumulated_bit = d_eval;
let accumulated_bit_next = d_next_eval;
let bit = extract_bit(accumulated_bit, accumulated_bit_next);
let bit_consistency = check_bit_consistency(bit);
let y_alpha = bit.square() * (y_beta_eval - Scalar::one()) + Scalar::one();
let x_alpha = x_beta_eval * bit;
let xy_consistency = ((bit * q_c_eval) - xy_alpha) * kappa;
let x_3 = acc_x_next;
let lhs = x_3 + (x_3 * xy_alpha * acc_x * acc_y * EDWARDS_D);
let rhs = (x_alpha * acc_y) + (y_alpha * acc_x);
let x_acc_consistency = (lhs - rhs) * kappa_sq;
let y_3 = acc_y_next;
let lhs = y_3 - (y_3 * xy_alpha * acc_x * acc_y * EDWARDS_D);
let rhs = (x_alpha * acc_x) + (y_alpha * acc_y);
let y_acc_consistency = (lhs - rhs) * kappa_cu;
let a = bit_consistency + x_acc_consistency + y_acc_consistency + xy_consistency;
q_ecc_poly * &(a * ecc_separation_challenge)
}
}