xref: /titanic_50/usr/src/lib/libc/i386/gen/lmul.s (revision b86efd96f8acd85ddaa930a2f0c1d664237e4aaf)
1/*
2 * CDDL HEADER START
3 *
4 * The contents of this file are subject to the terms of the
5 * Common Development and Distribution License, Version 1.0 only
6 * (the "License").  You may not use this file except in compliance
7 * with the License.
8 *
9 * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
10 * or http://www.opensolaris.org/os/licensing.
11 * See the License for the specific language governing permissions
12 * and limitations under the License.
13 *
14 * When distributing Covered Code, include this CDDL HEADER in each
15 * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
16 * If applicable, add the following below this CDDL HEADER, with the
17 * fields enclosed by brackets "[]" replaced with your own identifying
18 * information: Portions Copyright [yyyy] [name of copyright owner]
19 *
20 * CDDL HEADER END
21 */
22/*
23 * Copyright 2004 Sun Microsystems, Inc.  All rights reserved.
24 * Use is subject to license terms.
25 */
26
27	.ident	"%Z%%M%	%I%	%E% SMI"
28
29	.file	"%M%"
30
31/ Double long multiply routine.
32
33#include <sys/asm_linkage.h>
34
35	ANSI_PRAGMA_WEAK(lmul,function)
36
37#include "SYS.h"
38
39	.set	lop,8
40	.set	rop,16
41	.set	ans,0
42
43	ENTRY(lmul)
44	popl	%eax
45	xchgl	%eax,0(%esp)
46
47	pushl	%eax
48
49	movl	lop+4(%esp),%eax
50	mull	rop(%esp)	/ high(lop) * low(rop)
51	movl	%eax,%ecx	/ partial high(product)
52
53	movl	rop+4(%esp),%eax
54	mull	lop(%esp)	/ high(rop) * low(lop)
55	addl	%eax,%ecx	/ partial sum of high(product)
56
57	movl	rop(%esp),%eax
58	mull	lop(%esp)	/ low(rop) * low(lop)
59	addl	%edx,%ecx	/ final high(product)
60
61	movl	%eax,%edx
62	popl	%eax
63	movl	%edx,ans(%eax)
64	movl	%ecx,ans+4(%eax)
65
66	ret
67	SET_SIZE(lmul)
68