OpenSSL密码库算法笔记——第5.4.9章 椭圆曲线点的加法

先来回忆下前面提到的仿射点的加法。

曲线上的两个放射点P1＝(x1,y1)和P2＝(x2,y2)做加法P3＝P1＋P2＝(x3,y3)，当P1＝P2时，称此时的加法运算为二倍点运算。如果P1≠－P2，则P3＝P1＋P2＝(x3,y3)，可以由如下公式计算出来：

对不同的射影坐标系而言，只需要把射影点转化成对应的放射点再代入以上两个式子即可得到在对应射影坐标系下的射影点加法公式。

下面给出在各个不同的射影坐标系下做点的加法所需的运算量，其中Ａ表示放射坐标，Ｐ表示标准射影坐标，Ｊ表示Jacobi射影坐标，Ｃ表示Chudnovsky射影坐标，Ｍ表示乘法，Ｓ表示平方，Ｉ表示逆运算。

二倍点		普通加法		混合系数
2Ａ→Ａ 2Ｐ→Ｐ 2Ｊ→Ｊ 2Ｃ→Ｃ	1Ｉ＋2Ｍ＋2Ｓ 7Ｍ＋3Ｓ 4Ｍ＋4Ｓ 5Ｍ＋4Ｓ	Ａ＋Ａ→Ａ Ｐ＋Ｐ→Ｐ Ｊ＋Ｊ→Ｊ Ｃ＋Ｃ→Ｃ	1Ｉ＋2Ｍ＋1Ｓ 12Ｍ＋2Ｓ 12Ｍ＋4Ｓ 11Ｍ＋3Ｓ	Ｊ＋Ａ→Ｊ Ｊ＋Ｃ→Ｊ Ｃ＋Ａ→Ｃ	8Ｍ＋3Ｓ 11Ｍ＋3Ｓ 8Ｍ＋3Ｓ

表5.5 不同的射影坐标系下做点的加法所需的运算量

从上表可以看出，二倍点运算使用Jacobi射影坐标系最快，而加法则采用Jacobi－仿射混合系数最快。

       

如此就只需要1S＋1M，减少了二倍点的运算量。这就解释了在§5.1.1提出的问题：为什么在椭圆曲线点群的定义要加入一个参数int a_is_minus3。

下面来解释点的加法函数。

───────────────────────────────────────

int EC_POINT_add(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a, const EC_POINT *b)

功能：    点的加法

输入：    group，a，b

输出：    r ← a ＋ b

返回：    1【正常】or 0【出错】

出处：    ec_lib.c

调用：     ▼   int ec_GFp_simple_add(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a, const EC_POINT *b)

───────────────────────────────────────

接着来看看二倍点运算函数。

───────────────────────────────────────

int EC_POINT_dbl(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a)

功能：     二倍点运算

输入：     group，a

输出：     r ← 2a

返回：     1【正常】or 0【出错】

出处：     ec_lib.c

调用：     ▼   int ec_GFp_simple_dbl(const EC_GROUP *group, EC_POINT *r, const EC_POINT *a)

───────────────────────────────────────