Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
awilliam
GitHub Repository: awilliam/linux-vfio
 Path: blob/master/arch/x86/math-emu/mul_Xsig.S
10817 views
1
/*---------------------------------------------------------------------------+

2
 |  mul_Xsig.S                                                               |

3
 |                                                                           |

4
 | Multiply a 12 byte fixed point number by another fixed point number.      |

5
 |                                                                           |

6
 | Copyright (C) 1992,1994,1995                                              |

7
 |                       W. Metzenthen, 22 Parker St, Ormond, Vic 3163,      |

8
 |                       Australia.  E-mail [email protected] |

9
 |                                                                           |

10
 | Call from C as:                                                           |

11
 |   void mul32_Xsig(Xsig *x, unsigned b)                                    |

12
 |                                                                           |

13
 |   void mul64_Xsig(Xsig *x, unsigned long long *b)                         |

14
 |                                                                           |

15
 |   void mul_Xsig_Xsig(Xsig *x, unsigned *b)                                |

16
 |                                                                           |

17
 | The result is neither rounded nor normalized, and the ls bit or so may    |

18
 | be wrong.                                                                 |

19
 |                                                                           |

20
 +---------------------------------------------------------------------------*/

21
	.file	"mul_Xsig.S"

22


23


24
#include "fpu_emu.h"

25


26
.text

27
ENTRY(mul32_Xsig)

28
	pushl %ebp

29
	movl %esp,%ebp

30
	subl $16,%esp

31
	pushl %esi

32


33
	movl PARAM1,%esi

34
	movl PARAM2,%ecx

35


36
	xor %eax,%eax

37
	movl %eax,-4(%ebp)

38
	movl %eax,-8(%ebp)

39


40
	movl (%esi),%eax        /* lsl of Xsig */

41
	mull %ecx		/* msl of b */

42
	movl %edx,-12(%ebp)

43


44
	movl 4(%esi),%eax	/* midl of Xsig */

45
	mull %ecx		/* msl of b */

46
	addl %eax,-12(%ebp)

47
	adcl %edx,-8(%ebp)

48
	adcl $0,-4(%ebp)

49


50
	movl 8(%esi),%eax	/* msl of Xsig */

51
	mull %ecx		/* msl of b */

52
	addl %eax,-8(%ebp)

53
	adcl %edx,-4(%ebp)

54


55
	movl -12(%ebp),%eax

56
	movl %eax,(%esi)

57
	movl -8(%ebp),%eax

58
	movl %eax,4(%esi)

59
	movl -4(%ebp),%eax

60
	movl %eax,8(%esi)

61


62
	popl %esi

63
	leave

64
	ret

65


66


67
ENTRY(mul64_Xsig)

68
	pushl %ebp

69
	movl %esp,%ebp

70
	subl $16,%esp

71
	pushl %esi

72


73
	movl PARAM1,%esi

74
	movl PARAM2,%ecx

75


76
	xor %eax,%eax

77
	movl %eax,-4(%ebp)

78
	movl %eax,-8(%ebp)

79


80
	movl (%esi),%eax        /* lsl of Xsig */

81
	mull 4(%ecx)		/* msl of b */

82
	movl %edx,-12(%ebp)

83


84
	movl 4(%esi),%eax	/* midl of Xsig */

85
	mull (%ecx)		/* lsl of b */

86
	addl %edx,-12(%ebp)

87
	adcl $0,-8(%ebp)

88
	adcl $0,-4(%ebp)

89


90
	movl 4(%esi),%eax	/* midl of Xsig */

91
	mull 4(%ecx)		/* msl of b */

92
	addl %eax,-12(%ebp)

93
	adcl %edx,-8(%ebp)

94
	adcl $0,-4(%ebp)

95


96
	movl 8(%esi),%eax	/* msl of Xsig */

97
	mull (%ecx)		/* lsl of b */

98
	addl %eax,-12(%ebp)

99
	adcl %edx,-8(%ebp)

100
	adcl $0,-4(%ebp)

101


102
	movl 8(%esi),%eax	/* msl of Xsig */

103
	mull 4(%ecx)		/* msl of b */

104
	addl %eax,-8(%ebp)

105
	adcl %edx,-4(%ebp)

106


107
	movl -12(%ebp),%eax

108
	movl %eax,(%esi)

109
	movl -8(%ebp),%eax

110
	movl %eax,4(%esi)

111
	movl -4(%ebp),%eax

112
	movl %eax,8(%esi)

113


114
	popl %esi

115
	leave

116
	ret

117


118


119


120
ENTRY(mul_Xsig_Xsig)

121
	pushl %ebp

122
	movl %esp,%ebp

123
	subl $16,%esp

124
	pushl %esi

125


126
	movl PARAM1,%esi

127
	movl PARAM2,%ecx

128


129
	xor %eax,%eax

130
	movl %eax,-4(%ebp)

131
	movl %eax,-8(%ebp)

132


133
	movl (%esi),%eax        /* lsl of Xsig */

134
	mull 8(%ecx)		/* msl of b */

135
	movl %edx,-12(%ebp)

136


137
	movl 4(%esi),%eax	/* midl of Xsig */

138
	mull 4(%ecx)		/* midl of b */

139
	addl %edx,-12(%ebp)

140
	adcl $0,-8(%ebp)

141
	adcl $0,-4(%ebp)

142


143
	movl 8(%esi),%eax	/* msl of Xsig */

144
	mull (%ecx)		/* lsl of b */

145
	addl %edx,-12(%ebp)

146
	adcl $0,-8(%ebp)

147
	adcl $0,-4(%ebp)

148


149
	movl 4(%esi),%eax	/* midl of Xsig */

150
	mull 8(%ecx)		/* msl of b */

151
	addl %eax,-12(%ebp)

152
	adcl %edx,-8(%ebp)

153
	adcl $0,-4(%ebp)

154


155
	movl 8(%esi),%eax	/* msl of Xsig */

156
	mull 4(%ecx)		/* midl of b */

157
	addl %eax,-12(%ebp)

158
	adcl %edx,-8(%ebp)

159
	adcl $0,-4(%ebp)

160


161
	movl 8(%esi),%eax	/* msl of Xsig */

162
	mull 8(%ecx)		/* msl of b */

163
	addl %eax,-8(%ebp)

164
	adcl %edx,-4(%ebp)

165


166
	movl -12(%ebp),%edx

167
	movl %edx,(%esi)

168
	movl -8(%ebp),%edx

169
	movl %edx,4(%esi)

170
	movl -4(%ebp),%edx

171
	movl %edx,8(%esi)

172


173
	popl %esi

174
	leave

175
	ret

176


177


178