123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107 |
- #include "sqrt_l.h"
- const char sqrt_l_id[] = "@(#)$Id $" sqrt_l_h;
- #include "typedef.h"
- #include "basic_op.h"
- #include "count.h"
- #include "sqrt_l.tab"
- Word32 sqrt_l_exp (
- Word32 L_x,
- Word16 *exp
- )
- {
-
-
-
- Word16 e, i, a, tmp;
- Word32 L_y;
- test ();
- if (L_x <= (Word32) 0)
- {
- *exp = 0; move16 ();
- return (Word32) 0;
- }
- e = norm_l_ex (L_x) & 0xFFFE; logic16 ();
- L_x = L_shl_ex (L_x, e);
- *exp = e; move16 ();
- L_x = L_shr_ex (L_x, 9);
- i = extract_h_ex (L_x);
- L_x = L_shr_ex (L_x, 1);
- a = extract_l_ex (L_x);
- a = a & (Word16) 0x7fff; logic16 ();
- i = sub_ex (i, 16);
- L_y = L_deposit_h_ex (table[i]);
- tmp = sub_ex (table[i], table[i + 1]);
- L_y = L_msu_ex (L_y, tmp, a);
-
-
- return (L_y);
- }
|