diff options
| -rw-r--r-- | src/number.cpp | 43 | 
1 files changed, 26 insertions, 17 deletions
| diff --git a/src/number.cpp b/src/number.cpp index be9b8be..81aefbe 100644 --- a/src/number.cpp +++ b/src/number.cpp | |||
| @@ -439,7 +439,7 @@ void Number::divide( const Number &rhs, Number &q, Number &r ) const | |||
| 439 | nNum.aInt.append( digit( j ) ); | 439 | nNum.aInt.append( digit( j ) ); | 
| 440 | } | 440 | } | 
| 441 | // Bu::println("New numerator: %1").arg( nNum ); | 441 | // Bu::println("New numerator: %1").arg( nNum ); | 
| 442 | while( iAnchor > -iScale && !nNum.isZero()) // division loop | 442 | while( iAnchor > -iScale && (!nNum.isZero() || iAnchor > 0)) // division loop | 
| 443 | { | 443 | { | 
| 444 | // As long as the numerator is smaller than the divisor we keep | 444 | // As long as the numerator is smaller than the divisor we keep | 
| 445 | // including new digits. | 445 | // including new digits. | 
| @@ -449,10 +449,10 @@ void Number::divide( const Number &rhs, Number &q, Number &r ) const | |||
| 449 | nNum.aInt.insert( 0, digit(iAnchor-iNumShift) ); | 449 | nNum.aInt.insert( 0, digit(iAnchor-iNumShift) ); | 
| 450 | // Bu::println("Numerator too small, new numerator: %1").arg( nNum ); | 450 | // Bu::println("Numerator too small, new numerator: %1").arg( nNum ); | 
| 451 | } | 451 | } | 
| 452 | if( iAnchor < -iScale ) | ||
| 453 | break; | ||
| 454 | 452 | ||
| 455 | // Bu::println(" New anchor: %1").arg( iAnchor ); | 453 | // Bu::println(" New anchor: %1 vs iScale = %2").arg( iAnchor ).arg( iScale ); | 
| 454 | if( iAnchor < 0 && -1-iAnchor >= iScale ) | ||
| 455 | break; | ||
| 456 | 456 | ||
| 457 | // How many times does the divisor fit into the numerator | 457 | // How many times does the divisor fit into the numerator | 
| 458 | int iCount = 0; | 458 | int iCount = 0; | 
| @@ -477,23 +477,32 @@ void Number::divide( const Number &rhs, Number &q, Number &r ) const | |||
| 477 | q.aFrac.set( -1-iAnchor, iCount ); | 477 | q.aFrac.set( -1-iAnchor, iCount ); | 
| 478 | } | 478 | } | 
| 479 | } | 479 | } | 
| 480 | for( int j = -iAnchor; j < nNum.aInt.getSize(); j++ ) | 480 | // if( iAnchor <= 0 ) | 
| 481 | r.aInt.append( nNum.aInt[j] ); | ||
| 482 | int jtop = 0; | ||
| 483 | for(; nNum.aInt[jtop] == 0; jtop++ ) { } | ||
| 484 | // Bu::println("Computing remainder fraction from %1 to %2") | ||
| 485 | // .arg( -iAnchor-1 ).arg( jtop ); | ||
| 486 | r.aFrac.clear(); | ||
| 487 | for( int j = -iAnchor-1, k = 0; j >= jtop && k < r.iScale; j--, k++ ) | ||
| 488 | { | 481 | { | 
| 489 | if( j < nNum.aInt.getSize() ) | 482 | for( int j = -iAnchor; j < nNum.aInt.getSize(); j++ ) | 
| 483 | r.aInt.append( nNum.aInt[j] ); | ||
| 484 | int jtop = 0; | ||
| 485 | for(; nNum.aInt[jtop] == 0; jtop++ ) { } | ||
| 486 | // Bu::println("Computing remainder fraction from %1 to %2") | ||
| 487 | // .arg( -iAnchor-1 ).arg( jtop ); | ||
| 488 | r.aFrac.clear(); | ||
| 489 | for( int j = -iAnchor-1, k = 0; j >= jtop && k < r.iScale; j--, k++ ) | ||
| 490 | { | 490 | { | 
| 491 | r.aFrac.set( k, nNum.aInt[j] ); | 491 | if( j < nNum.aInt.getSize() ) | 
| 492 | // Bu::println("setting %1 to %2").arg( k ).arg( nNum.aInt[j] ); | 492 | { | 
| 493 | r.aFrac.set( k, nNum.aInt[j] ); | ||
| 494 | // Bu::println("setting %1 to %2").arg( k ).arg( nNum.aInt[j] ); | ||
| 495 | } | ||
| 493 | } | 496 | } | 
| 494 | } | 497 | } | 
| 495 | // Bu::println("Final numerator? %1").arg( nNum ); | 498 | /* else | 
| 496 | // Bu::println("Remainder? %1").arg( r ); | 499 | { | 
| 500 | Bu::println("iAnchor: %1").arg( iAnchor ); | ||
| 501 | } | ||
| 502 | Bu::println("Quotient: %1").arg( q ); | ||
| 503 | Bu::println("Final numerator? %1").arg( nNum ); | ||
| 504 | Bu::println("Remainder? %1").arg( r ); | ||
| 505 | */ | ||
| 497 | } | 506 | } | 
| 498 | 507 | ||
| 499 | bool Number::isZero() const | 508 | bool Number::isZero() const | 
