Remove more floating point constants

Author: jecisc

src/Math-Core-Process/PMIterativeProcess.class.st

@@ -99,10 +99,13 @@ PMIterativeProcess >> iterations [
 
 { #category : #private }
 PMIterativeProcess >> limitedSmallValue: aNumber [
-	"Private - prevent aNumber from being smaller in absolute value than a small number."
-	^aNumber abs < PMFloatingPointMachine new smallNumber
-			ifTrue: [ PMFloatingPointMachine new smallNumber]
-			ifFalse:[ aNumber]
+	"Private - prevent aNumber from being smaller in absolute value than a small."
+
+	| smallestAccepted |
+	smallestAccepted := Float fmin sqrt.
+	^ aNumber abs < smallestAccepted
+		  ifTrue: [ smallestAccepted ]
+		  ifFalse: [ aNumber ]
 ]
 
 { #category : #initialization }

src/Math-Distributions/PMFisherTippettDistribution.class.st

@@ -60,7 +60,7 @@ PMFisherTippettDistribution >> distributionValue: aNumber [
 		 the receiver with a value lower than or equal to aNumber."
 	| arg |
 	arg := ( aNumber - alpha) / beta.
-	arg := arg < PMFloatingPointMachine new largestExponentArgument negated
+	arg := arg < Float fmax ln negated
 					ifTrue: [ ^0]
 					ifFalse:[arg negated exp].
 	^arg negated exp
@@ -126,14 +126,17 @@ PMFisherTippettDistribution >> standardDeviation [
 
 { #category : #information }
 PMFisherTippettDistribution >> value: aNumber [
-		"Answers the probability that a random variable distributed according to the receiver
+	"Answers the probability that a random variable distributed according to the receiver
 		 gives a value between aNumber and aNumber + espilon (infinitesimal interval)."
+
 	| arg |
-	arg := ( aNumber - alpha) / beta.
-	arg := arg > PMFloatingPointMachine new largestExponentArgument ifTrue: [ ^0]
-						ifFalse:[arg negated exp + arg].
-	^arg > PMFloatingPointMachine new largestExponentArgument ifTrue: [ 0]
-						ifFalse:[ arg negated exp / beta]
+	arg := aNumber - alpha / beta.
+	arg := arg > Float fmax ln
+		       ifTrue: [ ^ 0 ]
+		       ifFalse: [ arg negated exp + arg ].
+	^ arg > Float fmax ln
+		  ifTrue: [ 0 ]
+		  ifFalse: [ arg negated exp / beta ]
 ]
 
 { #category : #information }

src/Math-Helpers/PMFloatingPointMachine.class.st

@@ -4,13 +4,8 @@ A `PMFloatingPointMachine` represents the numerical precision of this system.
 ##Instance Variables
 
 - `defaultNumericalPrecision`	The relative numerical precision that can be expected for a general numerical computation. One should consider to numbers a and b equal if the relative difference between them is less than the default machine precision,
-- `largestExponentArgument` Natural logarithm of largest number,
-- `largestNumber` The largest positive number that can be represented in the machine,
 - `machinePrecision` $r^{-(n+1)}$, with the largest n such that $(1 + r^{-n}) - 1$ != 0,
-- `negativeMachinePrecision` $r^{-(n+1)}$, with the largest n such that $(1 - r^{-n}) - 1$ != 0,
-- `radix` The radix of the floating point representation. This is often 2,
-- `smallNumber` A number that can be added to some value without noticeably changing the result of the computation,
-- `smallestNumber` The smallest positive number different from 0.
+- `negativeMachinePrecision` $r^{-(n+1)}$, with the largest n such that $(1 - r^{-n}) - 1$ != 0
 
 This class is detailed in Object Oriented Implementation of Numerical Methods, Section 1.4.1 and 1.4.2.
 
@@ -21,9 +16,7 @@ Class {
 	#instVars : [
 		'defaultNumericalPrecision',
 		'machinePrecision',
-		'negativeMachinePrecision',
-		'smallNumber',
-		'largestExponentArgument'
+		'negativeMachinePrecision'
 	],
 	#classVars : [
 		'UniqueInstance'
@@ -78,19 +71,6 @@ PMFloatingPointMachine >> defaultNumericalPrecision [
 	^ defaultNumericalPrecision
 ]
 
-{ #category : #information }
-PMFloatingPointMachine >> largestExponentArgument [
-
-	largestExponentArgument ifNil: [ largestExponentArgument := Float fmax ln ].
-	^ largestExponentArgument
-]
-
-{ #category : #information }
-PMFloatingPointMachine >> largestNumber [
-
-	^ Float fmax
-]
-
 { #category : #information }
 PMFloatingPointMachine >> machinePrecision [
 
@@ -104,10 +84,3 @@ PMFloatingPointMachine >> negativeMachinePrecision [
 	negativeMachinePrecision ifNil: [ self computeNegativeMachinePrecision ].
 	^ negativeMachinePrecision
 ]
-
-{ #category : #information }
-PMFloatingPointMachine >> smallNumber [
-
-	smallNumber ifNil: [ smallNumber := Float fmin sqrt ].
-	^ smallNumber
-]

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

@@ -42,21 +42,6 @@ PMFloatingPointMachineTestCase >> testMachinePrecisionLargestNumberIsLargest [
 	self assert: Float fmax * (1 + mach defaultNumericalPrecision) equals: Float infinity
 ]
 
-{ #category : #precision }
-PMFloatingPointMachineTestCase >> testMachinePrecisionSmallNumberLargerThanSmallestNumber [
-
-	| mach |
-	mach := PMFloatingPointMachine new.
-	self assert: Float fmin < mach smallNumber
-]
-
-{ #category : #precision }
-PMFloatingPointMachineTestCase >> testMachinePrecisionSmallNumberNotZero [
-	| mach |
-	mach := PMFloatingPointMachine new.
-	self assert: mach smallNumber > 0.0
-]
-
 { #category : #precision }
 PMFloatingPointMachineTestCase >> testUniqueInstance [