]>
Commit | Line | Data |
---|---|---|
11b7501a SB |
1 | /* Copyright 2013 Google Inc. All Rights Reserved.\r |
2 | \r | |
3 | Distributed under MIT license.\r | |
4 | See file LICENSE for detail or copy at https://opensource.org/licenses/MIT\r | |
5 | */\r | |
6 | \r | |
7 | /* Literal cost model to allow backward reference replacement to be efficient.\r | |
8 | */\r | |
9 | \r | |
10 | #include "./literal_cost.h"\r | |
11 | \r | |
dd4f667e LG |
12 | #include "../common/platform.h"\r |
13 | #include <brotli/types.h>\r | |
11b7501a | 14 | #include "./fast_log.h"\r |
11b7501a SB |
15 | #include "./utf8_util.h"\r |
16 | \r | |
17 | #if defined(__cplusplus) || defined(c_plusplus)\r | |
18 | extern "C" {\r | |
19 | #endif\r | |
20 | \r | |
21 | static size_t UTF8Position(size_t last, size_t c, size_t clamp) {\r | |
22 | if (c < 128) {\r | |
23 | return 0; /* Next one is the 'Byte 1' again. */\r | |
24 | } else if (c >= 192) { /* Next one is the 'Byte 2' of utf-8 encoding. */\r | |
25 | return BROTLI_MIN(size_t, 1, clamp);\r | |
26 | } else {\r | |
27 | /* Let's decide over the last byte if this ends the sequence. */\r | |
dd4f667e | 28 | if (last < 0xE0) {\r |
11b7501a SB |
29 | return 0; /* Completed two or three byte coding. */\r |
30 | } else { /* Next one is the 'Byte 3' of utf-8 encoding. */\r | |
31 | return BROTLI_MIN(size_t, 2, clamp);\r | |
32 | }\r | |
33 | }\r | |
34 | }\r | |
35 | \r | |
36 | static size_t DecideMultiByteStatsLevel(size_t pos, size_t len, size_t mask,\r | |
dd4f667e | 37 | const uint8_t* data) {\r |
11b7501a SB |
38 | size_t counts[3] = { 0 };\r |
39 | size_t max_utf8 = 1; /* should be 2, but 1 compresses better. */\r | |
40 | size_t last_c = 0;\r | |
11b7501a SB |
41 | size_t i;\r |
42 | for (i = 0; i < len; ++i) {\r | |
43 | size_t c = data[(pos + i) & mask];\r | |
dd4f667e | 44 | ++counts[UTF8Position(last_c, c, 2)];\r |
11b7501a SB |
45 | last_c = c;\r |
46 | }\r | |
47 | if (counts[2] < 500) {\r | |
48 | max_utf8 = 1;\r | |
49 | }\r | |
50 | if (counts[1] + counts[2] < 25) {\r | |
51 | max_utf8 = 0;\r | |
52 | }\r | |
53 | return max_utf8;\r | |
54 | }\r | |
55 | \r | |
56 | static void EstimateBitCostsForLiteralsUTF8(size_t pos, size_t len, size_t mask,\r | |
dd4f667e LG |
57 | const uint8_t* data, float* cost) {\r |
58 | /* max_utf8 is 0 (normal ASCII single byte modeling),\r | |
59 | 1 (for 2-byte UTF-8 modeling), or 2 (for 3-byte UTF-8 modeling). */\r | |
11b7501a SB |
60 | const size_t max_utf8 = DecideMultiByteStatsLevel(pos, len, mask, data);\r |
61 | size_t histogram[3][256] = { { 0 } };\r | |
62 | size_t window_half = 495;\r | |
63 | size_t in_window = BROTLI_MIN(size_t, window_half, len);\r | |
64 | size_t in_window_utf8[3] = { 0 };\r | |
65 | \r | |
11b7501a SB |
66 | size_t i;\r |
67 | { /* Bootstrap histograms. */\r | |
68 | size_t last_c = 0;\r | |
69 | size_t utf8_pos = 0;\r | |
70 | for (i = 0; i < in_window; ++i) {\r | |
71 | size_t c = data[(pos + i) & mask];\r | |
72 | ++histogram[utf8_pos][c];\r | |
73 | ++in_window_utf8[utf8_pos];\r | |
74 | utf8_pos = UTF8Position(last_c, c, max_utf8);\r | |
75 | last_c = c;\r | |
76 | }\r | |
77 | }\r | |
78 | \r | |
79 | /* Compute bit costs with sliding window. */\r | |
80 | for (i = 0; i < len; ++i) {\r | |
81 | if (i >= window_half) {\r | |
82 | /* Remove a byte in the past. */\r | |
83 | size_t c =\r | |
84 | i < window_half + 1 ? 0 : data[(pos + i - window_half - 1) & mask];\r | |
85 | size_t last_c =\r | |
86 | i < window_half + 2 ? 0 : data[(pos + i - window_half - 2) & mask];\r | |
87 | size_t utf8_pos2 = UTF8Position(last_c, c, max_utf8);\r | |
88 | --histogram[utf8_pos2][data[(pos + i - window_half) & mask]];\r | |
89 | --in_window_utf8[utf8_pos2];\r | |
90 | }\r | |
91 | if (i + window_half < len) {\r | |
92 | /* Add a byte in the future. */\r | |
93 | size_t c = data[(pos + i + window_half - 1) & mask];\r | |
94 | size_t last_c = data[(pos + i + window_half - 2) & mask];\r | |
95 | size_t utf8_pos2 = UTF8Position(last_c, c, max_utf8);\r | |
96 | ++histogram[utf8_pos2][data[(pos + i + window_half) & mask]];\r | |
97 | ++in_window_utf8[utf8_pos2];\r | |
98 | }\r | |
99 | {\r | |
100 | size_t c = i < 1 ? 0 : data[(pos + i - 1) & mask];\r | |
101 | size_t last_c = i < 2 ? 0 : data[(pos + i - 2) & mask];\r | |
102 | size_t utf8_pos = UTF8Position(last_c, c, max_utf8);\r | |
103 | size_t masked_pos = (pos + i) & mask;\r | |
104 | size_t histo = histogram[utf8_pos][data[masked_pos]];\r | |
105 | double lit_cost;\r | |
106 | if (histo == 0) {\r | |
107 | histo = 1;\r | |
108 | }\r | |
109 | lit_cost = FastLog2(in_window_utf8[utf8_pos]) - FastLog2(histo);\r | |
110 | lit_cost += 0.02905;\r | |
111 | if (lit_cost < 1.0) {\r | |
112 | lit_cost *= 0.5;\r | |
113 | lit_cost += 0.5;\r | |
114 | }\r | |
115 | /* Make the first bytes more expensive -- seems to help, not sure why.\r | |
116 | Perhaps because the entropy source is changing its properties\r | |
117 | rapidly in the beginning of the file, perhaps because the beginning\r | |
118 | of the data is a statistical "anomaly". */\r | |
119 | if (i < 2000) {\r | |
120 | lit_cost += 0.7 - ((double)(2000 - i) / 2000.0 * 0.35);\r | |
121 | }\r | |
122 | cost[i] = (float)lit_cost;\r | |
123 | }\r | |
124 | }\r | |
125 | }\r | |
126 | \r | |
127 | void BrotliEstimateBitCostsForLiterals(size_t pos, size_t len, size_t mask,\r | |
dd4f667e | 128 | const uint8_t* data, float* cost) {\r |
11b7501a SB |
129 | if (BrotliIsMostlyUTF8(data, pos, mask, len, kMinUTF8Ratio)) {\r |
130 | EstimateBitCostsForLiteralsUTF8(pos, len, mask, data, cost);\r | |
131 | return;\r | |
132 | } else {\r | |
133 | size_t histogram[256] = { 0 };\r | |
134 | size_t window_half = 2000;\r | |
135 | size_t in_window = BROTLI_MIN(size_t, window_half, len);\r | |
136 | \r | |
137 | /* Bootstrap histogram. */\r | |
138 | size_t i;\r | |
139 | for (i = 0; i < in_window; ++i) {\r | |
140 | ++histogram[data[(pos + i) & mask]];\r | |
141 | }\r | |
142 | \r | |
143 | /* Compute bit costs with sliding window. */\r | |
144 | for (i = 0; i < len; ++i) {\r | |
145 | size_t histo;\r | |
146 | if (i >= window_half) {\r | |
147 | /* Remove a byte in the past. */\r | |
148 | --histogram[data[(pos + i - window_half) & mask]];\r | |
149 | --in_window;\r | |
150 | }\r | |
151 | if (i + window_half < len) {\r | |
152 | /* Add a byte in the future. */\r | |
153 | ++histogram[data[(pos + i + window_half) & mask]];\r | |
154 | ++in_window;\r | |
155 | }\r | |
156 | histo = histogram[data[(pos + i) & mask]];\r | |
157 | if (histo == 0) {\r | |
158 | histo = 1;\r | |
159 | }\r | |
160 | {\r | |
161 | double lit_cost = FastLog2(in_window) - FastLog2(histo);\r | |
162 | lit_cost += 0.029;\r | |
163 | if (lit_cost < 1.0) {\r | |
164 | lit_cost *= 0.5;\r | |
165 | lit_cost += 0.5;\r | |
166 | }\r | |
167 | cost[i] = (float)lit_cost;\r | |
168 | }\r | |
169 | }\r | |
170 | }\r | |
171 | }\r | |
172 | \r | |
173 | #if defined(__cplusplus) || defined(c_plusplus)\r | |
174 | } /* extern "C" */\r | |
175 | #endif\r |