mpicoder.c 9.02 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
/* mpicoder.c  -  Coder for the external representation of MPIs
 * Copyright (C) 1998, 1999 Free Software Foundation, Inc.
 *
 * This file is part of GnuPG.
 *
 * GnuPG is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * GnuPG is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA
 */

21
#include <linux/bitops.h>
22
#include <linux/count_zeros.h>
23
#include <linux/byteorder/generic.h>
Herbert Xu's avatar
Herbert Xu committed
24
#include <linux/scatterlist.h>
25
#include <linux/string.h>
26 27 28 29
#include "mpi-internal.h"

#define MAX_EXTERN_MPI_BITS 16384

30 31 32 33 34 35 36 37 38 39 40 41 42
/**
 * mpi_read_raw_data - Read a raw byte stream as a positive integer
 * @xbuffer: The data to read
 * @nbytes: The amount of data to read
 */
MPI mpi_read_raw_data(const void *xbuffer, size_t nbytes)
{
	const uint8_t *buffer = xbuffer;
	int i, j;
	unsigned nbits, nlimbs;
	mpi_limb_t a;
	MPI val = NULL;

43
	while (nbytes > 0 && buffer[0] == 0) {
44 45 46 47 48 49 50 51 52 53
		buffer++;
		nbytes--;
	}

	nbits = nbytes * 8;
	if (nbits > MAX_EXTERN_MPI_BITS) {
		pr_info("MPI: mpi too large (%u bits)\n", nbits);
		return NULL;
	}
	if (nbytes > 0)
54
		nbits -= count_leading_zeros(buffer[0]) - (BITS_PER_LONG - 8);
55

56
	nlimbs = DIV_ROUND_UP(nbytes, BYTES_PER_MPI_LIMB);
57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80
	val = mpi_alloc(nlimbs);
	if (!val)
		return NULL;
	val->nbits = nbits;
	val->sign = 0;
	val->nlimbs = nlimbs;

	if (nbytes > 0) {
		i = BYTES_PER_MPI_LIMB - nbytes % BYTES_PER_MPI_LIMB;
		i %= BYTES_PER_MPI_LIMB;
		for (j = nlimbs; j > 0; j--) {
			a = 0;
			for (; i < BYTES_PER_MPI_LIMB; i++) {
				a <<= 8;
				a |= *buffer++;
			}
			i = 0;
			val->d[j - 1] = a;
		}
	}
	return val;
}
EXPORT_SYMBOL_GPL(mpi_read_raw_data);

81 82 83
MPI mpi_read_from_buffer(const void *xbuffer, unsigned *ret_nread)
{
	const uint8_t *buffer = xbuffer;
84 85
	unsigned int nbits, nbytes;
	MPI val;
86 87

	if (*ret_nread < 2)
88
		return ERR_PTR(-EINVAL);
89 90 91 92
	nbits = buffer[0] << 8 | buffer[1];

	if (nbits > MAX_EXTERN_MPI_BITS) {
		pr_info("MPI: mpi too large (%u bits)\n", nbits);
93
		return ERR_PTR(-EINVAL);
94 95
	}

96
	nbytes = DIV_ROUND_UP(nbits, 8);
97
	if (nbytes + 2 > *ret_nread) {
98 99
		pr_info("MPI: mpi larger than buffer nbytes=%u ret_nread=%u\n",
				nbytes, *ret_nread);
100 101 102
		return ERR_PTR(-EINVAL);
	}

103
	val = mpi_read_raw_data(buffer + 2, nbytes);
104
	if (!val)
105
		return ERR_PTR(-ENOMEM);
106

107
	*ret_nread = nbytes + 2;
108 109 110 111
	return val;
}
EXPORT_SYMBOL_GPL(mpi_read_from_buffer);

Michal Marek's avatar
Michal Marek committed
112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
static int count_lzeros(MPI a)
{
	mpi_limb_t alimb;
	int i, lzeros = 0;

	for (i = a->nlimbs - 1; i >= 0; i--) {
		alimb = a->d[i];
		if (alimb == 0) {
			lzeros += sizeof(mpi_limb_t);
		} else {
			lzeros += count_leading_zeros(alimb) / 8;
			break;
		}
	}
	return lzeros;
}

129 130 131 132 133 134 135
/**
 * mpi_read_buffer() - read MPI to a bufer provided by user (msb first)
 *
 * @a:		a multi precision integer
 * @buf:	bufer to which the output will be written to. Needs to be at
 *		leaset mpi_get_size(a) long.
 * @buf_len:	size of the buf.
136 137 138
 * @nbytes:	receives the actual length of the data written on success and
 *		the data to-be-written on -EOVERFLOW in case buf_len was too
 *		small.
139 140 141
 * @sign:	if not NULL, it will be set to the sign of a.
 *
 * Return:	0 on success or error code in case of error
142
 */
143 144
int mpi_read_buffer(MPI a, uint8_t *buf, unsigned buf_len, unsigned *nbytes,
		    int *sign)
