Line data Source code
1 : /* -*- Mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
2 : /*
3 : * This file is part of the LibreOffice project.
4 : *
5 : * This Source Code Form is subject to the terms of the Mozilla Public
6 : * License, v. 2.0. If a copy of the MPL was not distributed with this
7 : * file, You can obtain one at http://mozilla.org/MPL/2.0/.
8 : *
9 : * This file incorporates work covered by the following license notice:
10 : *
11 : * Licensed to the Apache Software Foundation (ASF) under one or more
12 : * contributor license agreements. See the NOTICE file distributed
13 : * with this work for additional information regarding copyright
14 : * ownership. The ASF licenses this file to you under the Apache
15 : * License, Version 2.0 (the "License"); you may not use this file
16 : * except in compliance with the License. You may obtain a copy of
17 : * the License at http://www.apache.org/licenses/LICENSE-2.0 .
18 : */
19 :
20 : #include <math.h>
21 :
22 : #include <rtl/ustrbuf.hxx>
23 : #include <tools/bigint.hxx>
24 :
25 :
26 : #include <string.h>
27 : #include <ctype.h>
28 :
29 : static const long MY_MAXLONG = 0x3fffffff;
30 : static const long MY_MINLONG = -MY_MAXLONG;
31 : static const long MY_MAXSHORT = 0x00007fff;
32 : static const long MY_MINSHORT = -MY_MAXSHORT;
33 :
34 : /*
35 : * The algorithms for Addition, Substraction, Multiplication and Divison
36 : * of large numbers originate from SEMINUMERICAL ALGORITHMS by
37 : * DONALD E. KNUTH in the series The Art of Computer Programming:
38 : * chapter 4.3.1. The Classical Algorithms.
39 : */
40 :
41 : // TODO: Needs conversion to sal_uInt16/INT16/sal_uInt32/sal_Int32
42 0 : void BigInt::MakeBigInt( const BigInt& rVal )
43 : {
44 0 : if ( rVal.bIsBig )
45 : {
46 0 : memcpy( (void*)this, (const void*)&rVal, sizeof( BigInt ) );
47 0 : while ( nLen > 1 && nNum[nLen-1] == 0 )
48 0 : nLen--;
49 : }
50 : else
51 : {
52 0 : long nTmp = rVal.nVal;
53 :
54 0 : nVal = rVal.nVal;
55 0 : bIsBig = true;
56 0 : if ( nTmp < 0 )
57 : {
58 0 : bIsNeg = true;
59 0 : nTmp = -nTmp;
60 : }
61 : else
62 0 : bIsNeg = false;
63 :
64 0 : nNum[0] = (sal_uInt16)(nTmp & 0xffffL);
65 0 : nNum[1] = (sal_uInt16)(nTmp >> 16);
66 0 : if ( nTmp & 0xffff0000L )
67 0 : nLen = 2;
68 : else
69 0 : nLen = 1;
70 : }
71 0 : }
72 :
73 0 : void BigInt::Normalize()
74 : {
75 0 : if ( bIsBig )
76 : {
77 0 : while ( nLen > 1 && nNum[nLen-1] == 0 )
78 0 : nLen--;
79 :
80 0 : if ( nLen < 3 )
81 : {
82 0 : if ( nLen < 2 )
83 0 : nVal = nNum[0];
84 0 : else if ( nNum[1] & 0x8000 )
85 0 : return;
86 : else
87 0 : nVal = ((long)nNum[1] << 16) + nNum[0];
88 :
89 0 : bIsBig = false;
90 :
91 0 : if ( bIsNeg )
92 0 : nVal = -nVal;
93 : }
94 : // else nVal is undefined !!! W.P.
95 : }
96 : // why? nVal is undefined ??? W.P.
97 0 : else if ( nVal & 0xFFFF0000L )
98 0 : nLen = 2;
99 : else
100 0 : nLen = 1;
101 : }
102 :
103 0 : void BigInt::Mult( const BigInt &rVal, sal_uInt16 nMul )
104 : {
105 0 : sal_uInt16 nK = 0;
106 0 : for ( int i = 0; i < rVal.nLen; i++ )
107 : {
108 0 : sal_uInt32 nTmp = (sal_uInt32)rVal.nNum[i] * (sal_uInt32)nMul + nK;
109 0 : nK = (sal_uInt16)(nTmp >> 16);
110 0 : nNum[i] = (sal_uInt16)nTmp;
111 : }
112 :
113 0 : if ( nK )
114 : {
115 0 : nNum[rVal.nLen] = nK;
116 0 : nLen = rVal.nLen + 1;
117 : }
118 : else
119 0 : nLen = rVal.nLen;
120 :
121 0 : bIsBig = true;
122 0 : bIsNeg = rVal.bIsNeg;
123 0 : }
124 :
125 0 : void BigInt::Div( sal_uInt16 nDiv, sal_uInt16& rRem )
126 : {
127 0 : sal_uInt32 nK = 0;
128 0 : for ( int i = nLen - 1; i >= 0; i-- )
129 : {
130 0 : sal_uInt32 nTmp = (sal_uInt32)nNum[i] + (nK << 16);
131 0 : nNum[i] = (sal_uInt16)(nTmp / nDiv);
132 0 : nK = nTmp % nDiv;
133 : }
134 0 : rRem = (sal_uInt16)nK;
135 :
136 0 : if ( nNum[nLen-1] == 0 )
137 0 : nLen -= 1;
138 0 : }
139 :
140 0 : bool BigInt::IsLess( const BigInt& rVal ) const
141 : {
142 0 : if ( rVal.nLen < nLen)
143 0 : return true;
144 0 : if ( rVal.nLen > nLen )
145 0 : return false;
146 :
147 : int i;
148 0 : for ( i = nLen - 1; i > 0 && nNum[i] == rVal.nNum[i]; i-- )
149 : {
150 : }
151 0 : return rVal.nNum[i] < nNum[i];
152 : }
153 :
154 0 : void BigInt::AddLong( BigInt& rB, BigInt& rErg )
155 : {
156 0 : if ( bIsNeg == rB.bIsNeg )
157 : {
158 : int i;
159 : char len;
160 :
161 : // if length of the two values differ, fill remaining positions
162 : // of the smaller value with zeros.
163 0 : if (nLen >= rB.nLen)
164 : {
165 0 : len = nLen;
166 0 : for (i = rB.nLen; i < len; i++)
167 0 : rB.nNum[i] = 0;
168 : }
169 : else
170 : {
171 0 : len = rB.nLen;
172 0 : for (i = nLen; i < len; i++)
173 0 : nNum[i] = 0;
174 : }
175 :
176 : // Add numerals, starting from the back
177 : long k;
178 0 : long nZ = 0;
179 0 : for (i = 0, k = 0; i < len; i++) {
180 0 : nZ = (long)nNum[i] + (long)rB.nNum[i] + k;
181 0 : if (nZ & 0xff0000L)
182 0 : k = 1;
183 : else
184 0 : k = 0;
185 0 : rErg.nNum[i] = (sal_uInt16)(nZ & 0xffffL);
186 : }
187 : // If an overflow occurred, add to solution
188 0 : if (nZ & 0xff0000L) // or if(k)
189 : {
190 0 : rErg.nNum[i] = 1;
191 0 : len++;
192 : }
193 : // Set length and sign
194 0 : rErg.nLen = len;
195 0 : rErg.bIsNeg = bIsNeg && rB.bIsNeg;
196 0 : rErg.bIsBig = true;
197 : }
198 : // If one of the values is negative, perform substraction instead
199 0 : else if (bIsNeg)
200 : {
201 0 : bIsNeg = false;
202 0 : rB.SubLong(*this, rErg);
203 0 : bIsNeg = true;
204 : }
205 : else
206 : {
207 0 : rB.bIsNeg = false;
208 0 : SubLong(rB, rErg);
209 0 : rB.bIsNeg = true;
210 : }
211 0 : }
212 :
213 0 : void BigInt::SubLong( BigInt& rB, BigInt& rErg )
214 : {
215 0 : if ( bIsNeg == rB.bIsNeg )
216 : {
217 : int i;
218 : char len;
219 : long nZ, k;
220 :
221 : // if length of the two values differ, fill remaining positions
222 : // of the smaller value with zeros.
223 0 : if (nLen >= rB.nLen)
224 : {
225 0 : len = nLen;
226 0 : for (i = rB.nLen; i < len; i++)
227 0 : rB.nNum[i] = 0;
228 : }
229 : else
230 : {
231 0 : len = rB.nLen;
232 0 : for (i = nLen; i < len; i++)
233 0 : nNum[i] = 0;
234 : }
235 :
236 0 : if ( IsLess(rB) )
237 : {
238 0 : for (i = 0, k = 0; i < len; i++)
239 : {
240 0 : nZ = (long)nNum[i] - (long)rB.nNum[i] + k;
241 0 : if (nZ < 0)
242 0 : k = -1;
243 : else
244 0 : k = 0;
245 0 : rErg.nNum[i] = (sal_uInt16)(nZ & 0xffffL);
246 : }
247 0 : rErg.bIsNeg = bIsNeg;
248 : }
249 : else
250 : {
251 0 : for (i = 0, k = 0; i < len; i++)
252 : {
253 0 : nZ = (long)rB.nNum[i] - (long)nNum[i] + k;
254 0 : if (nZ < 0)
255 0 : k = -1;
256 : else
257 0 : k = 0;
258 0 : rErg.nNum[i] = (sal_uInt16)(nZ & 0xffffL);
259 : }
260 : // if a < b, revert sign
261 0 : rErg.bIsNeg = !bIsNeg;
262 : }
263 0 : rErg.nLen = len;
264 0 : rErg.bIsBig = true;
265 : }
266 : // If one of the values is negative, perform addition instead
267 0 : else if (bIsNeg)
268 : {
269 0 : bIsNeg = false;
270 0 : AddLong(rB, rErg);
271 0 : bIsNeg = true;
272 0 : rErg.bIsNeg = true;
273 : }
274 : else
275 : {
276 0 : rB.bIsNeg = false;
277 0 : AddLong(rB, rErg);
278 0 : rB.bIsNeg = true;
279 0 : rErg.bIsNeg = false;
280 : }
281 0 : }
282 :
283 0 : void BigInt::MultLong( const BigInt& rB, BigInt& rErg ) const
284 : {
285 : int i, j;
286 : sal_uInt32 nZ, k;
287 :
288 0 : rErg.bIsNeg = bIsNeg != rB.bIsNeg;
289 0 : rErg.bIsBig = true;
290 0 : rErg.nLen = nLen + rB.nLen;
291 :
292 0 : for (i = 0; i < rErg.nLen; i++)
293 0 : rErg.nNum[i] = 0;
294 :
295 0 : for (j = 0; j < rB.nLen; j++)
296 : {
297 0 : for (i = 0, k = 0; i < nLen; i++)
298 : {
299 0 : nZ = (sal_uInt32)nNum[i] * (sal_uInt32)rB.nNum[j] +
300 0 : (sal_uInt32)rErg.nNum[i + j] + k;
301 0 : rErg.nNum[i + j] = (sal_uInt16)(nZ & 0xffffUL);
302 0 : k = nZ >> 16;
303 : }
304 0 : rErg.nNum[i + j] = (sal_uInt16)k;
305 : }
306 0 : }
307 :
308 0 : void BigInt::DivLong( const BigInt& rB, BigInt& rErg ) const
309 : {
310 : int i, j;
311 : long nTmp;
312 : sal_uInt16 nK, nQ, nMult;
313 0 : short nLenB = rB.nLen;
314 0 : short nLenB1 = rB.nLen - 1;
315 0 : BigInt aTmpA, aTmpB;
316 :
317 0 : nMult = (sal_uInt16)(0x10000L / ((long)rB.nNum[nLenB1] + 1));
318 :
319 0 : aTmpA.Mult( *this, nMult );
320 0 : if ( aTmpA.nLen == nLen )
321 : {
322 0 : aTmpA.nNum[aTmpA.nLen] = 0;
323 0 : aTmpA.nLen++;
324 : }
325 :
326 0 : aTmpB.Mult( rB, nMult );
327 :
328 0 : for (j = aTmpA.nLen - 1; j >= nLenB; j--)
329 : { // guess divisor
330 0 : nTmp = ( (long)aTmpA.nNum[j] << 16 ) + aTmpA.nNum[j - 1];
331 0 : if (aTmpA.nNum[j] == aTmpB.nNum[nLenB1])
332 0 : nQ = 0xFFFF;
333 : else
334 0 : nQ = (sal_uInt16)(((sal_uInt32)nTmp) / aTmpB.nNum[nLenB1]);
335 :
336 0 : if ( ((sal_uInt32)aTmpB.nNum[nLenB1 - 1] * nQ) >
337 0 : ((((sal_uInt32)nTmp) - aTmpB.nNum[nLenB1] * nQ) << 16) + aTmpA.nNum[j - 2])
338 0 : nQ--;
339 : // Start division
340 0 : nK = 0;
341 0 : nTmp = 0;
342 0 : for (i = 0; i < nLenB; i++)
343 : {
344 0 : nTmp = (long)aTmpA.nNum[j - nLenB + i]
345 0 : - ((long)aTmpB.nNum[i] * nQ)
346 0 : - nK;
347 0 : aTmpA.nNum[j - nLenB + i] = (sal_uInt16)nTmp;
348 0 : nK = (sal_uInt16) (nTmp >> 16);
349 0 : if ( nK )
350 0 : nK = (sal_uInt16)(0x10000UL - nK);
351 : }
352 0 : unsigned short& rNum( aTmpA.nNum[j - nLenB + i] );
353 0 : rNum = rNum - nK; // MSVC yields a warning on -= here, so don't use it
354 0 : if (aTmpA.nNum[j - nLenB + i] == 0)
355 0 : rErg.nNum[j - nLenB] = nQ;
356 : else
357 : {
358 0 : rErg.nNum[j - nLenB] = nQ - 1;
359 0 : nK = 0;
360 0 : for (i = 0; i < nLenB; i++)
361 : {
362 0 : nTmp = aTmpA.nNum[j - nLenB + i] + aTmpB.nNum[i] + nK;
363 0 : aTmpA.nNum[j - nLenB + i] = (sal_uInt16)(nTmp & 0xFFFFL);
364 0 : if (nTmp & 0xFFFF0000L)
365 0 : nK = 1;
366 : else
367 0 : nK = 0;
368 : }
369 : }
370 : }
371 :
372 0 : rErg.bIsNeg = bIsNeg != rB.bIsNeg;
373 0 : rErg.bIsBig = true;
374 0 : rErg.nLen = nLen - rB.nLen + 1;
375 0 : }
376 :
377 0 : void BigInt::ModLong( const BigInt& rB, BigInt& rErg ) const
378 : {
379 : short i, j;
380 : long nTmp;
381 : sal_uInt16 nK, nQ, nMult;
382 0 : short nLenB = rB.nLen;
383 0 : short nLenB1 = rB.nLen - 1;
384 0 : BigInt aTmpA, aTmpB;
385 :
386 0 : nMult = (sal_uInt16)(0x10000L / ((long)rB.nNum[nLenB1] + 1));
387 :
388 0 : aTmpA.Mult( *this, nMult);
389 0 : if ( aTmpA.nLen == nLen )
390 : {
391 0 : aTmpA.nNum[aTmpA.nLen] = 0;
392 0 : aTmpA.nLen++;
393 : }
394 :
395 0 : aTmpB.Mult( rB, nMult);
396 :
397 0 : for (j = aTmpA.nLen - 1; j >= nLenB; j--)
398 : { // Guess divisor
399 0 : nTmp = ( (long)aTmpA.nNum[j] << 16 ) + aTmpA.nNum[j - 1];
400 0 : if (aTmpA.nNum[j] == aTmpB.nNum[nLenB1])
401 0 : nQ = 0xFFFF;
402 : else
403 0 : nQ = (sal_uInt16)(((sal_uInt32)nTmp) / aTmpB.nNum[nLenB1]);
404 :
405 0 : if ( ((sal_uInt32)aTmpB.nNum[nLenB1 - 1] * nQ) >
406 0 : ((((sal_uInt32)nTmp) - aTmpB.nNum[nLenB1] * nQ) << 16) + aTmpA.nNum[j - 2])
407 0 : nQ--;
408 : // Start division
409 0 : nK = 0;
410 0 : nTmp = 0;
411 0 : for (i = 0; i < nLenB; i++)
412 : {
413 0 : nTmp = (long)aTmpA.nNum[j - nLenB + i]
414 0 : - ((long)aTmpB.nNum[i] * nQ)
415 0 : - nK;
416 0 : aTmpA.nNum[j - nLenB + i] = (sal_uInt16)nTmp;
417 0 : nK = (sal_uInt16) (nTmp >> 16);
418 0 : if ( nK )
419 0 : nK = (sal_uInt16)(0x10000UL - nK);
420 : }
421 0 : unsigned short& rNum( aTmpA.nNum[j - nLenB + i] );
422 0 : rNum = rNum - nK;
423 0 : if (aTmpA.nNum[j - nLenB + i] == 0)
424 0 : rErg.nNum[j - nLenB] = nQ;
425 : else
426 : {
427 0 : rErg.nNum[j - nLenB] = nQ - 1;
428 0 : nK = 0;
429 0 : for (i = 0; i < nLenB; i++) {
430 0 : nTmp = aTmpA.nNum[j - nLenB + i] + aTmpB.nNum[i] + nK;
431 0 : aTmpA.nNum[j - nLenB + i] = (sal_uInt16)(nTmp & 0xFFFFL);
432 0 : if (nTmp & 0xFFFF0000L)
433 0 : nK = 1;
434 : else
435 0 : nK = 0;
436 : }
437 : }
438 : }
439 :
440 0 : rErg = aTmpA;
441 0 : rErg.Div( nMult, nQ );
442 0 : }
443 :
444 0 : bool BigInt::ABS_IsLess( const BigInt& rB ) const
445 : {
446 0 : if (bIsBig || rB.bIsBig)
447 : {
448 0 : BigInt nA, nB;
449 0 : nA.MakeBigInt( *this );
450 0 : nB.MakeBigInt( rB );
451 0 : if (nA.nLen == nB.nLen)
452 : {
453 : int i;
454 0 : for (i = nA.nLen - 1; i > 0 && nA.nNum[i] == nB.nNum[i]; i--)
455 : {
456 : }
457 0 : return nA.nNum[i] < nB.nNum[i];
458 : }
459 : else
460 0 : return nA.nLen < nB.nLen;
461 : }
462 0 : if ( nVal < 0 )
463 0 : if ( rB.nVal < 0 )
464 0 : return nVal > rB.nVal;
465 : else
466 0 : return nVal > -rB.nVal;
467 : else
468 0 : if ( rB.nVal < 0 )
469 0 : return nVal < -rB.nVal;
470 : else
471 0 : return nVal < rB.nVal;
472 : }
473 :
474 0 : BigInt::BigInt( const BigInt& rBigInt )
475 : {
476 0 : if ( rBigInt.bIsBig )
477 0 : memcpy( (void*)this, (const void*)&rBigInt, sizeof( BigInt ) );
478 : else
479 : {
480 0 : bIsSet = rBigInt.bIsSet;
481 0 : bIsBig = false;
482 0 : nVal = rBigInt.nVal;
483 : }
484 0 : }
485 :
486 0 : BigInt::BigInt( const OUString& rString )
487 0 : : nLen(0)
488 : {
489 0 : bIsSet = true;
490 0 : bIsNeg = false;
491 0 : bIsBig = false;
492 0 : nVal = 0;
493 :
494 0 : bool bNeg = false;
495 0 : const sal_Unicode* p = rString.getStr();
496 0 : if ( *p == '-' )
497 : {
498 0 : bNeg = true;
499 0 : p++;
500 : }
501 0 : while( *p >= '0' && *p <= '9' )
502 : {
503 0 : *this *= 10;
504 0 : *this += *p - '0';
505 0 : p++;
506 : }
507 0 : if ( bIsBig )
508 0 : bIsNeg = bNeg ? sal_True : sal_False;
509 0 : else if( bNeg )
510 0 : nVal = -nVal;
511 0 : }
512 :
513 0 : BigInt::BigInt( double nValue )
514 : {
515 0 : bIsSet = true;
516 :
517 0 : if ( nValue < 0 )
518 : {
519 0 : nValue *= -1;
520 0 : bIsNeg = true;
521 : }
522 : else
523 : {
524 0 : bIsNeg = false;
525 : }
526 :
527 0 : if ( nValue < 1 )
528 : {
529 0 : bIsBig = false;
530 0 : nVal = 0;
531 : }
532 : else
533 : {
534 0 : bIsBig = true;
535 :
536 0 : int i=0;
537 :
538 0 : while ( ( nValue > 65536.0 ) && ( i < MAX_DIGITS ) )
539 : {
540 0 : nNum[i] = (sal_uInt16) fmod( nValue, 65536.0 );
541 0 : nValue -= nNum[i];
542 0 : nValue /= 65536.0;
543 0 : i++;
544 : }
545 0 : if ( i < MAX_DIGITS )
546 0 : nNum[i++] = (sal_uInt16) nValue;
547 :
548 0 : nLen = i;
549 :
550 0 : if ( i < 3 )
551 0 : Normalize();
552 : }
553 0 : }
554 :
555 0 : BigInt::BigInt( sal_uInt32 nValue )
556 : {
557 0 : bIsSet = true;
558 0 : if ( nValue & 0x80000000UL )
559 : {
560 0 : bIsBig = true;
561 0 : bIsNeg = false;
562 0 : nNum[0] = (sal_uInt16)(nValue & 0xffffUL);
563 0 : nNum[1] = (sal_uInt16)(nValue >> 16);
564 0 : nLen = 2;
565 : }
566 : else
567 : {
568 0 : bIsBig = false;
569 0 : nVal = nValue;
570 : }
571 0 : }
572 :
573 0 : BigInt::operator sal_uIntPtr() const
574 : {
575 0 : if ( !bIsBig )
576 0 : return (sal_uInt32)nVal;
577 0 : else if ( nLen == 2 )
578 : {
579 : sal_uInt32 nRet;
580 0 : nRet = ((sal_uInt32)nNum[1]) << 16;
581 0 : nRet += nNum[0];
582 0 : return nRet;
583 : }
584 0 : return 0;
585 : }
586 :
587 0 : BigInt::operator double() const
588 : {
589 0 : if ( !bIsBig )
590 0 : return (double) nVal;
591 : else
592 : {
593 0 : int i = nLen-1;
594 0 : double nRet = (double) ((sal_uInt32)nNum[i]);
595 :
596 0 : while ( i )
597 : {
598 0 : nRet *= 65536.0;
599 0 : i--;
600 0 : nRet += (double) ((sal_uInt32)nNum[i]);
601 : }
602 :
603 0 : if ( bIsNeg )
604 0 : nRet *= -1;
605 :
606 0 : return nRet;
607 : }
608 : }
609 :
610 0 : BigInt& BigInt::operator=( const BigInt& rBigInt )
611 : {
612 0 : if ( rBigInt.bIsBig )
613 0 : memcpy( (void*)this, (const void*)&rBigInt, sizeof( BigInt ) );
614 : else
615 : {
616 0 : bIsSet = rBigInt.bIsSet;
617 0 : bIsBig = false;
618 0 : nVal = rBigInt.nVal;
619 : }
620 0 : return *this;
621 : }
622 :
623 0 : BigInt& BigInt::operator+=( const BigInt& rVal )
624 : {
625 0 : if ( !bIsBig && !rVal.bIsBig )
626 : {
627 0 : if( nVal <= MY_MAXLONG && rVal.nVal <= MY_MAXLONG
628 0 : && nVal >= MY_MINLONG && rVal.nVal >= MY_MINLONG )
629 : { // No overflows may occur here
630 0 : nVal += rVal.nVal;
631 0 : return *this;
632 : }
633 :
634 0 : if( (nVal < 0) != (rVal.nVal < 0) )
635 : { // No overflows may occur here
636 0 : nVal += rVal.nVal;
637 0 : return *this;
638 : }
639 : }
640 :
641 0 : BigInt aTmp1, aTmp2;
642 0 : aTmp1.MakeBigInt( *this );
643 0 : aTmp2.MakeBigInt( rVal );
644 0 : aTmp1.AddLong( aTmp2, *this );
645 0 : Normalize();
646 0 : return *this;
647 : }
648 :
649 0 : BigInt& BigInt::operator-=( const BigInt& rVal )
650 : {
651 0 : if ( !bIsBig && !rVal.bIsBig )
652 : {
653 0 : if ( nVal <= MY_MAXLONG && rVal.nVal <= MY_MAXLONG &&
654 0 : nVal >= MY_MINLONG && rVal.nVal >= MY_MINLONG )
655 : { // No overflows may occur here
656 0 : nVal -= rVal.nVal;
657 0 : return *this;
658 : }
659 :
660 0 : if ( (nVal < 0) == (rVal.nVal < 0) )
661 : { // No overflows may occur here
662 0 : nVal -= rVal.nVal;
663 0 : return *this;
664 : }
665 : }
666 :
667 0 : BigInt aTmp1, aTmp2;
668 0 : aTmp1.MakeBigInt( *this );
669 0 : aTmp2.MakeBigInt( rVal );
670 0 : aTmp1.SubLong( aTmp2, *this );
671 0 : Normalize();
672 0 : return *this;
673 : }
674 :
675 0 : BigInt& BigInt::operator*=( const BigInt& rVal )
676 : {
677 0 : if ( !bIsBig && !rVal.bIsBig
678 0 : && nVal <= MY_MAXSHORT && rVal.nVal <= MY_MAXSHORT
679 0 : && nVal >= MY_MINSHORT && rVal.nVal >= MY_MINSHORT )
680 : // TODO: not optimal !!! W.P.
681 : { // No overflows may occur here
682 0 : nVal *= rVal.nVal;
683 : }
684 : else
685 : {
686 0 : BigInt aTmp1, aTmp2;
687 0 : aTmp1.MakeBigInt( rVal );
688 0 : aTmp2.MakeBigInt( *this );
689 0 : aTmp1.MultLong(aTmp2, *this);
690 0 : Normalize();
691 : }
692 0 : return *this;
693 : }
694 :
695 0 : BigInt& BigInt::operator/=( const BigInt& rVal )
696 : {
697 0 : if ( !rVal.bIsBig )
698 : {
699 0 : if ( rVal.nVal == 0 )
700 : {
701 : OSL_FAIL( "BigInt::operator/ --> divide by zero" );
702 0 : return *this;
703 : }
704 :
705 0 : if ( !bIsBig )
706 : {
707 : // No overflows may occur here
708 0 : nVal /= rVal.nVal;
709 0 : return *this;
710 : }
711 :
712 0 : if ( rVal.nVal == 1 )
713 0 : return *this;
714 :
715 0 : if ( rVal.nVal == -1 )
716 : {
717 0 : bIsNeg = !bIsNeg;
718 0 : return *this;
719 : }
720 :
721 0 : if ( rVal.nVal <= (long)0xFFFF && rVal.nVal >= -(long)0xFFFF )
722 : {
723 : // Divide BigInt with an sal_uInt16
724 : sal_uInt16 nTmp;
725 0 : if ( rVal.nVal < 0 )
726 : {
727 0 : nTmp = (sal_uInt16) -rVal.nVal;
728 0 : bIsNeg = !bIsNeg;
729 : }
730 : else
731 0 : nTmp = (sal_uInt16) rVal.nVal;
732 :
733 0 : Div( nTmp, nTmp );
734 0 : Normalize();
735 0 : return *this;
736 : }
737 : }
738 :
739 0 : if ( ABS_IsLess( rVal ) )
740 : {
741 0 : *this = BigInt( (long)0 );
742 0 : return *this;
743 : }
744 :
745 : // Divide BigInt with BigInt
746 0 : BigInt aTmp1, aTmp2;
747 0 : aTmp1.MakeBigInt( *this );
748 0 : aTmp2.MakeBigInt( rVal );
749 0 : aTmp1.DivLong(aTmp2, *this);
750 0 : Normalize();
751 0 : return *this;
752 : }
753 :
754 0 : BigInt& BigInt::operator%=( const BigInt& rVal )
755 : {
756 0 : if ( !rVal.bIsBig )
757 : {
758 0 : if ( rVal.nVal == 0 )
759 : {
760 : OSL_FAIL( "BigInt::operator/ --> divide by zero" );
761 0 : return *this;
762 : }
763 :
764 0 : if ( !bIsBig )
765 : {
766 : // No overflows may occur here
767 0 : nVal %= rVal.nVal;
768 0 : return *this;
769 : }
770 :
771 0 : if ( rVal.nVal <= (long)0xFFFF && rVal.nVal >= -(long)0xFFFF )
772 : {
773 : // Divide Bigint by short
774 : sal_uInt16 nTmp;
775 0 : if ( rVal.nVal < 0 )
776 : {
777 0 : nTmp = (sal_uInt16) -rVal.nVal;
778 0 : bIsNeg = !bIsNeg;
779 : }
780 : else
781 0 : nTmp = (sal_uInt16) rVal.nVal;
782 :
783 0 : Div( nTmp, nTmp );
784 0 : *this = BigInt( (long)nTmp );
785 0 : return *this;
786 : }
787 : }
788 :
789 0 : if ( ABS_IsLess( rVal ) )
790 0 : return *this;
791 :
792 : // Divide BigInt with BigInt
793 0 : BigInt aTmp1, aTmp2;
794 0 : aTmp1.MakeBigInt( *this );
795 0 : aTmp2.MakeBigInt( rVal );
796 0 : aTmp1.ModLong(aTmp2, *this);
797 0 : Normalize();
798 0 : return *this;
799 : }
800 :
801 0 : bool operator==( const BigInt& rVal1, const BigInt& rVal2 )
802 : {
803 0 : if ( rVal1.bIsBig || rVal2.bIsBig )
804 : {
805 0 : BigInt nA, nB;
806 0 : nA.MakeBigInt( rVal1 );
807 0 : nB.MakeBigInt( rVal2 );
808 0 : if ( nA.bIsNeg == nB.bIsNeg )
809 : {
810 0 : if ( nA.nLen == nB.nLen )
811 : {
812 : int i;
813 0 : for ( i = nA.nLen - 1; i > 0 && nA.nNum[i] == nB.nNum[i]; i-- )
814 : {
815 : }
816 :
817 0 : return nA.nNum[i] == nB.nNum[i];
818 : }
819 0 : return false;
820 : }
821 0 : return false;
822 : }
823 0 : return rVal1.nVal == rVal2.nVal;
824 : }
825 :
826 0 : bool operator<( const BigInt& rVal1, const BigInt& rVal2 )
827 : {
828 0 : if ( rVal1.bIsBig || rVal2.bIsBig )
829 : {
830 0 : BigInt nA, nB;
831 0 : nA.MakeBigInt( rVal1 );
832 0 : nB.MakeBigInt( rVal2 );
833 0 : if ( nA.bIsNeg == nB.bIsNeg )
834 : {
835 0 : if ( nA.nLen == nB.nLen )
836 : {
837 : int i;
838 0 : for ( i = nA.nLen - 1; i > 0 && nA.nNum[i] == nB.nNum[i]; i-- )
839 : {
840 : }
841 :
842 0 : if ( nA.bIsNeg )
843 0 : return nA.nNum[i] > nB.nNum[i];
844 : else
845 0 : return nA.nNum[i] < nB.nNum[i];
846 : }
847 0 : if ( nA.bIsNeg )
848 0 : return nA.nLen > nB.nLen;
849 : else
850 0 : return nA.nLen < nB.nLen;
851 : }
852 0 : return !nB.bIsNeg;
853 : }
854 0 : return rVal1.nVal < rVal2.nVal;
855 : }
856 :
857 0 : bool operator >(const BigInt& rVal1, const BigInt& rVal2 )
858 : {
859 0 : if ( rVal1.bIsBig || rVal2.bIsBig )
860 : {
861 0 : BigInt nA, nB;
862 0 : nA.MakeBigInt( rVal1 );
863 0 : nB.MakeBigInt( rVal2 );
864 0 : if ( nA.bIsNeg == nB.bIsNeg )
865 : {
866 0 : if ( nA.nLen == nB.nLen )
867 : {
868 : int i;
869 0 : for ( i = nA.nLen - 1; i > 0 && nA.nNum[i] == nB.nNum[i]; i-- )
870 : {
871 : }
872 :
873 0 : if ( nA.bIsNeg )
874 0 : return nA.nNum[i] < nB.nNum[i];
875 : else
876 0 : return nA.nNum[i] > nB.nNum[i];
877 : }
878 0 : if ( nA.bIsNeg )
879 0 : return nA.nLen < nB.nLen;
880 : else
881 0 : return nA.nLen > nB.nLen;
882 : }
883 0 : return !nA.bIsNeg;
884 : }
885 :
886 0 : return rVal1.nVal > rVal2.nVal;
887 : }
888 :
889 : /* vim:set shiftwidth=4 softtabstop=4 expandtab: */
|
