A
Size: a a a
2019 February 15
VG
Спасибо! Ну хоть бы в названии упомянули, что это Julia)
AK
KT
Начал читать, отметил для себя одну странность с identity matrix
2019 February 16
KT
Вот это место:
MatrixT(I, (n,n)) is the n×n identity matrix. The identity
operator I is available in the LinearAlgebra.jl package.
KT
Это конечно чистейшая правда. А ещё можно создать identity matrix так: задать сначала матрицу n x n из нулей, а потом в цикле заполнить диагональ единицами. Но зачем?
KT
I уже сам по себе является представителем типа UniformScaling. Это значит, что I во всех разумных ситуациях уже ведёт себя как единичная матрица (левая единица, правая и т.д.). При этом не выделяя памяти на бессмысленные нули.KT
x = rand(2, 4)
I * x # вернет х
x * I # вернет х
KT
В плане линейной алгебры Julia как раз очень продуманная, без всяких там
При этом, если система переопределенная, то он решит с помощью МНК, если матрица имеет специальный тип Sym что-то там, то скорее всего Холецким и т.д., т.е. от типа матрицы А зависит то, как Julia будет искать решение, поэтому если создаем 3-диагональную матрицу (к примеру) — нужно быть настолько специфичным, насколько это возможно. Типов разных там куча для матриц
eye(5) или типа того. К примеру, решение СЛУ всё время пишется одинаково (что-то типа
A = rand(3, 5)
b = rand(3)
A \ b
При этом, если система переопределенная, то он решит с помощью МНК, если матрица имеет специальный тип Sym что-то там, то скорее всего Холецким и т.д., т.е. от типа матрицы А зависит то, как Julia будет искать решение, поэтому если создаем 3-диагональную матрицу (к примеру) — нужно быть настолько специфичным, насколько это возможно. Типов разных там куча для матриц
KT
Ну и опечатки всякие попадаются, ну и странные штуки, типа:
Я не спорю, что так можно, но почему к примеру не
???
## initalizations
psi = fill!(Vector{Float64}(undef,2), 1e2)
Я не спорю, что так можно, но почему к примеру не
psi = 1e2 * ones(2)
???
2019 February 17
AZ
Спасибо, утащил в слак
Что за слак?
OB
Aleksey Zhdanov
Что за слак?
Скорее всего ODS, но это не точно
g
Крутой у Вас препод! Заливайте задания сразу на гит и пишите сюда - возможно, народ подключится и будет дополнять интересными примерами.
Я бы пожалел студентов. :)
g
В образовательном процессе лучше использовать популярные инструменты.
g
В плане линейной алгебры Julia как раз очень продуманная, без всяких там
При этом, если система переопределенная, то он решит с помощью МНК, если матрица имеет специальный тип Sym что-то там, то скорее всего Холецким и т.д., т.е. от типа матрицы А зависит то, как Julia будет искать решение, поэтому если создаем 3-диагональную матрицу (к примеру) — нужно быть настолько специфичным, насколько это возможно. Типов разных там куча для матриц
eye(5) или типа того. К примеру, решение СЛУ всё время пишется одинаково (что-то типа
A = rand(3, 5)
b = rand(3)
A \ b
При этом, если система переопределенная, то он решит с помощью МНК, если матрица имеет специальный тип Sym что-то там, то скорее всего Холецким и т.д., т.е. от типа матрицы А зависит то, как Julia будет искать решение, поэтому если создаем 3-диагональную матрицу (к примеру) — нужно быть настолько специфичным, насколько это возможно. Типов разных там куча для матриц
Не нашел для
Symmetric оператора (\). Видимо, всё-таки сначала нужно самому разложение сделать. :)g
А кто в julia слаке сидит? Там поживее, чем тут? Есть русскоязычная комната?
A
Скорее всего ODS, но это не точно
Да
g
меня несколько расстраивает переход к общему языку по той причине что это вызовет распыление сил и не приведет к никаким результатам, так как популярности набрать не выйдет
Поддержу. Есть множество отличных языков общего назначения. Зачем ещё один? Просто потому что julia? :)
g
Vladimir
Стоит упомянуть Байеса, как тут же разгорается холивар аж на 30 сообщений в канале, что есть невиданная активность 😄
Обычно активность возникает с подачи местных философов или фантазеров. А технические вопросы редко обсуждаются (потому что на julia почти никто не пишет?). Ну, это если сравнивать с активность в @rlang_ru, например.