What are some recommended ways to create networks of ensembles of neurons?

Hi! I’m very new to Julia and considering migrating my Izhikevich neuronal network simulations to Julia from MATLAB.

I was wondering about ways of implementing interactions among ensembles of units. I’d like to have populations of Izhikevich neurons (implemented in Julia here), each potentially with unique parameters, that form networks. So there might be hundreds of pyramidal cells and many types of interneurons.

What are the recommended ways to set this up? How might I create the topology? In MATLAB I was just using vectors and matrices and keeping track of the indices of each subpopulation, but this was getting clunky, and I had to create unique matrices for each connection between ensembles. How might I used modelingtoolkit for this?

Any suggestions welcomed!