Selected Papers from the 2011 International Conference on Reconfigurable Computing and FPGAs (ReconFig 2011)
View this Special IssueResearch Article  Open Access
Lyndon Judge, Suvarna Mane, Patrick Schaumont, "A HardwareAccelerated ECDLP with HighPerformance Modular Multiplication", International Journal of Reconfigurable Computing, vol. 2012, Article ID 439021, 14 pages, 2012. https://doi.org/10.1155/2012/439021
A HardwareAccelerated ECDLP with HighPerformance Modular Multiplication
Abstract
Elliptic curve cryptography (ECC) has become a popular public key cryptography standard. The security of ECC is due to the difficulty of solving the elliptic curve discrete logarithm problem (ECDLP). In this paper, we demonstrate a successful attack on ECC over prime field using the Pollard rho algorithm implemented on a hardwaresoftware cointegrated platform. We propose a highperformance architecture for multiplication over prime field using specialized DSP blocks in the FPGA. We characterize this architecture by exploring the design space to determine the optimal integer basis for polynomial representation and we demonstrate an efficient mapping of this design to multiple standard prime field elliptic curves. We use the resulting modular multiplier to demonstrate lowlatency multiplications for curves secp112r1 and P192. We apply our modular multiplier to implement a complete attack on secp112r1 using a Nallatech FSBCompute platform with Virtex5 FPGA. The measured performance of the resulting design is 114 cycles per Pollard rho step at 100 MHz, which gives 878 K iterations per second per ECC core. We extend this design to a multicore ECDLP implementation that achieves 14.05 M iterations per second with 16 parallel point addition cores.
1. Introduction
Elliptic curve cryptosystems (ECC), independently introduced by Miller [1] and Koblitz [2], have now found significant place in the academic literature, practical applications, and security standards. Their popularity is mainly because their shorter key sizes offer high levels of security relative to other public key cryptosystems, such as RSA. The security of ECC relies on the difficulty of elliptic curve discrete logarithmic problem (ECDLP) [3]. By definition, ECDLP is to find an integer for two points and on an elliptic curve defined over a finite field such that Here, denotes the scalar multiplication with .
The Pollard rho method [4] is the strongest known attack against ECC today. This method solves ECDLP by generating points on the curve iteratively using a pseudorandom iteration function such that . Since the elliptic curve is defined over a finite field, is finite and the walk will eventually encounter the same point twice resulting in a collision. When a collision occurs, the ECDLP is solved (for wellchosen form of ; see Section 3). Several optimizations to the Pollard rho method have been proposed to allow independent parallel walks [5], better iteration functions [6, 7], and more efficient collision detection [8].
There have been several different approaches to implement Pollard rho algorithm on software and hardware platforms. Most of the solutions are implemented on software platforms using general purpose workstations, such as clusters of PlayStation3 [9], Cell CPUs [10], and GPUs [11, 12]. These software approaches are inherently limited by the sequential nature of software on the target platform.
Programmable hardware platforms are an attractive alternative to the above because they efficiently support parallelization. However, most of the FPGAbased solutions that have been proposed do not deal well with the control complexity of ECDLP. Instead, they focus on the efficient implementation of datapath operations and ignore the system integration aspect of the solution.
There has been little work in the area of supporting or accelerating a full Pollard rho algorithm on a hardwaresoftware platform. Our solution, therefore, goes one step further as we demonstrate the parallelized Pollard rho algorithm on FPGA along with its integration to a software driver. We start from a reference software implementation, and demonstrate an efficient, parallel implementation of the prime field arithmetic for primes of the form . We also present a novel highperformance architecture for modular multiplication that can be applied to a variety of standard prime field elliptic curves specified by the National Institute of Standards and Technology (NIST) [13] and Standards for Efficient Cryptography Group (SECG) [14].
Most existing ECDLP solutions for hardware platforms target curves over binary fields because the arithmetic is more hardware friendly. However, our work shows that implementation of prime field arithmetic on hardware can be as feasible as binary field arithmetic.
We use a computing platform by Nallatech, which consists of a quadcore Xeon core (E7310, 1.6 GHz) with a tightly coupled Virtex5 (xq5vsx240t) FPGA. The hardware runs at 100MHz and uses 5229 slices per ECC core.
2. Related Work
The software solution proposed by Bernstein on Cell CPUs is the fastest existing software solution to the ECDLP for secp112r1 curve [15]. It uses the negation map and noninteger polynomialbasis arithmetic to report speedup over a similar solution by Bos et al. [9]. Both of these software solutions use prime field arithmetic in an affine coordinate system, and they exploit the SIMD architecture and rich instruction set of the Cell CPU. Another software solution by Bos [10] describes the implementation of parallel Pollard rho algorithm on synergistic processor units of cell broadband engine architecture to approach the ECC2K130 Certicom challenge. Highperformance ECDLP solutions for GPUs platforms have also been proposed in [12].
Fan et al. propose the use of a normalbasis, binary field multiplication to implement a highspeed attack on ECC2K130 [16] using Spartan3 FPGAs. Fan's solution outperforms attacks on the same curve using GPUs [12] and Cell CPUs [10], demonstrating the suitability of FPGA platforms for solving large ECDLPs in binary fields. Another binary field solution, for the COPACOBANA [17] FPGA cluster, targets a 160bit curve [18]. Güneysu et al. propose an architecture to solve ECDLP over prime fields using FPGAs and analyze its estimated performance for different ECC curves ranging from 64to 160bit fields [19].
Among hardwarebased solutions, Meurice de Dormale et al. propose an FPGA solution to attack the ECC Certicom challenge for [20]. Though it discusses the hardwaresoftware integration aspect of the solution, the authors did not confirm if their system was operational. The authors of [19] discuss some aspects of systemlevel integration for their primefield ECDLP system. A threelayer hybrid distributed system is described by Majkowski et al. to solve ECDLP over binary field [21]. It uses the general purpose computers with FPGAs and integrates them with a main server at the top level.
The outline of the paper is as follows. In the next section, we discuss the background of the parallel Pollard rho algorithm. We present our highperformance architecture for modular multiplication and demonstrate its applicability to primes of the form in Section 4. In Section 5, we discuss the additional modular arithmetic units and we describe the complete hardwaresoftware integrated system architecture in Section 6. Section 7 shows implementation results, including measured performance for secp112r1 ECDLP, and we conclude the paper in Section 8.
3. Pollard Rho Algorithm
Let be a prime and . Given the elliptic curve over of order , let be a point of order . Solving the ECDLP requires finding an integer given two points , such that .
The Pollard rho algorithm [4] uses a pseudorandom iteration function to solve the ECDLP. It conducts a pseudorandom walk by starting from a random seed point on the curve, for random , and generating subsequent points using the iteration function . Since the elliptic curve is defined over a finite field, the iteration function will eventually produce the same point twice, resulting in a cycle.
The name of the algorithm, rho, expresses the Greek letter , which shows a walk ending in a cycle. Cycles can be efficiently detected using Floyd's cyclefinding method or Brent's cyclefinding algorithm [8]. The collision point, which gives the solution to the ECDLP, is located at the starting point of the cycle. Therefore, the underlying idea of this algorithm is to search for two distinct points on the curve such that . Due to the birthday paradox, assuming a random iteration function, the expected number of iterations to find collision is [4].
The effectiveness of Pollard's rho method depends on the randomness of the iteration function. As such, studies have been conducted to evaluate the strength of various proposed iteration functions [6, 7, 22]. Teske proposes an additive iteration function: where is an index function and each element , for random values [7]. Based on analysis by Teske [6, 7] and Bos et al. [22], the additive iteration function is more similar to a truly random walk than Pollard's original function and other proposed variants. Furthermore, an additive walk with is very close to a true randomness and achieves speedup of 1.25× over Pollard's iteration function [6]. Thus, we perform an additive walk and choose .
Using the additive walk, a collision occurs when two points are found such that Based on the definition of the index function , the collision points can be rewritten as where and . The solution can then be obtained as
3.1. Point Addition
Each iteration of the Pollard's rho algorithm requires a point addition of the current point, , with a precomputed combination of and , . Addition of two distinct points for an elliptic curve over a finite field , is defined geometrically as shown in Figure 1. To compute the addition, a line is drawn through the two points and . The line intersects exactly one other point of the elliptic curve. The intersection point is , which can be reflected across the axis to find .
Algebraically, adding two points and gives the sum , where Each point addition requires 2 multiplications modulo , 1 modular squaring, 5 modular subtractions, 1 modular addition, and 1 modular inversion. The sequence of arithmetic operations for a point addition is shown in Table 1.

