New file |
0,0 → 1,294 |
if(!dojo._hasResource["dojox.math.matrix"]){ //_hasResource checks added by build. Do not use _hasResource directly in your code. |
dojo._hasResource["dojox.math.matrix"] = true; |
dojo.provide("dojox.math.matrix"); |
|
dojo.mixin(dojox.math.matrix, { |
iDF:0, |
ALMOST_ZERO: 1e-10, |
multiply: function(/* Array */a, /* Array */b){ |
// summary |
// Multiply matrix a by matrix b. |
var ay=a.length, ax=a[0].length, by=b.length, bx=b[0].length; |
if(ax!=by){ |
console.warn("Can't multiply matricies of sizes " + ax + "," + ay + " and " + bx + "," + by); |
return [[0]]; |
} |
var c=[]; |
for (var k=0; k<ay; k++) { |
c[k]=[]; |
for(var i=0; i<bx; i++){ |
c[k][i]=0; |
for(var m=0; m<ax; m++){ |
c[k][i]+=a[k][m]*b[m][i]; |
} |
} |
} |
return c; // Array |
}, |
product: function(/* Array... */){ |
// summary |
// Return the product of N matrices |
if (arguments.length==0){ |
console.warn("can't multiply 0 matrices!"); |
return 1; |
} |
var m=arguments[0]; |
for(var i=1; i<arguments.length; i++){ |
m=this.multiply(m, arguments[i]); |
} |
return m; // Array |
}, |
sum: function(/* Array... */){ |
// summary |
// Return the sum of N matrices |
if(arguments.length==0){ |
console.warn("can't sum 0 matrices!"); |
return 0; // Number |
} |
var m=this.copy(arguments[0]); |
var rows=m.length; |
if(rows==0){ |
console.warn("can't deal with matrices of 0 rows!"); |
return 0; |
} |
var cols=m[0].length; |
if(cols==0){ |
console.warn("can't deal with matrices of 0 cols!"); |
return 0; |
} |
for(var i=1; i<arguments.length; ++i){ |
var arg=arguments[i]; |
if(arg.length!=rows || arg[0].length!=cols){ |
console.warn("can't add matrices of different dimensions: first dimensions were " + rows + "x" + cols + ", current dimensions are " + arg.length + "x" + arg[0].length); |
return 0; |
} |
for(var r=0; r<rows; r++) { |
for(var c=0; c<cols; c++) { |
m[r][c]+=arg[r][c]; |
} |
} |
} |
return m; // Array |
}, |
inverse: function(/* Array */a){ |
// summary |
// Return the inversion of the passed matrix |
if(a.length==1 && a[0].length==1){ |
return [[1/a[0][0]]]; // Array |
} |
var tms=a.length, m=this.create(tms, tms), mm=this.adjoint(a), det=this.determinant(a), dd=0; |
if(det==0){ |
console.warn("Determinant Equals 0, Not Invertible."); |
return [[0]]; |
}else{ |
dd=1/det; |
} |
for(var i=0; i<tms; i++) { |
for (var j=0; j<tms; j++) { |
m[i][j]=dd*mm[i][j]; |
} |
} |
return m; // Array |
}, |
determinant: function(/* Array */a){ |
// summary |
// Calculate the determinant of the passed square matrix. |
if(a.length!=a[0].length){ |
console.warn("Can't calculate the determinant of a non-squre matrix!"); |
return 0; |
} |
var tms=a.length, det=1, b=this.upperTriangle(a); |
for (var i=0; i<tms; i++){ |
var bii=b[i][i]; |
if (Math.abs(bii)<this.ALMOST_ZERO) { |
return 0; // Number |
} |
det*=bii; |
} |
det*=this.iDF; |
return det; // Number |
}, |
upperTriangle: function(/* Array */m){ |
// Summary |
// Find the upper triangle of the passed matrix and return it. |
m=this.copy(m); |
var f1=0, temp=0, tms=m.length, v=1; |
this.iDF=1; |
for(var col=0; col<tms-1; col++){ |
if(typeof m[col][col]!="number") { |
console.warn("non-numeric entry found in a numeric matrix: m[" + col + "][" + col + "]=" + m[col][col]); |
} |
v=1; |
var stop_loop=0; |
while((m[col][col] == 0) && !stop_loop){ |
if (col+v>=tms){ |
this.iDF=0; |
stop_loop=1; |
}else{ |
for(var r=0; r<tms; r++){ |
temp=m[col][r]; |
m[col][r]=m[col+v][r]; |
m[col+v][r]=temp; |
} |
v++; |
this.iDF*=-1; |
} |
} |
for(var row=col+1; row<tms; row++){ |
if(typeof m[row][col]!="number"){ |
console.warn("non-numeric entry found in a numeric matrix: m[" + row + "][" + col + "]=" + m[row][col]); |
} |
if(typeof m[col][row]!="number"){ |
console.warn("non-numeric entry found in a numeric matrix: m[" + col + "][" + row + "]=" + m[col][row]); |
} |
if(m[col][col]!=0){ |
var f1=(-1)* m[row][col]/m[col][col]; |
for (var i=col; i<tms; i++){ |
m[row][i]=f1*m[col][i]+m[row][i]; |
} |
} |
} |
} |
return m; // Array |
}, |
create: function(/* Number */a, /* Number */b, /* Number? */value){ |
// summary |
// Create a new matrix with rows a and cols b, and pre-populate with value. |
value=value||0; |
var m=[]; |
for (var i=0; i<b; i++){ |
m[i]=[]; |
for(var j=0; j<a; j++) { |
m[i][j]=value; |
} |
} |
return m; // Array |
}, |
ones: function(/* Number */a, /* Number */b){ |
// summary |
// Create a matrix pre-populated with ones |
return this.create(a, b, 1); // Array |
}, |
zeros: function(/* Number */a, /* Number */b){ |
// summary |
// Create a matrix pre-populated with zeros |
return this.create(a, b); // Array |
}, |
identity: function(/* Number */size, /* Number? */scale){ |
// summary |
// Create an identity matrix based on the size and scale. |
scale=scale||1; |
var m=[]; |
for(var i=0; i<size; i++){ |
m[i]=[]; |
for(var j=0; j<size; j++){ |
m[i][j]=(i==j?scale:0); |
} |
} |
return m; // Array |
}, |
adjoint: function(/* Array */a){ |
// summary |
// Find the adjoint of the passed matrix |
var tms=a.length; |
if(tms<=1){ |
console.warn("Can't find the adjoint of a matrix with a dimension less than 2"); |
return [[0]]; |
} |
if(a.length!=a[0].length){ |
console.warn("Can't find the adjoint of a non-square matrix"); |
return [[0]]; |
} |
var m=this.create(tms, tms), ap=this.create(tms-1, tms-1); |
var ii=0, jj=0, ia=0, ja=0, det=0; |
for(var i=0; i<tms; i++){ |
for (var j=0; j<tms; j++){ |
ia=0; |
for(ii=0; ii<tms; ii++){ |
if(ii==i){ |
continue; |
} |
ja = 0; |
for(jj=0; jj<tms; jj++){ |
if(jj==j){ |
continue; |
} |
ap[ia][ja] = a[ii][jj]; |
ja++; |
} |
ia++; |
} |
det=this.determinant(ap); |
m[i][j]=Math.pow(-1, (i+j))*det; |
} |
} |
return this.transpose(m); // Array |
}, |
transpose: function(/* Array */a){ |
// summary |
// Transpose the passed matrix (i.e. rows to columns) |
var m=this.create(a.length, a[0].length); |
for(var i=0; i<a.length; i++){ |
for(var j=0; j<a[i].length; j++){ |
m[j][i]=a[i][j]; |
} |
} |
return m; // Array |
}, |
format: function(/* Array */a, /* Number? */points){ |
// summary |
// Return a string representation of the matrix, rounded to points (if needed) |
points=points||5; |
function format_int(x, dp){ |
var fac=Math.pow(10, dp); |
var a=Math.round(x*fac)/fac; |
var b=a.toString(); |
if(b.charAt(0)!="-"){ |
b=" "+b; |
} |
if(b.indexOf(".")>-1){ |
b+="."; |
} |
while(b.length<dp+3){ |
b+="0"; |
} |
return b; |
} |
var ya=a.length; |
var xa=ya>0?a[0].length:0; |
var buffer=""; |
for(var y=0; y<ya; y++){ |
buffer+="| "; |
for(var x=0; x<xa; x++){ |
buffer+=format_int(a[y][x], points)+" "; |
} |
buffer+="|\n"; |
} |
return buffer; // string |
}, |
copy: function(/* Array */a){ |
// summary |
// Create a copy of the passed matrix |
var ya=a.length, xa=a[0].length, m=this.create(xa, ya); |
for(var y=0; y<ya; y++){ |
for(var x=0; x<xa; x++){ |
m[y][x]=a[y][x]; |
} |
} |
return m; // Array |
}, |
scale: function(/* Array */a, /* Number */factor){ |
// summary |
// Create a copy of passed matrix and scale each member by factor. |
a=this.copy(a); |
var ya=a.length, xa=a[0].length; |
for(var y=0; y<ya; y++){ |
for(var x=0; x<xa; x++){ |
a[y][x]*=factor; |
} |
} |
return a; |
} |
}); |
|
} |