testHouseholder added

Author: mabdi

src/Math-Tests-Core/PMVectorTest.class.st

@@ -69,6 +69,23 @@ PMVectorTest >> testGreaterThan [
 	self assert: vec equals: vecCopy asPMVector.
 ]
 
+{ #category : #tests }
+PMVectorTest >> testHouseholder [
+	| u w |
+	u := #(-1 0 1) asPMVector. "`x <= 0` when x = -1"
+	w := u householder.
+	self
+		assert: w class equals: Array;
+		assert: (w at: 1) equals: 1.7071067811865475;
+		assert: (w at: 2) asArray equals: #(1.0 -0.0 -0.4142135623730951).
+	u := #(1.00001 2.00007) asPMVector. "`x <= 0` when x = 1.00001"
+	w := u householder.
+	self
+		assert: w class equals: Array;
+		assert: (w at: 1) equals: 0.5527953485259909;
+		assert: (w at: 2) asArray equals: #(1.0 -1.6180158992689828).
+]
+
 { #category : #tests }
 PMVectorTest >> testLessThan [
 	| vec vecCopy |
@@ -79,7 +96,7 @@ PMVectorTest >> testLessThan [
 	self assert: vec equals: vecCopy asPMVector.
 ]
 
-{ #category : #tests } 
+{ #category : #tests }
 PMVectorTest >> testMatrixConversionWithBothDims [
 	| vect result expected |
 	vect := #(1 0.5 0.2 3 1 -1 7 3 2 12 13 3) asPMVector .