rlab/example/matrixEx.js

var M = require("../rlab").M;

var c = console;
var v = [1,2,3];
c.log("v.sin()=", v.sin());
c.log("v.norm2()=", v.norm2());
c.log("v.norm2Squared()=", v.norm2Squared());

var A = [[1,2,3],[4,5,6],[7,3,9]];
var AiA = A.inv().dot(A);
c.log("AiA=\n", AiA.strM());
c.log("AiA.tr()=\n", AiA.tr().strM());
c.log("A=\n", A.str());
c.log("A.mul(0.1)=\n", A.mul(0.1).strM());
c.log("A.row(1)=", A.row(1));
c.log("A.col(1)=", A.col(1));
c.log("A.sumM()=", A.sumM());
c.log("A.rowSum(2)=", A.rowSum(2));
c.log("A.colSum(2)=", A.colSum(2));
c.log("A.mean(row)=", A.rowMean().str());
c.log("A.mean(col)=", A.colMean().str());

var D = M.diag(v);
c.log("D=", D);
/*
// console.log("M.eig=", M.eig);
var E = M.eig(A);
c.log("E=A.eig", E);
c.log("E.lambda=", E.lambda);
c.log("E.E", E.E);
// E*[lambda]*E-1
c.log("E=A.eig(); E*[lambda]*E-1.x-A =", E.E.dot(R.T.diag(E.lambda)).dot(E.E.inv()).x.sub(A).abs().sumM()<0.001);
*/

var Eλ = M.eigR(A);
var E = Eλ.E, λ=Eλ.lambda;
c.log("E*[λ]*E-1=", E.dot(λ.diag()).dot(E.inv()).strM());

// c.log("I=", M.identity(3).strM());