Added matrix utils

Author: perak

lib/quantum-circuit.js

@@ -2850,9 +2850,16 @@ QuantumCircuit.prototype.evalMathString = function(str, vars) {
 QuantumCircuit.prototype.parseMatrix = function(M, vars) {
 	var self = this;
 
+	var A = null;
+	if(typeof M == "string") {
+		A = self.evalMathString(M, vars);
+	} else {
+		A = M;
+	}
+
 	var res = [];
 
-	M.map(function(row, rowIndex) {
+	A.map(function(row, rowIndex) {
 		if(Array.isArray(row)) {
 			res.push([]);
 			row.map(function(val, colIndex) {
@@ -2982,6 +2989,95 @@ QuantumCircuit.prototype.stringifyMatrix = function(M, options) {
 };
 
 
+QuantumCircuit.prototype.matrixRe = function(U) {
+	var res = [];
+	for(var r = 0; r < U.length; r++) {
+		var row = U[r];
+		if(Array.isArray(row)) {
+			var newRow = [];
+			for(var c = 0; c < row.length; c++) {
+				if(typeof row[c] == "object") {
+					newRow.push(row[c].re);
+				} else {
+					newRow.push(row[c]);
+				}
+			}
+			res.push(newRow);
+		} else {
+			if(typeof row == "object") {
+				newRow.push(row.re);
+			} else {
+				newRow.push(row);
+			}
+		}
+	}
+
+	return res;
+};
+
+QuantumCircuit.prototype.matrixIm = function(U) {
+	var res = [];
+	for(var r = 0; r < U.length; r++) {
+		var row = U[r];
+		if(Array.isArray(row)) {
+			var newRow = [];
+			for(var c = 0; c < row.length; c++) {
+				if(typeof row[c] == "object") {
+					newRow.push(row[c].im);
+				} else {
+					newRow.push(0);
+				}
+			}
+			res.push(newRow);
+		} else {
+			if(typeof row == "object") {
+				newRow.push(row.im);
+			} else {
+				newRow.push(0);
+			}
+		}
+	}
+
+	return res;
+};
+
+QuantumCircuit.prototype.matrixAbs = function(U) {
+	var res = [];
+	for(var r = 0; r < U.length; r++) {
+		var row = U[r];
+		if(Array.isArray(row)) {
+			var newRow = [];
+			for(var c = 0; c < row.length; c++) {
+				newRow.push(math.abs(row[c]));
+			}
+			res.push(newRow);
+		} else {
+			res.push(math.abs(row));
+		}
+	}
+
+	return res;
+};
+
+QuantumCircuit.prototype.matrixArg = function(U) {
+	var res = [];
+	for(var r = 0; r < U.length; r++) {
+		var row = U[r];
+		if(Array.isArray(row)) {
+			var newRow = [];
+			for(var c = 0; c < row.length; c++) {
+				newRow.push(math.arg(row[c]));
+			}
+			res.push(newRow);
+		} else {
+			res.push(math.arg(row));
+		}
+	}
+
+	return res;
+};
+
+
 //
 // combineList is array of objects:
 //