blob: a708be6e9a120748bf56e3eeb9af32b77eb27271 [file] [log] [blame]
Gabor Mezeif049dbf2022-07-18 23:02:33 +02001/**
Janos Follath63184682022-08-11 17:42:59 +01002 * Modular bignum functions
Gilles Peskine7aab2fb2022-09-27 13:19:13 +02003 *
4 * This module implements operations on integers modulo some fixed modulus.
Gilles Peskine7f887bd2022-09-27 13:12:30 +02005 */
6
7/*
Gabor Mezeif049dbf2022-07-18 23:02:33 +02008 * Copyright The Mbed TLS Contributors
9 * SPDX-License-Identifier: Apache-2.0
10 *
11 * Licensed under the Apache License, Version 2.0 (the "License"); you may
12 * not use this file except in compliance with the License.
13 * You may obtain a copy of the License at
14 *
15 * http://www.apache.org/licenses/LICENSE-2.0
16 *
17 * Unless required by applicable law or agreed to in writing, software
18 * distributed under the License is distributed on an "AS IS" BASIS, WITHOUT
19 * WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
20 * See the License for the specific language governing permissions and
21 * limitations under the License.
22 */
23
24#ifndef MBEDTLS_BIGNUM_MOD_H
25#define MBEDTLS_BIGNUM_MOD_H
26
27#include "common.h"
28
29#if defined(MBEDTLS_BIGNUM_C)
30#include "mbedtls/bignum.h"
31#endif
32
Janos Follath296ea662022-08-11 14:58:29 +010033/* Skip 1 as it is slightly easier to accidentally pass to functions. */
34typedef enum
35{
36 MBEDTLS_MPI_MOD_REP_INVALID = 0,
37 MBEDTLS_MPI_MOD_REP_MONTGOMERY = 2,
38 MBEDTLS_MPI_MOD_REP_OPT_RED
39} mbedtls_mpi_mod_rep_selector;
40
41/* Make mbedtls_mpi_mod_rep_selector and mbedtls_mpi_mod_ext_rep disjoint to
42 * make it easier to catch when they are accidentally swapped. */
43typedef enum
44{
45 MBEDTLS_MPI_MOD_EXT_REP_INVALID = 0,
46 MBEDTLS_MPI_MOD_EXT_REP_LE = 8,
47 MBEDTLS_MPI_MOD_EXT_REP_BE
48} mbedtls_mpi_mod_ext_rep;
Gabor Mezeif049dbf2022-07-18 23:02:33 +020049
50typedef struct
51{
Gabor Mezeif049dbf2022-07-18 23:02:33 +020052 mbedtls_mpi_uint *p;
Gabor Mezeifd65e822022-08-12 18:09:12 +020053 size_t limbs;
Gabor Mezeif049dbf2022-07-18 23:02:33 +020054} mbedtls_mpi_mod_residue;
55
Hanno Beckercd860df2022-08-18 16:23:05 +010056typedef struct {
57 mbedtls_mpi_uint const *rr; /* The residue for 2^{2*n*biL} mod N */
58 mbedtls_mpi_uint mm; /* Montgomery const for -N^{-1} mod 2^{ciL} */
59} mbedtls_mpi_mont_struct;
60
Gabor Mezei89e31462022-08-12 15:36:56 +020061typedef void *mbedtls_mpi_opt_red_struct;
Gabor Mezeif049dbf2022-07-18 23:02:33 +020062
63typedef struct {
Janos Follathed5c8d32022-08-15 11:50:22 +010064 const mbedtls_mpi_uint *p;
Gabor Mezeifd65e822022-08-12 18:09:12 +020065 size_t limbs; // number of limbs
66 size_t bits; // bitlen of p
Janos Follath296ea662022-08-11 14:58:29 +010067 mbedtls_mpi_mod_rep_selector int_rep; // selector to signal the active member of the union
Gabor Mezeif049dbf2022-07-18 23:02:33 +020068 union rep
69 {
70 mbedtls_mpi_mont_struct mont;
71 mbedtls_mpi_opt_red_struct ored;
72 } rep;
73} mbedtls_mpi_mod_modulus;
74
Gabor Mezei37b06362022-08-02 17:22:18 +020075/** Setup a residue structure.
76 *
Janos Follathfc6fbb42022-11-25 15:43:17 +000077 * The residue will be set up with the buffer \p p and modulus \p m.
Janos Follath41427de2022-11-24 19:04:54 +000078 *
Janos Follathfc6fbb42022-11-25 15:43:17 +000079 * The memory pointed to by \p p will be used by the resulting residue structure.
80 * The value at the pointed-to memory will be the initial value of \p r and must
81 * hold a value that is less than the modulus. This value will be used as-is
Janos Follath41427de2022-11-24 19:04:54 +000082 * and interpreted according to the value of the `m->int_rep` field.
83 *
84 * The modulus \p m will be the modulus associated with \p r. The residue \p r
Janos Follath6eb92c02022-11-26 17:34:37 +000085 * should only be used in operations where the modulus is \p m.
Janos Follath41427de2022-11-24 19:04:54 +000086 *
Janos Follathee530cc2022-11-25 15:54:40 +000087 * \param[out] r The address of the residue to setup.
Janos Follath6b8a4ad2022-08-19 10:58:34 +010088 * \param[in] m The address of the modulus related to \p r.
Janos Follathfc6fbb42022-11-25 15:43:17 +000089 * \param[in] p The address of the limb array containing the value of \p r.
Janos Follath6b8a4ad2022-08-19 10:58:34 +010090 * The memory pointed to by \p p will be used by \p r and must
91 * not be modified in any way until after
Minos Galanakisaed832a2022-11-24 09:09:47 +000092 * mbedtls_mpi_mod_residue_release() is called. The data
Janos Follathfc6fbb42022-11-25 15:43:17 +000093 * pointed to by \p p must be less than the modulus (the value
94 * pointed to by `m->p`) and already in the representation
Janos Follath41427de2022-11-24 19:04:54 +000095 * indicated by `m->int_rep`.
Janos Follathee530cc2022-11-25 15:54:40 +000096 * \param p_limbs The number of limbs of \p p. Must be the same as the number
Janos Follath6eb92c02022-11-26 17:34:37 +000097 * of limbs in the modulus \p m.
Gabor Mezei37b06362022-08-02 17:22:18 +020098 *
99 * \return \c 0 if successful.
Janos Follathb7a88ec2022-08-19 12:24:40 +0100100 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p p_limbs is less than the
101 * limbs in \p m or if \p p is not less than \p m.
Gabor Mezei37b06362022-08-02 17:22:18 +0200102 */
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200103int mbedtls_mpi_mod_residue_setup( mbedtls_mpi_mod_residue *r,
Janos Follath6b8a4ad2022-08-19 10:58:34 +0100104 const mbedtls_mpi_mod_modulus *m,
Janos Follath8b718b52022-07-25 11:31:02 +0100105 mbedtls_mpi_uint *p,
Janos Follathb7a88ec2022-08-19 12:24:40 +0100106 size_t p_limbs );
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200107
Gabor Mezei37b06362022-08-02 17:22:18 +0200108/** Unbind elements of a residue structure.
109 *
Janos Follathdae11472022-08-08 11:50:02 +0100110 * This function removes the reference to the limb array that was passed to
111 * mbedtls_mpi_mod_residue_setup() to make it safe to free or use again.
112 *
113 * This function invalidates \p r and it must not be used until after
114 * mbedtls_mpi_mod_residue_setup() is called on it again.
115 *
Janos Follath6b8a4ad2022-08-19 10:58:34 +0100116 * \param[out] r The address of residue to release.
Gabor Mezei37b06362022-08-02 17:22:18 +0200117 */
118void mbedtls_mpi_mod_residue_release( mbedtls_mpi_mod_residue *r );
119
120/** Initialize a modulus structure.
121 *
Janos Follatha95f2042022-08-19 12:09:17 +0100122 * \param[out] m The address of the modulus structure to initialize.
Gabor Mezei37b06362022-08-02 17:22:18 +0200123 */
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200124void mbedtls_mpi_mod_modulus_init( mbedtls_mpi_mod_modulus *m );
125
Janos Follath63184682022-08-11 17:42:59 +0100126/** Setup a modulus structure.
Gabor Mezei37b06362022-08-02 17:22:18 +0200127 *
Janos Follath6b8a4ad2022-08-19 10:58:34 +0100128 * \param[out] m The address of the modulus structure to populate.
129 * \param[in] p The address of the limb array storing the value of \p m.
130 * The memory pointed to by \p p will be used by \p m and must
131 * not be modified in any way until after
Janos Follathdae11472022-08-08 11:50:02 +0100132 * mbedtls_mpi_mod_modulus_free() is called.
Janos Follathb7a88ec2022-08-19 12:24:40 +0100133 * \param p_limbs The number of limbs of \p p.
Janos Follathdae11472022-08-08 11:50:02 +0100134 * \param int_rep The internal representation to be used for residues
135 * associated with \p m (see #mbedtls_mpi_mod_rep_selector).
Gabor Mezei37b06362022-08-02 17:22:18 +0200136 *
137 * \return \c 0 if successful.
Janos Follathee530cc2022-11-25 15:54:40 +0000138 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p int_rep is invalid.
Gabor Mezei37b06362022-08-02 17:22:18 +0200139 */
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200140int mbedtls_mpi_mod_modulus_setup( mbedtls_mpi_mod_modulus *m,
Janos Follathed5c8d32022-08-15 11:50:22 +0100141 const mbedtls_mpi_uint *p,
Janos Follathb7a88ec2022-08-19 12:24:40 +0100142 size_t p_limbs,
Janos Follath296ea662022-08-11 14:58:29 +0100143 mbedtls_mpi_mod_rep_selector int_rep );
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200144
Janos Follathdae11472022-08-08 11:50:02 +0100145/** Free elements of a modulus structure.
146 *
147 * This function frees any memory allocated by mbedtls_mpi_mod_modulus_setup().
148 *
149 * \warning This function does not free the limb array passed to
150 * mbedtls_mpi_mod_modulus_setup() only removes the reference to it,
151 * making it safe to free or to use it again.
Gabor Mezei37b06362022-08-02 17:22:18 +0200152 *
Janos Follatha95f2042022-08-19 12:09:17 +0100153 * \param[in,out] m The address of the modulus structure to free.
Gabor Mezei37b06362022-08-02 17:22:18 +0200154 */
155void mbedtls_mpi_mod_modulus_free( mbedtls_mpi_mod_modulus *m );
156
Janos Follath5933f692022-11-02 14:35:17 +0000157/* BEGIN MERGE SLOT 1 */
158
159/* END MERGE SLOT 1 */
160
161/* BEGIN MERGE SLOT 2 */
162
163/* END MERGE SLOT 2 */
164
165/* BEGIN MERGE SLOT 3 */
Tom Cosgrove62b20482022-12-01 14:27:37 +0000166/**
167 * \brief Perform a fixed-size modular subtraction.
168 *
169 * Calculate `A - B modulo N`.
170 *
171 * \p A, \p B and \p X must all have the same number of limbs as \p N.
172 *
173 * \p X may be aliased to \p A or \p B, or even both, but may not overlap
174 * either otherwise.
175 *
176 * \note This function does not check that \p A or \p B are in canonical
177 * form (that is, are < \p N) - that will have been done by
178 * mbedtls_mpi_mod_residue_setup().
179 *
180 * \param[out] X The address of the result MPI. Must be initialized.
181 * Must have the same number of limbs as the modulus \p N.
182 * \param[in] A The address of the first MPI.
183 * \param[in] B The address of the second MPI.
184 * \param[in] N The address of the modulus. Used to perform a modulo
185 * operation on the result of the subtraction.
186 *
187 * \return \c 0 if successful.
188 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if the given MPIs do not
189 * have the correct number of limbs.
190 */
191int mbedtls_mpi_mod_sub( mbedtls_mpi_mod_residue *X,
192 const mbedtls_mpi_mod_residue *A,
193 const mbedtls_mpi_mod_residue *B,
194 const mbedtls_mpi_mod_modulus *N );
Tom Cosgrove4302d022022-12-13 10:46:39 +0000195
196/**
197 * \brief Perform modular inversion of an MPI with respect to a modulus \p N.
198 *
Tom Cosgroved692ba42022-12-14 09:53:45 +0000199 * \p A and \p X must be associated with the modulus \p N and will therefore
200 * have the same number of limbs as \p N.
201 *
Tom Cosgrove4302d022022-12-13 10:46:39 +0000202 * \p X may be aliased to \p A.
203 *
204 * \warning Currently only supports prime moduli, but does not check for them.
205 *
206 * \param[out] X The modular inverse of \p A with respect to \p N.
207 * \param[in] A The number to calculate the modular inverse of.
208 * Must not be 0.
209 * \param[in] N The modulus to use.
210 *
211 * \return \c 0 if successful.
212 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p A and \p N do not
213 * have the same number of limbs.
214 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p A is zero.
215 * \return #MBEDTLS_ERR_MPI_ALLOC_FAILED if couldn't allocate enough
216 * memory (needed for conversion to and from Mongtomery form
217 * when not in Montgomery form already, and for temporary use
218 * by the inversion calculation itself).
219 */
220
221int mbedtls_mpi_mod_inv( mbedtls_mpi_mod_residue *X,
222 const mbedtls_mpi_mod_residue *A,
223 const mbedtls_mpi_mod_modulus *N );
Janos Follath5933f692022-11-02 14:35:17 +0000224/* END MERGE SLOT 3 */
225
226/* BEGIN MERGE SLOT 4 */
227
228/* END MERGE SLOT 4 */
229
230/* BEGIN MERGE SLOT 5 */
Werner Lewise1b6b7c2022-11-29 12:25:05 +0000231/**
232 * \brief Perform a fixed-size modular addition.
233 *
234 * Calculate `A + B modulo N`.
235 *
Werner Lewiseed01aa2022-12-13 17:18:17 +0000236 * \p A, \p B and \p X must all be associated with the modulus \p N and must
237 * all have the same number of limbs as \p N.
Werner Lewise1b6b7c2022-11-29 12:25:05 +0000238 *
239 * \p X may be aliased to \p A or \p B, or even both, but may not overlap
240 * either otherwise.
241 *
242 * \note This function does not check that \p A or \p B are in canonical
243 * form (that is, are < \p N) - that will have been done by
244 * mbedtls_mpi_mod_residue_setup().
245 *
Werner Lewiseed01aa2022-12-13 17:18:17 +0000246 * \param[out] X The address of the result residue. Must be initialized.
Werner Lewise1b6b7c2022-11-29 12:25:05 +0000247 * Must have the same number of limbs as the modulus \p N.
Werner Lewiseed01aa2022-12-13 17:18:17 +0000248 * \param[in] A The address of the first input residue.
249 * \param[in] B The address of the second input residue.
Werner Lewise1b6b7c2022-11-29 12:25:05 +0000250 * \param[in] N The address of the modulus. Used to perform a modulo
251 * operation on the result of the addition.
252 *
253 * \return \c 0 if successful.
254 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if the given MPIs do not
255 * have the correct number of limbs.
256 */
257int mbedtls_mpi_mod_add( mbedtls_mpi_mod_residue *X,
258 const mbedtls_mpi_mod_residue *A,
259 const mbedtls_mpi_mod_residue *B,
260 const mbedtls_mpi_mod_modulus *N );
Janos Follath5933f692022-11-02 14:35:17 +0000261/* END MERGE SLOT 5 */
262
263/* BEGIN MERGE SLOT 6 */
264
265/* END MERGE SLOT 6 */
266
267/* BEGIN MERGE SLOT 7 */
Janos Follath41427de2022-11-24 19:04:54 +0000268/** Read a residue from a byte buffer.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000269 *
Janos Follath41427de2022-11-24 19:04:54 +0000270 * The residue will be automatically converted to the internal representation
Janos Follathfc6fbb42022-11-25 15:43:17 +0000271 * based on the value of the `m->int_rep` field.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000272 *
Janos Follath41427de2022-11-24 19:04:54 +0000273 * The modulus \p m will be the modulus associated with \p r. The residue \p r
274 * should only be used in operations where the modulus is \p m or a modulus
275 * equivalent to \p m (in the sense that all their fields or memory pointed by
276 * their fields hold the same value).
277 *
Janos Follath1f8afa22022-11-28 14:32:33 +0000278 * \param[out] r The address of the residue. It must have exactly the same
Janos Follathfc6fbb42022-11-25 15:43:17 +0000279 * number of limbs as the modulus \p m.
Janos Follath1f8afa22022-11-28 14:32:33 +0000280 * \param[in] m The address of the modulus.
281 * \param[in] buf The input buffer to import from.
Janos Follath3e3fc912022-11-24 18:02:46 +0000282 * \param buflen The length in bytes of \p buf.
283 * \param ext_rep The endianness of the number in the input buffer.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000284 *
285 * \return \c 0 if successful.
Janos Follath41427de2022-11-24 19:04:54 +0000286 * \return #MBEDTLS_ERR_MPI_BUFFER_TOO_SMALL if \p r isn't
Minos Galanakis81f4b112022-11-10 14:40:38 +0000287 * large enough to hold the value in \p buf.
Janos Follath41427de2022-11-24 19:04:54 +0000288 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p ext_rep
289 * is invalid or the value in the buffer is not less than \p m.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000290 */
291int mbedtls_mpi_mod_read( mbedtls_mpi_mod_residue *r,
Minos Galanakis8b375452022-11-24 11:04:11 +0000292 const mbedtls_mpi_mod_modulus *m,
293 const unsigned char *buf,
Janos Follath3e3fc912022-11-24 18:02:46 +0000294 size_t buflen,
295 mbedtls_mpi_mod_ext_rep ext_rep );
Janos Follath5933f692022-11-02 14:35:17 +0000296
Janos Follath41427de2022-11-24 19:04:54 +0000297/** Write a residue into a byte buffer.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000298 *
Janos Follath41427de2022-11-24 19:04:54 +0000299 * The modulus \p m must be the modulus associated with \p r (see
300 * mbedtls_mpi_mod_residue_setup() and mbedtls_mpi_mod_read()).
Minos Galanakis81f4b112022-11-10 14:40:38 +0000301 *
Janos Follath41427de2022-11-24 19:04:54 +0000302 * The residue will be automatically converted from the internal representation
303 * based on the value of `m->int_rep` field.
304 *
305 * \warning If the buffer is smaller than `m->bits`, the number of
Janos Follath6eb92c02022-11-26 17:34:37 +0000306 * leading zeroes is leaked through timing. If \p r is
Janos Follath41427de2022-11-24 19:04:54 +0000307 * secret, the caller must ensure that \p buflen is at least
308 * (`m->bits`+7)/8.
309 *
Janos Follath1f8afa22022-11-28 14:32:33 +0000310 * \param[in] r The address of the residue. It must have the same number of
311 * limbs as the modulus \p m. (\p r is an input parameter, but
312 * its value will be modified during execution and restored
313 * before the function returns.)
314 * \param[in] m The address of the modulus associated with \r.
315 * \param[out] buf The output buffer to export to.
Janos Follath3e3fc912022-11-24 18:02:46 +0000316 * \param buflen The length in bytes of \p buf.
Janos Follath41427de2022-11-24 19:04:54 +0000317 * \param ext_rep The endianness in which the number should be written into
318 * the output buffer.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000319 *
320 * \return \c 0 if successful.
321 * \return #MBEDTLS_ERR_MPI_BUFFER_TOO_SMALL if \p buf isn't
Janos Follath41427de2022-11-24 19:04:54 +0000322 * large enough to hold the value of \p r (without leading
323 * zeroes).
Janos Follathfc6fbb42022-11-25 15:43:17 +0000324 * \return #MBEDTLS_ERR_MPI_BAD_INPUT_DATA if \p ext_rep is invalid.
Janos Follath1f8afa22022-11-28 14:32:33 +0000325 * \return #MBEDTLS_ERR_MPI_ALLOC_FAILED if couldn't allocate enough
326 * memory for conversion. Can occur only for moduli with
327 * MBEDTLS_MPI_MOD_REP_MONTGOMERY.
Minos Galanakis81f4b112022-11-10 14:40:38 +0000328 */
Minos Galanakis8b375452022-11-24 11:04:11 +0000329int mbedtls_mpi_mod_write( const mbedtls_mpi_mod_residue *r,
330 const mbedtls_mpi_mod_modulus *m,
Minos Galanakis81f4b112022-11-10 14:40:38 +0000331 unsigned char *buf,
Janos Follath3e3fc912022-11-24 18:02:46 +0000332 size_t buflen,
333 mbedtls_mpi_mod_ext_rep ext_rep );
Janos Follath5933f692022-11-02 14:35:17 +0000334/* END MERGE SLOT 7 */
335
336/* BEGIN MERGE SLOT 8 */
337
338/* END MERGE SLOT 8 */
339
340/* BEGIN MERGE SLOT 9 */
341
342/* END MERGE SLOT 9 */
343
344/* BEGIN MERGE SLOT 10 */
345
346/* END MERGE SLOT 10 */
347
Gabor Mezeif049dbf2022-07-18 23:02:33 +0200348#endif /* MBEDTLS_BIGNUM_MOD_H */