145
{
146
	uint8_t *p;
147 148 149 150 151 152 153
#if BYTES_PER_MPI_LIMB == 4
	__be32 alimb;
#elif BYTES_PER_MPI_LIMB == 8
	__be64 alimb;
#else
#error please implement for this limb size.
#endif
154
	unsigned int n = mpi_get_size(a);
Michal Marek's avatar
Michal Marek committed
155
	int i, lzeros;
156

157
	if (!buf || !nbytes)
158
		return -EINVAL;
159 160 161

	if (sign)
		*sign = a->sign;
162

Michal Marek's avatar
Michal Marek committed
163
	lzeros = count_lzeros(a);
164

165 166 167 168 169
	if (buf_len < n - lzeros) {
		*nbytes = n - lzeros;
		return -EOVERFLOW;
	}

170
	p = buf;
Tadeusz Struk's avatar
Tadeusz Struk committed
171
	*nbytes = n - lzeros;
172

173 174 175
	for (i = a->nlimbs - 1 - lzeros / BYTES_PER_MPI_LIMB,
			lzeros %= BYTES_PER_MPI_LIMB;
		i >= 0; i--) {
176
#if BYTES_PER_MPI_LIMB == 4
177
		alimb = cpu_to_be32(a->d[i]);
178
#elif BYTES_PER_MPI_LIMB == 8
179
		alimb = cpu_to_be64(a->d[i]);
180 181 182
#else
#error please implement for this limb size.
#endif
183 184 185
		memcpy(p, (u8 *)&alimb + lzeros, BYTES_PER_MPI_LIMB - lzeros);
		p += BYTES_PER_MPI_LIMB - lzeros;
		lzeros = 0;
186
	}
187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204
	return 0;
}
EXPORT_SYMBOL_GPL(mpi_read_buffer);

/*
 * mpi_get_buffer() - Returns an allocated buffer with the MPI (msb first).
 * Caller must free the return string.
 * This function does return a 0 byte buffer with nbytes set to zero if the
 * value of A is zero.
 *
 * @a:		a multi precision integer.
 * @nbytes:	receives the length of this buffer.
 * @sign:	if not NULL, it will be set to the sign of the a.
 *
 * Return:	Pointer to MPI buffer or NULL on error
 */
void *mpi_get_buffer(MPI a, unsigned *nbytes, int *sign)
{
Tadeusz Struk's avatar
Tadeusz Struk committed
205
	uint8_t *buf;
206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228
	unsigned int n;
	int ret;

	if (!nbytes)
		return NULL;

	n = mpi_get_size(a);

	if (!n)
		n++;

	buf = kmalloc(n, GFP_KERNEL);

	if (!buf)
		return NULL;

	ret = mpi_read_buffer(a, buf, n, nbytes, sign);

	if (ret) {
		kfree(buf);
		return NULL;
	}
	return buf;
229 230 231
}
EXPORT_SYMBOL_GPL(mpi_get_buffer);

232 233 234 235 236 237 238 239 240
/**
 * mpi_write_to_sgl() - Funnction exports MPI to an sgl (msb first)
 *
 * This function works in the same way as the mpi_read_buffer, but it
 * takes an sgl instead of u8 * buf.
 *
 * @a:		a multi precision integer
 * @sgl:	scatterlist to write to. Needs to be at least
 *		mpi_get_size(a) long.
241 242
 * @nbytes:	the number of bytes to write.  Leading bytes will be
 *		filled with zero.
243 244 245 246
 * @sign:	if not NULL, it will be set to the sign of a.
 *
 * Return:	0 on success or error code in case of error
 */
247
int mpi_write_to_sgl(MPI a, struct scatterlist *sgl, unsigned nbytes,
248 249 250
		     int *sign)
{
	u8 *p, *p2;
251 252 253 254 255 256 257
#if BYTES_PER_MPI_LIMB == 4
	__be32 alimb;
#elif BYTES_PER_MPI_LIMB == 8
	__be64 alimb;
#else
#error please implement for this limb size.
#endif
258
	unsigned int n = mpi_get_size(a);
Herbert Xu's avatar
Herbert Xu committed
259
	struct sg_mapping_iter miter;
260
	int i, x, buf_len;
Herbert Xu's avatar
Herbert Xu committed
261
	int nents;
262 263 264 265

	if (sign)
		*sign = a->sign;

266
	if (nbytes < n)
267 268
		return -EOVERFLOW;

Herbert Xu's avatar
Herbert Xu committed
269 270 271
	nents = sg_nents_for_len(sgl, nbytes);
	if (nents < 0)
		return -EINVAL;
272

Herbert Xu's avatar
Herbert Xu committed
273 274 275 276
	sg_miter_start(&miter, sgl, nents, SG_MITER_ATOMIC | SG_MITER_TO_SG);
	sg_miter_next(&miter);
	buf_len = miter.length;
	p2 = miter.addr;
277

Herbert Xu's avatar
Herbert Xu committed
278
	while (nbytes > n) {
279 280 281 282
		i = min_t(unsigned, nbytes - n, buf_len);
		memset(p2, 0, i);
		p2 += i;
		nbytes -= i;
Herbert Xu's avatar
Herbert Xu committed
283 284 285 286 287 288 289

		buf_len -= i;
		if (!buf_len) {
			sg_miter_next(&miter);
			buf_len = miter.length;
			p2 = miter.addr;
		}
290 291 292
	}

	for (i = a->nlimbs - 1; i >= 0; i--) {
293
#if BYTES_PER_MPI_LIMB == 4
294
		alimb = a->d[i] ? cpu_to_be32(a->d[i]) : 0;
295
#elif BYTES_PER_MPI_LIMB == 8
296
		alimb = a->d[i] ? cpu_to_be64(a->d[i]) : 0;
297 298 299
#else
#error please implement for this limb size.
#endif
300
		p = (u8 *)&alimb;
301

302
		for (x = 0; x < sizeof(alimb); x++) {
303
			*p2++ = *p++;
Herbert Xu's avatar
Herbert Xu committed
304 305 306 307 308
			if (!--buf_len) {
				sg_miter_next(&miter);
				buf_len = miter.length;
				p2 = miter.addr;
			}
309 310
		}
	}
Herbert Xu's avatar
Herbert Xu committed
311 312

	sg_miter_stop(&miter);
313 314 315 316 317 318 319 320 321 322 323 324 325
	return 0;
}
EXPORT_SYMBOL_GPL(mpi_write_to_sgl);