By Fermat's Little Theorem, for . Thus, computing modular inversion requires modular exponentiation, making inversion the most computationally expensive operation of the point addition. Optimization at the system level can reduce the cost of inversion through the use of Montgomery's trick [23]. Montgomery's trick makes it possible to share the cost of an inversion among computations by performing inversions on vectors of points simultaneously. Since the computational cost of inversion is two orders of magnitude more than other arithmetic operations, this results in significant savings over computation of individual modular inversions. We apply Montgomery's trick by implementing a vectorized point addition datapath that computes random walks on a single ECC core and shares the inversion cost across all walks.
Despite optimizations, inversion remains an expensive operation. Since exponentiation is achieved using multiplication, the cost of inversion is directly tied to the cost of modular multiplication. Therefore, design of an efficient modular multiplier is an important aspect of the system to maximize overall performance.
3.2. Comparison of Field Arithmetic in and
Standards for elliptic curve cryptography have been defined for curves over both binary and prime fields [13, 14]. Binary fields can be represented with a polynomial basis or normal basis. In general, hardware implementations favor polynomial basis because it allows simplified reduction through the use of an irreducible polynomial of the form . Prior work has confirmed the idea that binary field curves are better suited for hardware platforms. A comparison of ECDLP engines by [20] shows a significant speedup for binary field ECDLP over prime field ECDLP implemented on FPGA platform. The performance discrepancy between ECDLPs for binary field and prime field curves is due to differences in the required arithmetic operations for each field and the use of search optimization techniques for selected curves.
The Pollard rho method is the best known attack on ECC for curves defined over both binary and prime fields, but the cost of the attack varies based on the properties of the finite field arithmetic. Elliptic curves over prime fields require conventional integer arithmetic operations followed by costly reduction modulo . However, the properties of curves over GF() allow optimizations that can reduce the cost of point addition. Additions and subtractions in the binary finite field are reduced to XOR operations with no carry. Furthermore, the cost of modular reduction is reduced to parallel XOR gates for binary fields. For comparison, reduction modulo requires sequential additions of carry bits and has a cycle cost approximately equal to that of multiplication. Although attacks on curves over can be considerably faster than attacks on curves over , NIST security standards recommend approximately equivalent field sizes for both binary and prime field curves at a given security level [13].
3.3. Parallelization
Van Oorschot and Wiener [5] described a parallelization technique that enables parallel walks on a single curve of Pollard rho algorithm to speed up the computation of ECDLP. The idea is to define a subset of as distinguished points (DPs), points which have a distinguishing characteristic. Each parallel walk starts from a distinct random seed point of the form and continues a Pollard rho walk until it encounters a point that satisfies the distinguishing property. Due to properties of the index function of the additive iteration function, a collision between two parallel walks at any point causes the walks to be identical from that point forward. Therefore, checking only DPs for collisions reduces communication and search overhead, while ensuring that any collision between walks is detected. Once a collision is detected, the server must derive the secret key using the collision point and the distinct seed points of the colliding walk. Therefore, whenever a DP is found, it is transmitted to the server along with the random seed point that generated the walk. The parallel Pollard rho method allows distribution of the random walks among multiple processing clients and sharing all DPs found with a central server that performs a collision search. This technique results in a linear speedup as the number of clients increases.
The expected number of DPs required to find a collision is a fraction of the expected path length. This depends on the density of DPs in a point set , which in turn depends on the chosen distinguishing property. Consider a distinguishing property defined as a point with bits fixed in its coordinate. This results in a probability of that a given point is a DP and total DPs on the curve. Each walk will require steps on average to find a DP. Based on the birthday paradox, the number DPs required to find a collision is Since each DP is generated from a walk with a unique random seed point, the number of parallel walks required to solve the ECDLP is also .
4. Modular Multiplication
We have designed a novel architecture for modular multiplication in a prime field. Typically, hardware solutions use binary field arithmetic, primarily due to the assumption that binary field avoids the costly carry propagation of prime field arithmetic. However, we demonstrate that prime field arithmetic can be efficiently implemented in hardware.
We target the 112bit secp112r1 elliptic curve in , but demonstrate that the proposed architecture can be generalized to perform multiplication modulo any prime of the form .
We perform prime field arithmetic using a redundant 128bit polynomial representation in an affine coordinate system, as proposed by Bernstein et al. in [15] and Bos et al. in [9]. We represent the integers in secp112r1 redundantly in the ring , where . This allows us to perform reduction modulo , rather than modulo . We perform all arithmetic over the 128bit field using following method.
Reduction for the unbalanced coefficient is a constant multiplication and an addition. Assuming , then mod . Integers represented redundantly mod can be converted into a canonical form in by multiplying with . Let and . Then, mod mod . Therefore, we start with a unique representation in GF(), perform arithmetic in , and canonicalize the results in to a unique representation by multiplying with . Similar redundant representation can be used to simplify reduction for any prime where by performing arithmetic modulo and canonicalizing results by multiplication with .
4.1. Modular Multiplication Algorithm
Bernstein et al. [15] use noninteger basis for polynomial representation of data to achieve an efficient software implementation. We choose integer basis representation to make the partial product computation uniform across all coefficients, which also makes the design scalable over larger fields. Each bit operand is represented as , where , is the length of the operand, and is the integer basis. For secp112r1 curve, and we discuss optimization of quantity in a subsequent section.
Our approach to modular multiplication is based on the design proposed by Güneysu and Paar in [25]. We use schoolbook multiplication to compute the product , for . Schoolbook multiplication requires computation of partial products in time, where . We parallelize computation of the partial products to perform multiplications in parallel by multiplying one field of one operand, , with all fields from the other operand, . As proposed by Güneysu and Paar in [25], placing in a shift register and in a rotating register ensures that each multiplier produces aligned results, that can be directly accumulated. Our design adds an important optimization of the reduction step. We integrate reduction modulo into the multiplication by multiplying the most significant field of by when rotating it to the least significant field. In order to accommodate the multiplication by of each field of , we represent each field of with bits.
Each partial product is an bit multiplication and produces a bit result. After accumulation, each field of the result can be up to bits. We reduce each field to bits using two parallel carry chains. The first carry chain is integrated into partial product accumulation and accumulates the lower bits of partial product with the upper bits of partial product . The second carry chain adds the upper bits of each accumulated partial product in column to the lower bits of the next accumulated partial product in column . In both carry chains, reduction modulo is achieve by multiplying the upper bits of the most significant field, , by before adding them to the least significant field, . The final bit result is obtained by concatenating the bit reduced outputs. The parallel carry chains allow us to overlap the multiplication and reduction stages of the the algorithm to achieve significantly lower latency relative to prior work [25].
4.2. Hardware Architecture
Figure 2 depicts our hardware architecture to implement modular multiplication for secp112r1 elliptic curve with general field size parameters ; , . This architecture targets a Xilinx Virtex 5 FPGA platform and makes use of dedicated DSP blocks for highperformance arithmetic. Each DSP block includes a bit multiplier and the FPGA fabric includes dedicated routing paths for highspeed connections between DSP blocks [26]. The DSP blocks allow each bit partial product multiplication and accumulation to be computed in a single clock cycle. The modular multiplication takes two 128bit inputs, and , divided into bit fields and produces a 128bit output of the product reduced modulo . There are DSP48 multipliers employed to compute partial products of bit coefficients. Since the partial products are required, it takes multiplication cycles to complete the full unreducedmultiplication result.
Reduction modulo adds to the cycle cost of a modular multiplication. By multiplying the shifting operand with 3, we perform the reduction in parallel with the multiplication. For the 128bit data field with and , it takes eight cycles of multiplication and 12 iterations of reduction. This results in a total cycle cost of 20 cycles per modular multiplication. However, our algorithm reduces this cost by overlapping reduction with partial product multiplication and accumulation using the previously described parallel carry adder chains. Thus, the cost has been reduced to 14 cycles, a significant improvement in latency over the architecture proposed by Güneysu and Paar in [25], which takes 70 clock cycles for 256bit modular multiplication.
4.3. Optimizations
To achieve optimal performance, we evaluate the impact of the integer basis for the polynomial representation on the latency of our design. The performance of the modular multiplier is heavily influenced by the value selected for , which determines the length of each field in the operands and the number of partial products to be computed. Our architecture relies on fast computation of partial products using the available DSP blocks on the FPGA; so performance is also impacted by the suitability of the selected field size to efficiently use these resources.
Modification of our architecture to support different field lengths for polynomial representation is straightforward requiring only that the number of multiplieradder columns be changed. The number of multiplieradder columns is given by , which is related to the field length . Increasing the field length requires fewer columns in the architecture and reduces the number of partial products computed. The results of this exploration will be discussed in Section 7.2.
4.4. Extension to Other Curves
Our modular multiplier is a generalized architecture that can be easily adapted to any overellipticcurve over , where has the form . We demonstrate this by adapting the architecture for NIST standard P192 curve, which uses the prime . The multiplier architecture is shown in Figure 3. The multiplication for this 192bit NIST curve is similar to that of secp112r1 curve, but the 192bit curve needs additional multiplicationadder columns. Additionally, 192 uses , rather than for secp112r1, which requires the following reduction operation. With polynomial field size , we define as the field containing , which is . For , the reduction modulo is mod .
We modify the modular multiplier architecture as follows to make it compatible with P192 curve.(i)Multiplication factor 3 is replaced by 1 in the rotation path of multiplication stage as well as in the addition reduction feedback path. This reduces the number of bits in each field of operand to bits and the length of partial product outputs to bits. (ii)In the rotation path, operand field is fed back to two and . An adder is introduced before for this purpose. Thus, . (iii)In addition reduction stage 1, carry bits of highest accumulator output are fed back to both and . Similarly in addition reduction stage 2, carry bits of are folded and added back to both and .
The adapted architecture has a latency of 17 clock cycles at 193 MHz for and has area of 364 Virtex 5 FPGA slices and 12 DSP cores. This is comparable to the multiplier performance for curve secp112r1 (given in Section 7.2) and demonstrates the flexibility of our design to accelerate modular multiplier for general standard prime field elliptic curve.
5. Additional Modular Arithmetic Units for ECDLP
The field operation of a point addition in corresponds to a Pollard rho step that consists of four subtractions, one addition, four modular multiplications, and one inversion. This requires design of modular adder/subtractor and modular inversion units. Subsequent sections explain the architecture of these arithmetic modules in detail.
5.1. Addition/Subtraction
We use an integer basis polynomial representation to perform addition and subtraction modulo using highspeed DSP blocks in the FPGA. We choose a field size of 32, which allows the use of four parallel DSP adder/subtractors to compute the sum of the two 128bit operands. Our design requires one cycle for the addition/subtraction and one additional cycle for reduction modulo .
5.2. Inversion
From Fermat's little theorem, it follows that the modular inverse of can be obtained by computing . Therefore, computation of the modular inverse requires exponentiation to the power, which is achieved through successive square and multiply operations. For secp112r1 curve, we start with and compute . We perform exponentiation using a variant of the left to right binary method, also known as square and multiply method [24]. This method scans a bitwise representation of the exponent left to right considering one bit at a time. A squaring is performed for each bit in the exponent, and an additional multiplication by the input operand is performed whenever the current bit is 1. The basic algorithm is shown as Algorithm 1. Due to dependencies between each operation, parallelization of the algorithm is not possible; it needs 112 squarings and 59 multiplications to compute an inversion for secp112r1 curve in 128bit arithmetic. We improve these figures by applying the following techniques to optimize the inversion operation.

5.2.1. Windowing Optimization
We apply the sliding window method [24] to reduce the number of squarings and multiplications required for inversion. The windowing method is an optimization of the square and multiply algorithm that considers multiple bits of the exponent simultaneously. Computing with window size requires precomputation of and the first odd powers of . This method considers up to bits of the exponent in each iteration and performs a multiplication with one of the precomputed powers of based on the value of those bits. The sliding window method for modular exponentiation is shown as Algorithm 2.

We achieve optimal performance for a window size of four, which allows inversion with only 108 squarings and 29 modular multiplications. This results in a total of 137 multiplication operations, rather than 171 operations without windowing, a speedup of 1.25.
5.2.2. Vector Inversion
To further reduce the cost of an inversion, we use Montgomery's trick [23], which allows multiple modular inverses to be computed together for significant latency savings. Montgomery's trick is based on the observation that given , the individual inverses can be easily computed. This allows computation of modular inverses using multiplications and one inversion. Although the cost of a single inversion is two orders of magnitude higher than multiplication, this approach allows the cost of inversion to be shared across operations so the marginal cost of an inversion becomes comparable to other arithmetic operations for large vector size . The algorithm for vector inversion using Montgomery's trick is shown as Algorithm 3.

We apply this method by implementing a vectorized inversion module that applies Algorithm 3 for any vector size . To take advantage of vectorized inversion, we design a vectorized point addition datapath that performs random walks on a vector of points simultaneously. We choose vector size , which yields speedup of 19× over 32 individual inversions, and leave optimization of this quantity for future work.
5.3. Squaring
An inversion involves a total of 137 modular multiplications out of which 75% are squaring operations. Having a dedicated squaring unit allows specific optimizations to take advantage of the properties of squaring and reduces the time required to compute an inversion. Consequently, a dedicated module optimized for squaring provides significant acceleration of the point addition operation.
Squaring is a special case of the general multiplication of two 128bit operands and . During squaring , all offdiagonal partial products, represented as for , are computed twice by the conventional schoolbook multiplication method used by our modular multiplier. To reduce computational cost of squaring, we avoid computing any partial product more than once by altering the operand control structure of our multiplier design as shown in Figure 4.
We represent the operand using the same integer basis for polynomial representation previously described for our modular multiplier. We copy into a shift register and a rotating register, each with bit fields for polynomial representation, to achieve aligned accumulation of partial products. As shown in Figure 4, the first cycle of multiplication computes the partial products on the diagonal, that is, where . All remaining partial products are offdiagonal. Since offdiagonal partial products are computed twice in schoolbook multiplication, we multiply each field of the rotating register by before computing these partial products. The multiplication by is implemented as a bit shift left at negligible cost. This reduces the number of partial product multiplications from to . The reduction stage for the square module is identical to that of multiplication.
6. ECDLP System Architecture
Our complete ECDLP system is implemented on a Nallatech cointegrated hardware software platform. Figure 5 depicts the architecture of the Nallatech system. It consists of one quadcore Xeon processor E7310 and three Virtex5 FPGAs (1 xc5vlx110, 2 xc5vsx240t). A fast North Bridge integrates highspeed components, including a Xeon, FPGA, and main memory. A slower South Bridge integrates peripherals into the system, including the hard disk. Both the Xeon and the FPGA can directly access system memory using a Front Side Bus (FSB). The FPGA performs the computationally expensive Pollard rho iterations, while the Xeon processor manages the central database of distinguished points and executes collision search. The communication between software and hardware is carried out only for the exchange of seed points and distinguished points, which minimizes the communication overhead.
6.1. Software Driver
The Xeon processor executes a software driver (in ) and manages software interface to FSB. The software driver mainly handles the communication interface with FPGA; seed point (SP) generation; storage and sorting of DPs. As shown in Figure 5, twoway communication between the Xeon and the FPGA takes place over the FSB.
When the program execution starts, the software calls APIs to configure FPGA card, to initialize the FSB link, and to allocate the workspace memory. It then generates random SPs on the curve and starts an attack by sending them to the FPGA over FSB. Every point has  and coordinates of 128bit length each.
The hardware finds a DP for each SP received and sends it to the software along with its corresponding SP. When the software receives SPDP pair from the FPGA, it performs a collision search among all the received DPs. Once a collision is detected, it computes the secret scalar and reports the solution to the ECDLP. As the software takes care of the central database of DPs, a collision search is conducted in parallel with hardware computations.
6.2. Hardware Accelerator
On the hardware side, as shown in Figure 6, the FPGA edge core provides an interface between the FSB and the ECC core. The edge core consists of control logic and two 256bit wide FIFOs. The RX FIFO buffers the incoming SPs received over FSB and the TX FIFO stores the DPs found for transmission to the Xeon.
The ECC core performs a random walk by computing the point addition operation iteratively until it finds a DP and stores that in the TX FIFO. We have defined DP as a point with 16 zeros in fixed positions such that and . The probability of a point being distinguished is almost exactly .
The distinguishing property of points allows to send only few points back to the Xeon, which reduces the communication overhead and minimizes the storage requirement in the hardware. The required bandwidth of communication bus is around 8 K bits/sec, which is well within the range of FSB. In the following are the details of key components in the design.
6.2.1. IO Controller
The IO controller manages the read/write interfaces of TXRX FIFOs and controls the ECC core operation. It loads the SPs from the RX FIFO to each ECC core and initiates a Pollard rho walk. When a DP is found, the IO controller halts the ECC core operation until a new SP is loaded from the RX FIFO. The computed DPs are buffered in the TX FIFO and then transferred to the Xeon along with corresponding SPs.
6.2.2. ECC Core
The ECC core consists of a microinstruction sequencer and the vectorized point addition (PA) datapath. The core operates on vectors of elements to perform independent random walks simultaneously. Each walk continues until it finds a DP or crosses the iteration limit, which is currently set to .
6.2.3. Vectorized Datapath
The datapath consists of modular arithmetic operators and memory. Each of these modules is designed to support vectorized point additions. A block diagram of the vectorized datapath is shown in Figure 7. As shown, each of the arithmetic modules interfaces directly to memory in order to load operands and store results. The modular add/sub, modular multiply, and distinguished point check modules are nonvectorized components and operate on inputs from one vector at a time. In order to complete each arithmetic operation over the vectorized walks, computation is serialized and repeated times for each of these operations. Conversely, the inversion module receives inputs from all vectors simultaneously and computes the inverse of all inputs together using the previously described vector inversion algorithm.
Based on the sequence of arithmetic operations required for point addition shown in Table 1, we need eight registers (,,,,,,,and ) to hold the intermediate results for a each point addition. For vector size of , we use entry register files to store these intermediate results. Each of the memories is 128 bits wide and has a depth equal to the vector size . The memory address corresponds to the vector index for each stored data value. We map these register files to the distributed RAMs available in the Xilinx Virtex 5 FPGA. We use an additional 256bit entry memory to store the original seed points for each parallel walk. When a distinguished point is found, that point, along with the original seed point stored in the SP RAM, is send to the IO controller for transmission to the software driver.
6.2.4. Microinstruction Sequencer
We implement a microinstruction sequencer that operates alongside the vectorized point addition datapath to provide control signals necessary to compute vectorized point addition. This allows control of the execution flow of the lowlevel arithmetic operations for a point addition without modification of the datapath. As shown in Figure 6, the sequencer consists of a microinstruction ROM, which stores instructions to implement point addition, and a next address logic unit (NALU), which determines instruction execution order based on outputs from the datapath. The interconnections between the point addition datapath and microinstruction sequencer are shown in Figure 7.
The vectorized point addition datapath requires a number of control signals to implement each arithmetic operation. A microcoded controller allows us to encode these various control signals into compact form, while maintaining maximum flexibility to modify design parameters including vector size and number of parallel ECC cores. It provides an efficient mechanism to separate our system's control and datapath. Since we use a very specific datapath to perform the particular sequence of operations for point addition, we include only microinstructions needed for the point addition. These correspond directly to the operations shown in Table 1. Common instructions required for general purpose microprocessors, such as jump, check flag and are not implemented for our design.
Our microinstruction format is shown in Figure 8. Each microinstruction includes three main segments, arithmetic operation selection, controlled by microinstruction bits 4–0, operand selection, controlled by bits 10–5, and register file write command, controlled by bits 20–11. Decoded control signals from each instruction are connected to datapath inputs as shown in Figure 7.
The next address logic unit (NALU) controls execution order of the microinstructions contained in the ROM. The NALU also handles vectorization of the point addition datapath by repeating each addition, subtraction, multiplication, and distinguished point check operation times for vector size . This ensures that the same microinstruction is applied sequentially to each element in the vector. Since the modular arithmetic operations have different latencies, the NALU relies on ready signals generated by datapath components to detect the completion of each instruction. This maximizes flexibility of the microinstruction sequencer and allows the NALU to be independent of the datapath implementation.
In a given cycle, the NALU reads a microinstruction and issues the corresponding control signals to point addition datapath. These control signals include a start pulse, operand selection, register file write, destination register select, and vector index. Unlike the other control signals, vector index is generated by a finitestate machine within the NALU, rather than decoded directly from microinstruction . To execute nonvectorized operations, that is all arithmetic operations other than inversion, the vector index is reset and control signals are applied to the datapath. The vector control finitestate machine increments the vector index when the datapath operation completes and asserts the control signals for the next vector element. This continues until the operation is completed for all vector elements. Then, the vector index resets and the NALU reads the next microinstruction from the ROM.
For an inversion operation, the NALU asserts control signals to load operands and store results from all vector elements. This allows loading of inputs into the inversion module and writing results into the corresponding register files. In this way, we efficiently implement vector inversion using Montgomery's trick within the point addition datapath. An inversion is performed only once per vectors, whereas the other instructions are repeated times.
6.2.5. Parallel Pollard Rho Walks
We have extended our design to support parallel Pollard rho walks by implementing multiple ECC cores on the FPGA. Each ECC core receives a separate set of SPs and performs a vectorized walk. We have extended the IO controller to support multiple ECC cores so that control signals from all cores are monitored. When a DP is found, the walk by the corresponding ECC core is halted while the IO controller retrieves the DP and loads a new SP. Since each ECC core contains its own microinstruction sequencer module, the multiple cores are not restricted to SIMD execution. When one core is in the process of sending DPs or receiving SPs, the other cores continue execution, avoiding costly latency associated with SIMD configuration.
Our design is scalable to support any number of parallel cores. The IO controller can be configured to support any number of cores. The total number of ECC cores in our design is only limited to the capacity of the FPGA. Our Nallatech computing platform includes two Xilinx Virtex 5 FPGA accelerators and both can be used simultaneously to implement ECC cores to maximize parallelization of Pollard rho walks to solve the ECDLP. Each core performs an independent walk and there is no additional overhead in the point addition iterations due to using multiple cores. Therefore, implementing multiple ECC cores produces a linear increase in performance.
7. Implementation Results
We have implemented our proposed system on the described Nallatech platform in both single and multicore variants. For demonstration purposes, the seed points that we generate are carefully chosen to be of order [15], which means we would need only steps to solve the ECDLP. This allows us to demonstrate collisions, proving that our solution works.
7.1. Overall Performance
The whole system runs at 100 MHz and uses 5229 slices which is 13% area of the Virtex5 device xq5vsx240t with a single ECC core. It takes 1.14 microseconds per Pollard rho step and can perform up to 878,000 iterations per second per ECC core.
Our design can easily be extended to include multiple ECC cores performing parallel walks. Each additional ECC core produces as linear increase in performance with negligible control overhead. With 16 ECC cores working in parallel, our system would achieve 14,050,000 iterations per second for secp112r1 ECDLP.
7.2. Evaluation of Modular Multiplier
We have modified the modular multiplier architecture shown in Figure 2 to vary the field length for the polynomial representation, , from 8 bits to 64 bits. Our results are given in Table 2. As shown, increasing the length of the field of the polynomial representation reduces the number of cycles required for the modular multiplication. However, using larger fields also degrades the maximum clock speed of the design. The 64bit field size requires the fewest cycles per modular multiplication, but achieves worse overall performance due to the low maximum clock frequency.

