Example #1
0
def interact(c1, c2, s):
    if c1.get_num_particles() <= s and \
       c2.get_num_particles() <= s:

        common.direct([c1.particles[i] for i in c1.indices],
                      [c2.particles[i] for i in c2.indices],
                      c1.region == c2.region)
    else:
        if c1.get_num_particles() > s:
            c1_children = c1.get_sub_cells()
        else:
            c1_children = [c1]
        if c2.get_num_particles() > s:
            c2_children = c2.get_sub_cells()
        else:
            c2_children = [c2]
        for i, j in taco.product(c1_children, c2_children):
            interact(i, j, s)
    return
Example #2
0
def interact(c1, c2, s):
    if c1.get_num_particles() <= s and \
       c2.get_num_particles() <= s:

        common.direct([c1.particles[i] for i in c1.indices],
                      [c2.particles[i] for i in c2.indices],
                      c1.region == c2.region)
    else:
        if c1.get_num_particles() > s:
            c1_children = c1.get_sub_cells()
        else:
            c1_children = [c1]
        if c2.get_num_particles() > s:
            c2_children = c2.get_sub_cells()
        else:
            c2_children = [c2]
        for i, j in taco.product(c1_children, c2_children):
            interact(i, j, s)
    return
Example #3
0
def interact(c1, c2, s):
    if c1.is_leaf() and c2.is_leaf():
        common.direct([c1.particles[i] for i in c1.indices],
                      [c2.particles[i] for i in c2.indices],
                      c1.region == c2.region)
    return