Julia Juno Jupyter part1 [無断転載禁止]©2ch.net

**デフォルトの名無しさん** · 2016/02/12(金) 19:05:13.83

LLVM で動く Julia のお勉強
Python に少し似てネイティブの速さで動く素敵な関数型言語
http://julialang.org/
http://junolab.org/
http://jupyter.org/

関連
Pythonista のための Julia 100 問 100 答
http://bicycle1885.hatenablog.com/entry/2014/12/23/170745

**デフォルトの名無しさん** · 2016/02/12(金) 19:06:27.32

みなさんの貢献が求められています
http://pkg.julialang.org/
http://docs.julialang.org/en/release-0.4/manual/packages/

**デフォルトの名無しさん** · 2016/02/12(金) 19:11:52.52

C の関数を呼ぶのも簡単
http://docs.julialang.org/en/release-0.4/manual/calling-c-and-fortran-code/
ツールが無くても可能ですが他にも便利なのがあるのでお好みで
Clang.jl
https://github.com/ihnorton/Clang.jl
Cpp.jl
https://github.com/timholy/Cpp.jl
Cxx.jl
https://github.com/Keno/Cxx.jl

**デフォルトの名無しさん** · 2016/02/12(金) 19:36:28.01

私が Julia を推す理由 (翻訳)
http://d.hatena.ne.jp/isobe1978/20140330/1396201066
原文
http://www.evanmiller.org/why-im-betting-on-julia.html
Julia の開発環境を構築してみる (IJulia)
http://qiita.com/miiton@github/items/33c2b62c4911216d5fdb
科学計算が得意な言語 Julia を試してみた
http://qiita.com/chezou/items/d090f26dcb31818d6964
Gadfly
http://gadflyjl.org/
Plotly
https://plot.ly/julia/

**デフォルトの名無しさん** · 2016/02/12(金) 23:29:57.13

実は数年前にここでjuliaスレ立てたんだが落ちたんだよなあ。。。。

**デフォルトの名無しさん** · 2016/02/13(土) 10:10:05.58

漏れもなんとなく見覚えがある気がするんだよなぁ・・・

落ちる条件ってなんだっけ

**デフォルトの名無しさん** · 2016/02/13(土) 10:11:00.83

という訳で保守がてら

juliaからCの関数を呼んでみた

julia> ccall(:printf, Void, (Ptr{Cchar}, Ptr{Cchar}), "%s", "abc")
abc
julia> ccall(:printf, Void, (Ptr{Cchar}, Ptr{Cchar}, Int), "%s%d", "abc", 9)
abc9

引数の数が自由なのがいいねっ！

**デフォルトの名無しさん** · 2016/02/13(土) 10:12:28.35

おれおれDLLだと呼び出し方がちょっと変わるみたい
(最近仕様が変わったらしいので古いドキュソの通りやると動かない
っていうか最新のドキュソもまだ動かないサンプルのままだぞω)

julia> libcurl = Base.Libdl.dlopen_e("./libcurl.dll")
Ptr{Void} @0x2c76e2c0
julia> curl_version = Base.Libdl.dlsym(libcurl, :curl_version)
Ptr{Void} @0x6b25c740
julia> curl_getenv = Base.Libdl.dlsym(libcurl, :curl_getenv)
Ptr{Void} @0x6b25c854
julia> curl_free = Base.Libdl.dlsym(libcurl, :curl_free)
Ptr{Void} @0x6b25cbf8
julia> p = ccall(curl_version, Ptr{Cchar}, ())
Ptr{Int8} @0x6b29c060
julia> ccall(:printf, Void, (Ptr{Cchar}, Ptr{Cchar}), "%s", p)
libcurl/7.30.0 OpenSSL/1.0.1e zlib/1.2.7
julia> q = ccall(curl_getenv, Ptr{Cchar}, (Ptr{Cchar}, ), "NUMBER_OF_PROCESSORS")
Ptr{Int8} @0x2c76df10
julia> ccall(:printf, Void, (Ptr{Cchar}, Ptr{Cchar}), "%s", q)
8
julia> ccall(curl_free, Cuint, (Ptr{Void}, ), q)
0x00000001
julia> Base.Libdl.dlclose(libcurl)

