Developer(s) | Cryptography LAB in Seoul National University |
---|---|
Initial release | May 15, 2016 |
Written in | C++ |
Type | Homomorphic encryption |
License | CC BY-NC 3.0 |
HEAAN (Homomorphic Encryption for Arithmetic of Approximate Numbers) is an open source homomorphic encryption (HE) library which implements an approximate HE scheme proposed by Cheon, Kim, Kim and Song (CKKS).[1] The first version of HEAAN was published on GitHub[2] on 15 May 2016, and later a new version of HEAAN with a bootstrapping algorithm[3] was released. Currently, the latest version is Version 2.1.[4][verification needed]
Unlike other HE schemes, the CKKS scheme supports approximate arithmetics over complex numbers (hence, real numbers). More precisely, the plaintext space of the CKKS scheme is [math]\displaystyle{ \mathbb{C}^{n/2} }[/math] for some power-of-two integer [math]\displaystyle{ n }[/math]. To deal with the complex plaintext vector efficiently, Cheon et al. proposed plaintext encoding/decoding methods which exploits a ring isomorphism [math]\displaystyle{ \phi: \mathbb{R}[X]/(X^n+1) \rightarrow \mathbb{C}^{n/2} }[/math].
Given a plaintext vector [math]\displaystyle{ \vec z = (z_1,z_2,...,z_{n/2}) \in \mathbb{C}^{n/2} }[/math] and a scaling factor [math]\displaystyle{ \Delta \gt 1 }[/math], the plaintext vector is encoded as a polynomial [math]\displaystyle{ m(X) \in R:= \mathbb{Z}[X]/(X^n+1) }[/math] by computing [math]\displaystyle{ m(X) = \lfloor \Delta \cdot \phi^{-1}(\vec z) \rceil \in R }[/math] where [math]\displaystyle{ \lfloor \cdot \rceil }[/math] denotes the coefficient-wise rounding function.
Given a message polynomial [math]\displaystyle{ m(X) \in R }[/math] and a scaling factor [math]\displaystyle{ \Delta \gt 1 }[/math], the message polynomial is decoded to a complex vector [math]\displaystyle{ \vec z \in \mathbb{C}^{n/2} }[/math] by computing [math]\displaystyle{ \vec z = \Delta^{-1}\cdot \phi(m(X)) \in \mathbb{C}^{n/2} }[/math].
Here the scaling factor [math]\displaystyle{ \Delta \gt 1 }[/math] enables us to control the encoding/decoding error which is occurred by the rounding process. Namely, one can obtain the approximate equation [math]\displaystyle{ \text{Dcd}(\text{Ecd}(\vec z; \Delta); \Delta) \approx \vec z }[/math] by controlling [math]\displaystyle{ \Delta }[/math] where [math]\displaystyle{ \text{Ecd} }[/math] and [math]\displaystyle{ \text{Dcd} }[/math] denote the encoding and decoding algorithm, respectively.
From the ring-isomorphic property of the mapping [math]\displaystyle{ \phi: \mathbb{R}[X]/(X^n+1) \rightarrow \mathbb{C}^{n/2} }[/math], for [math]\displaystyle{ m_1 = \text{Ecd}(\vec z_1;\Delta) }[/math] and [math]\displaystyle{ m_2 = \text{Ecd}(\vec z_2;\Delta) }[/math], the following hold:
where [math]\displaystyle{ \circ }[/math] denotes the Hadamard product of the same-length vectors. These properties guarantee the approximate correctness of the computations in the encoded state when the scaling factor [math]\displaystyle{ \Delta }[/math] is chosen appropriately.
The CKKS scheme basically consists of those algorithms: key Generation, encryption, decryption, homomorphic addition and multiplication, and rescaling. For a positive integer [math]\displaystyle{ q }[/math], let [math]\displaystyle{ R_q := R/qR }[/math] be the quotient ring of [math]\displaystyle{ R }[/math] modulo [math]\displaystyle{ q }[/math]. Let [math]\displaystyle{ \chi_s }[/math], [math]\displaystyle{ \chi_r }[/math] and [math]\displaystyle{ \chi_e }[/math] be distributions over [math]\displaystyle{ R }[/math] which output polynomials with small coefficients. These distributions, the initial modulus [math]\displaystyle{ Q }[/math], and the ring dimension [math]\displaystyle{ n }[/math] are predetermined before the key generation phase.
The key generation algorithm is following:
The encryption algorithm is following:
The decryption algorithm is following:
The decryption outputs an approximate value of the original message, i.e., [math]\displaystyle{ \text{Dec}(sk, \text{Enc}(pk, m)) \approx m }[/math], and the approximation error is determined by the choice of distributions [math]\displaystyle{ \chi_s, \chi_e, \chi_r }[/math]. When considering homomorphic operations, the evaluation errors are also included in the approximation error. Basic homomorphic operations, addition and multiplication, are done as follows.
The homomorphic addition algorithm is following:
The correctness holds as [math]\displaystyle{ \text{Dec}(sk, ct_\text{add}) \approx \text{Dec}(sk, ct) + \text{Dec}(sk, ct') }[/math].
The homomorphic multiplication algorithm is following:
The correctness holds as [math]\displaystyle{ \text{Dec}(sk, ct_\text{mult}) \approx \text{Dec}(sk, ct) \cdot \text{Dec}(sk, ct') }[/math].
Note that the approximation error (on the message) exponentially grows up on the number of homomorphic multiplications. To overcome this problem, most of HE schemes usually use a modulus-switching technique which was introduced by Brakerski, Gentry and Vaikuntanathan.[5] In case of HEAAN, the modulus-switching procedure is called rescaling. The Rescaling algorithm is very simple compared to Brakerski-Gentry-Vaikuntanathan's original algorithm. Applying the rescaling algorithm after a homomomorphic multiplication, the approximation error grows linearly, not exponentially.
The rescaling algorithm is following:
The total procedure of the CKKS scheme is as following: Each plaintext vector [math]\displaystyle{ \vec z }[/math] which consists of complex (or real) numbers is firstly encoded as a polynomial [math]\displaystyle{ m(X) \in R }[/math] by the encoding method, and then encrypted as a ciphertext [math]\displaystyle{ ct \in R_q^2 }[/math]. After several homomorphic operations, the resulting ciphertext is decrypted as a polynomial [math]\displaystyle{ m'(X) \in R }[/math] and then decoded as a plaintext vector [math]\displaystyle{ \vec z' }[/math] which is the final output.
The IND-CPA security of the CKKS scheme is based on the hardness assumption of the ring learning with errors (RLWE) problem, the ring variant of very promising lattice-based hard problem Learning with errors (LWE). Currently the best known attacks for RLWE over a power-of-two cyclotomic ring are general LWE attacks such as dual attack and primal attack. The bit security of the CKKS scheme based on known attacks was estimated by Albrecht's LWE estimator.[6]
Version 1.0, 1.1 and 2.1 have been released so far. Version 1.0 is the first implementation of the CKKS scheme without bootstrapping. In the second version, the bootstrapping algorithm was attached so that users are able to address large-scale homomorphic computations. In Version 2.1, currently the latest version, the multiplication of ring elements in [math]\displaystyle{ R_q }[/math] was accelerated by utilizing fast Fourier transform (FFT)-optimized number theoretic transform (NTT) implementation.