add checkBoundsToroidal

Author: aeroreyna

CS.m

@@ -31,7 +31,7 @@ function operators(obj)
                 s = s + stepsize .* randn(size(s));
                 newPop(i,:) = s;
             end
-            newPop = obj.checkBounds(newPop);
+            newPop = obj.checkBoundsToroidal(newPop);
             obj.remplacePopulation(newPop);
 
             % A fraction of worse nests are discovered with a probability pa
@@ -41,7 +41,7 @@ function operators(obj)
             stepsize = rand * (obj.population(randperm(obj.sizePopulation),:)-...
                                obj.population(randperm(obj.sizePopulation),:));
             newPop = obj.population + stepsize.*K;
-            newPop = obj.checkBounds(newPop);
+            newPop = obj.checkBoundsToroidal(newPop);
             obj.remplacePopulation(newPop);
         end
 

DE.m

@@ -40,8 +40,7 @@ function operators(obj)
                                obj.population(selectedPop(5),j));
                     end
                 end
-                solutionBase(solutionBase > 1) = 1;
-                solutionBase(solutionBase < 0) = 0;
+                solutionBase = obj.checkBoundsToroidal(solutionBase);
                 tempFit = obj.evalPopulation(solutionBase);
                 if tempFit < obj.fitness(selectedPop(1))
                     if obj.plotSpecial

PSO.m

@@ -50,7 +50,7 @@ function operators(obj)
             end
             % update solutions
             obj.population = obj.population + obj.velocity;
-            obj.checkBounds();
+            obj.checkBoundsToroidal();
 
             % eval fitness
             obj.evalPopulation();

SMS.m

@@ -69,7 +69,7 @@ function operators(obj)
             end
 
             obj.population = pobL;
-            obj.checkBounds();
+            obj.checkBoundsToroidal();
             obj.direction = dirL;
             obj.evalPopulation();
         end

WOA.m

@@ -48,7 +48,7 @@ function operators(obj)
                     end
                 end
             end
-            obj.checkBounds();
+            obj.checkBoundsToroidal();
             obj.evalPopulation();
             obj.updateBest();
         end

metaheuristic.m

@@ -161,6 +161,17 @@ function sortPopulation(obj)
             end
         end
 
+        function solutions = checkBoundsToroidal(obj, solutions)
+            if nargin == 1
+                solutions = obj.population;
+            end
+            solutions(solutions>1)=solutions(solutions>1)-1;
+            solutions(solutions<0)=solutions(solutions<0)+1;
+            if nargin == 1
+                obj.population = solutions;
+            end
+        end
+        
         function updateBest(obj)
             % update the best know so far solution and it's fitness.
             [bestFitTemp, bestIndex] = min(obj.fitness);