ref: ceab95a64994111974f983d23125916f3d8afc9d
parent: 3897556a13df0ceea0be6bc58820ee2610224238
author: Alexei Podtelezhnikov <[email protected]>
date: Fri Jul 4 19:01:32 EDT 2014
[base] Small optimization of the ancient code. * src/base/ftcalc.c (FT_MulDiv, FT_MulDiv_No_Round): Loosen up the condition for direct 32-bit calculations.
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,10 @@
+2014-07-04 Alexei Podtelezhnikov <[email protected]>
+
+ [base] Small optimization of the ancient code.
+
+ * src/base/ftcalc.c (FT_MulDiv, FT_MulDiv_No_Round): Loosen up the
+ condition for direct 32-bit calculations.
+
2014-06-27 Werner Lemberg <[email protected]>
Fix Apple standard glyph names.
--- a/src/base/ftcalc.c
+++ b/src/base/ftcalc.c
@@ -363,16 +363,18 @@
/* */
/* we compute 'a*b+c/2', then divide it by 'c'. (positive values) */
/* */
- /* 46340 is FLOOR(SQRT(2^31-1)). */
+ /* A sufficient condition to avoid overflow is as follows. */
/* */
- /* if ( a <= 46340 && b <= 46340 ) then ( a*b <= 0x7FFEA810 ) */
+ /* a + b <= 2 * sqrt( X - c/2 ) */
/* */
- /* 0x7FFFFFFF - 0x7FFEA810 = 0x157F0 */
+ /* where X = 2^31 - 1. After Taylor expansion, we make it stronger */
/* */
- /* if ( c < 0x157F0*2 ) then ( a*b+c/2 <= 0x7FFFFFFF ) */
+ /* a + b <= 92681.9 - c / 92681.9 */
/* */
- /* and 2*0x157F0 = 176096 */
+ /* with explicit 2*sqrt(X) = 92681.9. What we actually use is this */
/* */
+ /* a + b <= 92681 - (c >> 16) */
+ /* */
FT_EXPORT_DEF( FT_Long )
FT_MulDiv( FT_Long a,
@@ -390,7 +392,7 @@
s ^= b; b = FT_ABS( b );
s ^= c; c = FT_ABS( c );
- if ( a <= 46340L && b <= 46340L && c <= 176095L && c > 0 )
+ if ( (FT_ULong)a + (FT_ULong)b <= 92681UL - ( c >> 16 ) && c > 0 )
a = ( a * b + ( c >> 1 ) ) / c;
else if ( (FT_Int32)c > 0 )
@@ -427,7 +429,7 @@
s ^= b; b = FT_ABS( b );
s ^= c; c = FT_ABS( c );
- if ( a <= 46340L && b <= 46340L && c > 0 )
+ if ( (FT_ULong)a + (FT_ULong)b <= 92681UL && c > 0 )
a = a * b / c;
else if ( (FT_Int32)c > 0 )