% Esta funcion estandariza el problema inverso
%
% Sintaxis:
%           [Ui,si,Vis,R,ord,w]=stdbc(K,L)
% Donde:
%        K - Matriz del leadfield.
%        L - Laplaciano.
%
% Para mas informacion ver Readme.doc. 
%
% Autor: Nelson J. Trujillo Barreto.

function [Ui,si,Vis,R,ord,w]=stdbcm(K,L)
[ne,nggl]=size(K);

%H = eye(ne) - ones(ne,ne)/ne;
%K = H*K;

%Lr=L(1:3:nggl,1:3:nggl);
if L(1,1)<0 ;  L=-L; end;
[R,ord]=lchol(L);
clear L

%w = wiiv(K,1);
w=ones(nggl,1);
K = klinvr(R,K',ord,w);

[Ui,si,Vi] = svdrapid(K);

Vis = ljinvl(R,Vi,ord,w);
