Python\u306e\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u306b\u60a9\u3080\u30c7\u30d9\u30ed\u30c3\u30d1\u30fc\u306b\u6717\u5831\uff01numba\u3092\u4f7f\u3048\u3070\u3001\u7c21\u5358\u306b\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3067\u304d\u307e\u3059\u3002\u3053\u306e\u8a18\u4e8b\u3067\u306f\u3001numba\u306e\u57fa\u672c\u304b\u3089\u5b9f\u8df5\u7684\u306a\u30c6\u30af\u30cb\u30c3\u30af\u307e\u3067\u3001\u8a73\u3057\u304f\u89e3\u8aac\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u76ee\u6b21 <\/p>\n
Python\u306f\u3001\u30b7\u30f3\u30d7\u30eb\u3067\u66f8\u304d\u3084\u3059\u3044\u8a00\u8a9e\u3068\u3057\u3066\u77e5\u3089\u308c\u3066\u3044\u307e\u3059\u304c\u3001\u5b9f\u884c\u901f\u5ea6\u306e\u9762\u3067\u306fC\u8a00\u8a9e\u306a\u3069\u306b\u53ca\u3073\u307e\u305b\u3093\u3002\u305d\u3053\u3067\u767b\u5834\u3059\u308b\u306e\u304c\u3001Python\u306eJIT\u30b3\u30f3\u30d1\u30a4\u30e9\u300cnumba\u300d\u3067\u3059\u3002numba\u3092\u4f7f\u3048\u3070\u3001Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u306a\u30de\u30b7\u30f3\u30b3\u30fc\u30c9\u306b\u5909\u63db\u3057\u3001\u5b9f\u884c\u901f\u5ea6\u3092\u5927\u5e45\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u306f\u3001\u79d1\u5b66\u8a08\u7b97\u3084\u30c7\u30fc\u30bf\u5206\u6790\u3067\u3088\u304f\u4f7f\u308f\u308c\u308bNumPy\u3068\u3082\u7dca\u5bc6\u306b\u9023\u643a\u3057\u3066\u304a\u308a\u3001NumPy\u306e\u30c7\u30fc\u30bf\u578b\u3084\u95a2\u6570\u3092\u305d\u306e\u307e\u307e\u6d3b\u7528\u3067\u304d\u307e\u3059\u3002\u307e\u305f\u3001CPU\u3060\u3051\u3067\u306a\u304fGPU\u3067\u3082\u52d5\u4f5c\u3057\u3001\u4e26\u5217\u51e6\u7406\u306b\u3082\u5bfe\u5fdc\u3057\u3066\u3044\u308b\u305f\u3081\u3001\u5927\u898f\u6a21\u306a\u6570\u5024\u8a08\u7b97\u3092\u9ad8\u901f\u306b\u51e6\u7406\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u3092\u4f7f\u3046\u306b\u306f\u3001\u307e\u305apip\u3067\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
pip install numba<\/pre>\n\n\n\n\u6b21\u306b\u3001\u9ad8\u901f\u5316\u3057\u305f\u3044\u95a2\u6570\u306b@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4ed8\u3051\u308b\u3060\u3051\u3067\u3001JIT\u30b3\u30f3\u30d1\u30a4\u30eb\u304c\u6709\u52b9\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import jit\n\n@jit\ndef sum_squares(n):\n result = 0\n for i in range(n):\n result += i * i\n return result<\/pre>\n\n\n\n\u3053\u306e\u4f8b\u3067\u306f\u30011\u304b\u3089(n-1)\u307e\u3067\u306e\u6574\u6570\u306e\u4e8c\u4e57\u548c\u3092\u8a08\u7b97\u3059\u308b\u95a2\u6570sum_squares()\u3092\u5b9a\u7fa9\u3057\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3067JIT\u30b3\u30f3\u30d1\u30a4\u30eb\u3057\u3066\u3044\u307e\u3059\u3002<\/p>\n\n\n\n
\u3055\u3089\u306b\u3001nopython\u30e2\u30fc\u30c9\u3092\u4f7f\u3048\u3070\u3001Python\u30a4\u30f3\u30bf\u30d7\u30ea\u30bf\u3092\u4ecb\u3055\u305a\u306b\u76f4\u63a5\u30de\u30b7\u30f3\u30b3\u30fc\u30c9\u3092\u751f\u6210\u3067\u304d\u307e\u3059\u3002\u3053\u308c\u306b\u3088\u308a\u3001\u3088\u308a\u9ad8\u901f\u306a\u5b9f\u884c\u304c\u53ef\u80fd\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
@jit(nopython=True)\ndef sum_squares(n):\n result = 0\n for i in range(n):\n result += i * i\n return result<\/pre>\n\n\n\nnumba\u3092\u4f7f\u3048\u3070\u3001Python\u306e\u30b3\u30fc\u30c9\u3092\u307b\u3068\u3093\u3069\u5909\u66f4\u305b\u305a\u306b\u9ad8\u901f\u5316\u3067\u304d\u308b\u305f\u3081\u3001\u958b\u767a\u52b9\u7387\u3092\u7dad\u6301\u3057\u3064\u3064\u3001\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u3092\u5927\u5e45\u306b\u6539\u5584\u3067\u304d\u307e\u3059\u3002\u6b21\u9805\u3067\u306f\u3001numba\u3068Python\u306e\u4ee3\u8868\u7684\u306a\u6570\u5024\u8a08\u7b97\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3042\u308bnumpy\u306e\u9055\u3044\u306b\u3064\u3044\u3066\u8a73\u3057\u304f\u898b\u3066\u3044\u304d\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
numba\u3068numpy\u306e\u9055\u3044\u3092\u7406\u89e3\u3057\u3088\u3046<\/h2>\n\n\n\n
numba\u3068numpy\u306f\u3069\u3061\u3089\u3082Python\u306e\u6570\u5024\u8a08\u7b97\u3092\u9ad8\u901f\u5316\u3059\u308b\u305f\u3081\u306e\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3059\u304c\u3001\u3044\u304f\u3064\u304b\u306e\u91cd\u8981\u306a\u9055\u3044\u304c\u3042\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
numpy\u306f\u30e9\u30a4\u30d6\u30e9\u30ea\u3001numba\u306f\u30b3\u30f3\u30d1\u30a4\u30e9<\/h3>\n\n\n\n
numpy\u306fPython\u306e\u30e9\u30a4\u30d6\u30e9\u30ea\u306e\u4e00\u3064\u3067\u3001\u4e3b\u306b\u914d\u5217\u51e6\u7406\u306b\u7279\u5316\u3057\u3066\u3044\u307e\u3059\u3002\u4e00\u65b9\u3001numba\u306fPython\u306eJIT\u30b3\u30f3\u30d1\u30a4\u30e9\u3067\u3042\u308a\u3001Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u306a\u30de\u30b7\u30f3\u30b3\u30fc\u30c9\u306b\u5909\u63db\u3059\u308b\u3053\u3068\u3067\u3001\u6570\u5024\u8a08\u7b97\u5168\u822c\u3092\u9ad8\u901f\u5316\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u4f8b\u3048\u3070\u3001numpy\u3092\u4f7f\u3063\u3066\u914d\u5217\u306e\u52a0\u7b97\u3092\u884c\u3046\u5834\u5408\u306f\u6b21\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
import numpy as np\n\na = np.arange(10000)\nb = np.arange(10000)\nc = a + b<\/pre>\n\n\n\n\u3053\u308c\u306b\u5bfe\u3057\u3001numba\u3092\u4f7f\u3063\u305f\u5834\u5408\u306f\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u3063\u3066\u95a2\u6570\u3092\u30b3\u30f3\u30d1\u30a4\u30eb\u3057\u3001\u305d\u306e\u95a2\u6570\u5185\u3067\u914d\u5217\u3092\u52a0\u7b97\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import jit\nimport numpy as np\n\n@jit(nopython=True)\ndef add_arrays(x, y):\n return x + y\n\na = np.arange(10000)\nb = np.arange(10000)\nc = add_arrays(a, b)<\/pre>\n\n\n\n\u5f97\u610f\u3068\u3059\u308b\u51e6\u7406\u304c\u7570\u306a\u308b \u2013 numba\u306f\u6570\u5024\u8a08\u7b97\u306b\u7279\u5316<\/h3>\n\n\n\n
numpy\u306f\u914d\u5217\u51e6\u7406\u306b\u7279\u5316\u3057\u3066\u3044\u308b\u305f\u3081\u3001\u914d\u5217\u306e\u8981\u7d20\u3054\u3068\u306e\u6f14\u7b97\u3084\u884c\u5217\u8a08\u7b97\u306a\u3069\u3092\u9ad8\u901f\u306b\u51e6\u7406\u3067\u304d\u307e\u3059\u3002\u4e00\u65b9\u3001numba\u306f\u6570\u5024\u8a08\u7b97\u5168\u822c\u3092\u9ad8\u901f\u5316\u3067\u304d\u308b\u6c4e\u7528\u6027\u306e\u9ad8\u3044\u30c4\u30fc\u30eb\u3067\u3059\u3002<\/p>\n\n\n\n
\u307e\u305f\u3001numpy\u3067\u306f\u624b\u52d5\u3067\u306e\u30d9\u30af\u30c8\u30eb\u5316\u304c\u5fc5\u8981\u306a\u5834\u5408\u304c\u3042\u308a\u307e\u3059\u304c\u3001numba\u3067\u306f\u81ea\u52d5\u7684\u306b\u30d9\u30af\u30c8\u30eb\u5316\u3067\u304d\u307e\u3059\u3002\u3055\u3089\u306b\u3001numba\u306fGPU\u3067\u306e\u8a08\u7b97\u306b\u3082\u5bfe\u5fdc\u3057\u3066\u3044\u308b\u305f\u3081\u3001\u5927\u898f\u6a21\u306a\u6570\u5024\u8a08\u7b97\u3092\u975e\u5e38\u306b\u9ad8\u901f\u306b\u51e6\u7406\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u3001numba\u3068numpy\u306f\u6570\u5024\u8a08\u7b97\u3092\u9ad8\u901f\u5316\u3059\u308b\u3068\u3044\u3046\u70b9\u3067\u306f\u5171\u901a\u3057\u3066\u3044\u307e\u3059\u304c\u3001\u30e9\u30a4\u30d6\u30e9\u30ea\u3068\u30b3\u30f3\u30d1\u30a4\u30e9\u3068\u3044\u3046\u6839\u672c\u7684\u306a\u9055\u3044\u304c\u3042\u308a\u3001\u5f97\u610f\u3068\u3059\u308b\u51e6\u7406\u3082\u7570\u306a\u308a\u307e\u3059\u3002numba\u306e\u7279\u5fb4\u3092\u7406\u89e3\u3059\u308b\u3053\u3068\u3067\u3001Python\u3067\u306e\u6570\u5024\u8a08\u7b97\u3092\u3055\u3089\u306b\u52b9\u7387\u5316\u3067\u304d\u308b\u3067\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
\u3010\u30c6\u30af\u30cb\u30c3\u30af1\u3011@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3067\u95a2\u6570\u3092\u30b3\u30f3\u30d1\u30a4\u30eb<\/h2>\n\n\n\n
numba\u3092\u4f7f\u3063\u3066Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3059\u308b\u969b\u306b\u6b20\u304b\u305b\u306a\u3044\u306e\u304c\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3067\u3059\u3002@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u95a2\u6570\u306e\u524d\u306b\u4ed8\u3051\u308b\u3060\u3051\u3067\u3001\u305d\u306e\u95a2\u6570\u304cJIT\u30b3\u30f3\u30d1\u30a4\u30eb\u3055\u308c\u3001\u9ad8\u901f\u306b\u5b9f\u884c\u3055\u308c\u308b\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306e\u4f7f\u3044\u65b9<\/h3>\n\n\n\n
@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306e\u57fa\u672c\u7684\u306a\u4f7f\u3044\u65b9\u306f\u3001\u975e\u5e38\u306b\u30b7\u30f3\u30d7\u30eb\u3067\u3059\u3002\u9ad8\u901f\u5316\u3057\u305f\u3044\u95a2\u6570\u306e\u524d\u306b@jit\u3092\u4ed8\u3051\u308b\u3060\u3051\u3067\u3059\u3002<\/p>\n\n\n\n
from numba import jit\n\n@jit\ndef sum_squares(n):\n result = 0\n for i in range(n):\n result += i * i\n return result<\/pre>\n\n\n\n\u3053\u306e\u4f8b\u3067\u306f\u30011\u304b\u3089(n-1)\u307e\u3067\u306e\u6574\u6570\u306e\u4e8c\u4e57\u548c\u3092\u8a08\u7b97\u3059\u308b\u95a2\u6570sum_squares()\u3092\u5b9a\u7fa9\u3057\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4ed8\u3051\u3066\u3044\u307e\u3059\u3002\u3053\u308c\u306b\u3088\u308a\u3001sum_squares()\u95a2\u6570\u304cJIT\u30b3\u30f3\u30d1\u30a4\u30eb\u3055\u308c\u3001\u9ad8\u901f\u306b\u5b9f\u884c\u3055\u308c\u308b\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
\u30b3\u30f3\u30d1\u30a4\u30eb\u306e\u30aa\u30d7\u30b7\u30e7\u30f3\u5f15\u6570 \u2013 nopython\u3068cache<\/h3>\n\n\n\n
@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306b\u306f\u3001\u3044\u304f\u3064\u304b\u306e\u4fbf\u5229\u306a\u30aa\u30d7\u30b7\u30e7\u30f3\u5f15\u6570\u304c\u3042\u308a\u307e\u3059\u3002\u3088\u304f\u4f7f\u308f\u308c\u308b\u306e\u304c\u3001nopython\u3068cache\u3067\u3059\u3002<\/p>\n\n\n\n
nopython\u30e2\u30fc\u30c9\u3092\u4f7f\u3046\u3068\u3001Python\u30a4\u30f3\u30bf\u30d7\u30ea\u30bf\u3092\u4ecb\u3055\u305a\u306b\u30cd\u30a4\u30c6\u30a3\u30d6\u30b3\u30fc\u30c9\u3092\u751f\u6210\u3067\u304d\u307e\u3059\u3002\u3053\u308c\u306b\u3088\u308a\u3001\u6700\u9ad8\u306e\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u304c\u5f97\u3089\u308c\u307e\u3059\u3002<\/p>\n\n\n\n
@jit(nopython=True)\ndef sum_squares_nopython(n):\n result = 0\n for i in range(n):\n result += i * i\n return result<\/pre>\n\n\n\ncache\u30aa\u30d7\u30b7\u30e7\u30f3\u3092\u4f7f\u3046\u3068\u3001\u30b3\u30f3\u30d1\u30a4\u30eb\u7d50\u679c\u304c\u30c7\u30a3\u30b9\u30af\u306b\u30ad\u30e3\u30c3\u30b7\u30e5\u3055\u308c\u3001\u6b21\u56de\u4ee5\u964d\u306e\u5b9f\u884c\u6642\u306b\u518d\u5229\u7528\u3055\u308c\u307e\u3059\u3002\u3053\u308c\u306b\u3088\u308a\u3001\u30b3\u30f3\u30d1\u30a4\u30eb\u306e\u6642\u9593\u3092\u7bc0\u7d04\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
@jit(cache=True)\ndef sum_squares_cached(n):\n result = 0\n for i in range(n):\n result += i * i\n return result<\/pre>\n\n\n\n\u3053\u306e\u3088\u3046\u306b\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001\u7c21\u5358\u306bPython\u30b3\u30fc\u30c9\u3092JIT\u30b3\u30f3\u30d1\u30a4\u30eb\u3057\u3001\u9ad8\u901f\u5316\u3059\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002\u307e\u305f\u3001nopython\u3084cache\u306e\u30aa\u30d7\u30b7\u30e7\u30f3\u5f15\u6570\u3092\u4f7f\u3048\u3070\u3001\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316\u304c\u671f\u5f85\u3067\u304d\u308b\u3067\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
\u6b21\u9805\u3067\u306f\u3001numba\u306e\u578b\u6307\u5b9a\u3092\u4f7f\u3063\u3066\u3001\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316\u3092\u5b9f\u73fe\u3059\u308b\u65b9\u6cd5\u306b\u3064\u3044\u3066\u89e3\u8aac\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u3010\u30c6\u30af\u30cb\u30c3\u30af2\u3011numba\u306e\u578b\u6307\u5b9a\u3067\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316<\/h2>\n\n\n\n
numba\u3092\u4f7f\u3063\u3066Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3059\u308b\u969b\u306b\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3068\u4e26\u3093\u3067\u91cd\u8981\u306a\u306e\u304c\u578b\u6307\u5b9a\u3067\u3059\u3002\u95a2\u6570\u306e\u5f15\u6570\u3068\u623b\u308a\u5024\u306b\u578b\u3092\u6307\u5b9a\u3059\u308b\u3053\u3068\u3067\u3001\u30b3\u30f3\u30d1\u30a4\u30e9\u304c\u3088\u308a\u6700\u9069\u5316\u3055\u308c\u305f\u30b3\u30fc\u30c9\u3092\u751f\u6210\u3067\u304d\u3001\u3055\u3089\u306a\u308b\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u5411\u4e0a\u304c\u671f\u5f85\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
numpy\u306edtype\u306b\u5bfe\u5fdc\u3057\u305f\u578b\u6307\u5b9a\u306e\u65b9\u6cd5<\/h3>\n\n\n\n
numba\u306e\u578b\u6307\u5b9a\u306b\u306f\u3001numpy\u306edtype\u306b\u5bfe\u5fdc\u3057\u305f\u578b\u540d\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002\u4f8b\u3048\u3070\u3001float64\u306fnumpy.float64\u306b\u5bfe\u5fdc\u3057\u3066\u3044\u307e\u3059\u3002\u578b\u6307\u5b9a\u306f\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306e\u5f15\u6570\u3068\u3057\u3066\u6307\u5b9a\u3059\u308b\u65b9\u6cd5\u3068\u3001\u95a2\u6570\u306e\u30b7\u30b0\u30cd\u30c1\u30e3\u306b\u76f4\u63a5\u8a18\u8ff0\u3059\u308b\u65b9\u6cd5\u306e2\u901a\u308a\u304c\u3042\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306e\u5f15\u6570\u3068\u3057\u3066\u6307\u5b9a\u3059\u308b\u5834\u5408\u306f\u3001\u6b21\u306e\u3088\u3046\u306b\u8a18\u8ff0\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import jit, float64\n\n@jit(float64(float64, float64))\ndef distance(x, y):\n return ((x - y) ** 2) ** 0.5<\/pre>\n\n\n\n\u95a2\u6570\u306e\u30b7\u30b0\u30cd\u30c1\u30e3\u306b\u76f4\u63a5\u8a18\u8ff0\u3059\u308b\u5834\u5408\u306f\u3001\u6b21\u306e\u3088\u3046\u306b\u8a18\u8ff0\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import jit, float64\n\n@jit\ndef distance_typed(x: float64, y: float64) -> float64:\n return ((x - y) ** 2) ** 0.5<\/pre>\n\n\n\n\u578b\u6307\u5b9a\u3092\u3057\u306a\u3044\u5834\u5408\u3001numba\u306f\u5f15\u6570\u306e\u578b\u3092\u81ea\u52d5\u7684\u306b\u63a8\u6e2c\u3057\u307e\u3059\u304c\u3001\u578b\u6307\u5b9a\u3092\u3059\u308b\u3053\u3068\u3067\u3088\u308a\u6700\u9069\u5316\u3055\u308c\u305f\u30b3\u30fc\u30c9\u3092\u751f\u6210\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u578b\u6307\u5b9a\u306b\u3088\u308b\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u6539\u5584\u306e\u4f8b<\/h3>\n\n\n\n
\u578b\u6307\u5b9a\u306b\u3088\u308b\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u6539\u5584\u3092\u78ba\u8a8d\u3059\u308b\u305f\u3081\u3001\u578b\u6307\u5b9a\u306a\u3057\u306e\u95a2\u6570\u3068\u578b\u6307\u5b9a\u3042\u308a\u306e\u95a2\u6570\u306e\u5b9f\u884c\u6642\u9593\u3092\u6bd4\u8f03\u3057\u3066\u307f\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
from numba import jit, float64\nimport numpy as np\nimport time\n\n# \u578b\u6307\u5b9a\u306a\u3057\u306e\u95a2\u6570\n@jit\ndef sum_squares(n):\n result = 0\n for i in range(n):\n result += i * i\n return result\n\n# \u578b\u6307\u5b9a\u3042\u308a\u306e\u95a2\u6570\n@jit(float64(int64))\ndef sum_squares_typed(n):\n result = 0.0\n for i in range(n):\n result += i * i\n return result\n\n# \u5b9f\u884c\u6642\u9593\u306e\u8a08\u6e2c\nn = 10000000\nstart = time.time()\nsum_squares(n)\nend = time.time()\nprint(f\"\u578b\u6307\u5b9a\u306a\u3057\uff1a{end - start:.3f}\u79d2\")\n\nstart = time.time()\nsum_squares_typed(n)\nend = time.time()\nprint(f\"\u578b\u6307\u5b9a\u3042\u308a\uff1a{end - start:.3f}\u79d2\")<\/pre>\n\n\n\n\u50d5\u306e\u624b\u5143\u74b0\u5883\u3067\u306f\u5b9f\u884c\u7d50\u679c\u306f\u6b21\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
\u578b\u6307\u5b9a\u306a\u3057\uff1a0.034\u79d2\n\u578b\u6307\u5b9a\u3042\u308a\uff1a0.017\u79d2<\/pre>\n\n\n\n\u578b\u6307\u5b9a\u3042\u308a\u306e\u95a2\u6570\u306e\u65b9\u304c\u3001\u578b\u6307\u5b9a\u306a\u3057\u306e\u95a2\u6570\u3088\u308a\u3082\u9ad8\u901f\u3067\u3042\u308b\u3053\u3068\u304c\u308f\u304b\u308a\u307e\u3059\u3002\u3053\u306e\u3088\u3046\u306b\u3001\u578b\u6307\u5b9a\u306b\u3088\u3063\u3066\u5927\u5e45\u306a\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u6539\u5584\u304c\u671f\u5f85\u3067\u304d\u307e\u3059\u3002\u305f\u3060\u3057\u3001\u578b\u6307\u5b9a\u304c\u9069\u5207\u3067\u306a\u3044\u5834\u5408\u3001\u304b\u3048\u3063\u3066\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u304c\u4f4e\u4e0b\u3059\u308b\u3053\u3068\u304c\u3042\u308b\u306e\u3067\u6ce8\u610f\u304c\u5fc5\u8981\u3067\u3059\u3002<\/p>\n\n\n\n
\u6b21\u9805\u3067\u306f\u3001\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3092\u4f7f\u3063\u305f\u52b9\u7387\u7684\u306a\u8a08\u7b97\u65b9\u6cd5\u306b\u3064\u3044\u3066\u89e3\u8aac\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u3010\u30c6\u30af\u30cb\u30c3\u30af3\u3011\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3067\u52b9\u7387\u7684\u306a\u8a08\u7b97\u3092<\/h2>\n\n\n\n
numba\u3092\u4f7f\u3063\u3066Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3059\u308b\u969b\u306b\u3001\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u306f\u975e\u5e38\u306b\u91cd\u8981\u306a\u624b\u6cd5\u3067\u3059\u3002\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001\u914d\u5217\u306e\u8981\u7d20\u3054\u3068\u306e\u6f14\u7b97\u3092\u52b9\u7387\u7684\u306b\u884c\u3046\u3053\u3068\u304c\u3067\u304d\u3001\u51e6\u7406\u901f\u5ea6\u3092\u5927\u5e45\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
ufunc\u3092\u4f7f\u3063\u305f\u30d9\u30af\u30c8\u30eb\u5316\u306e\u65b9\u6cd5<\/h3>\n\n\n\n
numba\u3067\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3092\u884c\u3046\u306b\u306f\u3001numpy\u306eufunc\uff08\u30e6\u30cb\u30d0\u30fc\u30b5\u30eb\u95a2\u6570\uff09\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002ufunc\u306f\u3001\u914d\u5217\u306e\u8981\u7d20\u3054\u3068\u306b\u95a2\u6570\u3092\u9069\u7528\u3057\u3001\u7d50\u679c\u3092\u65b0\u3057\u3044\u914d\u5217\u3068\u3057\u3066\u8fd4\u3059\u6a5f\u80fd\u3092\u6301\u3063\u3066\u3044\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u3067ufunc\u3092\u4f7f\u3046\u306b\u306f\u3001@vectorize\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002@vectorize\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u95a2\u6570\u306b\u9069\u7528\u3059\u308b\u3053\u3068\u3067\u3001\u305d\u306e\u95a2\u6570\u304cufunc\u3068\u3057\u3066\u52d5\u4f5c\u3059\u308b\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import vectorize, float64\n\n@vectorize([float64(float64, float64)])\ndef add_arrays_vectorized(x, y):\n return x + y<\/pre>\n\n\n\n\u3053\u306e\u4f8b\u3067\u306f\u3001@vectorize\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u3063\u3066\u3001\u914d\u5217\u306e\u8981\u7d20\u3054\u3068\u306e\u52a0\u7b97\u3092\u884c\u3046ufunc\u3092\u5b9a\u7fa9\u3057\u3066\u3044\u307e\u3059\u3002@vectorize\u30c7\u30b3\u30ec\u30fc\u30bf\u306e\u5f15\u6570\u306b\u306f\u3001\u5165\u529b\u3068\u51fa\u529b\u306e\u578b\u3092\u6307\u5b9a\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u30d9\u30af\u30c8\u30eb\u5316\u306b\u3088\u308b\u51e6\u7406\u901f\u5ea6\u306e\u5411\u4e0a<\/h3>\n\n\n\n
\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001\u30eb\u30fc\u30d7\u3092\u4f7f\u3063\u305f\u6f14\u7b97\u3088\u308a\u3082\u5927\u5e45\u306b\u51e6\u7406\u901f\u5ea6\u3092\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002\u3053\u308c\u306f\u3001\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u304cCPU\u306eSIMD\uff08Single Instruction Multiple Data\uff09\u547d\u4ee4\u3092\u5229\u7528\u3057\u3066\u3001\u4e26\u5217\u51e6\u7406\u3092\u884c\u3046\u305f\u3081\u3067\u3059\u3002<\/p>\n\n\n\n
\u5b9f\u969b\u306b\u3001\u30d9\u30af\u30c8\u30eb\u5316\u306e\u52b9\u679c\u3092\u78ba\u8a8d\u3059\u308b\u305f\u3081\u306b\u3001\u30d9\u30af\u30c8\u30eb\u5316\u306a\u3057\u306e\u95a2\u6570\u3068\u30d9\u30af\u30c8\u30eb\u5316\u3042\u308a\u306e\u95a2\u6570\u3067\u3001\u5927\u898f\u6a21\u306a\u914d\u5217\u306e\u8981\u7d20\u3054\u3068\u306e\u52a0\u7b97\u3092\u884c\u3063\u3066\u307f\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
from numba import vectorize, float64\nimport numpy as np\nimport time\n\n# \u30d9\u30af\u30c8\u30eb\u5316\u306a\u3057\u306e\u95a2\u6570\ndef add_arrays(a, b):\n result = np.zeros_like(a)\n for i in range(len(a)):\n result[i] = a[i] + b[i]\n return result\n\n# \u30d9\u30af\u30c8\u30eb\u5316\u3042\u308a\u306e\u95a2\u6570\n@vectorize([float64(float64, float64)])\ndef add_arrays_vectorized(x, y):\n return x + y\n\n# \u5b9f\u884c\u6642\u9593\u306e\u8a08\u6e2c\nn = 10000000\na = np.random.rand(n)\nb = np.random.rand(n)\n\nstart = time.time()\nresult = add_arrays(a, b)\nend = time.time()\nprint(f\"\u30d9\u30af\u30c8\u30eb\u5316\u306a\u3057\uff1a{end - start:.3f}\u79d2\")\n\nstart = time.time()\nresult = add_arrays_vectorized(a, b)\nend = time.time()\nprint(f\"\u30d9\u30af\u30c8\u30eb\u5316\u3042\u308a\uff1a{end - start:.3f}\u79d2\")<\/pre>\n\n\n\n\u50d5\u306e\u624b\u5143\u74b0\u5883\u3067\u306f\u5b9f\u884c\u7d50\u679c\u306f\u6b21\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
\u30d9\u30af\u30c8\u30eb\u5316\u306a\u3057\uff1a1.425\u79d2\n\u30d9\u30af\u30c8\u30eb\u5316\u3042\u308a\uff1a0.009\u79d2<\/pre>\n\n\n\n\u30d9\u30af\u30c8\u30eb\u5316\u3042\u308a\u306e\u95a2\u6570\u306e\u65b9\u304c\u3001\u30d9\u30af\u30c8\u30eb\u5316\u306a\u3057\u306e\u95a2\u6570\u3088\u308a\u3082\u5727\u5012\u7684\u306b\u9ad8\u901f\u3067\u3042\u308b\u3053\u3068\u304c\u308f\u304b\u308a\u307e\u3059\u3002\u7279\u306b\u3001\u5927\u898f\u6a21\u306a\u914d\u5217\u306b\u5bfe\u3059\u308b\u6f14\u7b97\u3067\u306f\u3001\u30d9\u30af\u30c8\u30eb\u5316\u306b\u3088\u308b\u901f\u5ea6\u5411\u4e0a\u304c\u9855\u8457\u3067\u3059\u3002<\/p>\n\n\n\n
numba\u3067\u306f\u3001\u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3092\u6d3b\u7528\u3059\u308b\u3053\u3068\u3067\u3001Python\u30b3\u30fc\u30c9\u306e\u51e6\u7406\u901f\u5ea6\u3092\u5927\u5e45\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002\u6b21\u9805\u3067\u306f\u3001\u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3063\u3066\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316\u3092\u5b9f\u73fe\u3059\u308b\u65b9\u6cd5\u306b\u3064\u3044\u3066\u89e3\u8aac\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u3010\u30c6\u30af\u30cb\u30c3\u30af4\u3011\u4e26\u5217\u51e6\u7406\u3067\u30de\u30eb\u30c1\u30b3\u30a2CPU\u306e\u6027\u80fd\u3092\u5f15\u304d\u51fa\u3059<\/h2>\n\n\n\n
numba\u3092\u4f7f\u3063\u3066Python\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3059\u308b\u969b\u306b\u3001\u4e26\u5217\u51e6\u7406\u306f\u975e\u5e38\u306b\u5f37\u529b\u306a\u624b\u6cd5\u3067\u3059\u3002\u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001\u30de\u30eb\u30c1\u30b3\u30a2CPU\u3084\u30de\u30eb\u30c1\u30d7\u30ed\u30bb\u30c3\u30b5\u30b7\u30b9\u30c6\u30e0\u306e\u6027\u80fd\u3092\u6700\u5927\u9650\u306b\u5f15\u304d\u51fa\u3057\u3001\u30d7\u30ed\u30b0\u30e9\u30e0\u306e\u5b9f\u884c\u901f\u5ea6\u3092\u5927\u5e45\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u4e26\u5217\u51e6\u7406\u306e\u6982\u8981\u3068\u30e1\u30ea\u30c3\u30c8<\/h3>\n\n\n\n
\u4e26\u5217\u51e6\u7406\u3068\u306f\u3001\u8907\u6570\u306e\u30b3\u30a2\u3084\u30d7\u30ed\u30bb\u30c3\u30b5\u3092\u540c\u6642\u306b\u5229\u7528\u3057\u3066\u3001\u30bf\u30b9\u30af\u3092\u4e26\u884c\u3057\u3066\u5b9f\u884c\u3059\u308b\u624b\u6cd5\u3067\u3059\u3002numba\u3067\u306f\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306eparallel=True\u30aa\u30d7\u30b7\u30e7\u30f3\u3092\u4f7f\u3063\u3066\u3001\u4e26\u5217\u51e6\u7406\u3092\u6709\u52b9\u306b\u3059\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30e1\u30ea\u30c3\u30c8\u304c\u3042\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u306e\u4e26\u5217\u51e6\u7406\u3067\u306f\u3001prange()\u95a2\u6570\u3092\u4f7f\u3063\u3066\u30eb\u30fc\u30d7\u3092\u4e26\u5217\u5316\u3057\u307e\u3059\u3002prange()\u95a2\u6570\u306f\u3001range()\u95a2\u6570\u3068\u540c\u69d8\u306b\u4f7f\u7528\u3067\u304d\u307e\u3059\u304c\u3001\u30eb\u30fc\u30d7\u306e\u53cd\u5fa9\u51e6\u7406\u3092\u4e26\u5217\u306b\u5b9f\u884c\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u4ee5\u4e0b\u306f\u3001prange()\u95a2\u6570\u3092\u4f7f\u3063\u305f\u4e26\u5217\u51e6\u7406\u306e\u4f8b\u3067\u3059\u3002<\/p>\n\n\n\n
from numba import jit, prange\n\n@jit(nopython=True, parallel=True)\ndef sum_squares_parallel(n):\n result = 0\n for i in prange(n):\n result += i * i\n return result<\/pre>\n\n\n\n\u3053\u306e\u4f8b\u3067\u306f\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u306eparallel=True\u30aa\u30d7\u30b7\u30e7\u30f3\u3092\u6307\u5b9a\u3057\u3001prange()\u95a2\u6570\u3092\u4f7f\u3063\u3066\u30eb\u30fc\u30d7\u3092\u4e26\u5217\u5316\u3057\u3066\u3044\u307e\u3059\u3002<\/p>\n\n\n\n
\u5b9f\u969b\u306b\u3001\u4e26\u5217\u51e6\u7406\u306e\u52b9\u679c\u3092\u78ba\u8a8d\u3059\u308b\u305f\u3081\u306b\u3001\u901a\u5e38\u306e\u30eb\u30fc\u30d7\u3092\u4f7f\u3063\u305f\u95a2\u6570\u3068\u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3063\u305f\u95a2\u6570\u3067\u30011\u304b\u3089(n-1)\u307e\u3067\u306e\u6574\u6570\u306e\u4e8c\u4e57\u548c\u3092\u8a08\u7b97\u3057\u3066\u307f\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
from numba import jit, prange\nimport numpy as np\nimport time\n\n# \u901a\u5e38\u306e\u30eb\u30fc\u30d7\u3092\u4f7f\u3063\u305f\u95a2\u6570\ndef sum_squares(n):\n result = 0\n for i in range(n):\n result += i * i\n return result\n\n# \u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3063\u305f\u95a2\u6570\n@jit(nopython=True, parallel=True)\ndef sum_squares_parallel(n):\n result = 0\n for i in prange(n):\n result += i * i\n return result\n\n# \u5b9f\u884c\u6642\u9593\u306e\u8a08\u6e2c\nn = 100000000\n\nstart = time.time()\nresult = sum_squares(n)\nend = time.time()\nprint(f\"\u901a\u5e38\u306e\u30eb\u30fc\u30d7\uff1a{end - start:.3f}\u79d2\")\n\nstart = time.time()\nresult = sum_squares_parallel(n)\nend = time.time()\nprint(f\"\u4e26\u5217\u51e6\u7406\uff1a{end - start:.3f}\u79d2\")<\/pre>\n\n\n\n\u50d5\u306e\u624b\u5143\u74b0\u5883\u3067\u306f\u5b9f\u884c\u7d50\u679c\u306f\u6b21\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
\u901a\u5e38\u306e\u30eb\u30fc\u30d7\uff1a5.203\u79d2\n\u4e26\u5217\u51e6\u7406\uff1a0.659\u79d2<\/pre>\n\n\n\n\u4e26\u5217\u51e6\u7406\u3092\u4f7f\u3063\u305f\u95a2\u6570\u306e\u65b9\u304c\u3001\u901a\u5e38\u306e\u30eb\u30fc\u30d7\u3092\u4f7f\u3063\u305f\u95a2\u6570\u3088\u308a\u3082\u5927\u5e45\u306b\u9ad8\u901f\u3067\u3042\u308b\u3053\u3068\u304c\u308f\u304b\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
\u3053\u306e\u3088\u3046\u306b\u3001numba\u306e\u4e26\u5217\u51e6\u7406\u6a5f\u80fd\u3092\u6d3b\u7528\u3059\u308b\u3053\u3068\u3067\u3001Python\u30b3\u30fc\u30c9\u306e\u5b9f\u884c\u901f\u5ea6\u3092\u98db\u8e8d\u7684\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002\u6b21\u9805\u3067\u306f\u3001GPU\u3092\u4f7f\u3063\u305f\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316\u306b\u3064\u3044\u3066\u89e3\u8aac\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\u3010\u30c6\u30af\u30cb\u30c3\u30af5\u3011GPU\u3092\u6d3b\u7528\u3057\u3066\u5927\u898f\u6a21\u30c7\u30fc\u30bf\u3092\u9ad8\u901f\u51e6\u7406<\/h2>\n\n\n\n
\u3053\u308c\u307e\u3067\u7d39\u4ecb\u3057\u3066\u304d\u305f\u30c6\u30af\u30cb\u30c3\u30af\u306b\u52a0\u3048\u3066\u3001numba\u3067\u306fGPU\u3092\u6d3b\u7528\u3059\u308b\u3053\u3068\u3067\u3001\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316\u3092\u5b9f\u73fe\u3067\u304d\u307e\u3059\u3002GPU\u306f\u5927\u898f\u6a21\u306a\u4e26\u5217\u51e6\u7406\u306b\u9069\u3057\u3066\u3044\u308b\u305f\u3081\u3001\u7279\u306b\u884c\u5217\u7a4d\u306a\u3069\u306e\u6570\u5024\u8a08\u7b97\u3084\u3001\u6a5f\u68b0\u5b66\u7fd2\u306e\u30e2\u30c7\u30eb\u8a13\u7df4\u3001\u30b7\u30df\u30e5\u30ec\u30fc\u30b7\u30e7\u30f3\u3084\u79d1\u5b66\u8a08\u7b97\u306a\u3069\u3001\u5927\u91cf\u306e\u4e26\u5217\u8a08\u7b97\u304c\u5fc5\u8981\u306a\u5206\u91ce\u3067\u5a01\u529b\u3092\u767a\u63ee\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u306eCUDA\u30b5\u30dd\u30fc\u30c8 \u2013 GPU\u3078\u306e\u30aa\u30d5\u30ed\u30fc\u30c9<\/h3>\n\n\n\n
numba\u306f\u3001NVIDIA\u793e\u306eCUDA\u30d7\u30e9\u30c3\u30c8\u30d5\u30a9\u30fc\u30e0\u3092\u30b5\u30dd\u30fc\u30c8\u3057\u3066\u304a\u308a\u3001Python\u30b3\u30fc\u30c9\u3092GPU\u3067\u5b9f\u884c\u3059\u308b\u305f\u3081\u306e\u30b3\u30f3\u30d1\u30a4\u30e9\u6a5f\u80fd\u3092\u63d0\u4f9b\u3057\u3066\u3044\u307e\u3059\u3002\u3053\u306e\u6a5f\u80fd\u3092\u4f7f\u3046\u3053\u3068\u3067\u3001Python\u30b3\u30fc\u30c9\u3092GPU\u306b\u30aa\u30d5\u30ed\u30fc\u30c9\u3057\u3066\u5b9f\u884c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
GPU\u306b\u30aa\u30d5\u30ed\u30fc\u30c9\u3059\u308b\u306b\u306f\u3001@cuda.jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u7528\u3057\u307e\u3059\u3002@cuda.jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4ed8\u3051\u305f\u95a2\u6570\u306f\u3001GPU\u4e0a\u3067\u5b9f\u884c\u3055\u308c\u308b\u30ab\u30fc\u30cd\u30eb\u95a2\u6570\u3068\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
from numba import cuda\n\n@cuda.jit\ndef matmul_gpu(A, B, C):\n i, j = cuda.grid(2)\n if i < C.shape[0] and j < C.shape[1]:\n tmp = 0.\n for k in range(A.shape[1]):\n tmp += A[i, k] * B[k, j]\n C[i, j] = tmp<\/pre>\n\n\n\n\u3053\u306e\u4f8b\u3067\u306f\u3001@cuda.jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u3063\u3066GPU\u4e0a\u3067\u5b9f\u884c\u3059\u308b\u30ab\u30fc\u30cd\u30eb\u95a2\u6570\u3092\u5b9a\u7fa9\u3057\u3066\u3044\u307e\u3059\u3002\u30ab\u30fc\u30cd\u30eb\u95a2\u6570\u5185\u3067\u306f\u3001cuda.grid\u3092\u4f7f\u3063\u3066\u30b9\u30ec\u30c3\u30c9\u3068\u30d6\u30ed\u30c3\u30af\u306e\u30a4\u30f3\u30c7\u30c3\u30af\u30b9\u3092\u53d6\u5f97\u3057\u3001\u884c\u5217\u7a4d\u306e\u8a08\u7b97\u3092\u4e26\u5217\u306b\u5b9f\u884c\u3057\u3066\u3044\u307e\u3059\u3002<\/p>\n\n\n\n
GPU\u3092\u4f7f\u3063\u305f\u9ad8\u901f\u5316\u306e\u5b9f\u4f8b<\/h3>\n\n\n\n
\u5b9f\u969b\u306b\u3001GPU\u3092\u4f7f\u3063\u305f\u9ad8\u901f\u5316\u306e\u52b9\u679c\u3092\u78ba\u8a8d\u3059\u308b\u305f\u3081\u306b\u3001\u884c\u5217\u7a4d\u3092CPU\u3068GPU\u3067\u5b9f\u884c\u3057\u3001\u305d\u306e\u5b9f\u884c\u6642\u9593\u3092\u6bd4\u8f03\u3057\u3066\u307f\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
from numba import cuda\nimport numpy as np\nimport time\n\n# GPU\u3067\u5b9f\u884c\u3059\u308b\u95a2\u6570\n@cuda.jit\ndef matmul_gpu(A, B, C):\n i, j = cuda.grid(2)\n if i < C.shape[0] and j < C.shape[1]:\n tmp = 0.\n for k in range(A.shape[1]):\n tmp += A[i, k] * B[k, j]\n C[i, j] = tmp\n\n# CPU\u3067\u5b9f\u884c\u3059\u308b\u95a2\u6570\ndef matmul_cpu(A, B, C):\n for i in range(A.shape[0]):\n for j in range(B.shape[1]):\n tmp = 0.\n for k in range(A.shape[1]):\n tmp += A[i, k] * B[k, j]\n C[i, j] = tmp\n\n# \u884c\u5217\u306e\u30b5\u30a4\u30ba\nn = 1024\n\n# \u884c\u5217\u306e\u521d\u671f\u5316\nA = np.random.rand(n, n)\nB = np.random.rand(n, n)\nC_cpu = np.zeros((n, n))\nC_gpu = np.zeros((n, n))\n\n# CPU\u3067\u306e\u5b9f\u884c\u6642\u9593\u3092\u8a08\u6e2c\nstart = time.time()\nmatmul_cpu(A, B, C_cpu)\nend = time.time()\nprint(f\"CPU\u3067\u306e\u5b9f\u884c\u6642\u9593\uff1a{end - start:.3f}\u79d2\")\n\n# GPU\u3067\u306e\u5b9f\u884c\u6642\u9593\u3092\u8a08\u6e2c\nstart = time.time()\ndA = cuda.to_device(A)\ndB = cuda.to_device(B)\ndC = cuda.to_device(C_gpu)\nthreads_per_block = (16, 16)\nblocks_per_grid_x = (C_gpu.shape[0] + threads_per_block[0] - 1) \/\/ threads_per_block[0]\nblocks_per_grid_y = (C_gpu.shape[1] + threads_per_block[1] - 1) \/\/ threads_per_block[1]\nblocks_per_grid = (blocks_per_grid_x, blocks_per_grid_y)\nmatmul_gpu[blocks_per_grid, threads_per_block](dA, dB, dC)\ndC.copy_to_host(C_gpu)\nend = time.time()\nprint(f\"GPU\u3067\u306e\u5b9f\u884c\u6642\u9593\uff1a{end - start:.3f}\u79d2\")<\/pre>\n\n\n\n\u50d5\u306e\u624b\u5143\u74b0\u5883\u3067\u306f\u5b9f\u884c\u7d50\u679c\u306f\u6b21\u306e\u3088\u3046\u306b\u306a\u308a\u307e\u3059\u3002<\/p>\n\n\n\n
CPU\u3067\u306e\u5b9f\u884c\u6642\u9593\uff1a12.529\u79d2\nGPU\u3067\u306e\u5b9f\u884c\u6642\u9593\uff1a0.022\u79d2<\/pre>\n\n\n\nGPU\u3092\u4f7f\u3063\u305f\u884c\u5217\u7a4d\u306e\u65b9\u304c\u3001CPU\u3092\u4f7f\u3063\u305f\u5834\u5408\u3088\u308a\u3082\u5727\u5012\u7684\u306b\u9ad8\u901f\u3067\u3042\u308b\u3053\u3068\u304c\u308f\u304b\u308a\u307e\u3059\u3002\u3053\u306e\u3088\u3046\u306b\u3001numba\u306eGPU\u30b5\u30dd\u30fc\u30c8\u3092\u6d3b\u7528\u3059\u308b\u3053\u3068\u3067\u3001\u5927\u898f\u6a21\u306a\u6570\u5024\u8a08\u7b97\u3084\u30c7\u30fc\u30bf\u51e6\u7406\u3092\u9ad8\u901f\u306b\u884c\u3046\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
\u305f\u3060\u3057\u3001GPU\u3092\u4f7f\u3046\u305f\u3081\u306b\u306f\u3001NVIDIA\u793e\u88fd\u306eGPU\u3068CUDA\u30c4\u30fc\u30eb\u30ad\u30c3\u30c8\u304c\u5fc5\u8981\u3067\u3042\u308a\u3001\u74b0\u5883\u69cb\u7bc9\u306b\u4e00\u5b9a\u306e\u624b\u9593\u304c\u304b\u304b\u308a\u307e\u3059\u3002\u307e\u305f\u3001GPU\u306e\u30e1\u30e2\u30ea\u5bb9\u91cf\u3084\u30c7\u30fc\u30bf\u8ee2\u9001\u306e\u30aa\u30fc\u30d0\u30fc\u30d8\u30c3\u30c9\u306a\u3069\u306b\u3082\u6ce8\u610f\u304c\u5fc5\u8981\u3067\u3059\u3002GPU\u3092\u6d3b\u7528\u3059\u308b\u969b\u306f\u3001\u3053\u308c\u3089\u306e\u70b9\u3092\u8003\u616e\u3057\u3066\u3001\u9069\u5207\u306b\u8a2d\u8a08\u30fb\u5b9f\u88c5\u3059\u308b\u3053\u3068\u304c\u91cd\u8981\u3067\u3059\u3002<\/p>\n\n\n\n
\u307e\u3068\u3081 \u2013 numba\u3067\u7206\u901fPython\u30b3\u30fc\u30c7\u30a3\u30f3\u30b0\u3092\u4f53\u9a13\u3057\u3088\u3046\uff01<\/h2>\n\n\n\n
numba\u306f\u3001Python\u30b3\u30fc\u30c9\u306e\u9ad8\u901f\u5316\u306b\u7279\u5316\u3057\u305f\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3042\u308a\u3001JIT\u30b3\u30f3\u30d1\u30a4\u30e9\u3001\u30d9\u30af\u30c8\u30eb\u5316\u3001\u4e26\u5217\u51e6\u7406\u3001GPU\u306e\u6d3b\u7528\u306a\u3069\u3001\u69d8\u3005\u306a\u9ad8\u901f\u5316\u624b\u6cd5\u3092\u63d0\u4f9b\u3057\u3066\u3044\u307e\u3059\u3002\u3053\u306e\u8a18\u4e8b\u3067\u306f\u3001numba\u306e\u57fa\u672c\u7684\u306a\u4f7f\u3044\u65b9\u304b\u3089\u3001\u3088\u308a\u9ad8\u5ea6\u306a\u30c6\u30af\u30cb\u30c3\u30af\u307e\u3067\u3001\u5e45\u5e83\u304f\u89e3\u8aac\u3057\u3066\u304d\u307e\u3057\u305f\u3002<\/p>\n\n\n\n
5\u3064\u306e\u30c6\u30af\u30cb\u30c3\u30af\u306e\u632f\u308a\u8fd4\u308a<\/h3>\n\n\n\n
\n
- @jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u4f7f\u3063\u305f\u95a2\u6570\u306eJIT\u30b3\u30f3\u30d1\u30a4\u30eb<\/li>\n\n\n\n
- \u578b\u6307\u5b9a\u306b\u3088\u308b\u3055\u3089\u306a\u308b\u9ad8\u901f\u5316<\/li>\n\n\n\n
- \u30d9\u30af\u30c8\u30eb\u5316\u6f14\u7b97\u3067\u52b9\u7387\u7684\u306a\u8a08\u7b97<\/li>\n\n\n\n
- \u4e26\u5217\u51e6\u7406\u3067\u30de\u30eb\u30c1\u30b3\u30a2CPU\u306e\u6027\u80fd\u3092\u5f15\u304d\u51fa\u3059<\/li>\n\n\n\n
- GPU\u3092\u6d3b\u7528\u3057\u3066\u5927\u898f\u6a21\u30c7\u30fc\u30bf\u3092\u9ad8\u901f\u51e6\u7406<\/li>\n<\/ol>\n\n\n\n
\u3053\u308c\u3089\u306e\u30c6\u30af\u30cb\u30c3\u30af\u3092\u9069\u5207\u306b\u7d44\u307f\u5408\u308f\u305b\u308b\u3053\u3068\u3067\u3001Python\u30b3\u30fc\u30c9\u306e\u51e6\u7406\u901f\u5ea6\u3092\u5927\u5e45\u306b\u5411\u4e0a\u3055\u305b\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n\n\n\n
numba\u3092\u6d3b\u7528\u3057\u305f\u9ad8\u901f\u30b3\u30fc\u30c7\u30a3\u30f3\u30b0\u306e\u3059\u3059\u3081<\/h3>\n\n\n\n
numba\u3092\u52b9\u679c\u7684\u306b\u6d3b\u7528\u3059\u308b\u306b\u306f\u3001\u4ee5\u4e0b\u306e\u70b9\u306b\u7559\u610f\u3057\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
\n
- \u9ad8\u901f\u5316\u306e\u5bfe\u8c61\u3068\u306a\u308b\u95a2\u6570\u3092\u7279\u5b9a\u3057\u3001@jit\u30c7\u30b3\u30ec\u30fc\u30bf\u3092\u9069\u7528\u3059\u308b\u3002<\/li>\n\n\n\n
- \u53ef\u80fd\u306a\u9650\u308a\u578b\u6307\u5b9a\u3092\u884c\u3044\u3001\u3088\u308a\u6700\u9069\u5316\u3055\u308c\u305f\u30b3\u30fc\u30c9\u3092\u751f\u6210\u3059\u308b\u3002<\/li>\n\n\n\n
- \u30d9\u30af\u30c8\u30eb\u5316\u3084\u4e26\u5217\u51e6\u7406\u3092\u7a4d\u6975\u7684\u306b\u6d3b\u7528\u3057\u3001\u51e6\u7406\u901f\u5ea6\u3092\u5411\u4e0a\u3055\u305b\u308b\u3002<\/li>\n\n\n\n
- \u5fc5\u8981\u306b\u5fdc\u3058\u3066\u3001GPU\u306e\u6d3b\u7528\u3092\u691c\u8a0e\u3059\u308b\u3002<\/li>\n<\/ul>\n\n\n\n
\u307e\u305f\u3001numba\u3092\u6d3b\u7528\u3059\u308b\u305f\u3081\u306e\u5177\u4f53\u7684\u306a\u6b21\u306e\u30b9\u30c6\u30c3\u30d7\u3068\u3057\u3066\u3001\u4ee5\u4e0b\u3092\u63d0\u6848\u3057\u307e\u3059\u3002<\/p>\n\n\n\n
\n
- \u30b7\u30f3\u30d7\u30eb\u306a\u95a2\u6570\u3067JIT\u30b3\u30f3\u30d1\u30a4\u30eb\u3092\u8a66\u3057\u3001\u9ad8\u901f\u5316\u306e\u52b9\u679c\u3092\u78ba\u8a8d\u3059\u308b\u3002<\/li>\n\n\n\n
- \u578b\u6307\u5b9a\u3084\u30d9\u30af\u30c8\u30eb\u5316\u3001\u4e26\u5217\u51e6\u7406\u3092\u6bb5\u968e\u7684\u306b\u5c0e\u5165\u3057\u3001\u6700\u9069\u5316\u3092\u9032\u3081\u308b\u3002<\/li>\n\n\n\n
- \u30b3\u30fc\u30c9\u306e\u53ef\u8aad\u6027\u3068\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u306e\u30d0\u30e9\u30f3\u30b9\u3092\u8003\u616e\u3057\u306a\u304c\u3089\u3001\u6539\u5584\u3092\u7d9a\u3051\u308b\u3002<\/li>\n\n\n\n
- \u5927\u898f\u6a21\u306a\u6570\u5024\u8a08\u7b97\u3084\u30c7\u30fc\u30bf\u51e6\u7406\u304c\u5fc5\u8981\u306a\u5834\u5408\u306f\u3001GPU\u306e\u6d3b\u7528\u3092\u691c\u8a0e\u3059\u308b\u3002<\/li>\n<\/ol>\n\n\n\n
\u305f\u3060\u3057\u3001JIT\u30b3\u30f3\u30d1\u30a4\u30eb\u306e\u30aa\u30fc\u30d0\u30fc\u30d8\u30c3\u30c9\u3084\u3001\u578b\u6307\u5b9a\u306e\u9069\u5207\u3055\u3001GPU\u74b0\u5883\u306e\u69cb\u7bc9\u306a\u3069\u3001\u6ce8\u610f\u3059\u3079\u304d\u70b9\u3082\u3042\u308a\u307e\u3059\u3002\u3053\u308c\u3089\u306e\u70b9\u306b\u7559\u610f\u3057\u3064\u3064\u3001numba\u3092\u4e0a\u624b\u306b\u6d3b\u7528\u3059\u308b\u3053\u3068\u3067\u3001Python\u30b3\u30fc\u30c9\u306e\u9ad8\u901f\u5316\u3092\u5b9f\u73fe\u3057\u307e\u3057\u3087\u3046\u3002<\/p>\n\n\n\n
numba\u306f\u3001Python\u306e\u9ad8\u901f\u5316\u306b\u60a9\u3080\u3059\u3079\u3066\u306e\u958b\u767a\u8005\u306b\u3068\u3063\u3066\u3001\u5f37\u529b\u306a\u5473\u65b9\u3068\u306a\u308b\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3059\u3002\u3053\u306e\u8a18\u4e8b\u3067\u7d39\u4ecb\u3057\u305f\u30c6\u30af\u30cb\u30c3\u30af\u3092\u5b9f\u8df5\u3057\u3001numba\u3067\u7206\u901fPython\u30b3\u30fc\u30c7\u30a3\u30f3\u30b0\u3092\u4f53\u9a13\u3057\u3066\u307f\u3066\u304f\u3060\u3055\u3044\uff01<\/p>\n","protected":false},"excerpt":{"rendered":"
Python\u306e\u30d1\u30d5\u30a9\u30fc\u30de\u30f3\u30b9\u306b\u60a9\u3080\u30c7\u30d9\u30ed\u30c3\u30d1\u30fc\u306b\u6717\u5831\uff01numba\u3092\u4f7f\u3048\u3070\u3001\u7c21\u5358\u306b\u30b3\u30fc\u30c9\u3092\u9ad8\u901f\u5316\u3067\u304d\u307e\u3059\u3002\u3053\u306e\u8a18\u4e8b\u3067\u306f\u3001numba\u306e\u57fa\u672c\u304b\u3089\u5b9f\u8df5\u7684\u306a\u30c6\u30af\u30cb\u30c3\u30af\u307e\u3067\u3001\u8a73\u3057\u304f\u89e3\u8aac\u3057\u307e\u3059\u3002 Warning: Undefined … <\/p>\n","protected":false},"author":1,"featured_media":607,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[4],"tags":[],"class_list":{"0":"post-602","1":"post","2":"type-post","3":"status-publish","4":"format-standard","5":"has-post-thumbnail","7":"category-python"},"_links":{"self":[{"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/posts\/602","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/chocottopro.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=602"}],"version-history":[{"count":5,"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/posts\/602\/revisions"}],"predecessor-version":[{"id":608,"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/posts\/602\/revisions\/608"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/chocottopro.com\/index.php?rest_route=\/wp\/v2\/media\/607"}],"wp:attachment":[{"href":"https:\/\/chocottopro.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=602"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/chocottopro.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=602"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/chocottopro.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=602"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}