/*
 * mpi_read_raw_from_sgl() - Function allocates an MPI and populates it with
 *			     data from the sgl
 *
 * This function works in the same way as the mpi_read_raw_data, but it
 * takes an sgl instead of void * buffer. i.e. it allocates
 * a new MPI and reads the content of the sgl to the MPI.
 *
 * @sgl:	scatterlist to read from
326
 * @nbytes:	number of bytes to read
327 328 329
 *
 * Return:	Pointer to a new MPI or NULL on error
 */
330
MPI mpi_read_raw_from_sgl(struct scatterlist *sgl, unsigned int nbytes)
331
{
Herbert Xu's avatar
Herbert Xu committed
332
	struct sg_mapping_iter miter;
333
	unsigned int nbits, nlimbs;
Herbert Xu's avatar
Herbert Xu committed
334 335 336
	int x, j, z, lzeros, ents;
	unsigned int len;
	const u8 *buff;
337 338 339
	mpi_limb_t a;
	MPI val = NULL;

Herbert Xu's avatar
Herbert Xu committed
340 341 342
	ents = sg_nents_for_len(sgl, nbytes);
	if (ents < 0)
		return NULL;
343

Herbert Xu's avatar
Herbert Xu committed
344
	sg_miter_start(&miter, sgl, ents, SG_MITER_ATOMIC | SG_MITER_FROM_SG);
345

Herbert Xu's avatar
Herbert Xu committed
346 347 348
	lzeros = 0;
	len = 0;
	while (nbytes > 0) {
349
		while (len && !*buff) {
350
			lzeros++;
351 352 353
			len--;
			buff++;
		}
354 355 356 357

		if (len && *buff)
			break;

Herbert Xu's avatar
Herbert Xu committed
358 359 360 361
		sg_miter_next(&miter);
		buff = miter.addr;
		len = miter.length;

362
		nbytes -= lzeros;
363 364 365
		lzeros = 0;
	}

Herbert Xu's avatar
Herbert Xu committed
366 367
	miter.consumed = lzeros;

368
	nbytes -= lzeros;
369 370
	nbits = nbytes * 8;
	if (nbits > MAX_EXTERN_MPI_BITS) {
371
		sg_miter_stop(&miter);
372 373 374 375 376
		pr_info("MPI: mpi too large (%u bits)\n", nbits);
		return NULL;
	}

	if (nbytes > 0)
Herbert Xu's avatar
Herbert Xu committed
377
		nbits -= count_leading_zeros(*buff) - (BITS_PER_LONG - 8);
378

379 380
	sg_miter_stop(&miter);

381 382 383 384 385 386 387 388 389 390 391 392 393 394
	nlimbs = DIV_ROUND_UP(nbytes, BYTES_PER_MPI_LIMB);
	val = mpi_alloc(nlimbs);
	if (!val)
		return NULL;

	val->nbits = nbits;
	val->sign = 0;
	val->nlimbs = nlimbs;

	if (nbytes == 0)
		return val;

	j = nlimbs - 1;
	a = 0;
395 396
	z = BYTES_PER_MPI_LIMB - nbytes % BYTES_PER_MPI_LIMB;
	z %= BYTES_PER_MPI_LIMB;
397

Herbert Xu's avatar
Herbert Xu committed
398 399 400 401
	while (sg_miter_next(&miter)) {
		buff = miter.addr;
		len = miter.length;

402
		for (x = 0; x < len; x++) {
403
			a <<= 8;
Herbert Xu's avatar
Herbert Xu committed
404
			a |= *buff++;
405 406 407 408 409 410 411
			if (((z + x + 1) % BYTES_PER_MPI_LIMB) == 0) {
				val->d[j--] = a;
				a = 0;
			}
		}
		z += x;
	}
Herbert Xu's avatar
Herbert Xu committed
412

413 414 415
	return val;
}
EXPORT_SYMBOL_GPL(mpi_read_raw_from_sgl);