// genesis


int i, j, k, l, nsynapses12, nsynapses21, reciproc_synapses, reciproc_pairs, pairs
int target_index, sender_index
str name1, target1, sender1, name2, sender2, target2, target_path, sender_path, synsrc
str output, list
  
  sender_path = "/molecular_layer/Stellate"
  target_path = "/molecular_layer/Stellate"
  output = "number-reciproc-synapses.dat"
  list = ""
  
  reciproc_synapses = 0
  reciproc_pairs = 0
  pairs = 0

  for (k = 0; k < {number_stellate_cells - 1}; k = k + 1)

      l = {k + 1}

      while (l < {number_stellate_cells})

//        echo  k  {k}   l  {l}   i   {i}    j    {j}     sender   {sender_index}      target    {target_index}  

        //  sender = {source_path}[{sender_index}]/soma/spike
        sender1 = {sender_path} @ "[" @ {k} @ "]/soma/spike"
        target1 = {target_path} @ "[" @ {l} @ "]/##[][TYPE=compartment]/GABAA"

        sender2 = {sender_path} @ "[" @ {l} @ "]/soma/spike"
        target2 = {target_path} @ "[" @ {k} @ "]/##[][TYPE=compartment]/GABAA"

        nsynapses12 = 0

           foreach name1 ({el {target1}})
                   nsynapses12 = {{nsynapses12} + {getsyncount {sender1} {name1}}}
           end // foreach


        nsynapses21 = 0

           foreach name2 ({el {target2}})
                   nsynapses21 = {{nsynapses21} + {getsyncount {sender2} {name2}}}
           end // foreach

        if ({nsynapses12 > 0} && {nsynapses21 > 0})

           reciproc_synapses = {reciproc_synapses} + {nsynapses12} + {nsynapses21}
           reciproc_pairs = {reciproc_pairs} + 1
           pairs = {pairs} + 1
           echo {k} {l} {nsynapses12} {nsynapses21} {reciproc_synapses} {reciproc_pairs}  
           list = (list) @ " " @ {k} @ " " @ {l} @ " " @ {nsynapses12} @ " " @ {nsynapses21} @ {chr 10}

        elif ({nsynapses12 > 0} || {nsynapses21 > 0})
                pairs = {pairs} + 1
        end // if

        l = {l + 1}

    end // while 

end // for

list = (list) @ "number reciprocal synapses = " @ {reciproc_synapses} @ {chr 10}
list = (list) @ "number reciprocal pairs = " @ {reciproc_pairs} @ {chr 10}
list = (list) @ "number pairs = " @ {pairs} @ {chr 10}
echo {list} > {output}