170 похожих чатов

Понял, просто у меня функция scipy.sparse.linalg.spsolve() отрабатывает очень медленно.

Не знаю, что делать. Нужно, наверное, искать аналог на numpy, может быть видели?

19 ответов

22 просмотра

Думаешь она без нампая реализована?

Vadzim-_B Автор вопроса

Можешь попробовать итерационные методы

С учётом что он принимает ndarray вариантов мало

А размерность какая?

Vadzim-_B Автор вопроса
Tishka17
А размерность какая?

Всегда по-разному, но обычно очень большая. Обычно эта функция отрабатывает за 10-15 секунд

Vadzim-_B Автор вопроса
Vadzim _B
да вряд ли это поможет

Смотря какая точность решения нужна

Может лучше пример какой нибудь сгенерируешь с кодом чтобы вот тоже неоправданно долго работал? А то так только гадать остается

Vadzim-_B Автор вопроса
Владимир
Может лучше пример какой нибудь сгенерируешь с код...

Ab = scipy.sparse.coo_matrix((Sb, (Ib, Jb)), shape=(numRowsInA, alls)) Ag = scipy.sparse.coo_matrix((Sg, (Ig, Jg)), shape=(numRowsInA, alls)) Ar = scipy.sparse.coo_matrix((Sr, (Ir, Jr)), shape=(numRowsInA, alls)) Ab = Ab.tocsc() Ag = Ag.tocsc() Ar = Ar.tocsc() Rb = scipy.sparse.linalg.spsolve(Ab, Bb) Rg = scipy.sparse.linalg.spsolve(Ag, Bg) Rr = scipy.sparse.linalg.spsolve(Ar, Br)

Vadzim _B
Ab = scipy.sparse.coo_matrix((Sb, (Ib, Jb)), shape...

Хмм, погляжу, мб несколько позже что-то смогу предложить

А может можно параллелизовать?

Vadzim-_B Автор вопроса
Tishka17
А может можно параллелизовать?

думаю, что нужно это попробывать, спасибо за идею)

Vadzim-_B Автор вопроса
Tishka17
А может можно параллелизовать?

multiprocessing же ты имел ввиду?

Vadzim _B
multiprocessing же ты имел ввиду?

Нет, я имел ввиду встроенные возможности либы.

Vadzim-_B Автор вопроса
Tishka17
А 15 секунд - это долго?

да, желательно, чтобы раз в 40 быстрее отрабатывало

Похожие вопросы

Обсуждают сегодня

Господа, а что сейчас вообще с рынком труда на делфи происходит? Какова ситуация?
Rꙮman Yankꙮvsky
29
А вообще, что может смущать в самой Julia - бы сказал, что нет единого стандартного подхода по многим моментам, поэтому многое выглядит как "хаки" и произвол. Короче говоря, с...
Viktor G.
2
30500 за редактор? )
Владимир
47
а через ESC-код ?
Alexey Kulakov
29
Чёт не понял, я ж правильной функцией воспользовался чтобы вывести отладочную информацию? но что-то она не ловится
notme
18
У меня есть функция где происходит это: write_bit(buffer, 1); write_bit(buffer, 0); write_bit(buffer, 1); write_bit(buffer, 1); write_bit(buffer, 1); w...
~
14
Добрый день! Скажите пожалуйста, а какие программы вы бы рекомендовали написать для того, чтобы научиться управлять памятью? Можно написать динамический массив, можно связный ...
Филипп
7
Недавно Google Project Zero нашёл багу в SQLite с помощью LLM, о чём достаточно было шумно в определённых интернетах, которые сопровождались рассказами, что скоро всех "ибешни...
Alex Sherbakov
5
Ребят в СИ можно реализовать ООП?
Николай
33
https://github.com/erlang/otp/blob/OTP-27.1/lib/kernel/src/logger_h_common.erl#L174 https://github.com/erlang/otp/blob/OTP-27.1/lib/kernel/src/logger_olp.erl#L76 15 лет назад...
Maksim Lapshin
20
Карта сайта