うっかり解放忘れそう

**デフォルトの名無しさん** · 2016/02/13(土) 10:14:00.39

stdcallな関数だとさらにこんな感じだよ

julia> user32 = Base.Libdl.dlopen_e("user32.dll")
Ptr{Void} @0x2c767300
julia> MessageBoxA = Base.Libdl.dlsym(user32, :MessageBoxA)
Ptr{Void} @0x76c6cf50
julia> err = ccall(MessageBoxA, stdcall, UInt, (Ptr{UInt}, Ptr{Cchar}, Ptr{Cchar}, UInt), 0, "Hello, work!", "Message", 0)
0x00000001

wchar_tはArray{UInt16,1}らしい(UTF8を表示したければ先にwchar_tへの変換が必要)

julia> MessageBoxW = Base.Libdl.dlsym(user32, :MessageBoxW)
Ptr{Void} @0x76c6d280
julia> err = ccall(MessageBoxW, stdcall, UInt, (Ptr{UInt}, Ptr{Cwchar_t}, Ptr{Cwchar_t}, UInt), 0, Array{UInt16,1}([0x304C,0x3063,0xFF01,0]), Array{UInt16,1}([0x306C,0x308B,0x307D,0]), 0)
0x00000001
julia> Base.Libdl.dlclose(user32)

解放解放

引数や戻り値の定義を呼び出しのたびに何度も書くの面倒なのでマクロを書くのが良いらしいよ

**デフォルトの名無しさん** · 2016/02/13(土) 10:17:20.41

>っていうか最新のドキュソもまだ動かないサンプルのままだぞω

ああごめん動くわ

**デフォルトの名無しさん** · 2016/02/13(土) 10:22:27.58

ccall(:printf, Void, (Ptr{Cchar}, Ptr{Cchar}), "%s", ccall((:curl_version, :libcurl), Ptr{Cchar}, ()))

err = ccall((:MessageBoxW, :user32), stdcall, UInt, (Ptr{UInt}, Ptr{Cwchar_t}, Ptr{Cwchar_t}, UInt), 0, Array{UInt16,1}([0x304C,0x3063,0xFF01,0]), Array{UInt16,1}([0x306C,0x308B,0x307D,0]), 0)

Array{UInt16,1}のところはArray{Cwchar_t}とかVector{Cwchar_t}とも書ける

err = ccall((:MessageBoxW, :user32), stdcall, UInt, (Ptr{UInt}, Ptr{Cwchar_t}, Ptr{Cwchar_t}, UInt), 0, Vector{Cwchar_t}([0x304C,0x3063,0xFF01,0]), Array{Cwchar_t}([0x306C,0x308B,0x307D,0]), 0)

これでいいね
loadは省けるけど毎回loadと解放してくれてるんだろうか

**デフォルトの名無しさん** · 2016/02/13(土) 10:31:38.04

>>8
>OpenSSL/1.0.1e

ワロス

**デフォルトの名無しさん** · 2016/02/13(土) 11:59:14.88

もう少しちらうら
julia> bytes2hex([65,97,66,98,0x5c])
ERROR: MethodError: `bytes2hex` has no method matching bytes2hex(::Array{Int32,1})
julia> bytes2hex(Array{UInt8,1}([65,97,66,98,0x5c]))
"416142625c"
julia> bytestring(Array{UInt8,1}([65,97,66,98,0x5c]))
"AaBb\\"
julia> hex2bytes("416142625c")
5-element Array{UInt8,1}: 0x41 0x61 0x42 0x62 0x5c
julia> string2bytes("AaBb\\")
ERROR: UndefVarError: string2bytes not defined
julia> stringbytes("AaBb\\")
ERROR: UndefVarError: stringbytes not defined
julia> read(IOBuffer("AaBb\\"))
ERROR: MethodError: `read` has no method matching read(::Base.AbstractIOBuffer{Array{UInt8,1}})
julia> read(IOBuffer("AaBb\\"), UInt8)
0x41
julia> read(IOBuffer("AaBb\\"), UInt8, length("AaBb\\"))
5-element Array{UInt8,1}: 0x41 0x61 0x42 0x62 0x5c

