diff --git a/completeGlobal.m b/completeGlobal.m
index e09f5127c6b1c8e0c3cca1eb5d1b5b1c09e4c781..0e5a58d3d8622ee5f3b270bb456a95fe51cf2500 100644
--- a/completeGlobal.m
+++ b/completeGlobal.m
@@ -1,5 +1,5 @@
 function SRCompSet = completeGlobal(SRRefSet,varargin)
-%COMPLETE Summary of this function goes here
+%COMPLETEGLOBAL Summary of this function goes here
 %   Detailed explanation goes here
 SRSize  = [];
 maxSize = [];
diff --git a/completeLF.m b/completeLF.m
index f60eee2661cc7d0509f0ed81d72440511d2078e6..37dca6c9c94c31c4ced6fb98197b89cfebf738c3 100644
--- a/completeLF.m
+++ b/completeLF.m
@@ -1,5 +1,5 @@
 function [LFComp,LFMaskComp] = completeLF(LFRef,varargin)
-%COMPLETE Summary of this function goes here
+%COMPLETELF Summary of this function goes here
 %   Detailed explanation goes here
 
 validMethods = {'none','trace','rank'};
diff --git a/completeLocal.m b/completeLocal.m
index 54ecdb13d811ab72ea79e961090c888dbe2d0fb4..cb4864527db0e58545208eb6dcedccd6bc6b59fe 100644
--- a/completeLocal.m
+++ b/completeLocal.m
@@ -1,5 +1,5 @@
 function SRCompSet = completeLocal(SRRefSet,varargin)
-%COMPLETE Summary of this function goes here
+%COMPLETELOCAL Summary of this function goes here
 %   Detailed explanation goes here
 
 SRCompSet = struct(SRRefSet);
diff --git a/completeSR.m b/completeSR.m
index 5d5b443b5879dd00428221ec82ae5f44f1006355..060f5266e3781af3d8f3d882e0eb669554e176f1 100644
--- a/completeSR.m
+++ b/completeSR.m
@@ -1,5 +1,5 @@
 function [SRComp,SRMaskComp] = completeSR(SRRef,varargin)
