summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat (limited to '')
-rw-r--r--src/number.cpp43
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
499bool Number::isZero() const 508bool Number::isZero() const