Subversion Repositories Applications.papyrus

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;

	}

});

}