-%COMPLETE Summary of this function goes here
+%COMPLETESR Summary of this function goes here
 %   Detailed explanation goes here
 
 [SRComp,SRMaskComp] = cellfun(@(SRRef)(SR.completeLF(SRRef,varargin{:})),SRRef,...
diff --git a/factorize.m b/factorize.m
index b0adcc8942cd10a932f66bf2fcb0829e81411bc1..bded811ef5bd83d73d6fe2e9263cd7d1189610cd 100644
--- a/factorize.m
+++ b/factorize.m
@@ -1,5 +1,5 @@
 function [SRBSet,SRC,SRU,SRS,SRV] = factorize(SRRefSet)
-%SRTOBC Summary of this function goes here
+%FACTORIZE Summary of this function goes here
 %   Detailed explanation goes here
 
 SRBSet = struct(SRRefSet);
diff --git a/factorizeGlobal.m b/factorizeGlobal.m
index dcd46db03eb5189d336226bd54a07eb19c3911a8..f24af9e6ad385bd83d9c82b9f582ce08b9f3c31d 100644
--- a/factorizeGlobal.m
+++ b/factorizeGlobal.m
@@ -1,5 +1,5 @@
 function [SRBSet,LFC,LFU,LFS,LFV] = factorizeGlobal(SRRefSet)
-%SRTOBC Summary of this function goes here
+%FACTORIZEGLOBAL Summary of this function goes here
 %   Detailed explanation goes here
 SRSize  = [];
 maxSize = [];
@@ -11,6 +11,7 @@ LFRef = SRtoLF(SRRef);
 [LFB,LFC,LFU,LFS,LFV] = SR.factorizeLF(LFRef);
 
 SRB = LFtoSR(LFB);
+
 for lab=1:numel(SRRef)
     SRBSet(lab).Value = SRB{lab};
 end
diff --git a/factorizeLF.m b/factorizeLF.m
index 32f2a79833a8d7e545800fc3d631988b366b0f15..1dd184c3019d7bfe45e3b91cc51a7ee68dde9819 100644
--- a/factorizeLF.m
+++ b/factorizeLF.m
@@ -1,5 +1,5 @@
 function [B,C,U,S,V] = factorizeLF(LF)
-%LFTOBC Summary of this function goes here
+%FACTORIZELF Summary of this function goes here
 %   Detailed explanation goes here
 
 Mask = isnan(LF); LF(Mask) = 0;
diff --git a/factorizeLocal.m b/factorizeLocal.m
index 54fb31d8e99f13aebe26743c624d3e1d5438890d..1604c053d4683b6a7bbfc9205d41973c30b7b955 100644
--- a/factorizeLocal.m
+++ b/factorizeLocal.m
@@ -1,5 +1,5 @@
 function [SRBSet,SRC,SRU,SRS,SRV] = factorizeLocal(SRRefSet)
-%SRTOBC Summary of this function goes here
+%FACTORIZELOCAL Summary of this function goes here
 %   Detailed explanation goes here
 
 SRBSet = struct(SRRefSet);
diff --git a/factorizeSR.m b/factorizeSR.m
index 8b92c5ac03950e320a3cd59077a36215a3e8f5be..3105ebc097b49769993fa7c7b5715d72b12213a7 100644
--- a/factorizeSR.m
+++ b/factorizeSR.m
@@ -1,5 +1,5 @@
 function [SRB,SRC,SRU,SRS,SRV] = factorizeSR(SRRef)
-%SRTOBC Summary of this function goes here
+%FACTORIZESR Summary of this function goes here
 %   Detailed explanation goes here
 
 [SRB,SRC,SRU,SRS,SRV] = cellfun(@SR.factorizeLF,SRRef,'UniformOutput',false);
diff --git a/multiply.m b/multiply.m
index a0eb6f9dc09e52ceea28863f8f8417d850c46aa5..b8971a1cfad28e9a996af53797a8cf6f9e421c42 100644
--- a/multiply.m
+++ b/multiply.m
@@ -1,13 +1,13 @@
-function SRSet = multiply(BSet,SRC,varargin)
-%BCTOLF Summary of this function goes here
+function SRSet = multiply(SRBSet,SRC,varargin)
+%MULTIPLY Summary of this function goes here
 %   Detailed explanation goes here
 
-SRSet = struct(BSet);
-SRB = {BSet.Value};
+SRSet = struct(SRBSet);
+SRB = {SRBSet.Value};
 
 SRRef = SR.multiplySR(SRB,SRC,varargin{:});
 
-for lab=1:numel(BSet)
+for lab=1:numel(SRBSet)
     SRSet(lab).Value = SRRef{lab};
 end
 
diff --git a/multiplyLF.m b/multiplyLF.m
index 6c8fad01d50c4220891d0d5db779cf2d059ab832..a9abaf4cc87ae82cb651d78800be771de674688d 100644
--- a/multiplyLF.m
+++ b/multiplyLF.m
@@ -1,19 +1,19 @@
-function LF = multiplyLF(B,C,varargin)
-%BCTOLF Summary of this function goes here
+function LF = multiplyLF(LFB,LFC,varargin)
+%MULTIPLYLF Summary of this function goes here
 %   Detailed explanation goes here
 
 vecResh = @(x)   reshape(x,size(x,1)*size(x,2),[])';
-boxResh = @(x)   reshape(x',size(B));
+boxResh = @(x)   reshape(x',size(LFB));
 
 if nargin>2
     k = varargin{1};
-    [u,v] = ind2sub([size(B,1),size(B,2)],k);
-    B = B(u,v,:,:);
-    C = C(k,:);
+    [u,v] = ind2sub([size(LFB,1),size(LFB,2)],k);
+    LFB = LFB(u,v,:,:);
+    LFC = LFC(k,:);
 end
 
-Bvec = vecResh(B);
-LFvec = Bvec*C;
+Bvec = vecResh(LFB);
+LFvec = Bvec*LFC;
 LF = boxResh(LFvec);
 
 end
diff --git a/multiplyLocal.m b/multiplyLocal.m
new file mode 100644
index 0000000000000000000000000000000000000000..4d5af3945275b134265de17a0a00f40e6d506a2c
--- /dev/null
+++ b/multiplyLocal.m
@@ -0,0 +1,14 @@
+function SRSet = multiplyLocal(SRBSet,SRC,varargin)
+%MULTIPLYLOCAL Summary of this function goes here
+%   Detailed explanation goes here
+
+SRSet = struct(SRBSet);
+SRB = {SRBSet.Value};
+
+SRRef = SR.multiplySR(SRB,SRC,varargin{:});
+
+for lab=1:numel(SRBSet)
+    SRSet(lab).Value = SRRef{lab};
+end
+
+end
\ No newline at end of file
diff --git a/multiplySR.m b/multiplySR.m
index 5e2579bf75d5c5a53bfde654a436d7758d6ed88c..cd472ecbe33b68af98ce35e9d0d83fdd3e6b38ea 100644
--- a/multiplySR.m
+++ b/multiplySR.m
@@ -1,5 +1,5 @@
 function SRRef = multiplySR(SRB,SRC,varargin)
-%BCTOSR Summary of this function goes here
+%MULTIPLYSR Summary of this function goes here
 %   Detailed explanation goes here
 
 if nargin>2