adding all files from local

Author: proteekroy

BuildStochasticMatrix.m

@@ -0,0 +1,63 @@
+function M = BuildStochasticMatrix(opt, population)
+
+ 
+    switch(opt.UtilityOption)
+        case 1 
+            n = size(population,1);%size of population;
+            M = sparse(n, n);%large sqaure matrix
+            u = computeUtility1(population);%suppose this is L1-norm
+                        
+            for i = 1:n
+                for j = 1:n
+                    if(i~=j)
+                        if u(i)>u(j)
+                            M(j, i) = u(i) - u(j);%direction from (i) to (j)
+                        elseif u(i)<u(j)
+                            M(i, j) = u(j) - u(i);
+                        end
+                    end
+                end
+            end
+            
+            %normalize to make stochastic
+            for i = 1:n
+                M(:,i) = M(:,i)./sum(M(:,i));
+            end
+              
+        case 2
+            n = size(population,1);%size of population
+            M = sparse(n, n);%declare a sparse matrix
+            Mdl = KDTreeSearcher(population,'Distance','euclidean');
+            u = computeUtility1(population);%suppose this is L1-norm
+            d = size(population,2);
+            [Idx,~] = knnsearch(Mdl,population,'K',d+1);%k-d tree
+            for i = 1:n
+                for j = 2:d+1 %1st member is itself
+                    s = Idx(j);
+                    if(i~=s)
+                        if u(i)>u(s)
+                            %direction from (i) to (j), value iteration is reverse, it takes indegree
+                            %rather than outdegree
+                            M(s, i) = u(i) - u(s);
+                        elseif u(i)<u(s)
+                            M(i, s) = u(s) - u(i);
+                        end
+                    end
+                end
+            end
+            
+            %normalize M to make column stochastic
+            for i = 1:n
+                M(:,i) = M(:,i)./sum(M(:,i));
+            end
+            
+        otherwise
+            
+            
+    end
+
+end
+
+function u = computeUtility1(population)
+    u = sum(population,2); %should be normalized   ----   %%sum(population.*0.5,2);
+end

GD/.IGD.2D.txt

@@ -0,0 +1,21 @@
+136.0000000000 4.0000000000
+126.6304000000 4.0576000000
+117.0000000000 4.2500000000
+107.2336000000 4.6084000000
+96.8400000000 5.2100000000
+86.4100000000 6.1025000000
+75.4384000000 7.4596000000
+64.9800000000 9.2450000000
+55.3352000000 11.2338000000
+45.6968000000 13.6242000000
+36.6368000000 16.3592000000
+28.5768000000 19.3442000000
+21.2552000000 22.7138000000
+15.2352000000 26.2088000000
+10.5800000000 29.6450000000
+6.7712000000 33.2928000000
+4.0328000000 36.8082000000
+2.0808000000 40.3202000000
+0.8712000000 43.6178000000
+0.2048000000 46.8512000000
+0.0000000000 50.0000000000