Merge pull request #1 from gagern/realReciprocal

Fix reciprocal of a complex number with zero imaginary part

Author: orangeduck

documentation.html

@@ -1107,6 +1107,10 @@ <h1>Complex linear algebra</h1>
 OUT> {x:0.5588,y:-0.2353}
 IN> z.sub(w)
 OUT> {x:1, y:-4}
+IN> z.reciprocal()
+OUT> {x: 0.12, y: -0.16}
+IN> numeric.t(2, 0).reciprocal()
+OUT> {x: 0.5, y: 0}
 </pre>
 
 Complex vectors and matrices can also be handled:

src/documentation.html

@@ -1107,6 +1107,10 @@ <h1>Complex linear algebra</h1>
 OUT> {x:0.5588,y:-0.2353}
 IN> z.sub(w)
 OUT> {x:1, y:-4}
+IN> z.reciprocal()
+OUT> {x: 0.12, y: -0.16}
+IN> numeric.t(2, 0).reciprocal()
+OUT> {x: 0.5, y: 0}
 </pre>
 
 Complex vectors and matrices can also be handled:

src/numeric.js

@@ -1377,7 +1377,7 @@ numeric.T.prototype.reciprocal = function reciprocal() {
         var d = numeric.add(mul(this.x,this.x),mul(this.y,this.y));
         return new numeric.T(div(this.x,d),div(numeric.neg(this.y),d));
     }
-    return new T(div(1,this.x));
+    return new numeric.T(div(1,this.x), 0);
 }
 numeric.T.prototype.div = function div(y) {
     if(!(y instanceof numeric.T)) y = new numeric.T(y);