NetLogo Efficient way to create fixed number of links

后端 未结 2 926
别跟我提以往
别跟我提以往 2020-12-06 14:12

I have about 5000 agents (people) in my model. I want to give them an arbitrary number of friends and have reciprocal but random pairing. So if person A chooses person B the

2条回答
  •  醉话见心
    2020-12-06 14:54

    Great question. This is actually quite challenging to optimize. The problematic line is:

    let candidates other people with [ count my-links < friends ]

    This is slow because it has every agent checking with every other agent. With 5000 agents, that's 25,000,000 checks! Unfortunately, there isn't really a good way to optimize this particular line without some fancy data structures.

    Fortunately, there is a solution that generalizes really well to generating any degree distribution in the network (which it sounds like that's what you ultimately want). Unfortunately, the solution doesn't translate super well to NetLogo. Here it is though:

      let pairs [] ;; pairs will hold a pairs of turtles to be linked
      while [ pairs = [] ] [ ;; we might mess up creating these pairs (by making self loops), so we might need to try a couple of times
        let half-pairs reduce sentence [ n-values friends [ self ] ] of turtles ;; create a big list where each turtle appears once for each friend it wants to have
        set pairs (map list half-pairs shuffle half-pairs) ;; pair off the items of half-pairs with a randomized version of half-pairs, so we end up with a list like: [[ turtle 0 turtle 5 ] [ turtle 0 turtle 376 ] ... [ turtle 1 turtle 18 ]]
        ;; make sure that no turtle is paired with itself
        if not empty? filter [ first ? = last ? ] pairs [
          set pairs []
        ]
      ]
      ;; now that we have pairs that we know work, create the links
      foreach pairs [
        ask first ? [
          create-link-with last ?
        ]
      ]
    

    It doesn't matter if friends here is a global or a turtle variable. The amount of time this takes depends on the number of times that it needs to try making pairs, which is random. Experimenting, I found that it was usually about 3 seconds with 5000 agents, each with degree 5. This is compared to about 60 seconds on my machine with your original way of doing this (which, for what it's worth, is the way I would recommend when using fewer agents).

提交回复
热议问题