CN115062565B - 一种低时延椭圆曲线点乘电路设计方法 - Google Patents
一种低时延椭圆曲线点乘电路设计方法 Download PDFInfo
- Publication number
- CN115062565B CN115062565B CN202210713984.2A CN202210713984A CN115062565B CN 115062565 B CN115062565 B CN 115062565B CN 202210713984 A CN202210713984 A CN 202210713984A CN 115062565 B CN115062565 B CN 115062565B
- Authority
- CN
- China
- Prior art keywords
- multiplier
- multiplication
- modular
- karatuba
- multipliers
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000013461 design Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 230000009467 reduction Effects 0.000 claims description 11
- 239000000758 substrate Substances 0.000 claims description 6
- 238000007792 addition Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种低时延椭圆曲线点乘电路设计方法,包括:设计Karatsuba‑Ofman模乘法器;两个乘数输入乘法器,进行Karatsuba拆解,形成三个子乘法器;三个子乘法器的结果进行Karatsuba拼接,形成乘法结果;在Karatsuba‑Ofman模乘法器基础上,设计低时延点乘架构;点乘计算的迭代部分通过反复调用改进的Karatsuba‑Ofman模乘法器实现,点乘计算的模逆部分通过模平方器、模四次方器、一些寄存器与改进的Karatsuba‑Ofman模乘法器共同实现;基于低时延ECC点乘架构,设计相应的最优调度策略。本发明电路在小型有限域上和大型有限域上均拥有高时钟频率。
Description
技术领域
本发明属于数字集成电路技术领域,更具体的说是涉及一种低时延椭圆曲线点乘电路设计方法。
背景技术
椭圆曲线密码(Elliptic Curve Cryptography,ECC)由Neal Koblitz和VictorMiller于1985年分别独立的提出。相比于其他公钥密码算法,椭圆曲线密码由于其安全性高,计算量小,被广泛应用于互联网网络安全领域。ECC的密码体制可以分为三层,其中群运算层的点乘运算是目前所有ECC协议的基础,同时点乘运算相较于其他运算计算量大,时间、资源开销较高,因此点乘运算的性能决定了ECC整体的性能,ECC点乘成为学术界和工业界的热点研究方向。
ECC点乘的实现方法有软件实现和硬件实现。软件实现方案为利用编写的计算机、嵌入式软件程序,计算ECC点乘。软件实现ECC点乘虽然开发难度小、周期短,但由于软件平台中通用处理器针对ECC点乘此种特定运算的计算效率一般,同时通用处理器反复读写内存的时间开销大,因此使用软件实现ECC点乘的计算速度低、计算时延大(百毫秒级至秒级)。
随着人类社会互联网的不断发展,在保证网络安全的同时,对ECC点乘的计算速度提出了更高的要求。利用硬件手段实现ECC点乘,并针对ECC点乘的特殊性,优化硬件的设计,使得硬件实现ECC具备高性能的特征,单次ECC点乘的计算时延相较于软件实现具有显著降低(微秒级)。
硬件实现ECC点乘的平台一般为专用集成电路(Application SpecificIntegrated Circuit,ASIC)和现场可编程门阵列(Field Programmable Gate Array,FPGA)。其中利用ASIC实现ECC点乘通常可以获得最佳的计算速度,但ASIC的开发、生产周期长,往往达到数年,同时ASIC生产、测试的成本极高。相比于通过ASIC实现ECC点乘,利用FPGA实现具有成本低、开发周期短的优点,同时利用FPGA实现ECC点乘,也可以获得较佳的计算速度,足以应对绝大多数非严苛场景的需求。因此,利用FPGA硬件实现ECC点乘获得了工业界和学术界的广泛共识。
Khan提出了一种高速低延时的FPGA椭圆曲线点乘设计。该设计通过有限状态机调度一个乘法器、一个模平方器与若干加法器实现点乘。该设计虽然在小型有限域GF(2163)上表现较好,但在大型有限域上GF(2571)工作频率下降高达38%。
Salarifard设计了一种基于预计算技术的高性能椭圆曲线点乘电路,通过一系列预先计算,将固定基点不同倍数的特征值先行计算出,并存储在内存中,后续点乘计算过程中可大幅降低计算时钟周期数,进而有效降低计算时延。但此设计仅可在椭圆曲线签名场景中加以实际应用,在大量其他应用场景中,基点并不是固定不变的,因此本设计的优势也就难以发挥。
李丽娟设计了一种基于Koblitz曲线的快速椭圆曲线点乘电路,采用流水线的双字转换器将点乘算子转换至τ基表达形式,再基于τ基表达形式进行快速点乘。该设计虽然计算时延很低,但Koblitz曲线只是通用椭圆曲线的特例,故在大量曲线未指定的一般化应用场景无法工作。
在实际应用中,针对不同安全级别的应用场景,椭圆曲线点乘需要支持位宽不同的多种有限域。而现有的高速椭圆曲线点乘设计(如Khan设计)在小型有限域的工作时钟频率较高,然而当其支持大型有限域时,工作时钟频率会出现明显的陡降。同时,一些设计为了片面追求点乘的性能,其点乘架构仅支持特殊的曲线(李丽娟设计)、特殊的应用场景(Salarifard设计),并不具有普适性。
因此,如何提供一种低时延椭圆曲线点乘电路设计方法成为了本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种低时延椭圆曲线点乘电路设计方法,电路不仅在小型有限域GF(2163)上拥有高时钟频率,在大型有限域上GF(2571)依然可以维持较高的主频,在小型有限域上可以同时运行两次点乘运算,提高硬件资源利用率。
为了实现上述目的,本发明采用如下技术方案:
一种低时延椭圆曲线点乘电路设计方法,包括如下步骤:
步骤S1:设计Karatsuba-Ofman模乘法器;将两个乘数输入乘法器后,进行一次Karatsuba拆解,形成三个小位宽的子乘法器;三个子乘法器的结果进行一次Karatsuba拼接,形成最终的乘法结果;
步骤S2:在Karatsuba-Ofman模乘法器的基础上,设计低时延点乘架构;点乘计算的迭代部分通过反复调用改进的Karatsuba-Ofman模乘法器实现,点乘计算的模逆部分通过模平方器、模四次方器、一些寄存器与改进的Karatsuba-Ofman模乘法器共同实现;
步骤S3:基于低时延ECC点乘架构,设计相应的最优调度策略,点乘计算的迭代部分每轮迭代调度紧凑,模逆部分的计算时钟周期数取决于点乘架构当前计算所在的有限域大小。
进一步的,步骤S1中设计Karatsuba-Ofman模乘法器的具体方法为:将两个2w位宽的乘数A、B输入乘法器后,分别被拆分为高w位与低w位的小型操作数,依次记为Ah,Al,Bh与Bl,拆分后满足关系A=Ah×xw+Al,B=Bh×xw+Bl;其中Ah为乘数A中的高w位,Al为乘数A中的低w位,Bh为乘数B中的高w位与Bl为乘数A中的低w位,x为二进制域GF(2m)多项式基表达形式的基底;四个小型操作数经过拆分、相加处理,作为输入分别输入三个w位的乘法器中:乘法器1输入Ah与Bh并计算AhBh,乘法器2输入(Ah+Al)与(Bh+Bl)并计算(Ah+Al)×(Bh+Bl),乘法器3输入Al与Bl并计算AlBl;乘法器的结果输出后,需要根据当前进行点乘计算的有限域进行对应有限域的模约减;如果工作在小型有限域GF(2163),则利用A的高w位与B的高w位存储一组GF(2163)上的乘数,用A的低w位与B的低w位存储另一组GF(2163)上的乘数,乘法器1与乘法器3的结果分别为第一组乘数的乘积和第二组乘数的乘积;乘法器1和乘法器3的结果经过有限域GF(2163)的模约减,可同时得到两个GF(2163)乘法最终的结果;如果工作在大型有限域GF(2571),则A和B分别为两个GF(2571)上的乘数,乘法器1、乘法器2和乘法器3的结果需要经过Karatsuba拼接,以产生Karatsuba-Ofman乘法结果:
AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
其中x为二进制域GF(2m)多项式基表达形式的基底;乘法结果经过有限域GF(2571)的模约减,得到乘法的最终结果。
进一步的,乘法器的输出通过两个多路选择器控制,当乘法器工作在小型有限域GF(2163),两个MUX选通两个小型有限域GF(2163)的模约减结果;当乘法器工作在大型有限域GF(2571),两个MUX选通GF(2571)的模约减结果。
进一步的,A,B的乘积得到的方法为:设A,B为GF(2m)中两个m位宽的乘数,一般的,m为奇数,令2w=m+1,,则A,B表示为A=Ahxw+A1,B=Bhxw+B1,Karatsuba-Ofman乘法推导为:
A×B=(Ahxw+A1)×(Bhxw+B1)
=AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
只需要计算三次位宽为w的小型乘法分别计算AhBh与(Ah+Al)×(Bh×Bl)AlBl,再通过一系列加法将三个小型乘法的结果进行拼接,形成上式中xw的系数,最终得到A,B的乘积。
进一步的,步骤S2中,模平方、模四次方均通过纯组合逻辑电路实现;设置有限状态机用于调度模乘、模平方、模四次方运算,增加额外的寄存器用于缓存数据通路的数据,增加多路选择器MUX用于控制数据流;在每次点乘的后处理部分,基于Itoh Tsujii算法,通过有限状态机循环迭代乘法、平方、四次方运算以计算Z1与(xPZ1Z2)的模逆,进而最终计算(xQ,yQ);其中Z1、Z2分别为Montgomery Ladder的迭代两点结果(射影坐标),xP为点乘基点的横坐标(仿射坐标),xQ与yQ分别为点乘结果的横坐标(仿射坐标)与纵坐标(仿射坐标)。
进一步的,步骤S3中,电路上电后,复位为空闲状态,收到计算点乘的使能信号后,跳转至初始化状态执行初始化X1←xP,Z1←1,X2←xP 4+b,Z2←xP 2,初始化完成后进入主循环;主循环的两种情况,通过先默认计算ki=1,最后再判断是否交换Z1与Z2、X1与X2的值实现;计算公式T←Z1,Z1←(X1Z2+X2Z1)2,X1←xPZ1+X1X2TZ2与T←X2,X2←X2 4+bZ2 4,Z2←T2Z2 2的各个模运算拆分至6个时钟周期依次计算;其中,xP为点乘基点的横坐标(仿射坐标),(X1,Y1,Z1)、(X2,Y2,Z2)分别为Montgomery Ladder的迭代过程中两点结果(射影坐标),T(Temporary)为中间临时变量。
本发明的有益效果在于:
1、本发明的方法是面向全部GF(2m)椭圆曲线的,因此可以避免只局限于Koblitz曲线这种特殊情况。
2、本发明的方法是不用事先指定一个固定的基点P的,每次点乘的基点是可变的,因此可以避免只局限于固定基点P这种特殊情况。
3、本发明的硬件设计是经过充分优化后的电路结构图,电路的关键路径得到了充分优化。同时Karatsuba-Ofman模乘法器本身的性能较好,工作频率较高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本实用新式的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明Karatsuba-Ofman模乘法器架构图;
图2为本发明低时延点乘架构图;
图3为本发明低时延点乘架构的状态机跳转图;
图4为本发明ki=1,ki+1=1时的迭代调度图;
图5为本发明ki=1,ki+1=0时的迭代调度图。
具体实施方式
下面将结合本发明的实施例中,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种低时延椭圆曲线点乘电路设计方法,包括如下步骤:
步骤S1:设计Karatsuba-Ofman模乘法器;将两个乘数输入乘法器后,进行一次Karatsuba拆解,形成三个小位宽的子乘法器,子乘法器不再进行Karatsuba拆解,子乘法器采用复杂度为O(n2)的经典乘法器架构;三个子乘法器的结果进行一次Karatsuba拼接,形成最终的乘法结果。
Karatsuba-Ofman乘法是一种通用的,具有明显优化效果的乘法算法。然而Karatsuba-Ofman乘法算法并非拆解的次数越多优化效果越好,而是存在最优的拆解方法。经过多次实验,本发明提出一种单层Karatsuba拆解深度的改进的Karatsuba-Ofman模乘法器。对于两个乘数,进行一次Karatsuba拆解,形成三个较小位宽的子乘法器,子乘法器不再进行Karatsuba-Ofman拆解,采用经典乘法方法加以实现。三个子乘法器的结果进行一次Karatsuba拼接,即可形成最终的乘法结果。
步骤S2:在Karatsuba-Ofman模乘法器的基础上,设计低时延点乘架构;点乘计算的迭代部分通过反复调用改进的Karatsuba-Ofman模乘法器实现,点乘计算的模逆部分通过模平方器、模四次方器、一些寄存器与改进的Karatsuba-Ofman模乘法器共同实现。
基于改进的Karatsuba-Ofman模乘法器,利用有限状态机发明了低时延ECC点乘架构。点乘计算的迭代部分通过反复调用改进的Karatsuba-Ofman模乘法器实现,点乘计算的模逆部分通过模平方器、模四次方器、一些寄存器与改进的Karatsuba-Ofman模乘法器共同实现。
步骤S3:基于低时延ECC点乘架构,设计相应的最优调度策略,点乘计算的迭代部分每轮迭代调度紧凑,模逆部分的计算时钟周期数取决于点乘架构当前计算所在的有限域大小。
基于低时延ECC点乘架构,发明了相应的最优调度策略,各个时钟周期所需执行的操作如图4和图5所示。点乘计算的迭代部分每轮迭代调度紧凑,仅需6个时钟周期,模逆部分的计算时钟周期数取决于点乘架构当前计算所在的有限域大小。
(一)Karatsuba-Ofman模乘法器
A,B的乘积得到的方法为:设A,B为GF(2m)中两个m位宽的乘数,一般的,m为奇数,令2w=m+1,,则A,B表示为A=Ahxw+A1,B=Bhxw+B1,Karatsuba-Ofman乘法推导为:
A×B=(Ahxw+A1)×(Bhxw+B1)
=AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
只需要计算三次位宽为w的小型乘法分别计算AhBh与(Ah+Al)×(Bh+Bl)AlBl,再通过一系列加法将三个小型乘法的结果进行拼接,形成上式中xw的系数,最终得到A,B的乘积
根据以上推导,本发明设计Karatsuba-Ofman模乘法器如图1所示。将两个2w位宽的乘数A、B输入乘法器后,分别被拆分为高w位与低w位的小型操作数,依次记为Ah,Al,Bh与Bl,拆分后满足关系A=Ah×xw+Al,B=Bh×xw+Bl;其中Ah为乘数A中的高w位,Al为乘数A中的低w位,Bh为乘数B中的高w位与Bl为乘数A中的低w位,x为二进制域GF(2m)多项式基表达形式的基底;四个小型操作数经过拆分、相加处理,作为输入分别输入三个w位的乘法器中:乘法器1输入Ah与Bh并计算AhBh,乘法器2输入(Ah+Al)与(Bh+Bl)并计算(Ah+Al)×(Bh+Bl),乘法器3输入Al与Bl并计算AlBl;乘法器的结果输出后,需要根据当前进行点乘计算的有限域进行对应有限域的模约减;如果工作在小型有限域GF(2163),则利用A的高w位与B的高w位存储一组GF(2163)上的乘数,用A的低w位与B的低w位存储另一组GF(2163)上的乘数,乘法器1与乘法器3的结果分别为第一组乘数的乘积和第二组乘数的乘积;乘法器1和乘法器3的结果经过有限域GF(2163)的模约减,可同时得到两个GF(2163)乘法最终的结果;如果工作在大型有限域GF(2571),则A和B分别为两个GF(2571)上的乘数,乘法器1、乘法器2和乘法器3的结果需要经过Karatsuba拼接,以产生Karatsuba-Ofman乘法结果:
AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
其中x为二进制域GF(2m)多项式基表达形式的基底;乘法结果经过有限域GF(2571)的模约减,得到乘法的最终结果。乘法器的输出通过两个多路选择器(MUX)控制,当乘法器工作在小型有限域GF(2163),两个MUX选通两个GF(2163)的模约减结果;当乘法器工作在大型有限域GF(2571),两个MUX选通GF(2571)的模约减结果。
(二)低时延点乘架构
基于Montgomery Ladder与投影坐标系转化,在GF(2m)上通用椭圆曲线y2+xy=x3+ax2+b的点乘计算算法如下表所示。
表1
表1中,点乘就是在已知基点P的横纵坐标、标量k的基础上,计算Q=kP。P的坐标为(xP,yP),Q的坐标为(xQ,yQ);为了便于硬件电路执行点乘,本发明通常把标量k以二进制的形式写出来,就是k=(kt-1,kt-2,…,k0)2;对于k的二进制形式,我们一般默认最高位是有意义的1,即条件kt-1=1;k的脚坐标是从t-1一直到0的,由于k的最高位为1,所以下面的循环部分,只需要从t-2开始进行迭代,直到i=0;椭圆曲线的方程为y2+xy=x3+ax2+b,b即为方程中的参数b。
在Karatsuba-Ofman模乘法器的基础上,本发明设计低时延点乘架构如图2所示。其中模平方、模四次方均可通过纯组合逻辑电路实现。同时,设置有限状态机用于调度模乘、模平方、模四次方运算,增加额外的寄存器用于缓存数据通路的数据,增加多路选择器MUX用于控制数据流。在每次点乘的后处理部分,基于Itoh Tsujii算法,通过有限状态机循环迭代乘法、平方、四次方运算以计算Z1与(xPZ1Z2)的模逆,进而最终计算(xQ,yQ)。
(三)基于低时延点乘架构的调度策略
低时延点乘架构的状态机跳转图如图3所示。电路上电后,复位为空闲状态,收到计算点乘的使能信号后,跳转至初始化状态执行表1中的初始化X1←xP,Z1←1,X2←xP 4+b,Z2←xP 2,初始化完成后进入主循环。主循环的两种情况,可以通过先默认计算ki=1,最后再判断是否交换Z1与Z2、X1与X2的值实现。计算公式T←Z1,Z1←(X1Z2+X2Z1)2,X1←xPZ1+X1X2TZ2与T←X2,X2←X2 4+bZ2 4,Z2←T2Z2 2的各个模运算拆分至6个时钟周期依次计算,形成6个步骤,拆分方法见图4。每轮迭代均为6个步骤,但由于初次进入主循环与主循环迭代过程存在数据通路的不一致,故额外设置首轮步骤1状态,首轮步骤2状态,在非首轮迭代过程中,仅存在无需交换交换Z1与Z2、X1与X2的的迭代保持步骤1、迭代保持步骤2与需交换交换Z1与Z2、X1与X2的的迭代翻转步骤1、迭代翻转步骤2。
图4为ki=1,ki+1=1时的迭代调度图,图5为ki=1,ki+1=0时的迭代调度图,其中MUL为模乘,SQR为模平方,QUA为模四次方,ADD为模加(即异或),虚线表示时钟周期。由于k的当前执行位与上一轮执行位相同,故ki=1,ki+1=1时完成一轮迭代后,进入迭代保持步骤;ki=1,ki+1=0时,当前执行位与上一轮执行位不同,故进入迭代翻转步骤。
本发明的算法(表1)是面向全部GF(2m)椭圆曲线的,因此可以避免只局限于Koblitz曲线这种特殊情况。
本发明的算法(表1)是不用事先指定一个固定的基点P的,每次点乘的基点是可变的,因此可以避免只局限于固定基点P这种特殊情况。
本发明的硬件设计(图2)是经过充分优化后的电路结构图(通过开发工具分析电路的关键路径,关键路径就是一个电路结构中最耗时的那一条路径,可以类比于木桶原理,最差的这条关键路径决定了电路总体的性能。利用工具分析这条路径的起点与终点,移除一部分逻辑至其他的路径,可以类比于把木桶的相对较长的木条锯下来补在较短路径上。当然只是类比,实际优化过程还需要考虑电路的前后数据依赖关系,数据通路是否可以复用等等。最终,得到了优化后的电路),电路的关键路径得到了充分优化。同时Karatsuba-Ofman模乘法器本身的性能较好,工作频率较高。
本发明设计了一种二进制域支持通用椭圆曲线的低时延点乘电路,该电路不仅在小型有限域GF(2163)上拥有高时钟频率,在大型有限域上GF(2571)依然可以维持较高的主频。同时,在小型有限域上可以同时运行两次点乘运算,提高硬件资源利用率。本发明面向通用ECC曲线,普适性更好。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (1)
1.一种低时延椭圆曲线点乘电路设计方法,其特征在于,包括如下步骤:
步骤S1:设计Karatsuba-Ofman模乘法器;将两个乘数输入乘法器后,进行一次Karatsuba拆解,形成三个小位宽的子乘法器;三个子乘法器的结果进行一次Karatsuba拼接,形成最终的乘法结果;
步骤S2:在Karatsuba-Ofman模乘法器的基础上,设计低时延点乘架构;点乘计算的迭代部分通过反复调用改进的Karatsuba-Ofman模乘法器实现,点乘计算的模逆部分通过模平方器、模四次方器、一些寄存器与改进的Karatsuba-Ofman模乘法器共同实现;
步骤S3:基于低时延ECC点乘架构,设计相应的最优调度策略,点乘计算的迭代部分每轮迭代调度紧凑,模逆部分的计算时钟周期数取决于点乘架构当前计算所在的有限域大小;
步骤S1中设计Karatsuba-Ofman模乘法器的具体方法为:将两个2w位宽的乘数A、B输入乘法器后,分别被拆分为高w位与低w位的小型操作数,依次记为Ah,Al,Bh与Bl,拆分后满足关系A=Ah×xw+Al,B=Bh×xw+Bl;其中Ah为乘数A中的高w位,Al为乘数A中的低w位,Bh为乘数B中的高w位与Bl为乘数A中的低w位,x为二进制域GF(2m)多项式基表达形式的基底;四个小型操作数经过拆分、相加处理,作为输入分别输入三个w位的乘法器中:乘法器1输入Ah与Bh并计算AhBh,乘法器2输入(Ah+Al)与(Bh+Bl)并计算(Ah+Al)×(Bh+Bl),乘法器3输入Al与Bl并计算AlBl;乘法器的结果输出后,需要根据当前进行点乘计算的有限域进行对应有限域的模约减;如果工作在小型有限域GF(2163),则利用A的高w位与B的高w位存储一组GF(2163)上的乘数,用A的低w位与B的低w位存储另一组GF(2163)上的乘数,乘法器1与乘法器3的结果分别为第一组乘数的乘积和第二组乘数的乘积;乘法器1和乘法器3的结果经过有限域GF(2163)的模约减,可同时得到两个GF(2163)乘法最终的结果;如果工作在大型有限域GF(2571),则A和B分别为两个GF(2571)上的乘数,乘法器1、乘法器2和乘法器3的结果需要经过Karatsuba拼接,以产生Karatsuba-Ofman乘法结果:
AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
其中x为二进制域GF(2m)多项式基表达形式的基底;乘法结果经过有限域GF(2571)的模约减,得到乘法的最终结果;
乘法器的输出通过两个多路选择器控制,当乘法器工作在小型有限域GF(2163),两个MUX选通两个小型有限域GF(2163)的模约减结果;当乘法器工作在大型有限域GF(2571),两个MUX选通GF(2571)的模约减结果;
A,B的乘积得到的方法为:设A,B为GF(2m)中两个m位宽的乘数,一般的,m为奇数,令2w=m+1,则A,B表示为A=Ahxw+A1,B=Bhxw+B1,Karatsuba-Ofman乘法推导为:
A×B=(Ahxw+A1)×(Bhxw+B1)
=AhBhx2w+[(Ah+A1)×(Bh+B1)+AhBh+A1B1]xw+A1B1
只需要计算三次位宽为w的小型乘法分别计算AhBh与(Ah+Al)×(Bh+Bl)AlBl,再通过一系列加法将三个小型乘法的结果进行拼接,形成上式中xw的系数,最终得到A,B的乘积;
步骤S2中,模平方、模四次方均通过纯组合逻辑电路实现;设置有限状态机用于调度模乘、模平方、模四次方运算,增加额外的寄存器用于缓存数据通路的数据,增加多路选择器MUX用于控制数据流;在每次点乘的后处理部分,基于ltoh Tsujii算法,通过有限状态机循环迭代乘法、平方、四次方运算以计算Z1与(xpZ1Z2)的模逆,进而最终计算(xQ,yQ);其中Z1、Z2分别为Montgomery Ladder的迭代两点结果,xP为点乘基点的横坐标,xQ与yQ分别为点乘结果的横坐标与纵坐标;
步骤S3中,电路上电后,复位为空闲状态,收到计算点乘的使能信号后,跳转至初始化状态执行初始化X1←xP,Z1←1,X2←xP 4+b,Z2←xp 2,初始化完成后进入主循环;主循环的两种情况,通过先默认计算ki=1,最后再判断是否交换Z1与Z2、X1与X2的值实现;计算公式
T←Z1,Z1←(X1Z2+X2Z1)2,X1←xPZ1+X1X2TZ2与T←X2,X2←X2 4+bZ2 4,Z2←T2Z2 2的各个模运算拆分至6个时钟周期依次计算;其中,xP为点乘基点的横坐标,(X1,Y1,Z1)、(X2,Y2,Z2)分别为Montgomery Ladder的迭代过程中两点结果,T为中间临时变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210713984.2A CN115062565B (zh) | 2022-06-22 | 2022-06-22 | 一种低时延椭圆曲线点乘电路设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210713984.2A CN115062565B (zh) | 2022-06-22 | 2022-06-22 | 一种低时延椭圆曲线点乘电路设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115062565A CN115062565A (zh) | 2022-09-16 |
CN115062565B true CN115062565B (zh) | 2024-01-05 |
Family
ID=83201735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210713984.2A Active CN115062565B (zh) | 2022-06-22 | 2022-06-22 | 一种低时延椭圆曲线点乘电路设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115062565B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020094256A (ko) * | 2001-06-08 | 2002-12-18 | 이광엽 | 타원곡선암호화의 유한체 역원과 승산 장치 및 그 방법 |
DE102005028662A1 (de) * | 2005-03-04 | 2006-09-07 | IHP GmbH - Innovations for High Performance Microelectronics/Institut für innovative Mikroelektronik | Verfahren und Vorrichtung zum Berechnen einer Polynom-Multiplikation, insbesondere für die elliptische Kurven-Kryptographie |
CN102306091A (zh) * | 2011-07-08 | 2012-01-04 | 西安电子科技大学 | 椭圆曲线点乘硬件快速实现方法 |
CN104184578A (zh) * | 2014-07-30 | 2014-12-03 | 山东大学 | 一种基于fpga的椭圆曲线标量乘法加速电路及其算法 |
CN108306735A (zh) * | 2017-12-29 | 2018-07-20 | 成都锐成芯微科技股份有限公司 | 椭圆曲线点乘运算的硬件实现方法及其系统 |
CN109144472A (zh) * | 2018-07-17 | 2019-01-04 | 东南大学 | 一种二元扩域椭圆曲线的标量乘法及其实现电路 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7240084B2 (en) * | 2002-05-01 | 2007-07-03 | Sun Microsystems, Inc. | Generic implementations of elliptic curve cryptography using partial reduction |
-
2022
- 2022-06-22 CN CN202210713984.2A patent/CN115062565B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20020094256A (ko) * | 2001-06-08 | 2002-12-18 | 이광엽 | 타원곡선암호화의 유한체 역원과 승산 장치 및 그 방법 |
DE102005028662A1 (de) * | 2005-03-04 | 2006-09-07 | IHP GmbH - Innovations for High Performance Microelectronics/Institut für innovative Mikroelektronik | Verfahren und Vorrichtung zum Berechnen einer Polynom-Multiplikation, insbesondere für die elliptische Kurven-Kryptographie |
CN102306091A (zh) * | 2011-07-08 | 2012-01-04 | 西安电子科技大学 | 椭圆曲线点乘硬件快速实现方法 |
CN104184578A (zh) * | 2014-07-30 | 2014-12-03 | 山东大学 | 一种基于fpga的椭圆曲线标量乘法加速电路及其算法 |
CN108306735A (zh) * | 2017-12-29 | 2018-07-20 | 成都锐成芯微科技股份有限公司 | 椭圆曲线点乘运算的硬件实现方法及其系统 |
CN109144472A (zh) * | 2018-07-17 | 2019-01-04 | 东南大学 | 一种二元扩域椭圆曲线的标量乘法及其实现电路 |
Non-Patent Citations (5)
Title |
---|
D. McGrew ; Cisco Systems ; K. Igoe ; M. Salter ; National Security Agency ; .Fundamental Elliptic Curve Cryptography Algorithms.IETF rfc6090.2011,全文. * |
M. Lochter ; BSI ; J. Merkle ; secunet Security Networks ; .Elliptic Curve Cryptography (ECC) Brainpool StandardCurves and Curve Generation.IETF rfc5639.2010,全文. * |
刘金龙 ; 张玉婷 ; 王尧 ; .GF(2~m)域ECC点乘算法优化设计.通信技术.2020,(06),全文. * |
域椭圆曲线点乘的VLSI实现方法研究;李超;张强;曲英杰;;计算机测量与控制(12);全文 * |
罗鹏 ; 许应 ; 封君 ; 王新安 ; .基于分治算法的ECC乘法器结构及实现.计算机工程.2009,(13),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN115062565A (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Satoh et al. | A scalable dual-field elliptic curve cryptographic processor | |
Ansari et al. | High-performance architecture of elliptic curve scalar multiplication | |
CN110351087B (zh) | 流水线型的蒙哥马利模乘运算方法 | |
Li et al. | High-Performance Pipelined Architecture of Elliptic Curve Scalar Multiplication Over GF (${2}^{m} $) | |
Chung et al. | A high-performance elliptic curve cryptographic processor over GF (p) with SPA resistance | |
Liu et al. | High performance FPGA implementation of elliptic curve cryptography over binary fields | |
Tian et al. | High-speed FPGA implementation of SIKE based on an ultra-low-latency modular multiplier | |
CN104184578A (zh) | 一种基于fpga的椭圆曲线标量乘法加速电路及其算法 | |
JP3213628B2 (ja) | Mを法として長い整数を乗算するための算術ユニット及びそのような乗算デバイスを具えるr.s.a.変換器 | |
CN101630244B (zh) | 一种流水线型椭圆曲线双标量乘法系统及方法 | |
Niasar et al. | Optimized architectures for elliptic curve cryptography over Curve448 | |
CN109144472B (zh) | 一种二元扩域椭圆曲线的标量乘法及其实现电路 | |
Zeghid et al. | Speed/area-efficient ECC processor implementation over GF (2 m) on FPGA via novel algorithm-architecture co-design | |
CN103780381B (zh) | 基于高基脉动阵列的蒙哥马利算法的实现装置及方法 | |
Abdulrahman et al. | High-speed hybrid-double multiplication architectures using new serial-out bit-level mastrovito multipliers | |
Reyhani-Masoleh et al. | New multiplicative inverse architectures using Gaussian normal basis | |
CN115062565B (zh) | 一种低时延椭圆曲线点乘电路设计方法 | |
CN107092462B (zh) | 一种基于fpga的64位异步乘法器 | |
CN112099763A (zh) | 用于sm2的快速安全硬件乘法器及其应用 | |
CN109284085B (zh) | 一种基于fpga的高速模乘和模幂运算方法及装置 | |
Bos et al. | Topics in computational number theory inspired by Peter L. Montgomery | |
JP3659320B2 (ja) | 乗算モジュール、乗法逆元演算回路、乗法逆元演算制御方式、該乗法逆元演算を用いる装置、暗号装置、誤り訂正復号器 | |
CN103944714A (zh) | 基于ecc实现抗dpa攻击的标量乘算法的方法及装置 | |
Vangapandu et al. | FPGA Implementation of High-Performance Montgomery Modular Multiplication with Adaptive Hold Logic | |
CN114594925A (zh) | 适用于sm2加密运算的高效模乘电路及其运算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |