Revert "Redo range check in CPDF_SampledFunc::v_Call()."
[pdfium.git] / core / src / fpdfapi / fpdf_page / fpdf_page_func.cpp
1 // Copyright 2014 PDFium Authors. All rights reserved.
2 // Use of this source code is governed by a BSD-style license that can be
3 // found in the LICENSE file.
4
5 // Original code copyright 2014 Foxit Software Inc. http://www.foxitsoftware.com
6
7 #include <limits.h>
8
9 #include "../../../include/fpdfapi/fpdf_module.h"
10 #include "../../../include/fpdfapi/fpdf_page.h"
11 #include "../../../include/fxcrt/fx_safe_types.h"
12 #include "../../../third_party/base/numerics/safe_conversions_impl.h"
13 #include "pageint.h"
14
15 class CPDF_PSEngine;
16 typedef enum {PSOP_ADD, PSOP_SUB, PSOP_MUL, PSOP_DIV, PSOP_IDIV, PSOP_MOD,
17               PSOP_NEG, PSOP_ABS, PSOP_CEILING, PSOP_FLOOR, PSOP_ROUND, PSOP_TRUNCATE,
18               PSOP_SQRT, PSOP_SIN, PSOP_COS, PSOP_ATAN, PSOP_EXP, PSOP_LN, PSOP_LOG,
19               PSOP_CVI, PSOP_CVR, PSOP_EQ, PSOP_NE, PSOP_GT, PSOP_GE, PSOP_LT, PSOP_LE,
20               PSOP_AND, PSOP_OR, PSOP_XOR, PSOP_NOT, PSOP_BITSHIFT, PSOP_TRUE, PSOP_FALSE,
21               PSOP_IF, PSOP_IFELSE, PSOP_POP, PSOP_EXCH, PSOP_DUP, PSOP_COPY,
22               PSOP_INDEX, PSOP_ROLL, PSOP_PROC, PSOP_CONST
23              } PDF_PSOP;
24 class CPDF_PSProc
25 {
26 public:
27     ~CPDF_PSProc();
28     FX_BOOL     Parse(CPDF_SimpleParser& parser);
29     FX_BOOL     Execute(CPDF_PSEngine* pEngine);
30     CFX_PtrArray                m_Operators;
31 };
32 #define PSENGINE_STACKSIZE 100
33 class CPDF_PSEngine
34 {
35 public:
36     CPDF_PSEngine();
37     ~CPDF_PSEngine();
38     FX_BOOL     Parse(const FX_CHAR* string, int size);
39     FX_BOOL     Execute()
40     {
41         return m_MainProc.Execute(this);
42     }
43     FX_BOOL     DoOperator(PDF_PSOP op);
44     void        Reset()
45     {
46         m_StackCount = 0;
47     }
48     void        Push(FX_FLOAT value);
49     void        Push(int value)
50     {
51         Push((FX_FLOAT)value);
52     }
53     FX_FLOAT    Pop();
54     int         GetStackSize()
55     {
56         return m_StackCount;
57     }
58 private:
59     FX_FLOAT    m_Stack[PSENGINE_STACKSIZE];
60     int         m_StackCount;
61     CPDF_PSProc m_MainProc;
62 };
63 CPDF_PSProc::~CPDF_PSProc()
64 {
65     int size = m_Operators.GetSize();
66     for (int i = 0; i < size; i ++) {
67         if (m_Operators[i] == (void*)PSOP_PROC) {
68             delete (CPDF_PSProc*)m_Operators[i + 1];
69             i ++;
70         } else if (m_Operators[i] == (void*)PSOP_CONST) {
71             FX_Free((FX_FLOAT*)m_Operators[i + 1]);
72             i ++;
73         }
74     }
75 }
76 FX_BOOL CPDF_PSProc::Execute(CPDF_PSEngine* pEngine)
77 {
78     int size = m_Operators.GetSize();
79     for (int i = 0; i < size; i ++) {
80         PDF_PSOP op = (PDF_PSOP)(uintptr_t)m_Operators[i];
81         if (op == PSOP_PROC) {
82             i ++;
83         } else if (op == PSOP_CONST) {
84             pEngine->Push(*(FX_FLOAT*)m_Operators[i + 1]);
85             i ++;
86         } else if (op == PSOP_IF) {
87             if (i < 2 || m_Operators[i - 2] != (void*)PSOP_PROC) {
88                 return FALSE;
89             }
90             if ((int)pEngine->Pop()) {
91                 ((CPDF_PSProc*)m_Operators[i - 1])->Execute(pEngine);
92             }
93         } else if (op == PSOP_IFELSE) {
94             if (i < 4 || m_Operators[i - 2] != (void*)PSOP_PROC ||
95                     m_Operators[i - 4] != (void*)PSOP_PROC) {
96                 return FALSE;
97             }
98             if ((int)pEngine->Pop()) {
99                 ((CPDF_PSProc*)m_Operators[i - 3])->Execute(pEngine);
100             } else {
101                 ((CPDF_PSProc*)m_Operators[i - 1])->Execute(pEngine);
102             }
103         } else {
104             pEngine->DoOperator(op);
105         }
106     }
107     return TRUE;
108 }
109 CPDF_PSEngine::CPDF_PSEngine()
110 {
111     m_StackCount = 0;
112 }
113 CPDF_PSEngine::~CPDF_PSEngine()
114 {
115 }
116 void CPDF_PSEngine::Push(FX_FLOAT v)
117 {
118     if (m_StackCount == 100) {
119         return;
120     }
121     m_Stack[m_StackCount++] = v;
122 }
123 FX_FLOAT CPDF_PSEngine::Pop()
124 {
125     if (m_StackCount == 0) {
126         return 0;
127     }
128     return m_Stack[--m_StackCount];
129 }
130 const struct _PDF_PSOpName {
131     const FX_CHAR* name;
132     PDF_PSOP op;
133 } _PDF_PSOpNames[] = {
134     {"add", PSOP_ADD}, {"sub", PSOP_SUB}, {"mul", PSOP_MUL}, {"div", PSOP_DIV},
135     {"idiv", PSOP_IDIV}, {"mod", PSOP_MOD}, {"neg", PSOP_NEG}, {"abs", PSOP_ABS},
136     {"ceiling", PSOP_CEILING}, {"floor", PSOP_FLOOR}, {"round", PSOP_ROUND},
137     {"truncate", PSOP_TRUNCATE}, {"sqrt", PSOP_SQRT}, {"sin", PSOP_SIN},
138     {"cos", PSOP_COS}, {"atan", PSOP_ATAN}, {"exp", PSOP_EXP}, {"ln", PSOP_LN},
139     {"log", PSOP_LOG}, {"cvi", PSOP_CVI}, {"cvr", PSOP_CVR}, {"eq", PSOP_EQ},
140     {"ne", PSOP_NE}, {"gt", PSOP_GT}, {"ge", PSOP_GE}, {"lt", PSOP_LT},
141     {"le", PSOP_LE}, {"and", PSOP_AND}, {"or", PSOP_OR}, {"xor", PSOP_XOR},
142     {"not", PSOP_NOT}, {"bitshift", PSOP_BITSHIFT}, {"true", PSOP_TRUE},
143     {"false", PSOP_FALSE}, {"if", PSOP_IF}, {"ifelse", PSOP_IFELSE},
144     {"pop", PSOP_POP}, {"exch", PSOP_EXCH}, {"dup", PSOP_DUP},
145     {"copy", PSOP_COPY}, {"index", PSOP_INDEX}, {"roll", PSOP_ROLL},
146     {NULL, PSOP_PROC}
147 };
148 FX_BOOL CPDF_PSEngine::Parse(const FX_CHAR* string, int size)
149 {
150     CPDF_SimpleParser parser((uint8_t*)string, size);
151     CFX_ByteStringC word = parser.GetWord();
152     if (word != FX_BSTRC("{")) {
153         return FALSE;
154     }
155     return m_MainProc.Parse(parser);
156 }
157 FX_BOOL CPDF_PSProc::Parse(CPDF_SimpleParser& parser)
158 {
159     while (1) {
160         CFX_ByteStringC word = parser.GetWord();
161         if (word.IsEmpty()) {
162             return FALSE;
163         }
164         if (word == FX_BSTRC("}")) {
165             return TRUE;
166         }
167         if (word == FX_BSTRC("{")) {
168             CPDF_PSProc* pProc = new CPDF_PSProc;
169             m_Operators.Add((void*)PSOP_PROC);
170             m_Operators.Add(pProc);
171             if (!pProc->Parse(parser)) {
172                 return FALSE;
173             }
174         } else {
175             int i = 0;
176             while (_PDF_PSOpNames[i].name) {
177                 if (word == CFX_ByteStringC(_PDF_PSOpNames[i].name)) {
178                     m_Operators.Add((void*)_PDF_PSOpNames[i].op);
179                     break;
180                 }
181                 i ++;
182             }
183             if (_PDF_PSOpNames[i].name == NULL) {
184                 FX_FLOAT* pd = FX_Alloc(FX_FLOAT, 1);
185                 *pd = FX_atof(word);
186                 m_Operators.Add((void*)PSOP_CONST);
187                 m_Operators.Add(pd);
188             }
189         }
190     }
191 }
192 #define PI 3.1415926535897932384626433832795f
193 FX_BOOL CPDF_PSEngine::DoOperator(PDF_PSOP op)
194 {
195     int i1, i2;
196     FX_FLOAT d1, d2;
197     switch (op) {
198         case PSOP_ADD:
199             d1 = Pop();
200             d2 = Pop();
201             Push(d1 + d2);
202             break;
203         case PSOP_SUB:
204             d2 = Pop();
205             d1 = Pop();
206             Push(d1 - d2);
207             break;
208         case PSOP_MUL:
209             d1 = Pop();
210             d2 = Pop();
211             Push(d1 * d2);
212             break;
213         case PSOP_DIV:
214             d2 = Pop();
215             d1 = Pop();
216             Push(d1 / d2);
217             break;
218         case PSOP_IDIV:
219             i2 = (int)Pop();
220             i1 = (int)Pop();
221             Push(i1 / i2);
222             break;
223         case PSOP_MOD:
224             i2 = (int)Pop();
225             i1 = (int)Pop();
226             Push(i1 % i2);
227             break;
228         case PSOP_NEG:
229             d1 = Pop();
230             Push(-d1);
231             break;
232         case PSOP_ABS:
233             d1 = Pop();
234             Push((FX_FLOAT)FXSYS_fabs(d1));
235             break;
236         case PSOP_CEILING:
237             d1 = Pop();
238             Push((FX_FLOAT)FXSYS_ceil(d1));
239             break;
240         case PSOP_FLOOR:
241             d1 = Pop();
242             Push((FX_FLOAT)FXSYS_floor(d1));
243             break;
244         case PSOP_ROUND:
245             d1 = Pop();
246             Push(FXSYS_round(d1));
247             break;
248         case PSOP_TRUNCATE:
249             i1 = (int)Pop();
250             Push(i1);
251             break;
252         case PSOP_SQRT:
253             d1 = Pop();
254             Push((FX_FLOAT)FXSYS_sqrt(d1));
255             break;
256         case PSOP_SIN:
257             d1 = Pop();
258             Push((FX_FLOAT)FXSYS_sin(d1 * PI / 180.0f));
259             break;
260         case PSOP_COS:
261             d1 = Pop();
262             Push((FX_FLOAT)FXSYS_cos(d1 * PI / 180.0f));
263             break;
264         case PSOP_ATAN:
265             d2 = Pop();
266             d1 = Pop();
267             d1 = (FX_FLOAT)(FXSYS_atan2(d1, d2) * 180.0 / PI);
268             if (d1 < 0) {
269                 d1 += 360;
270             }
271             Push(d1);
272             break;
273         case PSOP_EXP:
274             d2 = Pop();
275             d1 = Pop();
276             Push((FX_FLOAT)FXSYS_pow(d1, d2));
277             break;
278         case PSOP_LN:
279             d1 = Pop();
280             Push((FX_FLOAT)FXSYS_log(d1));
281             break;
282         case PSOP_LOG:
283             d1 = Pop();
284             Push((FX_FLOAT)FXSYS_log10(d1));
285             break;
286         case PSOP_CVI:
287             i1 = (int)Pop();
288             Push(i1);
289             break;
290         case PSOP_CVR:
291             break;
292         case PSOP_EQ:
293             d2 = Pop();
294             d1 = Pop();
295             Push((int)(d1 == d2));
296             break;
297         case PSOP_NE:
298             d2 = Pop();
299             d1 = Pop();
300             Push((int)(d1 != d2));
301             break;
302         case PSOP_GT:
303             d2 = Pop();
304             d1 = Pop();
305             Push((int)(d1 > d2));
306             break;
307         case PSOP_GE:
308             d2 = Pop();
309             d1 = Pop();
310             Push((int)(d1 >= d2));
311             break;
312         case PSOP_LT:
313             d2 = Pop();
314             d1 = Pop();
315             Push((int)(d1 < d2));
316             break;
317         case PSOP_LE:
318             d2 = Pop();
319             d1 = Pop();
320             Push((int)(d1 <= d2));
321             break;
322         case PSOP_AND:
323             i1 = (int)Pop();
324             i2 = (int)Pop();
325             Push(i1 & i2);
326             break;
327         case PSOP_OR:
328             i1 = (int)Pop();
329             i2 = (int)Pop();
330             Push(i1 | i2);
331             break;
332         case PSOP_XOR:
333             i1 = (int)Pop();
334             i2 = (int)Pop();
335             Push(i1 ^ i2);
336             break;
337         case PSOP_NOT:
338             i1 = (int)Pop();
339             Push((int)!i1);
340             break;
341         case PSOP_BITSHIFT: {
342                 int shift = (int)Pop();
343                 int i = (int)Pop();
344                 if (shift > 0) {
345                     Push(i << shift);
346                 } else {
347                     Push(i >> -shift);
348                 }
349                 break;
350             }
351         case PSOP_TRUE:
352             Push(1);
353             break;
354         case PSOP_FALSE:
355             Push(0);
356             break;
357         case PSOP_POP:
358             Pop();
359             break;
360         case PSOP_EXCH:
361             d2 = Pop();
362             d1 = Pop();
363             Push(d2);
364             Push(d1);
365             break;
366         case PSOP_DUP:
367             d1 = Pop();
368             Push(d1);
369             Push(d1);
370             break;
371         case PSOP_COPY: {
372                 int n = (int)Pop();
373                 if (n < 0 || n > PSENGINE_STACKSIZE || m_StackCount + n > PSENGINE_STACKSIZE || n > m_StackCount) {
374                     break;
375                 }
376                 for (int i = 0; i < n; i ++) {
377                     m_Stack[m_StackCount + i] = m_Stack[m_StackCount + i - n];
378                 }
379                 m_StackCount += n;
380                 break;
381             }
382         case PSOP_INDEX: {
383                 int n = (int)Pop();
384                 if (n < 0 || n >= m_StackCount) {
385                     break;
386                 }
387                 Push(m_Stack[m_StackCount - n - 1]);
388                 break;
389             }
390         case PSOP_ROLL: {
391                 int j = (int)Pop();
392                 int n = (int)Pop();
393                 if (m_StackCount == 0) {
394                     break;
395                 }
396                 if (n < 0 || n > m_StackCount) {
397                     break;
398                 }
399                 if (j < 0)
400                     for (int i = 0; i < -j; i ++) {
401                         FX_FLOAT first = m_Stack[m_StackCount - n];
402                         for (int ii = 0; ii < n - 1; ii ++) {
403                             m_Stack[m_StackCount - n + ii] = m_Stack[m_StackCount - n + ii + 1];
404                         }
405                         m_Stack[m_StackCount - 1] = first;
406                     }
407                 else
408                     for (int i = 0; i < j; i ++) {
409                         FX_FLOAT last = m_Stack[m_StackCount - 1];
410                         int ii;
411                         for (ii = 0; ii < n - 1; ii ++) {
412                             m_Stack[m_StackCount - ii - 1] = m_Stack[m_StackCount - ii - 2];
413                         }
414                         m_Stack[m_StackCount - ii - 1] = last;
415                     }
416                 break;
417             }
418         default:
419             break;
420     }
421     return TRUE;
422 }
423 static FX_FLOAT PDF_Interpolate(FX_FLOAT x, FX_FLOAT xmin, FX_FLOAT xmax, FX_FLOAT ymin, FX_FLOAT ymax)
424 {
425     return ((x - xmin) * (ymax - ymin) / (xmax - xmin)) + ymin;
426 }
427 static FX_DWORD _GetBits32(const uint8_t* pData, int bitpos, int nbits)
428 {
429     int result = 0;
430     for (int i = 0; i < nbits; i ++)
431         if (pData[(bitpos + i) / 8] & (1 << (7 - (bitpos + i) % 8))) {
432             result |= 1 << (nbits - i - 1);
433         }
434     return result;
435 }
436 typedef struct {
437     FX_FLOAT    encode_max, encode_min;
438     int                 sizes;
439 } SampleEncodeInfo;
440 typedef struct {
441     FX_FLOAT    decode_max, decode_min;
442 } SampleDecodeInfo;
443 class CPDF_SampledFunc : public CPDF_Function
444 {
445 public:
446     CPDF_SampledFunc();
447     virtual ~CPDF_SampledFunc();
448     virtual FX_BOOL             v_Init(CPDF_Object* pObj);
449     virtual FX_BOOL             v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const;
450     SampleEncodeInfo*   m_pEncodeInfo;
451     SampleDecodeInfo*   m_pDecodeInfo;
452     FX_DWORD    m_nBitsPerSample, m_SampleMax;
453     CPDF_StreamAcc*     m_pSampleStream;
454 };
455 CPDF_SampledFunc::CPDF_SampledFunc()
456 {
457     m_pSampleStream = NULL;
458     m_pEncodeInfo = NULL;
459     m_pDecodeInfo = NULL;
460 }
461 CPDF_SampledFunc::~CPDF_SampledFunc()
462 {
463     delete m_pSampleStream;
464     if (m_pEncodeInfo) {
465         FX_Free(m_pEncodeInfo);
466     }
467     if (m_pDecodeInfo) {
468         FX_Free(m_pDecodeInfo);
469     }
470 }
471 FX_BOOL CPDF_SampledFunc::v_Init(CPDF_Object* pObj)
472 {
473     if (pObj->GetType() != PDFOBJ_STREAM) {
474         return FALSE;
475     }
476     CPDF_Stream* pStream = (CPDF_Stream*)pObj;
477     CPDF_Dictionary* pDict = pStream->GetDict();
478     CPDF_Array* pSize = pDict->GetArray(FX_BSTRC("Size"));
479     CPDF_Array* pEncode = pDict->GetArray(FX_BSTRC("Encode"));
480     CPDF_Array* pDecode = pDict->GetArray(FX_BSTRC("Decode"));
481     m_nBitsPerSample = pDict->GetInteger(FX_BSTRC("BitsPerSample"));
482     m_SampleMax = 0xffffffff >> (32 - m_nBitsPerSample);
483     m_pSampleStream = new CPDF_StreamAcc;
484     m_pSampleStream->LoadAllData(pStream, FALSE);
485     m_pEncodeInfo = FX_Alloc(SampleEncodeInfo, m_nInputs);
486     int i;
487     FX_DWORD nTotalSamples = 1;
488     for (i = 0; i < m_nInputs; i ++) {
489         m_pEncodeInfo[i].sizes = pSize ? pSize->GetInteger(i) : 0;
490         if (!pSize && i == 0) {
491             m_pEncodeInfo[i].sizes = pDict->GetInteger(FX_BSTRC("Size"));
492         }
493         if (nTotalSamples > 0 && (uint32_t)(m_pEncodeInfo[i].sizes) > UINT_MAX / nTotalSamples) {
494             return FALSE;
495         }
496         nTotalSamples *= m_pEncodeInfo[i].sizes;
497         if (pEncode) {
498             m_pEncodeInfo[i].encode_min = pEncode->GetFloat(i * 2);
499             m_pEncodeInfo[i].encode_max = pEncode->GetFloat(i * 2 + 1);
500         } else {
501             m_pEncodeInfo[i].encode_min = 0;
502             if (m_pEncodeInfo[i].sizes == 1) {
503                 m_pEncodeInfo[i].encode_max = 1;
504             } else {
505                 m_pEncodeInfo[i].encode_max = (FX_FLOAT)m_pEncodeInfo[i].sizes - 1;
506             }
507         }
508     }
509     if (nTotalSamples > 0 && m_nBitsPerSample > UINT_MAX / nTotalSamples) {
510         return FALSE;
511     }
512     nTotalSamples *= m_nBitsPerSample;
513     if (nTotalSamples > 0 && ((uint32_t)m_nOutputs) > UINT_MAX / nTotalSamples) {
514         return FALSE;
515     }
516     nTotalSamples *= m_nOutputs;
517     if (nTotalSamples == 0 || m_pSampleStream->GetSize() * 8 < nTotalSamples) {
518         return FALSE;
519     }
520     m_pDecodeInfo = FX_Alloc(SampleDecodeInfo, m_nOutputs);
521     for (i = 0; i < m_nOutputs; i ++) {
522         if (pDecode) {
523             m_pDecodeInfo[i].decode_min = pDecode->GetFloat(2 * i);
524             m_pDecodeInfo[i].decode_max = pDecode->GetFloat(2 * i + 1);
525         } else {
526             m_pDecodeInfo[i].decode_min = m_pRanges[i * 2];
527             m_pDecodeInfo[i].decode_max = m_pRanges[i * 2 + 1];
528         }
529     }
530     return TRUE;
531 }
532 FX_BOOL CPDF_SampledFunc::v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const
533 {
534     int pos = 0;
535     CFX_FixedBufGrow<FX_FLOAT, 16> encoded_input_buf(m_nInputs);
536     FX_FLOAT* encoded_input = encoded_input_buf;
537     CFX_FixedBufGrow<int, 32> int_buf(m_nInputs * 2);
538     int* index = int_buf;
539     int* blocksize = index + m_nInputs;
540     for (int i = 0; i < m_nInputs; i ++) {
541         if (i == 0) {
542             blocksize[i] = 1;
543         } else {
544             blocksize[i] = blocksize[i - 1] * m_pEncodeInfo[i - 1].sizes;
545         }
546         encoded_input[i] = PDF_Interpolate(inputs[i], m_pDomains[i * 2], m_pDomains[i * 2 + 1],
547                                            m_pEncodeInfo[i].encode_min, m_pEncodeInfo[i].encode_max);
548         index[i] = (int)encoded_input[i];
549         if (index[i] < 0) {
550             index[i] = 0;
551         } else if (index[i] > m_pEncodeInfo[i].sizes - 1) {
552             index[i] = m_pEncodeInfo[i].sizes - 1;
553         }
554         pos += index[i] * blocksize[i];
555     }
556     FX_SAFE_INT32 bitpos = pos;
557     bitpos *= m_nBitsPerSample;
558     bitpos *= m_nOutputs;
559     if (!bitpos.IsValid()) {
560         return FALSE;
561     }
562     const uint8_t* pSampleData = m_pSampleStream->GetData();
563     if (pSampleData == NULL) {
564         return FALSE;
565     }
566     FX_SAFE_INT32 bitpos1 = m_nOutputs - 1 > 0 ? m_nOutputs - 1 : 0;
567     bitpos1 *= m_nBitsPerSample;
568     bitpos1 += bitpos.ValueOrDie();
569     if (!bitpos1.IsValid()) {
570         return FALSE;
571     }
572     for (int j = 0; j < m_nOutputs; j ++) {
573         FX_DWORD sample = _GetBits32(pSampleData, bitpos.ValueOrDie() + j * m_nBitsPerSample, m_nBitsPerSample);
574         FX_FLOAT encoded = (FX_FLOAT)sample;
575         for (int i = 0; i < m_nInputs; i ++) {
576             if (index[i] == m_pEncodeInfo[i].sizes - 1) {
577                 if (index[i] == 0) {
578                     encoded = encoded_input[i] * (FX_FLOAT)sample;
579                 }
580             } else {
581                 FX_SAFE_INT32 bitpos2 = blocksize[i];
582                 bitpos2 += pos;
583                 bitpos2 *= m_nOutputs;
584                 bitpos2 += j;
585                 bitpos2 *= m_nBitsPerSample;
586                 if (!bitpos2.IsValid()) {
587                     return FALSE;
588                 }
589                 FX_DWORD sample1 = _GetBits32(pSampleData, bitpos2.ValueOrDie(), m_nBitsPerSample);
590                 encoded += (encoded_input[i] - index[i]) * ((FX_FLOAT)sample1 - (FX_FLOAT)sample);
591             }
592         }
593         results[j] = PDF_Interpolate(encoded, 0, (FX_FLOAT)m_SampleMax,
594                                      m_pDecodeInfo[j].decode_min, m_pDecodeInfo[j].decode_max);
595     }
596     return TRUE;
597 }
598 class CPDF_PSFunc : public CPDF_Function
599 {
600 public:
601     virtual FX_BOOL             v_Init(CPDF_Object* pObj);
602     virtual FX_BOOL             v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const;
603     CPDF_PSEngine m_PS;
604 };
605 FX_BOOL CPDF_PSFunc::v_Init(CPDF_Object* pObj)
606 {
607     CPDF_Stream* pStream = (CPDF_Stream*)pObj;
608     CPDF_StreamAcc acc;
609     acc.LoadAllData(pStream, FALSE);
610     return m_PS.Parse((const FX_CHAR*)acc.GetData(), acc.GetSize());
611 }
612 FX_BOOL CPDF_PSFunc::v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const
613 {
614     CPDF_PSEngine& PS = (CPDF_PSEngine&)m_PS;
615     PS.Reset();
616     int i;
617     for (i = 0; i < m_nInputs; i ++) {
618         PS.Push(inputs[i]);
619     }
620     PS.Execute();
621     if (PS.GetStackSize() < m_nOutputs) {
622         return FALSE;
623     }
624     for (i = 0; i < m_nOutputs; i ++) {
625         results[m_nOutputs - i - 1] = PS.Pop();
626     }
627     return TRUE;
628 }
629 class CPDF_ExpIntFunc : public CPDF_Function
630 {
631 public:
632     CPDF_ExpIntFunc();
633     virtual ~CPDF_ExpIntFunc();
634     virtual FX_BOOL             v_Init(CPDF_Object* pObj);
635     virtual FX_BOOL             v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const;
636     FX_FLOAT    m_Exponent;
637     FX_FLOAT*   m_pBeginValues;
638     FX_FLOAT*   m_pEndValues;
639     int         m_nOrigOutputs;
640 };
641 CPDF_ExpIntFunc::CPDF_ExpIntFunc()
642 {
643     m_pBeginValues = NULL;
644     m_pEndValues = NULL;
645 }
646 CPDF_ExpIntFunc::~CPDF_ExpIntFunc()
647 {
648     if (m_pBeginValues) {
649         FX_Free(m_pBeginValues);
650     }
651     if (m_pEndValues) {
652         FX_Free(m_pEndValues);
653     }
654 }
655 FX_BOOL CPDF_ExpIntFunc::v_Init(CPDF_Object* pObj)
656 {
657     CPDF_Dictionary* pDict = pObj->GetDict();
658     if (pDict == NULL) {
659         return FALSE;
660     }
661     CPDF_Array* pArray0 = pDict->GetArray(FX_BSTRC("C0"));
662     if (m_nOutputs == 0) {
663         m_nOutputs = 1;
664         if (pArray0) {
665             m_nOutputs = pArray0->GetCount();
666         }
667     }
668     CPDF_Array* pArray1 = pDict->GetArray(FX_BSTRC("C1"));
669     m_pBeginValues = FX_Alloc2D(FX_FLOAT, m_nOutputs, 2);
670     m_pEndValues = FX_Alloc2D(FX_FLOAT, m_nOutputs, 2);
671     for (int i = 0; i < m_nOutputs; i ++) {
672         m_pBeginValues[i] = pArray0 ? pArray0->GetFloat(i) : 0.0f;
673         m_pEndValues[i] = pArray1 ? pArray1->GetFloat(i) : 1.0f;
674     }
675     m_Exponent = pDict->GetFloat(FX_BSTRC("N"));
676     m_nOrigOutputs = m_nOutputs;
677     if (m_nOutputs && m_nInputs > INT_MAX / m_nOutputs) {
678         return FALSE;
679     }
680     m_nOutputs *= m_nInputs;
681     return TRUE;
682 }
683 FX_BOOL CPDF_ExpIntFunc::v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const
684 {
685     for (int i = 0; i < m_nInputs; i ++)
686         for (int j = 0; j < m_nOrigOutputs; j ++) {
687             results[i * m_nOrigOutputs + j] = m_pBeginValues[j] + (FX_FLOAT)FXSYS_pow(inputs[i], m_Exponent) *
688                                               (m_pEndValues[j] - m_pBeginValues[j]);
689         }
690     return TRUE;
691 }
692 class CPDF_StitchFunc : public CPDF_Function
693 {
694 public:
695     CPDF_StitchFunc();
696     virtual ~CPDF_StitchFunc();
697     virtual FX_BOOL             v_Init(CPDF_Object* pObj);
698     virtual FX_BOOL             v_Call(FX_FLOAT* inputs, FX_FLOAT* results) const;
699     int                 m_nSubs;
700     CPDF_Function** m_pSubFunctions;
701     FX_FLOAT*   m_pBounds;
702     FX_FLOAT*   m_pEncode;
703 };
704 CPDF_StitchFunc::CPDF_StitchFunc()
705 {
706     m_nSubs = 0;
707     m_pSubFunctions = NULL;
708     m_pBounds = NULL;
709     m_pEncode = NULL;
710 }
711 CPDF_StitchFunc::~CPDF_StitchFunc()
712 {
713     for (int i = 0; i < m_nSubs; i ++)
714         delete m_pSubFunctions[i];
715     if (m_pSubFunctions) {
716         FX_Free(m_pSubFunctions);
717     }
718     if (m_pBounds) {
719         FX_Free(m_pBounds);
720     }
721     if (m_pEncode) {
722         FX_Free(m_pEncode);
723     }
724 }
725 FX_BOOL CPDF_StitchFunc::v_Init(CPDF_Object* pObj)
726 {
727     CPDF_Dictionary* pDict = pObj->GetDict();
728     if (pDict == NULL) {
729         return FALSE;
730     }
731     CPDF_Array* pArray = pDict->GetArray(FX_BSTRC("Functions"));
732     if (pArray == NULL) {
733         return FALSE;
734     }
735     m_nSubs = pArray->GetCount();
736     if (m_nSubs == 0) {
737         return FALSE;
738     }
739     m_pSubFunctions = FX_Alloc(CPDF_Function*, m_nSubs);
740     m_nOutputs = 0;
741     int i;
742     for (i = 0; i < m_nSubs; i ++) {
743         CPDF_Object* pSub = pArray->GetElementValue(i);
744         if (pSub == pObj) {
745             return FALSE;
746         }
747         m_pSubFunctions[i] = CPDF_Function::Load(pSub);
748         if (m_pSubFunctions[i] == NULL) {
749             return FALSE;
750         }
751         if (m_pSubFunctions[i]->CountOutputs() > m_nOutputs) {
752             m_nOutputs = m_pSubFunctions[i]->CountOutputs();
753         }
754     }
755     m_pBounds = FX_Alloc(FX_FLOAT, m_nSubs + 1);
756     m_pBounds[0] = m_pDomains[0];
757     pArray = pDict->GetArray(FX_BSTRC("Bounds"));
758     if (pArray == NULL) {
759         return FALSE;
760     }
761     for (i = 0; i < m_nSubs - 1; i ++) {
762         m_pBounds[i + 1] = pArray->GetFloat(i);
763     }
764     m_pBounds[m_nSubs] = m_pDomains[1];
765     m_pEncode = FX_Alloc2D(FX_FLOAT, m_nSubs, 2);
766     pArray = pDict->GetArray(FX_BSTRC("Encode"));
767     if (pArray == NULL) {
768         return FALSE;
769     }
770     for (i = 0; i < m_nSubs * 2; i ++) {
771         m_pEncode[i] = pArray->GetFloat(i);
772     }
773     return TRUE;
774 }
775 FX_BOOL CPDF_StitchFunc::v_Call(FX_FLOAT* inputs, FX_FLOAT* outputs) const
776 {
777     FX_FLOAT input = inputs[0];
778     int i;
779     for (i = 0; i < m_nSubs - 1; i ++)
780         if (input < m_pBounds[i + 1]) {
781             break;
782         }
783     if (m_pSubFunctions[i] == NULL) {
784         return FALSE;
785     }
786     input = PDF_Interpolate(input, m_pBounds[i], m_pBounds[i + 1], m_pEncode[i * 2], m_pEncode[i * 2 + 1]);
787     int nresults;
788     m_pSubFunctions[i]->Call(&input, m_nInputs, outputs, nresults);
789     return TRUE;
790 }
791 CPDF_Function* CPDF_Function::Load(CPDF_Object* pFuncObj)
792 {
793     if (pFuncObj == NULL) {
794         return NULL;
795     }
796     CPDF_Function* pFunc = NULL;
797     int type;
798     if (pFuncObj->GetType() == PDFOBJ_STREAM) {
799         type = ((CPDF_Stream*)pFuncObj)->GetDict()->GetInteger(FX_BSTRC("FunctionType"));
800     } else if (pFuncObj->GetType() == PDFOBJ_DICTIONARY) {
801         type = ((CPDF_Dictionary*)pFuncObj)->GetInteger(FX_BSTRC("FunctionType"));
802     } else {
803         return NULL;
804     }
805     if (type == 0) {
806         pFunc = new CPDF_SampledFunc;
807     } else if (type == 2) {
808         pFunc = new CPDF_ExpIntFunc;
809     } else if (type == 3) {
810         pFunc = new CPDF_StitchFunc;
811     } else if (type == 4) {
812         pFunc = new CPDF_PSFunc;
813     } else {
814         return NULL;
815     }
816     if (!pFunc->Init(pFuncObj)) {
817         delete pFunc;
818         return NULL;
819     }
820     return pFunc;
821 }
822 CPDF_Function::CPDF_Function()
823 {
824     m_pDomains = NULL;
825     m_pRanges = NULL;
826 }
827 CPDF_Function::~CPDF_Function()
828 {
829     if (m_pDomains) {
830         FX_Free(m_pDomains);
831         m_pDomains = NULL;
832     }
833     if (m_pRanges) {
834         FX_Free(m_pRanges);
835         m_pRanges = NULL;
836     }
837 }
838 FX_BOOL CPDF_Function::Init(CPDF_Object* pObj)
839 {
840     CPDF_Dictionary* pDict;
841     if (pObj->GetType() == PDFOBJ_STREAM) {
842         pDict = ((CPDF_Stream*)pObj)->GetDict();
843     } else {
844         pDict = (CPDF_Dictionary*)pObj;
845     }
846     CPDF_Array* pDomains = pDict->GetArray(FX_BSTRC("Domain"));
847     if (pDomains == NULL) {
848         return FALSE;
849     }
850     m_nInputs = pDomains->GetCount() / 2;
851     if (m_nInputs == 0) {
852         return FALSE;
853     }
854     m_pDomains = FX_Alloc2D(FX_FLOAT, m_nInputs, 2);
855     for (int i = 0; i < m_nInputs * 2; i ++) {
856         m_pDomains[i] = pDomains->GetFloat(i);
857     }
858     CPDF_Array* pRanges = pDict->GetArray(FX_BSTRC("Range"));
859     m_nOutputs = 0;
860     if (pRanges) {
861         m_nOutputs = pRanges->GetCount() / 2;
862         m_pRanges = FX_Alloc2D(FX_FLOAT, m_nOutputs, 2);
863         for (int i = 0; i < m_nOutputs * 2; i ++) {
864             m_pRanges[i] = pRanges->GetFloat(i);
865         }
866     }
867     FX_DWORD old_outputs = m_nOutputs;
868     FX_BOOL ret = v_Init(pObj);
869     if (m_pRanges && m_nOutputs > (int)old_outputs) {
870         m_pRanges = FX_Realloc(FX_FLOAT, m_pRanges, m_nOutputs * 2);
871         if (m_pRanges) {
872             FXSYS_memset(m_pRanges + (old_outputs * 2), 0, sizeof(FX_FLOAT) * (m_nOutputs - old_outputs) * 2);
873         }
874     }
875     return ret;
876 }
877 FX_BOOL CPDF_Function::Call(FX_FLOAT* inputs, int ninputs, FX_FLOAT* results, int& nresults) const
878 {
879     if (m_nInputs != ninputs) {
880         return FALSE;
881     }
882     nresults = m_nOutputs;
883     for (int i = 0; i < m_nInputs; i ++) {
884         if (inputs[i] < m_pDomains[i * 2]) {
885             inputs[i] = m_pDomains[i * 2];
886         } else if (inputs[i] > m_pDomains[i * 2 + 1]) {
887             inputs[i] = m_pDomains[i * 2] + 1;
888         }
889     }
890     v_Call(inputs, results);
891     if (m_pRanges) {
892         for (int i = 0; i < m_nOutputs; i ++) {
893             if (results[i] < m_pRanges[i * 2]) {
894                 results[i] = m_pRanges[i * 2];
895             } else if (results[i] > m_pRanges[i * 2 + 1]) {
896                 results[i] = m_pRanges[i * 2 + 1];
897             }
898         }
899     }
900     return TRUE;
901 }