んー
マニュアルにはreaddim()があるって書いてあるのに・・・
何かimport忘れてるとか？
length省略できんのもなんだかな
勘違いならいいけど

**デフォルトの名無しさん** · 2016/02/13(土) 17:12:04.67

ちょっと前から同名のAV女優が出てきて検索がちょっとアレになったよね
運が悪い

**デフォルトの名無しさん** · 2016/02/13(土) 17:45:27.64

"どうすればPythonをJuliaと同じくらい速く動かせるのか？ : 様々なやり方で計算の高速化を図る"
http://postd.cc/python_meets_julia_micro_performance/

**デフォルトの名無しさん** · 2016/02/13(土) 22:04:46.56

フォートランっぽい文法狂おしいほど嫌い

**デフォルトの名無しさん** · 2016/02/14(日) 11:00:37.07

パッケージがおれおれやってみたばっかりで
まともなものが少ねーじゃねーか
これ PyPI よりひでーぞ
github っつーかオプソは害だな

**デフォルトの名無しさん** · 2016/02/14(日) 11:01:38.90

>>16
21世紀にもなって添え字が 1 から始まるとかもうね

**デフォルトの名無しさん** · 2016/02/14(日) 11:11:08.67

おまいらなんだかんだ言って
みんなそれなりにやってんだな
ﾆﾔﾆﾔ

**デフォルトの名無しさん** · 2016/02/14(日) 19:24:28.30

>1 乙
>14 前スレで見た覚えがある
>16 それな

**デフォルトの名無しさん** · 2016/02/15(月) 10:16:09.06

コンパイル中のエラーよりも実行時の型エラーの方が厄介な言語だなこいつは。
慣れれば問題ないが。

**デフォルトの名無しさん** · 2016/02/15(月) 10:20:41.18

>>13
非推奨っぽいが、
これでいいらしい。

"AaBb\\".data

代入も出来るようだ。

s = "AaBb\\"
s.data[2] = UInt8(64)
s

**デフォルトの名無しさん** · 2016/02/15(月) 10:21:13.05

これは割りと嬉しいかも。

julia> @sprintf "ab [%7s]" "xあx"
"ab [ xあx]"

julia> @sprintf "ab [%-7s]" "xあx"
"ab [xあx ]"

julia> @sprintf "ab [%7s]" "xxx"
"ab [ xxx]"

julia> @sprintf "ab [%-7s]" "xxx"
"ab [xxx ]"

**デフォルトの名無しさん** · 2016/02/15(月) 10:31:07.93

うむむ
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes(""))
"ab [ ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141E381824141"))
"ab [AAあAA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141E981824141"))
"ab [AA遂AA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141Ee81824141"))
"ab [AA\ue042AA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141Ef81824141"))
"ab [AA\uf042AA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141E881824141"))
"ab [AA聂AA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141E8e8824141"))
"ab [AA騂AA ]"
julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("4141E9e8824141"))
"ab [AAꨂAA ]"

julia> @sprintf "ab [%-7s]" bytestring(hex2bytes("414155aa554141"))
ERROR: UnicodeError: invalid character index
in next at unicode/utf8.jl:65
in strwidth at strings/basic.jl:205
[inlined code] from printf.jl:159
in anonymous at no file:0

julia> @sprintf "ab [%-7s]" hex2bytes("414155aa554141")
"ab [UInt8[0x41,0x41,0x55,0xaa,0x55,0x41,0x41]]"

