function [u,d,v]=SVDBiDiaJacob(B)
[m,n]=size(B);
B1=B(1:n,1:n);
[u,d,v]=jacobi_svd(B1);
d(n+1:m,1:n)=zeros(m-n,n);
u(n+1:m,n+1:m)=eye(m-n);

