RS
Size: a a a
2021 May 09
VA
Подскажите, пожалуйста, как сделать, чтобы после каждого этапа у решателя обновлялось Р с помощью функции из 147 строчки? Функция возвращает массив значений по смыслу равный Р
АО
Но ведь
P это же скаляр, а не вектор?АО
А можно посмотреть на систему уравнений, которая решается? Не в виде Джулия кода, а в виде математических выражений.
🐙
Спасибо, гляну 👌🏽
RS
В 148 строке повторяющаяся сумма внутри цикла. К вопросу о производительности
VA
P - это вектор
VA
Могу скинуть код, там сложные выводы функций через другие функции
VA
dF/dz - du
dz - tspan
R1,R2 .. - это подфункции, прописанные ранее
dz - tspan
R1,R2 .. - это подфункции, прописанные ранее
P
Если я правильно понял, функцию с парциальными давлениями можно векторизовать
(P_total / sum(u)) .* u (не в тему коммент)VA
примерно так выглядят R1
VA
функция get_pressure возвращает вектор и так
P
Я чисто из стилистики. comprehension кмк лишний
P
но это так, до скрепок
P
а нельзя вычислять парциальные давления прямо в
molstreams, раз они есть функция от u?VA
Я пробовал, но видимо как-то не так
Я в самом конце функции вызывал get_pressure
А где вы предлагаете обозначить вычисление парциальных?
Я в самом конце функции вызывал get_pressure
А где вы предлагаете обозначить вычисление парциальных?
P
строки 153-159 это тоже парциальные?
VA
нет, это другие параметры
VA
они задаются в самом начале, а парциальные рассчитываются из мольных потоков u