prepare structure for R package

[valse.git] / OLD_MATLAB / InputParameters / grillelambda.m

diff --git a/OLD_MATLAB/InputParameters/grillelambda.m b/OLD_MATLAB/InputParameters/grillelambda.m
new file mode 100644 (file)
index 0000000..8d90665
--- /dev/null
+++ b/OLD_MATLAB/InputParameters/grillelambda.m
@@ -0,0 +1,17 @@
+function[gridLambda] = grillelambda(phiInit,rhoInit,piInit,gamInit,X,Y,gamma,mini,maxi,tau)
+
+       n = size(X, 1);
+       [p,m,k] = size(phiInit);
+       [phi,rho,pi,~,S] = EMGLLF(phiInit,rhoInit,piInit,gamInit,mini,maxi,1,0,X,Y,tau);
+
+       gridLambda = zeros(p,m,k);
+       for j=1:p
+               for mm=1:m
+                       gridLambda(j,mm,:) = abs(reshape(S(j,mm,:),[k,1])) ./ (n*pi.^gamma);
+               end
+       end
+       
+       gridLambda = unique(gridLambda);
+       gridLambda(gridLambda()>1) = [];
+
+end