We also evaluate the impact of the modular multiplier design on overall point addition performance. These results are also shown in Table 2. Modular multiplication is the dominant operation in each point addition, with approximately 85% of iteration cycle count used for multiplication. Thus, we show that decreasing the cycle count of the modular multiplication operation has significant impact on the overall performance of our design. Our results show that best overall performance is achieved for field length of 16bits when the full system runs at 180 MHz. However, since maximum clock frequency of our complete design is 100 MHz, we select field length of 32 bits, which gives 1.3X increase in performance relative to 16 bit field size.
Güneysu and Paar's architecture described in [25] targets 256bit prime arithmetic over two fixed NIST primes. Our solution shows an improvement in terms of latency for the important ECC arithmetic operations. Assuming the cycle cost for 256bit arithmetic as twice of that for 128bit arithmetic (worstcase scenario), our architecture has cycle cost of 14 for a modular multiplication and 228 for the point addition. This is 5X and 3.5X times lower latency for a modular multiplication and point addition operation respectively, than those of the design in [25].
7.3. Comparison with Prior Work
Due to inherent differences between elliptic curves defined over binary and prime fields, performance comparison among various ECDLP implementations is not straightforward. Performance figures are also influenced by a variety of other factors, including target platform, curve size, and coordinate system. To minimize the impact of the target platform and coordinate system, we compare performance of different ECDLP implementations based on the total number of point addition operations performed per second.
Table 3 compares our solution with other software implementations. It shows that a multicore implementation of our design using 16 ECC cores achieves 1.6X improvement over the fastest existing software attack on secp112r1 curve.

