Skip to content
Snippets Groups Projects
Commit f67ed9e8 authored by Lars Wiebach's avatar Lars Wiebach
Browse files

mlapfds

parents 5bd80581 74814eb1
No related branches found
No related tags found
No related merge requests found
...@@ -20,15 +20,16 @@ def pso (n, pop, sample, c1=1/3, c2=1/3, c3=1/3, abort = 50): ...@@ -20,15 +20,16 @@ def pso (n, pop, sample, c1=1/3, c2=1/3, c3=1/3, abort = 50):
population = list() population = list()
# create pop-many particle # create pop-many particle
for i in range(pop): for i in range(pop):
new_particle = part.Particle(n) new_particle = part.Particle(n)
population.append(new_particle) population.append(new_particle)
for p in population: for p in population:
p.local_attractor=p.position p.local_attractor=p.position
p.local_attractor_cost = pl.cost_function(sample,p.local_attractor) p.local_attractor_cost = pl.cost_function(sample,p.local_attractor)
if(p.local_attractor_cost < global_attractor_cost): if(p.local_attractor_cost < global_attractor_cost):
global_attractor = p.position global_attractor = p.position
global_attractor_cost = p.local_attractor_cost global_attractor_cost = p.local_attractor_cost
...@@ -58,6 +59,6 @@ def pso (n, pop, sample, c1=1/3, c2=1/3, c3=1/3, abort = 50): ...@@ -58,6 +59,6 @@ def pso (n, pop, sample, c1=1/3, c2=1/3, c3=1/3, abort = 50):
if(p.cost<global_attractor_cost): if(p.cost<global_attractor_cost):
global_attractor=p.position global_attractor=p.position
global_attractor_cost = p.cost global_attractor_cost = p.cost
return test_array return test_array
File added
File added
File added
No preview for this file type
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment