1 /*-------------------------------------------------------------------------
2 _mullong.c - routine for multiplication of 32 bit (unsigned) long
4 Written By - Jean Louis VERN jlvern@writeme.com (1999)
5 Written By - Sandeep Dutta . sandeep.dutta@usa.net (1999)
7 This library is free software; you can redistribute it and/or modify it
8 under the terms of the GNU Library General Public License as published by the
9 Free Software Foundation; either version 2, or (at your option) any
12 This library is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 GNU Library General Public License for more details.
17 You should have received a copy of the GNU Library General Public License
18 along with this program; if not, write to the Free Software
19 Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
21 In other words, you are welcome to use, share and improve this program.
22 You are forbidden to forbid anyone else to use, share and improve
23 what you give them. Help stamp out software-hoarding!
24 -------------------------------------------------------------------------*/
26 /* Signed and unsigned multiplication are the same - as long as the output
27 has the same precision as the input.
29 To do: _mululong and _mulslong should be replaced by _mullong.
31 bernhard@bernhardheld.de
33 Assembler-functions are provided for:
35 mcs51 small stack-auto
38 #if !defined(SDCC_USE_XSTACK) && !defined(_SDCC_NO_ASM_LIB_FUNCS)
39 # if defined(SDCC_mcs51)
40 # if defined(SDCC_MODEL_SMALL)
41 # if defined(SDCC_STACK_AUTO) && !defined(SDCC_PARMS_IN_BANK1)
42 # define _MULLONG_ASM_SMALL_AUTO
44 # define _MULLONG_ASM_SMALL
46 # elif defined(SDCC_MODEL_LARGE)
47 # if !defined(SDCC_STACK_AUTO)
48 # define _MULLONG_ASM_LARGE
54 #if defined _MULLONG_ASM_SMALL || defined _MULLONG_ASM_SMALL_AUTO
57 _mullong_dummy (void) _naked
62 __mululong: ; obsolete
63 __mulslong: ; obsolete
66 .globl __mululong ; obsolete
67 .globl __mulslong ; obsolete
69 ; the result c will be stored in r4...r7
81 ; c1 a1 * b0 + a0 * b1
82 ; c2 a2 * b0 + a1 * b1 + a0 * b2
83 ; c3 a3 * b0 + a2 * b1 + a1 * b2 + a0 * b3
85 #if !defined(SDCC_STACK_AUTO) || defined(SDCC_PARMS_IN_BANK1)
86 #if defined(SDCC_PARMS_IN_BANK1)
92 #if defined(SDCC_NOOVERLAY)
99 __mululong_PARM_2: ; obsolete
100 __mulslong_PARM_2: ; obsolete
102 .globl __mullong_PARM_2
103 .globl __mululong_PARM_2 ; obsolete
104 .globl __mulslong_PARM_2 ; obsolete
108 b0 = __mullong_PARM_2
109 b1 = (__mullong_PARM_2+1)
110 b2 = (__mullong_PARM_2+2)
111 b3 = (__mullong_PARM_2+3)
116 ; parameter a comes in a, b, dph, dpl
117 mov r2,b ; save parameter a
207 #else // SDCC_STACK_AUTO
209 ; parameter a comes in a, b, dph, dpl
210 mov r2,b ; save parameter a
220 mov a,#-2-3 ; 1 return address 2 bytes, b 4 bytes
222 mov r0,a ; 1 r0 points to b0
227 mov b0,b ; we need b0 several times
228 inc r0 ; r0 points to b1
318 #endif // SDCC_STACK_AUTO
324 #elif defined _MULLONG_ASM_LARGE
327 _mullong_dummy (void) _naked
332 __mululong: ; obsolete
333 __mulslong: ; obsolete
336 .globl __mululong ; obsolete
337 .globl __mulslong ; obsolete
339 ; the result c will be stored in r4...r7
346 ; c1 a1 * b0 + a0 * b1
347 ; c2 a2 * b0 + a1 * b1 + a0 * b2
348 ; c3 a3 * b0 + a2 * b1 + a1 * b2 + a0 * b3
350 #if !defined(SDCC_PARMS_IN_BANK1)
354 __mululong_PARM_2: ; obsolete
355 __mulslong_PARM_2: ; obsolete
357 .globl __mullong_PARM_2
358 .globl __mululong_PARM_2 ; obsolete
359 .globl __mulslong_PARM_2 ; obsolete
365 ; parameter a comes in a, b, dph, dpl
366 mov r0,dpl ; save parameter a
378 #if defined(SDCC_PARMS_IN_BANK1)
381 mov dptr,#__mullong_PARM_2
390 #if defined(SDCC_PARMS_IN_BANK1)
404 #if defined(SDCC_PARMS_IN_BANK1)
422 #if defined(SDCC_PARMS_IN_BANK1)
435 #if defined(SDCC_PARMS_IN_BANK1)
449 #if defined(SDCC_PARMS_IN_BANK1)
452 mov dptr,#__mullong_PARM_2
464 #if defined(SDCC_PARMS_IN_BANK1)
474 #if defined(SDCC_PARMS_IN_BANK1)
485 #if defined(SDCC_PARMS_IN_BANK1)
496 #if defined(SDCC_PARMS_IN_BANK1)
513 #else // _MULLONG_ASM
520 struct {unsigned char b0,b1,b2,b3 ;} b;
521 struct {unsigned int lo,hi ;} i;
523 struct { unsigned char b0; unsigned int i12; unsigned char b3;} bi;
525 #if defined(SDCC_MODEL_LARGE) || defined (SDCC_ds390)
526 #define bcast(x) ((union bil xdata *)&(x))
527 #elif defined(__z80) || defined(__gbz80)
528 #define bcast(x) ((union bil *)&(x))
530 #define bcast(x) ((union bil near *)&(x))
536 ----------------------------
541 ----------------------------
549 |-------> only this side 32 x 32 -> 32
552 _mululong (unsigned long a, unsigned long b) // in future: _mullong
556 t.i.hi = bcast(a)->b.b0 * bcast(b)->b.b2; // A
557 t.i.lo = bcast(a)->b.b0 * bcast(b)->b.b0; // A
559 t.b.b3 += bcast(a)->b.b3 *
561 t.b.b3 += bcast(a)->b.b2 *
563 t.i.hi += bcast(a)->b.b2 * bcast(b)->b.b0; // E <- b lost in .lst
564 // bcast(a)->i.hi is free !
565 t.i.hi += bcast(a)->b.b1 * bcast(b)->b.b1; // D <- b lost in .lst
567 bcast(a)->bi.b3 = bcast(a)->b.b1 *
569 bcast(a)->bi.i12 = bcast(a)->b.b1 *
572 bcast(b)->bi.b3 = bcast(a)->b.b0 *
574 bcast(b)->bi.i12 = bcast(a)->b.b0 *
576 bcast(b)->bi.b0 = 0; // B
577 bcast(a)->bi.b0 = 0; // C
584 _mulslong (long a, long b) // obsolete
586 return _mululong (a, b);
589 #endif // _MULLONG_ASM