mirror of
https://github.com/logos-storage/plonky2.git
synced 2026-02-26 16:53:12 +00:00
power function complete
This commit is contained in:
Dmitry Vagner
parent
779a1a3f98
commit
77ec96f6f4
@ -70,7 +70,7 @@ miller_end:
|
||||
|
||||
miller_one:
|
||||
// stack: 0xnm, times, O, P, Q, out, retdest
|
||||
DUP1 %gt_const(0x10)
|
||||
DUP1 %lt_const(0x10)
|
||||
// stack: skip?, 0xnm, times, O, P, Q, out, retdest
|
||||
%jumpi(miller_zero)
|
||||
// stack: 0xnm, times, O, P, Q, out, retdest
|
||||
|
||||
@ -16,15 +16,17 @@
|
||||
/// return y2 * y1 * y0
|
||||
|
||||
global power:
|
||||
// stack: sqr, out, retdest
|
||||
// stack: sqr, out, retdest
|
||||
PUSH 1 DUP1 DUP1
|
||||
// stack: 1, 1, 1, sqr, out, retdest
|
||||
// stack: 1, 1, 1, sqr, out, retdest
|
||||
%mstore_kernel_general(200) %mstore_kernel_general(212) %mstore_kernel_general(224)
|
||||
// stack: sqr, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
// stack: sqr, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
PUSH power_loop_2 PUSH power_loop_1 PUSH power_return
|
||||
// stack: power_return, power_loop_1, power_loop_2, sqr, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
// stack: power_return, power_loop_1, power_loop_2, sqr, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
SWAP3
|
||||
// stack: sqr, power_loop_1, power_loop_2, power_return, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
// stack: sqr, power_loop_1, power_loop_2, power_return, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
PUSH 65 PUSH 62 PUSH 65
|
||||
// stack: 65, 62, 65, sqr, power_loop_1, power_loop_2, power_return, out, retdest {200: y0, 212: y1, 224: y2}
|
||||
%jump(power_loop_0)
|
||||
|
||||
power_return:
|
||||
@ -67,37 +69,167 @@ power_return_4:
|
||||
// stack: 212, 236, out, power_return_5, retdest {236: y0, 212: y1, 224: y2}
|
||||
%jump(mul_fp12)
|
||||
power_return_5:
|
||||
// stack: out, retdest {236: y0, 212: y1, 224: y2}
|
||||
// stack: out, retdest {236: y0, 212: y1, 224: y2}
|
||||
PUSH 224 DUP2
|
||||
// stack: out, 224, out, retdest {236: y0, 212: y1, 224: y2}
|
||||
// stack: out, 224, out, retdest {236: y0, 212: y1, 224: y2}
|
||||
%jump(mul_fp12)
|
||||
|
||||
/// def power_loop_0():
|
||||
/// for i in range(1, len4):
|
||||
/// if EXP4[-i]:
|
||||
/// abc = load(power_data_0)
|
||||
/// if a:
|
||||
/// y1 *= square
|
||||
/// if EXP2[-i]:
|
||||
/// if b:
|
||||
/// y2 *= square
|
||||
/// if EXP0[-i]:
|
||||
/// if c:
|
||||
/// y0 *= square
|
||||
/// square = square_fp12(square)
|
||||
/// y1 *= square
|
||||
///
|
||||
/// def power_loop_1():
|
||||
/// for i in range(len4, len2):
|
||||
/// if EXP2[-i]:
|
||||
/// ab = load(power_data_1)
|
||||
/// if a:
|
||||
/// y2 *= square
|
||||
/// if EXP0[-i]:
|
||||
/// if b:
|
||||
/// y0 *= square
|
||||
/// square = square_fp12(square)
|
||||
/// y2 *= square
|
||||
///
|
||||
/// def power_loop_2():
|
||||
/// for i in range(len2, len0):
|
||||
/// if EXP0[-i]:
|
||||
/// a = load(power_data_1)
|
||||
/// if a:
|
||||
/// y0 *= square
|
||||
/// square = square_fp12(square)
|
||||
/// y0 *= square
|
||||
|
||||
power_loop_0:
|
||||
// stack: i , j, k, sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: break?, i , j, k, sqr, retdest
|
||||
%jumpi(power_loop_0_end)
|
||||
// stack: i , j, k, sqr, retdest
|
||||
%sub_const(1)
|
||||
// stack: i-1, j, k, sqr, retdest
|
||||
DUP1 %mload_kernel_code(power_data_0)
|
||||
// stack: abc, i-1, j, k, sqr, retdest
|
||||
DUP1 %lt_const(100)
|
||||
// stack: skip?, abc, i-1, j, k, sqr, retdest
|
||||
%jumpi(power_loop_0_b)
|
||||
// stack: abc, i-1, j, k, sqr, retdest
|
||||
%sub_const(100)
|
||||
// stack: bc, i-1, j, k, sqr, retdest
|
||||
PUSH power_loop_0_b PUSH 212 DUP1 DUP8
|
||||
// stack: sqr, 212, 212, power_loop_0_b, bc, i-1, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_0_b:
|
||||
// stack: bc, i, j, k, sqr, retdest
|
||||
DUP1 %lt_const(10)
|
||||
// stack: skip?, bc, i, j, k, sqr, retdest
|
||||
%jumpi(power_loop_0_c)
|
||||
// stack: bc, i, j, k, sqr, retdest
|
||||
%sub_const(10)
|
||||
// stack: c, i, j, k, sqr, retdest
|
||||
PUSH power_loop_0_c PUSH 224 DUP1 DUP8
|
||||
// stack: sqr, 224, 224, power_loop_0_c, c, i, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_0_c:
|
||||
// stack: c, i, j, k, sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: skip?, c, i, j, k, sqr, retdest
|
||||
%jumpi(power_loop_0_sq)
|
||||
// stack: c, i, j, k, sqr, retdest
|
||||
POP
|
||||
// stack: i, j, k, sqr, retdest
|
||||
PUSH power_loop_0_sq PUSH 200 DUP1 DUP7
|
||||
// stack: sqr, 200, 200, power_loop_0_sq, i, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_0_sq:
|
||||
// stack: i, j, k, sqr, retdest
|
||||
PUSH power_loop_0 DUP5 DUP1
|
||||
// stack: sqr, sqr, power_loop_0, i, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_0_end:
|
||||
// stack: 0, j, k, sqr, retdest
|
||||
POP
|
||||
// stack: j, k, sqr, retdest
|
||||
PUSH power_loop_1 PUSH 212 DUP1 DUP6
|
||||
// stack: sqr, 212, 212, power_loop_1, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
|
||||
power_loop_1:
|
||||
// stack: j , k, sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: break?, j , k, sqr, retdest
|
||||
%jumpi(power_loop_1_end)
|
||||
// stack: j , k, sqr, retdest
|
||||
%sub_const(1)
|
||||
// stack: j-1, k, sqr, retdest
|
||||
DUP1 %mload_kernel_code(power_data_1)
|
||||
// stack: ab, j-1, k, sqr, retdest
|
||||
DUP1 %lt_const(10)
|
||||
// stack: skip?, ab, j-1, k, sqr, retdest
|
||||
%jumpi(power_loop_1_b)
|
||||
// stack: ab, j-1, k, sqr, retdest
|
||||
%sub_const(10)
|
||||
// stack: b, j-1, k, sqr, retdest
|
||||
PUSH power_loop_1_b PUSH 224 DUP1 DUP7
|
||||
// stack: sqr, 224, 224, power_loop_1_b, b, j-1, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_1_b:
|
||||
// stack: b, j, k, sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: skip?, b, j, k, sqr, retdest
|
||||
%jumpi(power_loop_1_sq)
|
||||
// stack: b, j, k, sqr, retdest
|
||||
POP
|
||||
// stack: j, k, sqr, retdest
|
||||
PUSH power_loop_1_sq PUSH 200 DUP1 DUP6
|
||||
// stack: sqr, 200, 200, power_loop_1_sq, j, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_1_sq:
|
||||
// stack: j, k, sqr, retdest
|
||||
PUSH power_loop_1 DUP4 DUP1
|
||||
// stack: sqr, sqr, power_loop_1, j, k, sqr, retdest
|
||||
%jump(square_fp12)
|
||||
power_loop_1_end:
|
||||
// stack: 0, k, sqr, retdest
|
||||
POP
|
||||
// stack: k, sqr, retdest
|
||||
PUSH power_loop_2 PUSH 224 DUP1 DUP6
|
||||
// stack: sqr, 224, 224, power_loop_2, k, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
|
||||
|
||||
power_loop_2:
|
||||
// stack: k , sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: break?, k , sqr, retdest
|
||||
%jumpi(power_loop_2_end)
|
||||
// stack: k , sqr, retdest
|
||||
%sub_const(1)
|
||||
// stack: k-1, sqr, retdest
|
||||
DUP1 %mload_kernel_code(power_data_2)
|
||||
// stack: a, k-1, sqr, retdest
|
||||
DUP1 ISZERO
|
||||
// stack: skip?, a, k-1, sqr, retdest
|
||||
%jumpi(power_loop_2_sq)
|
||||
// stack: a, k-1, sqr, retdest
|
||||
POP
|
||||
// stack: k-1, sqr, retdest
|
||||
PUSH power_loop_2_sq PUSH 200 DUP1 DUP5
|
||||
// stack: sqr, 200, 200, power_loop_2_sq, k-1, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
power_loop_2_sq:
|
||||
// stack: k, sqr, retdest
|
||||
PUSH power_loop_2 DUP3 DUP1
|
||||
// stack: sqr, sqr, power_loop_2, k, sqr, retdest
|
||||
%jump(square_fp12)
|
||||
power_loop_2_end:
|
||||
// stack: 0, sqr, retdest
|
||||
POP
|
||||
// stack: sqr, retdest
|
||||
PUSH power_return PUSH 200 DUP1 DUP4
|
||||
// stack: sqr, 200, 200, power_return, sqr, retdest
|
||||
%jump(mul_fp12)
|
||||
|
||||
Loading…
x
Reference in New Issue
Block a user