Julia is a high-level, high-performance dynamic programming language for technical computing, with syntax that is familiar to users of other technical

15794 ワード

http://julialang.org/
julia  |  source  |  downloads  |  docs  |  blog  |  community  |  teaching  |  publications  |  gsoc  |  juliacon  |  rss
Julia is a high-level, high-performance dynamic programming language for technical computing, with syntax that is familiar to users of other technical computing environments. It provides a sophisticated compiler,  distributed parallel execution , numerical accuracy, and an  extensive mathematical function library . The library, largely written in Julia itself, also integrates mature, best-of-breed C and Fortran libraries for  linear algebrarandom number generationsignal processing , and  string processing . In addition, the Julia developer community is contributing a number of  external packages  through Julia’s built-in package manager at a rapid pace.  IJulia , a collaboration between the  IPython  and Julia communities, provides a powerful browser-based graphical notebook interface to Julia.
Julia programs are organized around  multiple dispatch ; by defining functions and overloading them for different combinations of argument types, which can also be user-defined. For a more in-depth discussion of the rationale and advantages of Julia over other systems, see the following highlights or read the  introduction  in the  online manual .
A Summary of Features
  • Multiple dispatch : providing ability to define function behavior across many combinations of argument types
  • Dynamic type system: types for documentation, optimization, and dispatch
  • Good performance, approaching that of statically-compiled languages like C
  • Built-in package manager
  • Lisp-like macros  and other  metaprogramming facilities
  • Call Python functions: use the  PyCall  package
  • Call C functions  directly: no wrappers or special APIs
  • Powerful shell-like capabilities for  managing other processes
  • Designed for  parallelism and distributed computation
  • Coroutines : lightweight “green” threading
  • User-defined types are as fast and compact as built-ins
  • Automatic generation of efficient, specialized code for different argument types
  • Elegant and extensible conversions and promotions for numeric and other types
  • Efficient support for  Unicode , including but not limited to  UTF-8
  • MIT licensed : free and open source

  • High-Performance JIT Compiler
    Julia’s LLVM-based just-in-time (JIT) compiler combined with the language’s design allow it to approach and often match the performance of C. To get a sense of relative performance of Julia compared to other languages that can or could be used for numerical and scientific computing, we’ve written a small set of micro-benchmarks in a variety of languages:  CFortran , JuliaPythonMatlab/OctaveRJavaScriptGo , and  Mathematica . We encourage you to skim the code to get a sense for how easy or difficult numerical programming in each language is. The following micro-benchmark results were obtained on a single core (serial execution) on an Intel® Xeon® CPU E7-8850 2.00GHz CPU with 1TB of 1067MHz DDR3 RAM, running Linux:
     
    Fortran
    Julia
    Python
    R
    Matlab
    Octave
    Mathe-matica
    JavaScript
    Go
     
    gcc 4.8.1
    0.2
    2.7.3
    3.0.2
    R2012a
    3.6.4
    8.0
    V8 3.7.12.22
    go1
    fib
    0.26
    0.91
    30.37
    411.36
    1992.00
    3211.81
    64.46
    2.18
    1.03
    parse_int
    5.03
    1.60
    13.95
    59.40
    1463.16
    7109.85
    29.54
    2.43
    4.79
    quicksort
    1.11
    1.14
    31.98
    524.29
    101.84
    1132.04
    35.74
    3.51
    1.25
    mandel
    0.86
    0.85
    14.19
    106.97
    64.58
    316.95
    6.07
    3.49
    2.36
    pi_sum
    0.80
    1.00
    16.33
    15.42
    1.29
    237.41
    1.32
    0.84
    1.41
    rand_mat_stat
    0.64
    1.66
    13.52
    10.84
    6.61
    14.98
    4.52
    3.28
    8.12
    rand_mat_mul
    0.96
    1.01
    3.41
    3.98
    1.10
    3.41
    1.16
    14.60
    8.51
    Figure: benchmark times relative to C (smaller is better, C performance = 1.0).
    C compiled by gcc 4.8.1, taking best timing from all optimization levels (-O0 through -O3). C, Fortran and Julia use  OpenBLAS  v0.2.8. The Python implementations of rand_mat_stat and rand_mat_mul use NumPy (v1.6.1) functions; the rest are pure Python implementations.Benchmarks can also be seen  here as a plot  created with  Gadfly .
    These benchmarks, while not comprehensive, do test compiler performance on a range of common code patterns, such as function calls, string parsing, sorting, numerical loops, random number generation, and array operations. It is important to note that these benchmark implementations are not written for absolute maximal performance (the fastest code to compute fib(20)  is the constant literal  6765 ). Rather, all of the benchmarks are written to test the performance of specific algorithms, expressed in a reasonable idiom in each language. In particular, all languages use the same algorithm: the Fibonacci benchmarks are all recursive while the pi summation benchmarks are all iterative; the “algorithm” for random matrix multiplication is to call LAPACK, except where that’s not possible, such as in JavaScript. The point of these benchmarks is to compare the performance of specific algorithms across language implementations, not to compare the fastest means of computing a result, which in most high-level languages relies on calling C code. Raw benchmark numbers in CSV format are available  here . Various other Julia benchmarks, tracked over time can be found at  speed.julialang.org .
    To give a quick taste of what Julia looks like, here is the code used in the Mandelbrot and random matrix statistics benchmarks:
    function mandel(z)
        c = z
        maxiter = 80
        for n = 1:maxiter
            if abs(z) > 2
                return n-1
            end
            z = z^2 + c
        end
        return maxiter
    end
    
    function randmatstat(t)
        n = 5
        v = zeros(t)
        w = zeros(t)
        for i = 1:t
            a = randn(n,n)
            b = randn(n,n)
            c = randn(n,n)
            d = randn(n,n)
            P = [a b c d]
            Q = [a b; c d]
            v[i] = trace((P.'*P)^4)
            w[i] = trace((Q.'*Q)^4)
        end
        std(v)/mean(v), std(w)/mean(w)
    end
    

    The code above is quite clear, and should feel familiar to anyone who has programmed in other mathematical languages. The Julia implementation of  randmatstat  is considerably simpler than the equivalent  C implementation , without giving up much performance. Planned compiler optimizations will close this performance gap in the future. By design, Julia allows you to range from tight low-level loops, up to a high-level programming style, while sacrificing some performance, but gaining the ability to express complex algorithms easily. This continuous spectrum of programming levels is a hallmark of the Julia approach to programming and is very much an intentional feature of the language.
    Designed for Parallelism and Cloud Computing
    Julia does not impose any particular style of parallelism on the user. Instead, it provides a number of  key building blocks for distributed computation , making it flexible enough to support a number of styles of parallelism, and allowing users to add more. The following simple example demonstrates how to count the number of heads in a large number of coin tosses in parallel.
    nheads = @parallel (+) for i=1:100000000
      int(randbool())
    end
    

    This computation is automatically distributed across all available compute nodes, and the result, reduced by summation ( + ), is returned at the calling node.
    Here is a screenshot of a web-based interactive  IJulia  session, using  Gadfly  to produce various plots with  D3  as a rendering backend in the browser (SVG, PDF, PNG and various other backends are also supported):
    Julia is a high-level, high-performance dynamic programming language for technical computing, with syntax that is familiar to users of other technical_第1张图片
    This paves the way for fully cloud-based operation, including data management, code editing and sharing, execution, debugging, collaboration, analysis, data exploration, and visualization. The eventual goal is to let people stop worrying about administering machines and managing data and get straight to the real problem.
    Free, Open Source and Library-Friendly
    The core of the Julia implementation is licensed under the  MIT license . Various libraries used by the Julia environment include their own licenses such as the  GPLLGPL , and  BSD  (therefore the environment, which consists of the language, user interfaces, and libraries, is under the GPL). The language can be built as a shared library, so users can combine Julia with their own C/Fortran code or proprietary third-party libraries. Furthermore, Julia makes it  simple to call external functions  in C and Fortran shared libraries, without writing any wrapper code or even recompiling existing code. You can try calling external library functions directly from Julia’s interactive prompt, getting immediate feedback. See  LICENSE  for the full terms of Julia’s licensing.