login_cap.c revision 1.18 1 /* $NetBSD: login_cap.c,v 1.18 2004/05/09 23:28:11 lukem Exp $ */
2
3 /*-
4 * Copyright (c) 1995,1997 Berkeley Software Design, Inc. All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 * 1. Redistributions of source code must retain the above copyright
10 * notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 * notice, this list of conditions and the following disclaimer in the
13 * documentation and/or other materials provided with the distribution.
14 * 3. All advertising materials mentioning features or use of this software
15 * must display the following acknowledgement:
16 * This product includes software developed by Berkeley Software Design,
17 * Inc.
18 * 4. The name of Berkeley Software Design, Inc. may not be used to endorse
19 * or promote products derived from this software without specific prior
20 * written permission.
21 *
22 * THIS SOFTWARE IS PROVIDED BY BERKELEY SOFTWARE DESIGN, INC. ``AS IS'' AND
23 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
24 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
25 * ARE DISCLAIMED. IN NO EVENT SHALL BERKELEY SOFTWARE DESIGN, INC. BE LIABLE
26 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
27 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
28 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
29 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
30 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
31 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
32 * SUCH DAMAGE.
33 *
34 * BSDI login_cap.c,v 2.13 1998/02/07 03:17:05 prb Exp
35 */
36
37 #include <sys/cdefs.h>
38 #if defined(LIBC_SCCS) && !defined(lint)
39 __RCSID("$NetBSD: login_cap.c,v 1.18 2004/05/09 23:28:11 lukem Exp $");
40 #endif /* LIBC_SCCS and not lint */
41
42 #include <sys/types.h>
43 #include <sys/stat.h>
44 #include <sys/time.h>
45 #include <sys/resource.h>
46
47 #include <assert.h>
48 #include <ctype.h>
49 #include <err.h>
50 #include <errno.h>
51 #include <fcntl.h>
52 #include <limits.h>
53 #include <login_cap.h>
54 #include <paths.h>
55 #include <pwd.h>
56 #include <stdio.h>
57 #include <stdlib.h>
58 #include <string.h>
59 #include <syslog.h>
60 #include <unistd.h>
61 #include <util.h>
62
63 static void setuserpath(login_cap_t *, char *);
64 static u_quad_t multiply(u_quad_t, u_quad_t);
65 static u_quad_t strtolimit(char *, char **, int);
66 static u_quad_t strtosize(char *, char **, int);
67 static int gsetrl(login_cap_t *, int, char *, int type);
68 static int setuserenv(login_cap_t *);
69 static int isinfinite(const char *);
70
71 login_cap_t *
72 login_getclass(char *class)
73 {
74 const char *classfiles[2];
75 login_cap_t *lc;
76 int res;
77
78 /* class may be NULL */
79
80 if (secure_path(_PATH_LOGIN_CONF) == 0) {
81 classfiles[0] = _PATH_LOGIN_CONF;
82 classfiles[1] = NULL;
83 } else {
84 classfiles[0] = NULL;
85 }
86
87 if ((lc = malloc(sizeof(login_cap_t))) == NULL) {
88 syslog(LOG_ERR, "%s:%d malloc: %m", __FILE__, __LINE__);
89 return (0);
90 }
91
92 lc->lc_cap = 0;
93 lc->lc_style = 0;
94
95 if (class == NULL || class[0] == '\0')
96 class = LOGIN_DEFCLASS;
97
98 if ((lc->lc_class = strdup(class)) == NULL) {
99 syslog(LOG_ERR, "%s:%d strdup: %m", __FILE__, __LINE__);
100 free(lc);
101 return (0);
102 }
103
104 /*
105 * Not having a login.conf file is not an error condition.
106 * The individual routines deal reasonably with missing
107 * capabilities and use default values.
108 */
109 if (classfiles[0] == NULL)
110 return(lc);
111
112 if ((res = cgetent(&lc->lc_cap, classfiles, lc->lc_class)) != 0) {
113 lc->lc_cap = 0;
114 switch (res) {
115 case 1:
116 syslog(LOG_ERR, "%s: couldn't resolve 'tc'",
117 lc->lc_class);
118 break;
119 case -1:
120 if ((res = open(classfiles[0], 0)) >= 0)
121 close(res);
122 if (strcmp(lc->lc_class, LOGIN_DEFCLASS) == 0 &&
123 res < 0)
124 return (lc);
125 syslog(LOG_ERR, "%s: unknown class", lc->lc_class);
126 break;
127 case -2:
128 syslog(LOG_ERR, "%s: getting class information: %m",
129 lc->lc_class);
130 break;
131 case -3:
132 syslog(LOG_ERR, "%s: 'tc' reference loop",
133 lc->lc_class);
134 break;
135 default:
136 syslog(LOG_ERR, "%s: unexpected cgetent error",
137 lc->lc_class);
138 break;
139 }
140 free(lc->lc_class);
141 free(lc);
142 return (0);
143 }
144 return (lc);
145 }
146
147 login_cap_t *
148 login_getpwclass(const struct passwd *pwd)
149 {
150
151 /* pwd may be NULL */
152
153 return login_getclass(pwd ? pwd->pw_class : NULL);
154 }
155
156 char *
157 login_getcapstr(login_cap_t *lc, char *cap, char *def, char *e)
158 {
159 char *res = NULL;
160 int status;
161
162 errno = 0;
163
164 _DIAGASSERT(cap != NULL);
165
166 if (!lc || !lc->lc_cap)
167 return (def);
168
169 switch (status = cgetstr(lc->lc_cap, cap, &res)) {
170 case -1:
171 if (res)
172 free(res);
173 return (def);
174 case -2:
175 syslog(LOG_ERR, "%s: getting capability %s: %m",
176 lc->lc_class, cap);
177 if (res)
178 free(res);
179 return (e);
180 default:
181 if (status >= 0)
182 return (res);
183 syslog(LOG_ERR, "%s: unexpected error with capability %s",
184 lc->lc_class, cap);
185 if (res)
186 free(res);
187 return (e);
188 }
189 }
190
191 quad_t
192 login_getcaptime(login_cap_t *lc, char *cap, quad_t def, quad_t e)
193 {
194 char *ep;
195 char *res = NULL, *sres;
196 int status;
197 quad_t q, r;
198
199 _DIAGASSERT(cap != NULL);
200
201 errno = 0;
202 if (!lc || !lc->lc_cap)
203 return (def);
204
205 switch (status = cgetstr(lc->lc_cap, cap, &res)) {
206 case -1:
207 if (res)
208 free(res);
209 return (def);
210 case -2:
211 syslog(LOG_ERR, "%s: getting capability %s: %m",
212 lc->lc_class, cap);
213 errno = ERANGE;
214 if (res)
215 free(res);
216 return (e);
217 default:
218 if (status >= 0)
219 break;
220 syslog(LOG_ERR, "%s: unexpected error with capability %s",
221 lc->lc_class, cap);
222 errno = ERANGE;
223 if (res)
224 free(res);
225 return (e);
226 }
227
228 if (isinfinite(res))
229 return (RLIM_INFINITY);
230
231 errno = 0;
232
233 q = 0;
234 sres = res;
235 while (*res) {
236 r = strtoq(res, &ep, 0);
237 if (!ep || ep == res ||
238 ((r == QUAD_MIN || r == QUAD_MAX) && errno == ERANGE)) {
239 invalid:
240 syslog(LOG_ERR, "%s:%s=%s: invalid time",
241 lc->lc_class, cap, sres);
242 errno = ERANGE;
243 free(sres);
244 return (e);
245 }
246 switch (*ep++) {
247 case '\0':
248 --ep;
249 break;
250 case 's': case 'S':
251 break;
252 case 'm': case 'M':
253 r *= 60;
254 break;
255 case 'h': case 'H':
256 r *= 60 * 60;
257 break;
258 case 'd': case 'D':
259 r *= 60 * 60 * 24;
260 break;
261 case 'w': case 'W':
262 r *= 60 * 60 * 24 * 7;
263 break;
264 case 'y': case 'Y': /* Pretty absurd */
265 r *= 60 * 60 * 24 * 365;
266 break;
267 default:
268 goto invalid;
269 }
270 res = ep;
271 q += r;
272 }
273 free(sres);
274 return (q);
275 }
276
277 quad_t
278 login_getcapnum(login_cap_t *lc, char *cap, quad_t def, quad_t e)
279 {
280 char *ep;
281 char *res = NULL;
282 int status;
283 quad_t q;
284
285 _DIAGASSERT(cap != NULL);
286
287 errno = 0;
288 if (!lc || !lc->lc_cap)
289 return (def);
290
291 switch (status = cgetstr(lc->lc_cap, cap, &res)) {
292 case -1:
293 if (res)
294 free(res);
295 return (def);
296 case -2:
297 syslog(LOG_ERR, "%s: getting capability %s: %m",
298 lc->lc_class, cap);
299 errno = ERANGE;
300 if (res)
301 free(res);
302 return (e);
303 default:
304 if (status >= 0)
305 break;
306 syslog(LOG_ERR, "%s: unexpected error with capability %s",
307 lc->lc_class, cap);
308 errno = ERANGE;
309 if (res)
310 free(res);
311 return (e);
312 }
313
314 if (isinfinite(res))
315 return (RLIM_INFINITY);
316
317 errno = 0;
318 q = strtoq(res, &ep, 0);
319 if (!ep || ep == res || ep[0] ||
320 ((q == QUAD_MIN || q == QUAD_MAX) && errno == ERANGE)) {
321 syslog(LOG_ERR, "%s:%s=%s: invalid number",
322 lc->lc_class, cap, res);
323 errno = ERANGE;
324 free(res);
325 return (e);
326 }
327 free(res);
328 return (q);
329 }
330
331 quad_t
332 login_getcapsize(login_cap_t *lc, char *cap, quad_t def, quad_t e)
333 {
334 char *ep;
335 char *res = NULL;
336 int status;
337 quad_t q;
338
339 _DIAGASSERT(cap != NULL);
340
341 errno = 0;
342
343 if (!lc || !lc->lc_cap)
344 return (def);
345
346 switch (status = cgetstr(lc->lc_cap, cap, &res)) {
347 case -1:
348 if (res)
349 free(res);
350 return (def);
351 case -2:
352 syslog(LOG_ERR, "%s: getting capability %s: %m",
353 lc->lc_class, cap);
354 errno = ERANGE;
355 if (res)
356 free(res);
357 return (e);
358 default:
359 if (status >= 0)
360 break;
361 syslog(LOG_ERR, "%s: unexpected error with capability %s",
362 lc->lc_class, cap);
363 errno = ERANGE;
364 if (res)
365 free(res);
366 return (e);
367 }
368
369 errno = 0;
370 q = strtolimit(res, &ep, 0);
371 if (!ep || ep == res || (ep[0] && ep[1]) ||
372 ((q == QUAD_MIN || q == QUAD_MAX) && errno == ERANGE)) {
373 syslog(LOG_ERR, "%s:%s=%s: invalid size",
374 lc->lc_class, cap, res);
375 errno = ERANGE;
376 free(res);
377 return (e);
378 }
379 free(res);
380 return (q);
381 }
382
383 int
384 login_getcapbool(login_cap_t *lc, char *cap, u_int def)
385 {
386
387 _DIAGASSERT(cap != NULL);
388
389 if (!lc || !lc->lc_cap)
390 return (def);
391
392 return (cgetcap(lc->lc_cap, cap, ':') != NULL);
393 }
394
395 void
396 login_close(login_cap_t *lc)
397 {
398
399 if (lc) {
400 if (lc->lc_class)
401 free(lc->lc_class);
402 if (lc->lc_cap)
403 free(lc->lc_cap);
404 if (lc->lc_style)
405 free(lc->lc_style);
406 free(lc);
407 }
408 }
409
410 #define R_CTIME 1
411 #define R_CSIZE 2
412 #define R_CNUMB 3
413
414 static struct {
415 int what;
416 int type;
417 char * name;
418 } r_list[] = {
419 { RLIMIT_CPU, R_CTIME, "cputime", },
420 { RLIMIT_FSIZE, R_CSIZE, "filesize", },
421 { RLIMIT_DATA, R_CSIZE, "datasize", },
422 { RLIMIT_STACK, R_CSIZE, "stacksize", },
423 { RLIMIT_RSS, R_CSIZE, "memoryuse", },
424 { RLIMIT_MEMLOCK, R_CSIZE, "memorylocked", },
425 { RLIMIT_NPROC, R_CNUMB, "maxproc", },
426 { RLIMIT_NOFILE, R_CNUMB, "openfiles", },
427 { RLIMIT_CORE, R_CSIZE, "coredumpsize", },
428 { RLIMIT_SBSIZE, R_CSIZE, "sbsize", },
429 { -1, 0, 0 }
430 };
431
432 static int
433 gsetrl(login_cap_t *lc, int what, char *name, int type)
434 {
435 struct rlimit rl;
436 struct rlimit r;
437 char name_cur[32];
438 char name_max[32];
439
440 _DIAGASSERT(name != NULL);
441
442 sprintf(name_cur, "%s-cur", name);
443 sprintf(name_max, "%s-max", name);
444
445 if (getrlimit(what, &r)) {
446 syslog(LOG_ERR, "getting resource limit: %m");
447 return (-1);
448 }
449
450 #define RCUR r.rlim_cur
451 #define RMAX r.rlim_max
452
453 switch (type) {
454 case R_CTIME:
455 RCUR = login_getcaptime(lc, name, RCUR, RCUR);
456 RMAX = login_getcaptime(lc, name, RMAX, RMAX);
457 rl.rlim_cur = login_getcaptime(lc, name_cur, RCUR, RCUR);
458 rl.rlim_max = login_getcaptime(lc, name_max, RMAX, RMAX);
459 break;
460 case R_CSIZE:
461 RCUR = login_getcapsize(lc, name, RCUR, RCUR);
462 RMAX = login_getcapsize(lc, name, RMAX, RMAX);
463 rl.rlim_cur = login_getcapsize(lc, name_cur, RCUR, RCUR);
464 rl.rlim_max = login_getcapsize(lc, name_max, RMAX, RMAX);
465 break;
466 case R_CNUMB:
467 RCUR = login_getcapnum(lc, name, RCUR, RCUR);
468 RMAX = login_getcapnum(lc, name, RMAX, RMAX);
469 rl.rlim_cur = login_getcapnum(lc, name_cur, RCUR, RCUR);
470 rl.rlim_max = login_getcapnum(lc, name_max, RMAX, RMAX);
471 break;
472 default:
473 return (-1);
474 }
475
476 if (setrlimit(what, &rl)) {
477 syslog(LOG_ERR, "%s: setting resource limit %s: %m",
478 lc->lc_class, name);
479 return (-1);
480 }
481 #undef RCUR
482 #undef RMAX
483 return (0);
484 }
485
486 static int
487 setuserenv(login_cap_t *lc)
488 {
489 char *stop = ", \t";
490 int i, count;
491 char *ptr;
492 char **res;
493 char *str = login_getcapstr(lc, "setenv", NULL, NULL);
494
495 if (str == NULL || *str == '\0')
496 return 0;
497
498 /* count the sub-strings */
499 for (i = 1, ptr = str; *ptr; i++) {
500 ptr += strcspn(ptr, stop);
501 if (*ptr)
502 ptr++;
503 }
504
505 /* allocate ptr array and string */
506 count = i;
507 res = malloc(count * sizeof(char *) + strlen(str) + 1);
508
509 if (!res)
510 return -1;
511
512 ptr = (char *)(void *)res + count * sizeof(char *);
513 strcpy(ptr, str);
514
515 /* split string */
516 for (i = 0; *ptr && i < count; i++) {
517 res[i] = ptr;
518 ptr += strcspn(ptr, stop);
519 if (*ptr)
520 *ptr++ = '\0';
521 }
522
523 res[i] = NULL;
524
525 for (i = 0; i < count && res[i]; i++) {
526 if (*res[i] != '\0') {
527 if ((ptr = strchr(res[i], '=')) != NULL)
528 *ptr++ = '\0';
529 else
530 ptr = "";
531 setenv(res[i], ptr, 1);
532 }
533 }
534
535 free(res);
536 return 0;
537 }
538
539 int
540 setclasscontext(char *class, u_int flags)
541 {
542 int ret;
543 login_cap_t *lc;
544
545 flags &= LOGIN_SETRESOURCES | LOGIN_SETPRIORITY | LOGIN_SETUMASK |
546 LOGIN_SETPATH;
547
548 lc = login_getclass(class);
549 ret = lc ? setusercontext(lc, NULL, 0, flags) : -1;
550 login_close(lc);
551 return (ret);
552 }
553
554 int
555 setusercontext(login_cap_t *lc, struct passwd *pwd, uid_t uid, u_int flags)
556 {
557 login_cap_t *flc;
558 quad_t p;
559 int i;
560
561 flc = NULL;
562
563 if (!lc)
564 flc = lc = login_getclass(pwd ? pwd->pw_class : NULL);
565
566 /*
567 * Without the pwd entry being passed we cannot set either
568 * the group or the login. We could complain about it.
569 */
570 if (pwd == NULL)
571 flags &= ~(LOGIN_SETGROUP|LOGIN_SETLOGIN);
572
573 if (flags & LOGIN_SETRESOURCES)
574 for (i = 0; r_list[i].name; ++i)
575 if (gsetrl(lc, r_list[i].what, r_list[i].name,
576 r_list[i].type))
577 /* XXX - call syslog()? */;
578
579 if (flags & LOGIN_SETPRIORITY) {
580 p = login_getcapnum(lc, "priority", 0LL, 0LL);
581
582 if (setpriority(PRIO_PROCESS, 0, (int)p) < 0)
583 syslog(LOG_ERR, "%s: setpriority: %m", lc->lc_class);
584 }
585
586 if (flags & LOGIN_SETUMASK) {
587 p = login_getcapnum(lc, "umask", (quad_t) LOGIN_DEFUMASK,
588 (quad_t) LOGIN_DEFUMASK);
589 umask((mode_t)p);
590 }
591
592 if (flags & LOGIN_SETGROUP) {
593 if (setgid(pwd->pw_gid) < 0) {
594 syslog(LOG_ERR, "setgid(%d): %m", pwd->pw_gid);
595 login_close(flc);
596 return (-1);
597 }
598
599 if (initgroups(pwd->pw_name, pwd->pw_gid) < 0) {
600 syslog(LOG_ERR, "initgroups(%s,%d): %m",
601 pwd->pw_name, pwd->pw_gid);
602 login_close(flc);
603 return (-1);
604 }
605 }
606
607 if (flags & LOGIN_SETLOGIN)
608 if (setlogin(pwd->pw_name) < 0) {
609 syslog(LOG_ERR, "setlogin(%s) failure: %m",
610 pwd->pw_name);
611 login_close(flc);
612 return (-1);
613 }
614
615 if (flags & LOGIN_SETUSER)
616 if (setuid(uid) < 0) {
617 syslog(LOG_ERR, "setuid(%d): %m", uid);
618 login_close(flc);
619 return (-1);
620 }
621
622 if (flags & LOGIN_SETENV)
623 setuserenv(lc);
624
625 if (flags & LOGIN_SETPATH)
626 setuserpath(lc, pwd ? pwd->pw_dir : "");
627
628 login_close(flc);
629 return (0);
630 }
631
632 static void
633 setuserpath(login_cap_t *lc, char *home)
634 {
635 size_t hlen, plen;
636 int cnt = 0;
637 char *path;
638 char *p, *q;
639
640 _DIAGASSERT(home != NULL);
641
642 hlen = strlen(home);
643
644 p = path = login_getcapstr(lc, "path", NULL, NULL);
645 if (p) {
646 while (*p)
647 if (*p++ == '~')
648 ++cnt;
649 plen = (p - path) + cnt * (hlen + 1) + 1;
650 p = path;
651 q = path = malloc(plen);
652 if (q) {
653 while (*p) {
654 p += strspn(p, " \t");
655 if (*p == '\0')
656 break;
657 plen = strcspn(p, " \t");
658 if (hlen == 0 && *p == '~') {
659 p += plen;
660 continue;
661 }
662 if (q != path)
663 *q++ = ':';
664 if (*p == '~') {
665 strcpy(q, home);
666 q += hlen;
667 ++p;
668 --plen;
669 }
670 memcpy(q, p, plen);
671 p += plen;
672 q += plen;
673 }
674 *q = '\0';
675 } else
676 path = _PATH_DEFPATH;
677 } else
678 path = _PATH_DEFPATH;
679 if (setenv("PATH", path, 1))
680 warn("could not set PATH");
681 }
682
683 /*
684 * Convert an expression of the following forms
685 * 1) A number.
686 * 2) A number followed by a b (mult by 512).
687 * 3) A number followed by a k (mult by 1024).
688 * 5) A number followed by a m (mult by 1024 * 1024).
689 * 6) A number followed by a g (mult by 1024 * 1024 * 1024).
690 * 7) A number followed by a t (mult by 1024 * 1024 * 1024 * 1024).
691 * 8) Two or more numbers (with/without k,b,m,g, or t).
692 * separated by x (also * for backwards compatibility), specifying
693 * the product of the indicated values.
694 */
695 static u_quad_t
696 strtosize(char *str, char **endptr, int radix)
697 {
698 u_quad_t num, num2;
699 char *expr, *expr2;
700
701 _DIAGASSERT(str != NULL);
702 /* endptr may be NULL */
703
704 errno = 0;
705 num = strtouq(str, &expr, radix);
706 if (errno || expr == str) {
707 if (endptr)
708 *endptr = expr;
709 return (num);
710 }
711
712 switch(*expr) {
713 case 'b': case 'B':
714 num = multiply(num, (u_quad_t)512);
715 ++expr;
716 break;
717 case 'k': case 'K':
718 num = multiply(num, (u_quad_t)1024);
719 ++expr;
720 break;
721 case 'm': case 'M':
722 num = multiply(num, (u_quad_t)1024 * 1024);
723 ++expr;
724 break;
725 case 'g': case 'G':
726 num = multiply(num, (u_quad_t)1024 * 1024 * 1024);
727 ++expr;
728 break;
729 case 't': case 'T':
730 num = multiply(num, (u_quad_t)1024 * 1024);
731 num = multiply(num, (u_quad_t)1024 * 1024);
732 ++expr;
733 break;
734 }
735
736 if (errno)
737 goto erange;
738
739 switch(*expr) {
740 case '*': /* Backward compatible. */
741 case 'x':
742 num2 = strtosize(expr+1, &expr2, radix);
743 if (errno) {
744 expr = expr2;
745 goto erange;
746 }
747
748 if (expr2 == expr + 1) {
749 if (endptr)
750 *endptr = expr;
751 return (num);
752 }
753 expr = expr2;
754 num = multiply(num, num2);
755 if (errno)
756 goto erange;
757 break;
758 }
759 if (endptr)
760 *endptr = expr;
761 return (num);
762 erange:
763 if (endptr)
764 *endptr = expr;
765 errno = ERANGE;
766 return (UQUAD_MAX);
767 }
768
769 static u_quad_t
770 strtolimit(char *str, char **endptr, int radix)
771 {
772
773 _DIAGASSERT(str != NULL);
774 /* endptr may be NULL */
775
776 if (isinfinite(str)) {
777 if (endptr)
778 *endptr = str + strlen(str);
779 return ((u_quad_t)RLIM_INFINITY);
780 }
781 return (strtosize(str, endptr, radix));
782 }
783
784 static int
785 isinfinite(const char *s)
786 {
787 static const char *infs[] = {
788 "infinity",
789 "inf",
790 "unlimited",
791 "unlimit",
792 NULL
793 };
794 const char **i;
795
796 _DIAGASSERT(s != NULL);
797
798 for (i = infs; *i; i++) {
799 if (!strcasecmp(s, *i))
800 return 1;
801 }
802 return 0;
803 }
804
805 static u_quad_t
806 multiply(u_quad_t n1, u_quad_t n2)
807 {
808 static int bpw = 0;
809 u_quad_t m;
810 u_quad_t r;
811 int b1, b2;
812
813 /*
814 * Get rid of the simple cases
815 */
816 if (n1 == 0 || n2 == 0)
817 return (0);
818 if (n1 == 1)
819 return (n2);
820 if (n2 == 1)
821 return (n1);
822
823 /*
824 * sizeof() returns number of bytes needed for storage.
825 * This may be different from the actual number of useful bits.
826 */
827 if (!bpw) {
828 bpw = sizeof(u_quad_t) * 8;
829 while (((u_quad_t)1 << (bpw-1)) == 0)
830 --bpw;
831 }
832
833 /*
834 * First check the magnitude of each number. If the sum of the
835 * magnatude is way to high, reject the number. (If this test
836 * is not done then the first multiply below may overflow.)
837 */
838 for (b1 = bpw; (((u_quad_t)1 << (b1-1)) & n1) == 0; --b1)
839 ;
840 for (b2 = bpw; (((u_quad_t)1 << (b2-1)) & n2) == 0; --b2)
841 ;
842 if (b1 + b2 - 2 > bpw) {
843 errno = ERANGE;
844 return (UQUAD_MAX);
845 }
846
847 /*
848 * Decompose the multiplication to be:
849 * h1 = n1 & ~1
850 * h2 = n2 & ~1
851 * l1 = n1 & 1
852 * l2 = n2 & 1
853 * (h1 + l1) * (h2 + l2)
854 * (h1 * h2) + (h1 * l2) + (l1 * h2) + (l1 * l2)
855 *
856 * Since h1 && h2 do not have the low bit set, we can then say:
857 *
858 * (h1>>1 * h2>>1 * 4) + ...
859 *
860 * So if (h1>>1 * h2>>1) > (1<<(bpw - 2)) then the result will
861 * overflow.
862 *
863 * Finally, if MAX - ((h1 * l2) + (l1 * h2) + (l1 * l2)) < (h1*h2)
864 * then adding in residual amout will cause an overflow.
865 */
866
867 m = (n1 >> 1) * (n2 >> 1);
868
869 if (m >= ((u_quad_t)1 << (bpw-2))) {
870 errno = ERANGE;
871 return (UQUAD_MAX);
872 }
873
874 m *= 4;
875
876 r = (n1 & n2 & 1)
877 + (n2 & 1) * (n1 & ~(u_quad_t)1)
878 + (n1 & 1) * (n2 & ~(u_quad_t)1);
879
880 if ((u_quad_t)(m + r) < m) {
881 errno = ERANGE;
882 return (UQUAD_MAX);
883 }
884 m += r;
885
886 return (m);
887 }
888