Linear least square solver using svd, written for numeric js

gistfile1.txt

function svd_solve(A,b){
	var A_svd,r,s,ut,d,x;

	A_svd = numeric.svd(A);
	s = numeric.transpose( [ A_svd.S ] ); 
	r = 1;


	while( r < numeric.dim(A)[1] & s[r] >= Math.max.apply( Math, numeric.dim(A) )*numeric.epsilon*s[0] ) r++; 
  	d = numeric.dot( numeric.transpose( A_svd.U ), b);
	x = numeric.dot(A_svd.V, numeric.div(d,s) );
	return x;
}

r, which seems to calculate the rank(A), is calculated but not used.

Damn, you're right. That's annoying. So long ago now, I can't event remember my source.