Как я могу использовать мод в массиве с помощью numpy

Проблема связана с тем, что numpy не может применить операцию np.mod к массивам с заданной формой, одним из решений является изменение формы, например:

import numpy as np

a = np.array([6, 8, 10, 13, 15, 18, 21]).reshape((-1, 1))
b = np.array([2, 5])

print(a[np.mod(a, b).all(1)].reshape(-1))

Выход

[13 21]

Обратите внимание, что вам нужно изменить форму, чтобы получить запрошенный результат. Лучшее решение - предложенное @PaulPanzer:

import numpy as np

a = np.array([6, 8, 10, 13, 15, 18, 21])
b = np.array([2, 5])

print(a[np.mod.outer(a, b).all(1)])

Выход

[13 21]

Дальше

1. О numpy-вещании см. 1 и 2.
2. На внешний.