1package gates
 2
 3import (
 4	"fmt"
 5	"regexp"
 6	"strconv"
 7
 8	"github.com/consensys/gnark/frontend"
 9	gl "github.com/succinctlabs/gnark-plonky2-verifier/goldilocks"
10)
11
12var mulExtensionGateRegex = regexp.MustCompile("MulExtensionGate { num_ops: (?P<numOps>[0-9]+) }")
13
14func deserializeMulExtensionGate(parameters map[string]string) Gate {
15	// Has the format "MulExtensionGate { num_ops: 13 }"
16	numOps, hasNumOps := parameters["numOps"]
17	if !hasNumOps {
18		panic("Missing field num_ops in MulExtensionGate")
19	}
20
21	numOpsInt, err := strconv.Atoi(numOps)
22	if err != nil {
23		panic("Invalid num_ops field in MulExtensionGate")
24	}
25
26	return NewMultiplicationExtensionGate(uint64(numOpsInt))
27}
28
29type MultiplicationExtensionGate struct {
30	numOps uint64
31}
32
33func NewMultiplicationExtensionGate(numOps uint64) *MultiplicationExtensionGate {
34	return &MultiplicationExtensionGate{
35		numOps: numOps,
36	}
37}
38
39func (g *MultiplicationExtensionGate) Id() string {
40	return fmt.Sprintf("MulExtensionGate { num_ops: %d }", g.numOps)
41}
42
43func (g *MultiplicationExtensionGate) wiresIthMultiplicand0(i uint64) Range {
44	return Range{3 * gl.D * i, 3*gl.D*i + gl.D}
45}
46
47func (g *MultiplicationExtensionGate) wiresIthMultiplicand1(i uint64) Range {
48	return Range{3*gl.D*i + gl.D, 3*gl.D*i + 2*gl.D}
49}
50
51func (g *MultiplicationExtensionGate) wiresIthOutput(i uint64) Range {
52	return Range{3*gl.D*i + 2*gl.D, 3*gl.D*i + 3*gl.D}
53}
54
55func (g *MultiplicationExtensionGate) EvalUnfiltered(
56	api frontend.API,
57	glApi *gl.Chip,
58	vars EvaluationVars,
59) []gl.QuadraticExtensionVariable {
60	const0 := vars.localConstants[0]
61	constraints := []gl.QuadraticExtensionVariable{}
62	for i := uint64(0); i < g.numOps; i++ {
63		multiplicand0 := vars.GetLocalExtAlgebra(g.wiresIthMultiplicand0(i))
64		multiplicand1 := vars.GetLocalExtAlgebra(g.wiresIthMultiplicand1(i))
65		output := vars.GetLocalExtAlgebra(g.wiresIthOutput(i))
66
67		mul := glApi.MulExtensionAlgebra(multiplicand0, multiplicand1)
68		computed_output := glApi.ScalarMulExtensionAlgebra(const0, mul)
69
70		diff := glApi.SubExtensionAlgebra(output, computed_output)
71		for j := 0; j < gl.D; j++ {
72			constraints = append(constraints, diff[j])
73		}
74	}
75	return constraints
76}