Merge branch 'master' into master

Author: olekscode

src/Math-Core-Distribution/PMNormalDistribution.class.st

@@ -8,9 +8,22 @@ Class {
 	#classVars : [
 		'NextRandom'
 	],
-	#category : 'Math-Core-Distribution'
+	#category : #'Math-Core-Distribution'
 }
 
+{ #category : #information }
+PMNormalDistribution class >> boxMullerTransform [
+|v1 v2 w y|
+			[ v1 := Number random * 2 - 1.
+			v2 := Number random * 2 - 1.
+			w := v1 squared + v2 squared.
+			w > 1 ] whileTrue.
+			y := (w ln * 2 negated / w) sqrt.
+			v1 := y * v1.
+			NextRandom := y * v2.
+			^v1.
+]
+
 { #category : #information }
 PMNormalDistribution class >> distributionName [
 
@@ -40,22 +53,16 @@ PMNormalDistribution class >> new: aNumber1 sigma: aNumber2 [
 
 { #category : #information }
 PMNormalDistribution class >> random [
-	"A random value generated from the standard normal distribution N(0,1). Can be used to generate values from any other univariate normal distribution N(mu, sigma): first we generate Z ~ N(0,1) and then, to get X ~ N(mu, sigma), we calculate X = sigma * Z + mu (see #random method)"
-	
-	| v1 v2 w y |
-	NextRandom isNil
-		ifTrue: [ [ v1 := Number random * 2 - 1.
-					v2 := Number random * 2 - 1.
-					w := v1 squared + v2 squared.
-					w > 1 ] whileTrue: [].
-				  y := ( ( w ln * 2 negated) / w) sqrt.
-				v1 := y * v1.
-				NextRandom := y * v2.
-				]
-		ifFalse:[ v1 :=NextRandom.
-				  NextRandom := nil.
-				].
-	^v1
+	"Answer a random number distributed according to a (0,1) normal distribution."
+
+	| v1  |
+	NextRandom
+		ifNil: [
+			
+			v1 := self boxMullerTransform ]
+		ifNotNil: [ v1 := NextRandom.
+			NextRandom := nil ].
+	^ v1
 ]
 
 { #category : #information }
@@ -79,10 +86,8 @@ PMNormalDistribution >> distributionValue: aNumber [
 
 { #category : #initialization }
 PMNormalDistribution >> initialize: aNumber1 sigma: aNumber2 [
-
 	mu := aNumber1.
-	sigma := aNumber2.
-	^self
+	sigma := aNumber2
 ]
 
 { #category : #information }
@@ -99,8 +104,9 @@ PMNormalDistribution >> parameters [
 
 { #category : #information }
 PMNormalDistribution >> random [
-	"Answer a random number distributed accroding to the receiver."
-	^self class random * sigma + mu
+	"Answer a random number distributed according to the receiver."
+
+	^ self class random * sigma + mu
 ]
 
 { #category : #information }

src/Math-Numerical/Number.extension.st

@@ -70,15 +70,16 @@ Number >> productWithVector: aVector [
 ]
 
 { #category : #'*Math-Numerical' }
-Number >> random [
-		"Answers a random number distributed between 0 and the receiver."
-	^self class random * self
+Number class >> random [
+	"Answers a floating random number between 0 and 1 excluded"
+
+	^ PMMitchellMooreGenerator new floatValue
 ]
 
 { #category : #'*Math-Numerical' }
-Number class >> random [
-		"Answers a random number between 0 and the receiver."
-	^ PMMitchellMooreGenerator new floatValue
+Number >> random [
+		"Answers a random number distributed between 0 and the receiver."
+	^self class random * self
 ]
 
 { #category : #'*Math-Numerical' }

src/Math-Numerical/PMCauchyDistribution.class.st

@@ -59,10 +59,10 @@ PMCauchyDistribution >> distributionValue: aNumber [
 
 { #category : #initialization }
 PMCauchyDistribution >> initialize: aNumber1 scale: aNumber2 [
-	"Private - Initialize the parameters of the receiver. "
+	"Initialize the parameters of the receiver. "
+
 	mu := aNumber1.
-	beta := aNumber2.
-	^self
+	beta := aNumber2
 ]
 
 { #category : #information }

src/Math-Numerical/PMExponentialDistribution.class.st

@@ -66,8 +66,7 @@ PMExponentialDistribution >> initialize: aNumber [
 
 	aNumber > 0
 		ifFalse: [ self error: 'Illegal distribution parameters' ].
-	beta := aNumber.
-	^ self
+	beta := aNumber
 ]
 
 { #category : #information }

src/Math-Numerical/PMFisherTippettDistribution.class.st

@@ -68,12 +68,12 @@ PMFisherTippettDistribution >> distributionValue: aNumber [
 
 { #category : #initialization }
 PMFisherTippettDistribution >> initialize: aNumber1 scale: aNumber2 [
-		"Private - Initialize the parameters of the receiver."
-	aNumber2 > 0
-		ifFalse: [ self error: 'Illegal distribution parameters'].
+	"Initialize the parameters of the receiver."
+
+	aNumber2 <= 0
+		ifTrue: [ self error: 'Illegal distribution parameters' ].
 	alpha := aNumber1.
-	beta := aNumber2.
-	^self
+	beta := aNumber2
 ]
 
 { #category : #information }
@@ -104,11 +104,12 @@ PMFisherTippettDistribution >> parameters [
 
 { #category : #information }
 PMFisherTippettDistribution >> random [
-		"Answer a random number distributed according to the receiver."
+	"Answer a random number distributed according to the receiver."
+
 	| t |
 	[ t := flatGenerator floatValue ln negated.
-	  t > 0] whileFalse: [].
-	^t ln negated * beta + alpha
+	t > 0 ] whileFalse.
+	^ t ln negated * beta + alpha
 ]
 
 { #category : #information }