**デフォルトの名無しさん** · 2016/02/15(月) 11:19:37.91

>>11
julia> ccall((:curl_version, :libcurl), Ptr{Cchar}, ())
Ptr{Int8} @0x6b29c060

julia> @sprintf "%s" ccall((:curl_version, :libcurl), Ptr{Cchar}, ())
"Ptr{Int8} @0x6b29c060"

んーこれも期待した動作と違うなぁ。

**デフォルトの名無しさん** · 2016/02/15(月) 11:28:27.71

>>16
ifに括弧なしは一見便利だけどend強制されるのはPythonに慣れてるとうざく感じる。

julia> Array{UInt8}([1:n;]) == map(a -> UInt8(a), 1:n)
true

ちょっと気持ち悪いところもある。

**デフォルトの名無しさん** · 2016/02/15(月) 11:32:12.26

julia> Array{UInt8}([1:n;]) == map(a -> UInt8(a), (1:n))
true

julia> Array{UInt8}([1:n]) == map(a -> UInt8(a), (1:n))
警告出るが
true
(二回目からは警告出なくなる？)

julia> Array{UInt8}((1:n)) == map(a -> UInt8(a), (1:n))
エラー

**デフォルトの名無しさん** · 2016/02/15(月) 11:45:55.52

>>25
bytestring(ccall((:curl_version, :libcurl), Ptr{Cchar}, ()))

**デフォルトの名無しさん** · 2016/02/15(月) 11:49:21.03

>>28
有賀㌧。
>>24でも書いたけどbytestringってこんなデータ来ると死ぬから使いたくないなぁ。
julia> bytestring(hex2bytes("414155aa554141"))

**デフォルトの名無しさん** · 2016/02/15(月) 12:16:36.12

julia> pointer_to_array(ccall((:curl_version, :libcurl), Ptr{Cchar}, ()), n)
でnバイト分のArray{Int8}になるみたいだけどなぜかデフォでUInt8じゃないのが気に入らないなぁ。

julia> Array{UInt8}(pointer_to_array(ccall((:curl_version, :libcurl), Ptr{Cchar}, ()), n))
とかいちいち書くの面倒だなぁ。

**デフォルトの名無しさん** · 2016/02/15(月) 12:19:08.24

>>30
Cchar == UInt8 # false
Cchar == Int8 # true

**デフォルトの名無しさん** · 2016/02/15(月) 12:21:09.84

>>31
julia> pointer_to_array(ccall((:curl_version, :libcurl), Ptr{UInt8}, ()), n)
できました！！
ありがとうございます！！！

**デフォルトの名無しさん** · 2016/02/15(月) 18:04:10.07

意味わからずbytestring多用するのは危険

**デフォルトの名無しさん** · 2016/02/16(火) 11:03:22.23

今はまだ julia 0.3.2 対応みたいだけど web から試せるよ
https://try.jupyter.org/

IJulia をインストールすれば自分のサイトでもできるよ
https://github.com/JuliaLang/IJulia.jl

**デフォルトの名無しさん** · 2016/02/16(火) 11:49:26.30

>>14
ちょっと前から声が大きいだけで連投するうざいのが住み着いて2ちゃんはちょっとアレになったよね
筋が悪い

**デフォルトの名無しさん** · 2016/02/16(火) 12:08:19.40

日本語サイトしか見ない病だ

**デフォルトの名無しさん** · 2016/02/16(火) 17:07:24.65

ツイッターでjulia検索してもひどい
最悪だ

**デフォルトの名無しさん** · 2016/02/16(火) 17:38:07.56

英語だが Julia をつかった科学計算を丁寧に解説しているサイト
ちなみに著者の一人 Sargent はノーベル経済学賞をとっている

http://quant-econ.net/jl/

**デフォルトの名無しさん** · 2016/02/17(水) 08:55:28.12

Python版はこちら
http://quant-econ.net/py/