Замена для бесшумного вещания с использованием scipy.sparse.csc_matrix
В моем коде есть следующее выражение:
a = (b / x[:, np.newaxis]).sum(axis=1)
гдеb
ndarray формы(M, N)
, а такжеx
ndarray формы(M,)
, Сейчас,b
на самом деле редко, поэтому для эффективности памяти я хотел бы заменить вscipy.sparse.csc_matrix
или жеcsr_matrix
, Однако вещание таким способом не реализовано (даже если деление или умножение гарантированно сохраняют разреженность) (записиx
ненулевые), и поднимаетNotImplementedError
, Есть лиsparse
функция яЯ не знаю, что будет делать то, что я хочу? (dot()
суммирует по неправильной оси.)