Замена для бесшумного вещания с использованием 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() суммирует по неправильной оси.)

Ответы на вопрос(2)

Ваш ответ на вопрос