| 12345678 |
- /*
- Copyright (c) 2004-2009, The Dojo Foundation All Rights Reserved.
- Available via Academic Free License >= 2.1 OR the modified BSD license.
- see: http://dojotoolkit.org/license for details
- */
- if(!dojo._hasResource["dojox.math.matrix"]){dojo._hasResource["dojox.math.matrix"]=true;dojo.provide("dojox.math.matrix");dojo.mixin(dojox.math.matrix,{iDF:0,ALMOST_ZERO:1e-10,multiply:function(a,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;},product:function(){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;},sum:function(){if(arguments.length==0){console.warn("can't sum 0 matrices!");return 0;}var m=this.copy(arguments[0]);var _1=m.length;if(_1==0){console.warn("can't deal with matrices of 0 rows!");return 0;}var _2=m[0].length;if(_2==0){console.warn("can't deal with matrices of 0 cols!");return 0;}for(var i=1;i<arguments.length;++i){var _3=arguments[i];if(_3.length!=_1||_3[0].length!=_2){console.warn("can't add matrices of different dimensions: first dimensions were "+_1+"x"+_2+", current dimensions are "+_3.length+"x"+_3[0].length);return 0;}for(var r=0;r<_1;r++){for(var c=0;c<_2;c++){m[r][c]+=_3[r][c];}}}return m;},inverse:function(a){if(a.length==1&&a[0].length==1){return [[1/a[0][0]]];}var _4=a.length,m=this.create(_4,_4),mm=this.adjoint(a),_5=this.determinant(a),dd=0;if(_5==0){console.warn("Determinant Equals 0, Not Invertible.");return [[0]];}else{dd=1/_5;}for(var i=0;i<_4;i++){for(var j=0;j<_4;j++){m[i][j]=dd*mm[i][j];}}return m;},determinant:function(a){if(a.length!=a[0].length){console.warn("Can't calculate the determinant of a non-squre matrix!");return 0;}var _6=a.length,_7=1,b=this.upperTriangle(a);for(var i=0;i<_6;i++){var _8=b[i][i];if(Math.abs(_8)<this.ALMOST_ZERO){return 0;}_7*=_8;}_7*=this.iDF;return _7;},upperTriangle:function(m){m=this.copy(m);var f1=0,_9=0,_a=m.length,v=1;this.iDF=1;for(var _b=0;_b<_a-1;_b++){if(typeof m[_b][_b]!="number"){console.warn("non-numeric entry found in a numeric matrix: m["+_b+"]["+_b+"]="+m[_b][_b]);}v=1;var _c=0;while((m[_b][_b]==0)&&!_c){if(_b+v>=_a){this.iDF=0;_c=1;}else{for(var r=0;r<_a;r++){_9=m[_b][r];m[_b][r]=m[_b+v][r];m[_b+v][r]=_9;}v++;this.iDF*=-1;}}for(var _d=_b+1;_d<_a;_d++){if(typeof m[_d][_b]!="number"){console.warn("non-numeric entry found in a numeric matrix: m["+_d+"]["+_b+"]="+m[_d][_b]);}if(typeof m[_b][_d]!="number"){console.warn("non-numeric entry found in a numeric matrix: m["+_b+"]["+_d+"]="+m[_b][_d]);}if(m[_b][_b]!=0){var f1=(-1)*m[_d][_b]/m[_b][_b];for(var i=_b;i<_a;i++){m[_d][i]=f1*m[_b][i]+m[_d][i];}}}}return m;},create:function(a,b,_e){_e=_e||0;var m=[];for(var i=0;i<b;i++){m[i]=[];for(var j=0;j<a;j++){m[i][j]=_e;}}return m;},ones:function(a,b){return this.create(a,b,1);},zeros:function(a,b){return this.create(a,b);},identity:function(_f,_10){_10=_10||1;var m=[];for(var i=0;i<_f;i++){m[i]=[];for(var j=0;j<_f;j++){m[i][j]=(i==j?_10:0);}}return m;},adjoint:function(a){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);},transpose:function(a){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;},format:function(a,_11){_11=_11||5;function _12(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 _13="";for(var y=0;y<ya;y++){_13+="| ";for(var x=0;x<xa;x++){_13+=_12(a[y][x],_11)+" ";}_13+="|\n";}return _13;},copy:function(a){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;},scale:function(a,_14){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]*=_14;}}return a;}});}
|