Comparison with prior hardware ECDLP implementations is shown in Table 4. Our design achieves highperformance relative to other prime field solutions. When adjusting for curve size, we demonstrate significant speedup over prior designs by [18, 19].

We also show comparison with existing binary field solutions. However, it is difficult to make a fair comparison due to significant differences between binary and prime field arithmetic. In particular, properties of binary field arithmetic and reduction are wellsuited for hardware designs, which allows them to achieve much higher performance.
The solution proposed in [16] reports 111 M point addition iterations per second. The authors of [16] claim their system is capable of solving ECC2K130 within a year using five COPACOBANA machines. Similarly, the solution reported in [20] achieves 20 M iterations per second. We assume that the difference of performance figures exists due to factors including binary field arithmetic, different curve sizes, and use of pipelined architectures.
8. Conclusion
We successfully demonstrate a complete multicore ECC cryptanalytic machine to solve ECDLP on a hardwaresoftware cointegrated platform. We also implement a novel architecture on hardware to perform modular multiplication over prime field and this is the most efficient implementation reported at present for prime field multiplication. We then present a generalized modular multiplication architecture for primes of the form and demonstrate its application to NIST standard P192 curve. This work also demonstrates the use of microinstructionbased sequencing logic to support a vectorized point addition datapath with variable vector sizes. We compare our performance results with the previous implementations and show that a multicore implementation of our solution has competitive performance relative to existing hardware and software solutions.
Acknowledgments
An earlier version of this paper appeared as “An Integrated primefield ECDLP hardware accelerator with highperformance modular arithmetic units” in the Proceedings of the 2011 International Conference on Reconfigurable Computing and FPGAs (RECONFIG) [27]. This research was supported in part by the National Science Foundation Grant no. 0644070.
References
 V. Miller, “Use of elliptic curves in cryptography,” in Advances in Cryptology CRYPTO 85 Proceedings, H. Williams, Ed., vol. 218 of Lecture Notes in Computer Science, pp. 417–426, Springer, Berlin, Germany, 1986. View at: Publisher Site  Google Scholar
 N. Koblitz, “Elliptic curve cryptosystems,” Mathematics of Computation, vol. 48, no. 177, pp. 203–209, 1987. View at: Google Scholar
 I. Blake, G. Seroussi, N. Smart, and J. W. S. Cassels, Advances in Elliptic Curve Cryptography, London Mathematical Society Lecture Note Series, Cambridge University Press, New York, NY, USA, 2005.
 J. M. Pollard, “Monte carlo methods of index computation (mod p),” Mathematics of Computation, vol. 32, no. 143, pp. 918–924, 1978. View at: Google Scholar
 P. C. Van Oorschot and M. J. Wiener, “Parallel collision search with cryptanalytic applications,” Journal of Cryptology, vol. 12, no. 1, pp. 1–28, 1999. View at: Google Scholar
 E. Teske, “On random walks for Pollard's rho method,” Mathematics of Computation, vol. 70, no. 234, pp. 809–825, 2001. View at: Publisher Site  Google Scholar
 E. Teske, “Speeding up Pollard’s rho method for computing discrete logarithms,” in Algorithmic Number Theory, J. Buhler, Ed., vol. 1423 of Lecture Notes in Computer Science, pp. 541–554, Springer, Berlin, Germany, 1998. View at: Google Scholar
 R. P. Brent, “An improved Monte Carlo factorization algorithm,” BIT Numerical Mathematics, vol. 20, pp. 176–184, 1980. View at: Publisher Site  Google Scholar
 J. W. Bos, M. E. Kaihara, T. Kleinjung, A. K. Lenstra, and P. L. Montgomery, “On the security of 1024bit RSA and 160bit elliptic curve cryptography,” Report 2009/389, IACR Cryptology ePrint Archive, 2009, http://eprint.iacr.org/2009/389. View at: Google Scholar
 J. Bos, T. Kleinjung, R. Niederhagen, and P. Schwabe, “ECC2K130 on cell CPUs,” in Progress in Cryptology—AFRICACRYPT 2010, D. Bernstein and T. Lange, Eds., vol. 6055 of Lecture Notes in Computer Science, pp. 225–242, Springer, Berlin, Germany, 2010. View at: Google Scholar
 D. V. Bailey, L. Batina, D. J. Bernstein et al., “Breaking ECC2K130,” Report 2009/541, IACR Cryptology ePrint Archive, 2009, http://eprint.iacr.org/2009/541. View at: Google Scholar
 D. Bernstein, H.C. Chen, C.M. Cheng et al., “ECC2K130 on Nvidia GPUs,” in Progress in Cryptology—INDOCRYPT 2010, G. Gong and K. Gupta, Eds., vol. 6498 of Lecture Notes in Computer Science, pp. 328–346, Springer, Berlin, Germany, 2010. View at: Publisher Site  Google Scholar
 FIPS 1863: Digital Signature Standard (DSS), National Institute of Standards and Technology (NIST), June 2009, http://csrc.nist.gov/publications/fips/fips1863/fips_1863.pdf.
 SEC 2: Recommended Elliptic Curve Domain Parameters, Standards for Efficient Cryptography Group (SECG), January 2010, http://www.secg.org/download/aid784/sec2v2.pdf.
 D. Bernstein, T. Lange, and P. Schwabe, “On the correct use of the negation map in the Pollard rho method,” in Public Key Cryptography—PKC 2011, D. Catalano, N. Fazio, R. Gennaro, and A. Nicolosi, Eds., vol. 6571 of Lecture Notes in Computer Science, pp. 128–146, Springer, Berlin, Germany, 2011. View at: Publisher Site  Google Scholar
 J. Fan, D. V. Bailey, L. Batina, T. Güneysu, C. Paar, and I. Verbauwhede, “Breaking elliptic curve cryptosystems using reconfigurable hardware,” in Proceedings of the 20th International Conference on Field Programmable Logic and Applications (FPL '10), pp. 133–138, IEEE Computer Society, September 2010. View at: Publisher Site  Google Scholar
 S. Kumar, C. Paar, J. Pelzl, G. Pfeiffer, and M. Schimmler, “Breaking ciphers with COPACOBANA a costoptimized parallel code breaker,” in Cryptographic Hardware and Embedded Systems—CHES 2006, L. Goubin and M. Matsui, Eds., vol. 4249 of Lecture Notes in Computer Science, pp. 101–118, Springer, Berlin, Germany, 2006. View at: Publisher Site  Google Scholar
 T. Gueneysu, C. Paar, and J. Pelzl, “Attacking elliptic curve cryptosystems with specialpurpose hardware,” in Proceedings of the 15th ACM/SIGDA International Symposium on FieldProgrammable Gate Arrays (FPGA '07), pp. 207–215, ACM, New York, NY, USA, February 2007. View at: Publisher Site  Google Scholar
 T. Güneysu, C. Paar, and J. Pelzl, “Specialpurpose hardware for solving the elliptic curve discrete logarithm problem,” ACM Transactions on Reconfigurable Technology and Systems, vol. 1, no. 2, pp. 8:1–8:21, 2008. View at: Google Scholar
 G. Meurice de Dormale, P. Bulens, and J.J. Quisquater, “Collision search for elliptic curve discrete logarithm over GF(${2}^{m}$) with FPGA,” in Cryptographic Hardware and Embedded Systems—CHES 2007, P. Paillier and I. Verbauwhede, Eds., vol. 4727 of Lecture Notes in Computer Science, pp. 378–393, Springer, Berlin, Germany, 2007. View at: Publisher Site  Google Scholar
 P. Majkowski, T. Wojciechowski, M. Wojtyński, M. Rawski, and Z. Kotulski, “Heterogenic distributed system for cryptanalysis of elliptic curve based cryptosystems,” in Proceedings of the 19th International Conference on Systems Engineering (ICSEng '08), pp. 300–305, August 2008. View at: Publisher Site  Google Scholar
 J. W. Bos, M. E. Kaihara, T. Kleinjung, A. K. Lenstra, and P. L. Montgomery, “Solving a 112bit prime elliptic curve discrete logarithm problem on game consoles using sloppy reduction,” International Journal of Applied Cryptography, vol. 2, no. 3, pp. 212–228, 2012. View at: Google Scholar
 P. L. Montgomery, “Speeding the Pollard and elliptic curve methods of factorization,” Mathematics of Computation, vol. 48, no. 177, pp. 243–264, 1987. View at: Google Scholar
 D. R. Hankerson, S. A. Vanstone, and A. J. Menezes, Guide to Elliptic Curve Cryptography, Springer, New York, NY, USA, 2004.
 T. Güneysu and C. Paar, “Ultra high performance ECC over NIST primes on commercial FPGAs,” in Cryptographic Hardware and Embedded Systems—CHES 2008, E. Oswald and P. Rohatgi, Eds., vol. 5154 of Lecture Notes in Computer Science, pp. 62–78, Springer, Berlin, Germany, 2008. View at: Publisher Site  Google Scholar
 Virtex5 FPGA XtremeDSP Design Considerations, Xilinx, Inc., January 2010, http://www.xilinx.com/support/documentation/ user guides/ug193.pdf.
 S. Mane, L. Judge, and P. Schaumont, “An integrated primefield ECDLP hardware accelerator with highperformance modular arithmetic units,” in Proceedings of the International Conference on Reconfigurable Computing and FPGAs (ReConFig '11), P. M. Athanas, J. Becker, and R. Cumplido, Eds., pp. 198–203, IEEE Computer Society, December 2011. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 Lyndon Judge et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.