Use the smallest number provided by Pharo

Author: jecisc

src/Math-Helpers/PMFloatingPointMachine.class.st

@@ -22,7 +22,6 @@ Class {
 		'defaultNumericalPrecision',
 		'machinePrecision',
 		'negativeMachinePrecision',
-		'smallestNumber',
 		'smallNumber',
 		'largestExponentArgument'
 	],
@@ -72,25 +71,6 @@ PMFloatingPointMachine >> computeNegativeMachinePrecision [
 	tmp - one = zero ] whileFalse: [ negativeMachinePrecision := negativeMachinePrecision * inverseRadix ]
 ]
 
-{ #category : #information }
-PMFloatingPointMachine >> computeSmallestNumber [
-
-	| one floatingRadix inverseRadix fullMantissaNumber |
-	one := 1 asFloat.
-	floatingRadix := Float radix asFloat.
-	inverseRadix := one / floatingRadix.
-	fullMantissaNumber := one - (floatingRadix * self negativeMachinePrecision).
-	smallestNumber := fullMantissaNumber.
-	[
-	[
-	fullMantissaNumber := fullMantissaNumber * inverseRadix.
-	fullMantissaNumber = 0.0 ifTrue: [ Error signal ].
-	smallestNumber := fullMantissaNumber.
-	true ] whileTrue: [  ] ]
-		on: Error
-		do: [ :signal | signal return: nil ]
-]
-
 { #category : #information }
 PMFloatingPointMachine >> defaultNumericalPrecision [
 
@@ -125,40 +105,9 @@ PMFloatingPointMachine >> negativeMachinePrecision [
 	^ negativeMachinePrecision
 ]
 
-{ #category : #display }
-PMFloatingPointMachine >> showParameters [
-
-	Transcript
-		cr;
-		cr;
-		nextPutAll: 'Floating-point machine parameters';
-		cr;
-		nextPutAll: '---------------------------------'.
-	Transcript
-		cr;
-		nextPutAll: 'Machine precision: '.
-	self machinePrecision printOn: Transcript.
-	Transcript
-		cr;
-		nextPutAll: 'Negative machine precision: '.
-	self negativeMachinePrecision printOn: Transcript.
-	Transcript
-		cr;
-		nextPutAll: 'Smallest number: '.
-	self smallestNumber printOn: Transcript.
-	Transcript flush
-]
-
 { #category : #information }
 PMFloatingPointMachine >> smallNumber [
 
-	smallNumber ifNil: [ smallNumber := self smallestNumber sqrt ].
+	smallNumber ifNil: [ smallNumber := Float fmin sqrt ].
 	^ smallNumber
 ]
-
-{ #category : #information }
-PMFloatingPointMachine >> smallestNumber [
-
-	smallestNumber ifNil: [ self computeSmallestNumber ].
-	^ smallestNumber
-]

src/Math-Tests-Numerical/PMFloatingPointMachineTestCase.class.st

@@ -44,9 +44,10 @@ PMFloatingPointMachineTestCase >> testMachinePrecisionLargestNumberIsLargest [
 
 { #category : #precision }
 PMFloatingPointMachineTestCase >> testMachinePrecisionSmallNumberLargerThanSmallestNumber [
+
 	| mach |
 	mach := PMFloatingPointMachine new.
-	self assert: mach smallestNumber < mach smallNumber
+	self assert: Float fmin < mach smallNumber
 ]
 
 { #category : #precision }
@@ -56,13 +57,6 @@ PMFloatingPointMachineTestCase >> testMachinePrecisionSmallNumberNotZero [
 	self assert: mach smallNumber > 0.0
 ]
 
-{ #category : #precision }
-PMFloatingPointMachineTestCase >> testMachinePrecisionSmallestNumberNotZero [
-	| mach |
-	mach := PMFloatingPointMachine new.
-	self assert: mach smallestNumber > 0.0
-]
-
 { #category : #precision }
 PMFloatingPointMachineTestCase >> testUniqueInstance [
 

src/Math-Tests-Numerical/PMNumericalMethodsTestCase.class.st

@@ -98,7 +98,7 @@ PMNumericalMethodsTestCase >> testErrorFunction [
 	self assert: Float fmax errorFunction > (1 - PMFloatingPointMachine new machinePrecision).
 	"add some code to require initialize to run"
 	PMErfApproximation reset.
-	self assert: Float fmax negated errorFunction < PMFloatingPointMachine new smallestNumber
+	self assert: Float fmax negated errorFunction < Float fmin
 ]
 
 { #category : #estimation }