1
/*___________________________________________________________________________
3
| Function Name : extract_h |
7
| Return the 16 MSB of L_var1. |
9
| Complexity weight : 1 |
14
| 32 bit long signed integer (Word32 ) whose value falls in the |
15
| range : 0x8000 0000 <= L_var1 <= 0x7fff ffff. |
24
| 16 bit short signed integer (Word16) whose value falls in the |
25
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
26
|___________________________________________________________________________|
28
PJ_INLINE(Word16) extract_h (Word32 L_var1)
30
return (Word16) (L_var1 >> 16);
32
/* ------------------------- End of extract_h() ------------------------- */
35
/*___________________________________________________________________________
37
| Function Name : extract_l |
41
| Return the 16 LSB of L_var1. |
43
| Complexity weight : 1 |
48
| 32 bit long signed integer (Word32 ) whose value falls in the |
49
| range : 0x8000 0000 <= L_var1 <= 0x7fff ffff. |
58
| 16 bit short signed integer (Word16) whose value falls in the |
59
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
60
|___________________________________________________________________________|
62
PJ_INLINE(Word16) extract_l (Word32 L_var1)
64
return (Word16) L_var1;
67
/* ------------------------- End of extract_l() ------------------------- */
70
/*___________________________________________________________________________
72
| Function Name : saturate |
76
| Limit the 32 bit input to the range of a 16 bit word. |
81
| 32 bit long signed integer (Word32) whose value falls in the |
82
| range : 0x8000 0000 <= L_var1 <= 0x7fff ffff. |
91
| 16 bit short signed integer (Word16) whose value falls in the |
92
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
93
|___________________________________________________________________________|
95
PJ_INLINE(Word16) saturate (Word32 L_var1)
97
Word16 val16 = (Word16) L_var1;
106
/* ------------------------- End of saturate() ------------------------- */
110
PJ_INLINE(Word32) L_saturate (pj_int64_t LL_var1)
112
pj_int32_t L_var1 = (pj_int32_t)LL_var1;
113
if (LL_var1 == L_var1)
115
else if (LL_var1 > MAX_32)
123
/*___________________________________________________________________________
125
| Function Name : add |
129
| Performs the addition (var1+var2) with overflow control and saturation;|
130
| the 16 bit result is set at +32767 when overflow occurs or at -32768 |
131
| when underflow occurs. |
133
| Complexity weight : 1 |
138
| 16 bit short signed integer (Word16) whose value falls in the |
139
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
142
| 16 bit short signed integer (Word16) whose value falls in the |
143
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
152
| 16 bit short signed integer (Word16) whose value falls in the |
153
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
154
|___________________________________________________________________________|
156
PJ_INLINE(Word16) add (Word16 var1, Word16 var2)
158
return saturate (var1 + var2);
160
/* ------------------------- End of add() ------------------------- */
163
/*___________________________________________________________________________
165
| Function Name : sub |
169
| Performs the subtraction (var1+var2) with overflow control and satu- |
170
| ration; the 16 bit result is set at +32767 when overflow occurs or at |
171
| -32768 when underflow occurs. |
173
| Complexity weight : 1 |
178
| 16 bit short signed integer (Word16) whose value falls in the |
179
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
182
| 16 bit short signed integer (Word16) whose value falls in the |
183
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
192
| 16 bit short signed integer (Word16) whose value falls in the |
193
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
194
|___________________________________________________________________________|
196
PJ_INLINE(Word16) sub (Word16 var1, Word16 var2)
198
return saturate ((Word32) var1 - var2);
200
/* ------------------------- End of sub() ------------------------- */
203
/*___________________________________________________________________________
205
| Function Name : negate |
209
| Negate var1 with saturation, saturate in the case where input is -32768:|
210
| negate(var1) = sub(0,var1). |
212
| Complexity weight : 1 |
217
| 16 bit short signed integer (Word16) whose value falls in the |
218
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
227
| 16 bit short signed integer (Word16) whose value falls in the |
228
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
229
|___________________________________________________________________________|
231
PJ_INLINE(Word16) negate (Word16 var1)
233
return (Word16)((var1 == MIN_16) ? MAX_16 : -var1);
235
/* ------------------------- End of negate() ------------------------- */
240
/*___________________________________________________________________________
242
| Function Name : L_add |
246
| 32 bits addition of the two 32 bits variables (L_var1+L_var2) with |
247
| overflow control and saturation; the result is set at +2147483647 when |
248
| overflow occurs or at -2147483648 when underflow occurs. |
250
| Complexity weight : 2 |
254
| L_var1 32 bit long signed integer (Word32) whose value falls in the |
255
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
257
| L_var2 32 bit long signed integer (Word32) whose value falls in the |
258
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
267
| 32 bit long signed integer (Word32) whose value falls in the |
268
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
269
|___________________________________________________________________________|
271
PJ_INLINE(Word32) L_add (Word32 L_var1, Word32 L_var2)
274
return L_saturate(((pj_int64_t)L_var1) + L_var2);
278
L_var_out = L_var1 + L_var2;
280
if (((L_var1 ^ L_var2) & MIN_32) == 0)
282
if ((L_var_out ^ L_var1) & MIN_32)
285
L_var_out = (L_var1 < 0) ? MIN_32 : MAX_32;
292
/* ------------------------- End of L_add() ------------------------- */
297
/*___________________________________________________________________________
299
| Function Name : L_sub |
303
| 32 bits subtraction of the two 32 bits variables (L_var1-L_var2) with |
304
| overflow control and saturation; the result is set at +2147483647 when |
305
| overflow occurs or at -2147483648 when underflow occurs. |
307
| Complexity weight : 2 |
311
| L_var1 32 bit long signed integer (Word32) whose value falls in the |
312
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
314
| L_var2 32 bit long signed integer (Word32) whose value falls in the |
315
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
324
| 32 bit long signed integer (Word32) whose value falls in the |
325
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
326
|___________________________________________________________________________|
328
PJ_INLINE(Word32) L_sub (Word32 L_var1, Word32 L_var2)
331
return L_saturate((pj_int64_t)L_var1 - L_var2);
335
L_var_out = L_var1 - L_var2;
337
if (((L_var1 ^ L_var2) & MIN_32) != 0)
339
if ((L_var_out ^ L_var1) & MIN_32)
341
L_var_out = (L_var1 < 0L) ? MIN_32 : MAX_32;
348
/* ------------------------- End of L_sub() ------------------------- */
350
/*___________________________________________________________________________
352
| Function Name : L_mult |
356
| L_mult is the 32 bit result of the multiplication of var1 times var2 |
357
| with one shift left i.e.: |
358
| L_mult(var1,var2) = L_shl((var1 times var2),1) and |
359
| L_mult(-32768,-32768) = 2147483647. |
361
| Complexity weight : 1 |
366
| 16 bit short signed integer (Word16) whose value falls in the |
367
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
370
| 16 bit short signed integer (Word16) whose value falls in the |
371
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
380
| 32 bit long signed integer (Word32) whose value falls in the |
381
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
382
|___________________________________________________________________________|
384
PJ_INLINE(Word32) L_mult (Word16 var1, Word16 var2)
388
L_var_out = (Word32) var1 *(Word32) var2;
390
if (L_var_out != (Word32) 0x40000000L)
392
return L_var_out << 1;
400
/* ------------------------- End of L_mult() ------------------------- */
403
/*___________________________________________________________________________
405
| Function Name : L_mac |
409
| Multiply var1 by var2 and shift the result left by 1. Add the 32 bit |
410
| result to L_var3 with saturation, return a 32 bit result: |
411
| L_mac(L_var3,var1,var2) = L_add(L_var3,L_mult(var1,var2)). |
413
| Complexity weight : 1 |
417
| L_var3 32 bit long signed integer (Word32) whose value falls in the |
418
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
421
| 16 bit short signed integer (Word16) whose value falls in the |
422
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
425
| 16 bit short signed integer (Word16) whose value falls in the |
426
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
435
| 32 bit long signed integer (Word32) whose value falls in the |
436
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
437
|___________________________________________________________________________|
439
PJ_INLINE(Word32) L_mac (Word32 L_var3, Word16 var1, Word16 var2)
441
return L_add (L_var3, L_mult (var1, var2));
443
/* ------------------------- End of L_mac() ------------------------- */
446
/*___________________________________________________________________________
448
| Function Name : round |
452
| Round the lower 16 bits of the 32 bit input number into the MS 16 bits |
453
| with saturation. Shift the resulting bits right by 16 and return the 16 |
455
| round(L_var1) = extract_h(L_add(L_var1,32768)) |
457
| Complexity weight : 1 |
462
| 32 bit long signed integer (Word32 ) whose value falls in the |
463
| range : 0x8000 0000 <= L_var1 <= 0x7fff ffff. |
472
| 16 bit short signed integer (Word16) whose value falls in the |
473
| range : 0xffff 8000 <= var_out <= 0x0000 7fff. |
474
|___________________________________________________________________________|
476
PJ_INLINE(Word16) itu_round (Word32 L_var1)
478
return extract_h (L_add (L_var1, (Word32) 0x00008000L));
480
/* ------------------------- End of round() ------------------------- */
483
/*___________________________________________________________________________
485
| Function Name : L_shr |
489
| Arithmetically shift the 32 bit input L_var1 right var2 positions with |
490
| sign extension. If var2 is negative, arithmetically shift L_var1 left |
491
| by -var2 and zero fill the -var2 LSB of the result. Saturate the result |
492
| in case of underflows or overflows. |
494
| Complexity weight : 2 |
498
| L_var1 32 bit long signed integer (Word32) whose value falls in the |
499
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
502
| 16 bit short signed integer (Word16) whose value falls in the |
503
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
512
| 32 bit long signed integer (Word32) whose value falls in the |
513
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
514
|___________________________________________________________________________|
516
PJ_INLINE(Word32) L_shr_nocheck(Word32 L_var1, Word16 var2)
519
return L_var1 >> var2;
523
return (L_var1 < 0L) ? -1 : 0;
529
return ~((~L_var1) >> var2);
533
return L_var1 >> var2;
539
PJ_INLINE(Word32) L_shl_nocheck (Word32 L_var1, Word16 var2)
542
return L_saturate( ((pj_int64_t)L_var1) << var2 );
544
for (; var2 > 0; var2--)
546
if (L_var1 > (Word32) 0X3fffffffL)
553
if (L_var1 < (Word32) 0xc0000000L)
565
PJ_INLINE(Word32) L_shr (Word32 L_var1, Word16 var2)
571
return L_shl_nocheck (L_var1, (Word16) -var2);
575
return L_shr_nocheck(L_var1, var2);
578
/* ------------------------- End of L_shr() ------------------------- */
580
/*___________________________________________________________________________
582
| Function Name : L_shl |
586
| Arithmetically shift the 32 bit input L_var1 left var2 positions. Zero |
587
| fill the var2 LSB of the result. If var2 is negative, arithmetically |
588
| shift L_var1 right by -var2 with sign extension. Saturate the result in |
589
| case of underflows or overflows. |
591
| Complexity weight : 2 |
595
| L_var1 32 bit long signed integer (Word32) whose value falls in the |
596
| range : 0x8000 0000 <= L_var3 <= 0x7fff ffff. |
599
| 16 bit short signed integer (Word16) whose value falls in the |
600
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
609
| 32 bit long signed integer (Word32) whose value falls in the |
610
| range : 0x8000 0000 <= L_var_out <= 0x7fff ffff. |
611
|___________________________________________________________________________|
613
PJ_INLINE(Word32) L_shl (Word32 L_var1, Word16 var2)
619
return L_shr_nocheck(L_var1, (Word16) -var2);
623
return L_shl_nocheck(L_var1, var2);
626
/* ------------------------- End of L_shl() ------------------------- */
629
/*___________________________________________________________________________
631
| Function Name : abs_s |
635
| Absolute value of var1; abs_s(-32768) = 32767. |
637
| Complexity weight : 1 |
642
| 16 bit short signed integer (Word16) whose value falls in the |
643
| range : 0xffff 8000 <= var1 <= 0x0000 7fff. |
652
| 16 bit short signed integer (Word16) whose value falls in the |
653
| range : 0x0000 0000 <= var_out <= 0x0000 7fff. |
654
|___________________________________________________________________________|
656
PJ_INLINE(Word16) abs_s (Word16 var1)
661
else if (var1 == MIN_16)
664
return (Word16)-var1;
674
return (Word16)-var1;
683
/* ------------------------- End of abs_s() ------------------------- */
686
PJ_INLINE(Word16) shl_nocheck(Word16 var1, Word16 var2)
689
/* blp: this should be more optimized */
690
return saturate (((Word32)var1) << var2);
692
/* Original algorithm */
693
Word32 result = (Word32) var1 *((Word32) 1 << var2);
695
if ((var2 > 15 && var1 != 0) || (result != (Word32) ((Word16) result)))
698
return (Word16) ((var1 > 0) ? MAX_16 : MIN_16);
702
return extract_l (result);
707
PJ_INLINE(Word16) shr_nocheck(Word16 var1, Word16 var2)
710
/* blp: this should yield the same value */
711
return (Word16) (var1 >> var2);
713
/* Original algorithm */
716
return (Word16)((var1 < 0) ? -1 : 0);
722
return (Word16) (~((~var1) >> var2));
726
return (Word16)(var1 >> var2);