Update EIP-2565 to clarify specification

EIPS/eip-2565.md

@@ -4,33 +4,36 @@ title: Repricing of the EIP-198 ModExp precompile
 author: Kelly Olson (@ineffectualproperty), Sean Gulley (@sean-sn), Simon Peffers (@simonatsn), Justin Drake (@justindrake), Dankrad Feist (@dankrad)
 discussions-to: https://ethereum-magicians.org/t/big-integer-modular-exponentiation-eip-198-gas-cost/4150
 status: Last Call
-review-period-end: 2020-07-31
+review-period-end: 2020-10-09
 type: Standards Track
 category: Core
 created: 2020-03-20
 requires: 198
 ---
 
 ## Simple Summary
-The [EIP-198](https://eips.ethereum.org/EIPS/eip-198) ‘big integer modular exponentiation’, or `ModExp`, precompile is currently overpriced. Re-pricing this precompile will enable more cost efficient verification of RSA signatures, verifiable delay functions (VDFs), primality checks, and more.
+The [EIP-198](./eip-198.md) ‘big integer modular exponentiation’, or `ModExp`, precompile is currently overpriced. Re-pricing this precompile will enable more cost efficient verification of RSA signatures, verifiable delay functions (VDFs), primality checks, and more.
 
 ## Abstract
-After benchmarking the ModExp precompile, we discovered that it is ‘overpriced’ relative to other precompiles. We also discovered that the current gas pricing formula could be improved to better estimate the computational complexity of various ModExp input variables. To improve the gas cost pricing for this precompile the following options are available:
+After benchmarking the ModExp precompile, we discovered that it is ‘overpriced’ relative to other precompiles. We also discovered that the current gas pricing formula could be improved to better estimate the computational complexity of various ModExp input variables. To improve the gas cost pricing for this precompile, this EIP specifies:
 
-1. Modifying the gas pricing formula to better reflect the computational complexity of ModExp operations
-2. Changing the value of the `GQUADDIVISOR` parameter in the ModExp pricing formula to bring its costs more in-line with other precompiles
-3. Improving the underlying libraries beneath the ModExp Precompile
-4. Any combination of (1), (2), and (3)
-
-We recommend **Option (1) and (2)** which provides a large practical improvement to gas estimation while keeping implementation complexity low. Option (3) can be implemented at a later point to improve the gas pricing even further.
+1. A change to the `mult_complexity` formula to better reflect the computational complexity of ModExp operations
+2. A change to the value of the `GQUADDIVISOR` parameter in the ModExp pricing formula to bring its costs more in-line with other precompiles
+3. A minimum cost to call the precompile to prevent underpricing for small inputs
 
 ## Motivation
 Modular exponentiation is a foundational arithmetic operation for many cryptographic functions including signatures, VDFs, SNARKs, accumulators, and more. Unfortunately, the ModExp precompile is currently over-priced, making these operations inefficient and expensive. By reducing the cost of this precompile, these cryptographic functions become more practical, enabling improved security, stronger randomness (VDFs), and more.
 
 ## Specification
-The current gas pricing formula is defined in EIP-198. This formula divides a ‘computational complexity’ function by a ‘gas conversion’ parameter called `GQUADDIVISOR` to arrive at a gas cost: `floor(mult_complexity(x)/GQUADDIVISOR)`
+The current gas pricing formula is defined in [EIP-198](./eip-198.md) as follows: 
+
+```
+floor(mult_complexity(max(length_of_MODULUS, length_of_BASE)) * max(ADJUSTED_EXPONENT_LENGTH, 1) / GQUADDIVISOR)
+```
+
+As of `FORK_BLOCK_NUMBER` make the following changes to the pricing formula for the ModExp precompile:
 
-### **Recommended** Option (1): Modify ‘computational complexity’ function and add minimum gas cost
+### 1: Modify the `mult_complexity` function
 The current complexity function, as defined in EIP-198 is as follow:
 
 ```
@@ -45,50 +48,56 @@ This complexity formula was meant to approximate the difficulty of Karatsuba mul
 
 ```
 def mult_complexity(x):
-    ceiling(x/64)^2
+    ceiling(x/8)^2
 ```
-where is `x` is `max(length_of_MODULUS, length_of_BASE)`. `x` is divided by 64 to account for the number of limbs in multiprecision arithmetic.
-
-You can find comparison of these two complexity fomulas for the current test vectors as the following [spreadsheet](https://docs.google.com/spreadsheets/d/1-xBzA-2-l2ZQDQ1eh3XXGZjcRSBQ_Hnp7NubXpbiSUY/edit?usp=sharing).
+where is `x` is `max(length_of_MODULUS, length_of_BASE)`. `x` is divided by 8 to account for the number of limbs in multiprecision arithmetic.
 
-In addition to modifying the `mult_complexity` formula as above, we also recommend wrapping the entire function with a minimum gas price of 100 to ensure that a minimum amount of gas is used when the precompile is called e.g. `max(100,floor(mult_complexity(x)/GQUADDIVISOR))`
-
-### **Recommended** Option (2): Change value of GQUADDIVISOR
+### 2. Change value of `GQUADDIVISOR`
 `GQUADDIVISOR` is set to `20` per EIP-198. We recommend changing the value of this parameter to `3` to account for the changes in the recommended 'computational complexity' formula above.
 
-### Option (3): Replace libraries used by ModExp precompiles
-ModExp benchmarks for different libraries can be found at the following [spreadsheet](https://docs.google.com/spreadsheets/d/1Fq3d3wUjGN0R_FX-VPj7TKhCK33ac--P4QXB9MPQ8iw/edit?usp=sharing).
-
-While alternative libraries can provide improved a further 2-5x improvement in performance, this option is not recommended at this time.
+### 3. Set a minimum price for calling the precompile
+We recommend wrapping the entire function with a minimum gas price of 200 to ensure that a minimum amount of gas is used when the precompile is called e.g. `max(200,floor(mult_complexity(max(length_of_MODULUS, length_of_BASE)) * max(ADJUSTED_EXPONENT_LENGTH, 1) / GQUADDIVISOR))`
 
 ## Rationale
 
-### **Recommended** Option (1): Modify ‘computational complexity’ formula
-A comparison of the current ‘complexity’ function and the proposed function described above can be found at the following [spreadsheet](https://docs.google.com/spreadsheets/d/1-xBzA-2-l2ZQDQ1eh3XXGZjcRSBQ_Hnp7NubXpbiSUY/edit?usp=sharing).
+### 1. Modify ‘computational complexity’ formula to better reflect the computational complexity
+A comparison of the current ‘complexity’ function and the proposed function against the execution time can be seen below:
 
 ![Option 1 Graph](../assets/eip-2565/Complexity_Regression.png)
 
 The new complexity function has a better fit vs. the execution time when compared to the current complexity function. This better fit is because the new complexity formula accounts for the use of binary exponentiation algorithms that are used by ‘bigint’ libraries for large exponents. You may also notice the regression line of the proposed complexity function bisects the test vector data points. This is because the run time varies depending on if the modulus is even or odd.
 
-### **Recommended** Option (2): Change value of GQUADDIVISOR:
+### 2. Change the value of GQUADDIVISOR
 After changing the 'computational complexity' formula it is necessary to change `QGUADDIVSOR` to bring the gas costs inline with their runtime. We recommend changing the value from '20' to '3'. With this change, the cost of the ModExp precompile will have a higher cost (gas/second) than other precompiles such as ECRecover.
 
 ![Option 2 Graph](../assets/eip-2565/GQuad_Change.png)
 
-### Option (3): Improving the ModExp precompile implementations
-
-![Option 3 Graph](../assets/eip-2565/Library_Benchmarks.png)
-
-Replacing the underlying library can improve the performance of the ModExp precompile by 2x-4x for large exponents, but comes at a high implementation cost. We do not recommend this option at this time.
+### 3. Set a minimum gas cost to prevent abuse
+This prevents the precompile from underpricing small input values.
 
 ## Test Cases
-There are no changes to the underlying interface or arithmetic algorithms, so the existing test vectors can be reused. Gas values will need to be updated in the existing unit tests based on the final pricing decided in this EIP. This will ensure that the updated gas calculations are done correctly.
+There are no changes to the underlying interface or arithmetic algorithms, so the existing test vectors can be reused. Below is a table with the updated test vectors:
+
+| Test Case  | EIP-198 Pricing | New Pricing |
+| ------------- | ------------- | ------------- |
+| modexp_nagydani_1_square | 204  | 200  |
+| modexp_nagydani_1_qube | 204  | 200  |
+| modexp_nagydani_1_pow0x10001 | 3276  | 341  |
+| modexp_nagydani_2_square  | 665  | 200  |
+| modexp_nagydani_2_qube  | 665  | 200  |
+| modexp_nagydani_2_pow0x10001  | 10649  | 1365  |
+| modexp_nagydani_3_square  | 1894  | 341  |
+| modexp_nagydani_3_qube  | 1894  | 341  |
+| modexp_nagydani_3_pow0x10001  | 30310  | 5461  |
+| modexp_nagydani_4_square  | 5580  | 1365  |
+| modexp_nagydani_4_qube  | 5580  | 1365  |
+| modexp_nagydani_4_pow0x10001  | 89292  | 21845  |
+| modexp_nagydani_5_square  | 17868  | 5461  |
+| modexp_nagydani_5_qube  | 17868  | 5461  |
+| modexp_nagydani_5_pow0x10001  | 285900 | 87381  |
 
 ## Security Considerations
 The biggest security consideration for this EIP is creating a potential DoS vector by making ModExp operations too inexpensive relative to their computation time.
 
-## References
-[EIP-198](https://eips.ethereum.org/EIPS/eip-198) 
-
 ## Copyright
 Copyright and related rights waived via [CC0](https://creativecommons.org/publicdomain/zero/1.0/).