{"id":681,"date":"2011-12-15T01:50:22","date_gmt":"2011-12-14T16:50:22","guid":{"rendered":"http:\/\/daily.belltail.jp\/?p=681"},"modified":"2011-12-15T12:29:50","modified_gmt":"2011-12-15T03:29:50","slug":"boostublasbounded_vector%e3%82%92%e4%bd%bf%e3%81%86%e3%81%a8%e3%81%8d%e3%81%af%e5%b0%91%e3%81%97%e6%b0%97%e3%82%92%e3%81%a4%e3%81%91%e3%81%a6%e3%81%ad%e3%81%a8%e3%81%84%e3%81%86%e8%a9%b1","status":"publish","type":"post","link":"https:\/\/daily.belltail.jp\/?p=681","title":{"rendered":"boost::ublas::bounded_vector\u3092\u4f7f\u3046\u3068\u304d\u306f\u5c11\u3057\u6c17\u3092\u3064\u3051\u3066\u306d\u3068\u3044\u3046\u8a71"},"content":{"rendered":"

\u3053\u306e\u8a18\u4e8b\u306f\uff0c[http:\/\/partake.in\/events\/597a0fc3-0e3a-47a3-8fc3-4f32ad846a3d:title=Boost Advent Calendar 2011]\u306e14\u65e5\u76ee\u306e\u8a18\u4e8b\u3067\u3059\uff0e
\n(2\u6642\u9593\u8fd1\u304f\u3082\u9045\u523b\u3057\u3066\u3059\u307f\u307e\u305b\u3093\uff1e\uff1c)<\/span><\/p>\n

\u30b9\u30ad\u30eb\u306f\u307e\u3060\u307e\u3060\u3067\u3059\u304c\uff0c\u305d\u308c\u3067\u3082\u3082\u3046boost\u304c\u7121\u3044\u3068\u751f\u304d\u3066\u3044\u3051\u306a\u3044\u4f53\u306b\u306a\u3063\u3066\u3057\u307e\u3044\u307e\u3057\u305f\uff0e
\n\u305d\u3046\u3044\u3048\u3070[http:\/\/longgate.co.jp\/books\/grimoire-vol2.html:title=\u9b54\u5c0e\u66f8Vol2]\u306e[http:\/\/d.hatena.ne.jp\/kikairoya\/:title=id:kikairoya]\u5148\u751f\u306e\u7ae0(\u4e00\u756a\u6700\u5f8c\u3067\u3059\u306d)\u3042\u305f\u308a\u306f\uff0c\u30db\u30f3\u30c8\u306b\u300cboost\u3092\u4f7f\u3044\u5012\u3059\uff01\u300d\u304b\u3093\u3058\u304c\u3057\u3066\u6700\u9ad8\u306b\u8aad\u307f\u5fdc\u3048\u304c\u3042\u308a\u307e\u3057\u305f\uff0e\u305c\u3072\u304a\u8aad\u307f\u304f\u3060\u3055\u3044\uff0e<\/p>\n

\u3055\u3066\uff0c\u672c\u6765\u306f\u3053\u306e\u30a2\u30c9\u30d9\u3067\u306f\u7121\u96e3\u306b\u300c\u3042\u307e\u308a\u65e5\u672c\u8a9e\u60c5\u5831\u306e\u51fa\u56de\u3063\u3066\u306a\u3044\u30e9\u30a4\u30d6\u30e9\u30ea\u3092(\u307b\u3068\u3093\u3069\u7ffb\u8a33\u3067)\u7d39\u4ecb\u3059\u308b\u300d\u3064\u3082\u308a\u3060\u3063\u305f\u3093\u3067\u3059\u304c\uff0c\u610f\u5916\u3068\u30cd\u30bf\u63a2\u3057\u306b\u3064\u307e\u3065\u3044\u305f\u306e\u3067\uff0c\u65b9\u91dd\u8ee2\u63db\u3057\u3066\u50d5\u304c\u5f15\u3063\u304b\u304b\u3063\u3066\u3057\u307e\u3063\u305f\u843d\u3068\u3057\u7a74<\/strong>\u3092\u7d39\u4ecb\u3059\u308b\u3053\u3068\u306b\uff0e
\n\u7b4b\u9055\u3044\u306a\u3053\u3068\u3092\u8a00\u3063\u3066\u305f\u3089\uff0c\u305f\u3076\u3093\u8a00\u3063\u3066\u307e\u3059\u304c\uff0c\u30c4\u30c3\u30b3\u30df\u3092\u304a\u9858\u3044\u3057\u307e\u3059\uff0e
\n\u30db\u30f3\u30c8\u306b\u4e0d\u5f53\u306adis\u304c\u3042\u3063\u305f\u3089\u5168\u529b\u3067\u8b1d\u7f6a\u8a18\u4e8b\u3092<\/strong>\u66f8\u304b\u305b\u3066\u9802\u304d\u307e\u3059\u2026\uff0e<\/p>\n

(2011.12.15 12:30) \u7de8\u96c6\u4e0a\u306e\u30df\u30b9\u306b\u3064\u3044\u3066\u6307\u6458\u3092\u3044\u305f\u3060\u3044\u305f\u306e\u3067\u8a02\u6b63<\/span><\/p>\n

[:contents]<\/p>\n

<\/p>\n

\u203b\u6b63\u78ba\u306b\u306fboost::numeric::ublas\u3067\u3059\u304c\uff0c\u4fbf\u5b9c\u7684\u306bboost::ublas\u3068\u8868\u8a18\u3057\u307e\u3059\uff0e<\/p>\n

*boost::ublas\u3068\u306f<\/p>\n

\u7279\u306b\u7d30\u304b\u3044\u8aac\u660e\u306e\u5fc5\u8981\u306f\u7121\u3044\u3068\u306f\u601d\u3044\u307e\u3059\u304c\uff0e
\n1.29.0\u3067\u6b63\u5f0f\u306bboost\u5165\u308a\u3057\u305f\uff0c\u7dda\u5f62\u4ee3\u6570\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3059\uff0e
\nBLAS\u3068\u3044\u3046\u30e1\u30b8\u30e3\u30fc\u306a\u7dda\u5f62\u4ee3\u6570API\u304c\u3042\u3063\u3066(\u30e9\u30a4\u30d6\u30e9\u30ea\u3068\u3044\u3046\u3088\u308a\u306f\u898f\u683c\u3068\u8a00\u3063\u305f\u307b\u3046\u304c\u6b63\u78ba\u304b\u3082\u3057\u308c\u307e\u305b\u3093)\uff0c\u3053\u306e\u30b5\u30d6\u30bb\u30c3\u30c8\u3068\u3057\u3066uBLAS\u304c\u3042\u308a\u307e\u3059\uff0e\u8aad\u307f\u306f\u300c\u3086\u30fc\u3076\u3089\u3059\u300d\uff0e<\/p>\n

\r\n#include <iostream>\r\n#include <algorithm>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/matrix.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\tvector<double> vec(3);\r\n\tvec[0] = 1;\r\n\tvec[1] = 2;\r\n\tvec[2] = 3;\r\n\r\n\tmatrix<double> mat(3, 2);\r\n\tmat(0, 0) = 1; mat(0, 1) = 2;\r\n\tmat(1, 0) = 3; mat(1, 1) = 4;\r\n\tmat(2, 0) = 5; mat(2, 1) = 6;\r\n\r\n\tstd::cout << vec << std::endl;\r\n\tstd::cout << mat << std::endl;\r\n\tstd::cout << (3 * vec) << std::endl;\r\n\tstd::cout << prod(vec, mat) << std::endl;\r\n}\r\n<\/pre>\n
\u306b\u5bfe\u3057\u3066<\/p>\n
\r\n [][3](1,2,3)[]\r\n [][3,2]((1,2),(3,4),(5,6))[]\r\n [][3](3,6,9)[]\r\n [][2](22,28)[]\r\n<\/pre>\n
\u306a\u304b\u3093\u3058\uff0e<\/p>\n
uBLAS\u3092\u7d39\u4ecb\u3059\u308b\u3053\u3068\u306f\u672c\u7b4b\u3058\u3083\u306a\u3044\u306e\u3067\uff0c\u3053\u306e\u3042\u305f\u308a\u306f\u3055\u3089\u3063\u3068\uff0e<\/p>\n
*ublas::bounded_vector\u306b\u3064\u3044\u3066<\/p>\n
\u3055\u3066\u3055\u3066\uff0e
\n\u307f\u306a\u3055\u3093\u306f\u81ea\u5206\u306e\u30d7\u30ed\u30b0\u30e9\u30e0\u3067\u3061\u3087\u3063\u3068\u3057\u305f\u7dda\u5f62\u4ee3\u6570\u306e\u6570\u5024\u8a08\u7b97\u3092\u3059\u308b\u3068\u304d\uff0c\u3060\u3044\u305f\u3044\u4f55\u6b21\u5143\u306e\u884c\u5217\u306a\u308a\u30d9\u30af\u30c8\u30eb\u3092\u5229\u7528\u3055\u308c\u307e\u3059\u304b\uff1f
\n\u30c7\u30fc\u30bf\u51e6\u7406\u3084\u6a5f\u68b0\u5b66\u7fd2\u306e\u5206\u91ce\u3067\u306f\uff0c\u6570\u4e07\u304b\u3089\u6570\u767e\u4e07\u6b21\u5143\u3068\u3044\u3046\u30ec\u30d9\u30eb\u304c\u308f\u308a\u3068\u30b6\u30e9\u3067\u3059\uff0e\u5929\u6c17\u4e88\u5831\u3084X\u7ddaCT\u306a\u3093\u304b\u3082\u3042\u308c\u306f\u4e2d\u3067\u8d85\u9ad8\u6b21\u5143\u306a\u4e00\u6b21\u9023\u7acb\u65b9\u7a0b\u5f0f\u3092\u89e3\u3044\u3066\u308b\u3060\u3051\u3067\u3059\uff0e
\n\u4e16\u306b\u51fa\u56de\u3063\u3066\u308b\u7dda\u5f62\u4ee3\u6570\u30e9\u30a4\u30d6\u30e9\u30ea\u306f\u57fa\u672c\u7684\u306b\u306f\u3053\u306e\u3088\u3046\u306a\u30bf\u30b9\u30af\u306b\u6700\u9069\u5316\u3055\u308c\u3066\u3044\u307e\u3059(\u7279\u306bC\u7cfb\u306e\u3082\u306e\u306f)\uff0e<\/p>\n
\u3067\u3082\u305f\u3076\u3093\u666e\u901a\u306e(\u666e\u901a\u306e\u3063\u3066\u4f55\u3088\u3063\u3066\u8a71\u3067\u3059\u304c)\u30d7\u30ed\u30b0\u30e9\u30e0\u66f8\u304f\u4e0a\u3067\u4f7f\u3046\u3068\u3057\u305f\u3089\uff0c\u5e73\u9762\u30d9\u30af\u30c8\u30eb\u3068\u3057\u3066\u306e2\u6b21\u5143\u304b\u3089\u7a7a\u9593\u306e\u30a2\u30d5\u30a3\u30f3\u5909\u63db\u3068\u3057\u3066\u306e4×4\u6b21\u5143\u3068\u3044\u3046\u3054\u304f\u4f4e\u6b21\u5143<\/strong>\u304b\u3064\u56fa\u5b9a\u6b21\u5143<\/strong>\u306a\u5834\u5408\u304c\u307b\u3068\u3093\u3069\u3067\u3057\u3087\u3046\uff0e\u3053\u3046\u4f4e\u6b21\u5143\u3068\u3082\u306a\u308b\u3068\uff0c\u30c7\u30fc\u30bf\u69cb\u9020\u3084\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u3082\u9ad8\u6b21\u5143\u7528\u306e\u3082\u306e\u3067\u306f\u975e\u5e38\u306b\u975e\u52b9\u7387\u306b\u306a\u308a\u307e\u3059\uff0e\u3061\u3063\u3055\u3044\u9818\u57df\u304c\u5185\u90e8\u3067\u3046\u3058\u3083\u3046\u3058\u3083\u3068\u52d5\u7684\u30e1\u30e2\u30ea\u78ba\u4fdd\u306a\u3093\u304b\u3055\u308c\u308b\u3068\u306a\u3093\u304b\u3082\u3046\u306d\uff0e
\n\u305d\u3093\u306a\u3068\u304d\u6027\u80fd\u9762\u3084\u6271\u3044\u3084\u3059\u3055\u306e\u9762\u3067\u4fbf\u5229\u306a\u306e\u304c\uff0c\u56fa\u5b9a\u6b21\u5143\u30d9\u30af\u30c8\u30eb\u578bublas::bounded_vector\u3067\u3059(\u6b21\u5143\u3054\u3068\u306e\u6700\u9069\u306a\u6570\u5024\u8a08\u7b97\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u306e\u8a71\u306f\u7f6e\u3044\u3068\u304f\u3068\u3057\u3066)\uff0e<\/p>\n
-ublas::bounded_vector
\n-\u6b21\u5143\u306f\u30c6\u30f3\u30d7\u30ec\u30fc\u30c8\u5f15\u6570N\u3067\u4e0e\u3048\u308b(\u30b3\u30f3\u30d1\u30a4\u30eb\u6642\u306b\u56fa\u5b9a\uff01)<\/strong>
\n-\u5185\u90e8\u3067\u56fa\u5b9a\u9577\u914d\u5217<\/strong>\u3092\u7528\u3044\u308b(\u52d5\u7684\u30e1\u30e2\u30ea\u78ba\u4fdd\u3092\u884c\u308f\u306a\u3044\uff01)<\/strong><\/p>\n
\u3068\u3044\u3046\u7279\u5fb4\u3092\u6301\u3063\u3066\u3044\u3066\uff0c<\/p>\n
-\u4f4e\u6b21\u5143\u304b\u3064\u56fa\u5b9a\u6b21\u5143\u3068\u3044\u3046\u76ee\u7684\u306b\u3074\u3063\u305f\u308a
\n-\u5341\u5206\u6b21\u5143\u304c\u4f4e\u3051\u308c\u3070\uff0c\u6027\u80fd\u304a\u3088\u3073\u6271\u3044\u306e\u9762\u3067\u30d7\u30ea\u30df\u30c6\u30a3\u30d6\u578b\u306e\u611f\u899a\u3067\u4f7f\u3048\u308b<\/p>\n
\u3068\u3044\u3046\u5229\u70b9\u304c\u3042\u308a\uff0c\u304b\u306a\u308a\u30ac\u30f3\u30ac\u30f3\u4f7f\u3048\u307e\u3059\uff0e
\n\u3082\u3061\u308d\u3093\uff0c\u901a\u5e38\u306eublas::vector\u3068\u6df7\u5728\u3057\u3066\u4f7f\u3046\u3053\u3068\u304c\u3067\u304d\u307e\u3059\uff0e<\/p>\n
\u5b9f\u969b\u306b\u5c11\u3057\u3060\u3051\u4f7f\u3063\u3066\u307f\u305f\u4f8b\uff0e<\/p>\n
\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\t\/\/\u5185\u90e8\u3067\u56fa\u5b9a\u9577\u914d\u5217\u306a\u306e\u3067\uff0c\u4f7f\u3046\u9818\u57df\u306f\u6b21\u5143\u6570\u00d7sizeof(T)+4(\u6b21\u5143\u6570\u3082\u4fdd\u7ba1\u3059\u308b\u305f\u3081)\r\n\tstd::cout << sizeof(bounded_vector<int, 4>) << std::endl;\r\n\tstd::cout << sizeof(bounded_vector<float, 4>) << std::endl;\r\n\tstd::cout << sizeof(bounded_vector<double, 4>) << std::endl;\r\n\tstd::cout << sizeof(bounded_vector<double, 8>) << std::endl;\r\n\r\n\t\/\/\u5185\u90e8\u3067\u56fa\u5b9a\u9577\u914d\u5217\u306a\u306e\u3067\uff0c\r\n\t\/\/vector<char> huge1(16777216);\t\t\t\/\/\u30d2\u30fc\u30d7\u3067\u306a\u3089\u306a\u3093\u3067\u3082\u306a\u3044\u91cf(16MB)\u3067\u3082\r\n\t\/\/bounded_vector<char, 16777216> huge2;\t\t\/\/\u30b9\u30bf\u30c3\u30af\u306b\u914d\u7f6e\u3055\u308c\u3066\u308b\u3068segfault\u203b\u30b9\u30bf\u30c3\u30af\u30b5\u30a4\u30ba\u4f9d\u5b58\r\n\r\n\t\/\/bounded_vector\u3068\r\n\tbounded_vector<double, 3> vec1;\r\n\tvec1[0] = 1; vec1[1] = 2; vec1[2] = 3;\r\n\r\n\t\/\/\u3075\u3064\u3046\u306evector\u306e\r\n\tvector<double> vec2(3);\r\n\tvec2[0] = 1; vec2[1] = 2; vec2[2] = 3;\r\n\r\n\t\/\/\u6df7\u5728\u3057\u305f\u6f14\u7b97\u3082\u4f59\u88d5\r\n\tstd::cout << (vec1 + vec2) << std::endl;\r\n\r\n\t\/\/\u3068\u3044\u3046\u304b\u30af\u30e9\u30b9\u5b9a\u7fa9\u3092\u898b\u305f\u3089\u306a\u30fc\u3093\u3060\u3063\u3066\u306a\u308a\u307e\u3059\r\n}\r\n<\/pre>\n
\u3067<\/p>\n
\r\n 20\r\n 20\r\n 36\r\n 68\r\n [][3](2,4,6)[]\r\n<\/pre>\n
\u3044\u3044\u304b\u3093\u3058\u3067\u3059\u306d\uff0e\u4f4e\u6b21\u5143\u56fa\u5b9a\u6b21\u5143\u3067\u4f7f\u3044\u3084\u3059\u305d\u3046\uff0e
\n\u3067\u3082\u306d\uff0c\u4f7f\u3046\u3068\u304d\u306f\u3061\u3087\u3063\u3068\u6c17\u3092\u3064\u3051\u3066\u6b32\u3057\u3044\u3093\u3067\u3059\uff0e<\/p>\n
*\u3061\u3087\u3063\u3068\u8a2d\u8a08\u304c\u4e2d\u9014\u534a\u7aef\u306a\u6c17\u3082\u3059\u308b\uff01<\/p>\n
ublas::bounded_vector\u306f\uff0c\u6b21\u306e\u30b3\u30fc\u30c9\u306e\u30b3\u30f3\u30d1\u30a4\u30eb\u306b\u6210\u529f<\/strong>\u3057\uff0c\u5b9f\u884c\u6642\u306b\u4f8b\u5916<\/strong>\u3092\u5410\u304d\u307e\u3059\uff0e<\/p>\n
\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\tbounded_vector<double, 3> vec3;\r\n\tbounded_vector<double, 4> vec4;\r\n\r\n\tstd::cout << (vec3 + vec4) << std::endl;\r\n}\r\n<\/pre>\n
\r\nCheck failed in file \/usr\/local\/include\/boost\/numeric\/ublas\/vector_expression.hpp at line 548:\r\nsize1 == size2\r\nterminate called after throwing an instance of 'boost::numeric::ublas::bad_argument'\r\n  what():  bad argument\r\nzsh: abort      .\/a.out\r\n<\/pre>\n
\u304a\u3044\u3061\u3087\u3063\u3068\u5f85\u3066\u3068\uff0e
\n\u306a\u3093\u3067\uff0e\u30b3\u30f3\u30d1\u30a4\u30eb\u6642\u306b\u30a8\u30e9\u30fc\u306b\u3055\u305b\u308c\u3070\u3044\u3044\u3058\u3083\u3093\uff01<\/strong>\u3067\u304d\u308b\u3058\u3083\u3093\uff01<\/p>\n
\u307e\u3041\u3053\u308c\u3060\u3051\u306a\u3089\u6c17\u3092\u3064\u3051\u3066\u308c\u3070\u3067\u7247\u4ed8\u3051\u3089\u308c\u308b\u304b\u306a\u3063\u3066\u6c17\u3082\u3057\u307e\u3059\uff0e<\/p>\n
\u3067\uff0c\u6b21\u306e\u30b3\u30fc\u30c9\u306f\u30b3\u30f3\u30d1\u30a4\u30eb\u304c\u901a\u308a\u4f8b\u5916\u3082\u51fa\u307e\u305b\u3093\uff0e<\/p>\n
\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\tbounded_vector<double, 3> vec3;\r\n\tbounded_vector<double, 4> vec4;\r\n\r\n\tvec4 = vec3;\r\n\r\n\tstd::cout << vec3 << std::endl;\r\n\tstd::cout << vec4 << std::endl;\r\n}\r\n<\/pre>\n
\r\n[3](-1.08785e-41,-4.66095e-42,-1.0874e-41)\r\n[3](-1.08785e-41,-4.66095e-42,-1.0874e-41)\r\n<\/pre>\n
\u3084\u30fc\u3093\uff1e\uff1c
\n\u6b21\u6570\u304c\u9055\u3046\u306e\u306b\u4ee3\u5165\u3067\u304d\u308b\u306e\uff1e\uff1c
\nvec4\u306e\u6b21\u6570\u304c3\u306b\u306a\u3063\u3066\u308b\uff1e\uff1c<\/p>\n
\u4e00\u65b9\u3067\uff0c<\/p>\n
\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\tbounded_vector<double, 3> vec3;\r\n\tbounded_vector<double, 4> vec4;\r\n\r\n\tvec3 = vec4;\t\/\/\u4ee3\u5165\u3092\u9006\u306b\u3057\u3066\u307f\u305f\r\n\r\n\tstd::cout << vec3 << std::endl;\r\n\tstd::cout << vec4 << std::endl;\r\n}\r\n<\/pre>\n
\r\nCheck failed in file \/usr\/local\/include\/boost\/numeric\/ublas\/storage.hpp at line 343:\r\nsize &lt;= N\r\nterminate called after throwing an instance of 'boost::numeric::ublas::bad_size'\r\n  what():  bad size\r\nzsh: abort      .\/a.out\r\n<\/pre>\n
\u3053\u308c\u306f\u843d\u3061\u307e\u3059\u304b\u2026<\/p>\n
\u3067\u3082\u2193\u3053\u3046\u66f8\u3044\u3061\u3083\u3048\u3070\uff0c\u9ad8\u6b21\u5143\u2192\u4f4e\u6b21\u5143\u3082\u3044\u3051\u308b\u3057\uff0c\u7570\u306a\u308b\u6b21\u5143\u306e\u6f14\u7b97\u3082\u3044\u3051\u307e\u3059\uff0e<\/p>\n
\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\n\r\nusing namespace boost::numeric::ublas;\r\n\r\nint main()\r\n{\r\n\tbounded_vector<double, 3> vec3;\r\n\tbounded_vector<double, 4> vec4;\r\n\r\n\tvec4.resize(3);\r\n\r\n\tvec3 = vec4;\r\n\tstd::cout << (vec3 + vec4) << std::endl;\r\n\tstd::cout << vec3 << std::endl;\r\n\tstd::cout << vec4 << std::endl;\r\n}\r\n<\/pre>\n
\r\n[][3](-2.38555e-41,-0.578328,4.37691e-314)[]\r\n[][3](-1.19277e-41,-0.289164,2.18846e-314)[]\r\n[][3](-1.19277e-41,-0.289164,2.18846e-314)[]\r\n<\/pre>\n
\u3048\u30fc\u3067\u3082\u306a\u3093\u304b\u306a\u3041\u306a\u3093\u304b\u306a\u3041\uff1e\uff1c\uff1b\uff1b<\/p>\n
\u3053\u306e\u8fba\u308a\u3063\u3066\u306e\u306f\u8a2d\u8a08\u601d\u60f3\u304b\u3089\u3059\u308b\u3068\u5f53\u7136\u3068\u8a00\u3048\u3070\u5f53\u7136\u3067\uff0c\u8981\u306fbounded_vector\u306f\u300c\u30b9\u30c8\u30ec\u30fc\u30b8\u3092\u9577\u3055N\u306e\u56fa\u5b9a\u9577\u914d\u5217\u3068\u3059\u308b\u3053\u3068\u3067\uff0c\u6700\u5927N\u6b21\u5143\u306e\u30d9\u30af\u30c8\u30eb\u3092\u8868\u73fe\u3067\u304d\u308b\u3088\u3046\u306b\u3057\u305f\u300d\u3082\u306e\u3067\u3042\u3063\u3066\uff0c\u300c\u6b21\u5143\u3092\u56fa\u5b9a\u5024N\u306b\u3057\u305f\u300d\u3068\u3044\u3046\u308f\u3051\u3067\u306f\u306a\u3044<\/strong>\u306e\u3067\u3059\u306d\u30fc\uff0e<\/p>\n
\u3046\u30fc\u3093\uff0c\u3067\u3082\u3084\u3063\u3071\u308a<\/p>\n
\r\nbounded_vector<double, 4> vec;\r\n<\/pre>\n
\u3066\u66f8\u3044\u3066\u3042\u3063\u305f\u3089\uff0cvec\u306f4\u6b21\u5143\u306e\u56fa\u5b9a\u9577\u30d9\u30af\u30c8\u30eb\u3060\u3057\u305d\u308c\u306b\u5bfe\u3057\u30663\u6b21\u5143\u30d9\u30af\u30c8\u30eb\u3092\u4ee3\u5165\u3057\u305f\u308a\u3057\u305f\u3089\u30b3\u30f3\u30d1\u30a4\u30eb\u6642\u3067\u30a8\u30e9\u30fc\u306b\u306a\u308b\u3093\u3060\u308d\u3046\u3068\u8ab0\u3067\u3082\u671f\u5f85\u3057\u3061\u3083\u3046\u6c17\u304c\u3057\u307e\u3059\uff0e
\n\u7406\u89e3\u304c\u6d45\u3044\u306a\u3089\u4f7f\u3046\u306a\u3063\u3066\u3044\u3046\u8a71\u3067\u578b\u4ed8\u3051<\/del>\u7247\u4ed8\u3051\u3061\u3083\u3063\u3066\u3082\u3044\u3044\u3093\u3067\u3059\u304c\uff0e<\/p>\n
\u3068\u3044\u3046\u304b\uff0c\u4e00\u3064\u306e(\u5b9f\u8cea\u30d7\u30ea\u30df\u30c6\u30a3\u30d6\u306a)\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u30923\u6b21\u5143\u306b\u306f\u4f7f\u3046\u308f4\u6b21\u5143\u306b\u306f\u4f7f\u3046\u308f\uff0c\u305d\u306e\u6642\u70b9\u3067\u7fd2\u6163\u3068\u3057\u3066\u306f\u826f\u304f\u306a\u3044\u3057\uff0c\u305d\u306e\u3046\u3048(\u4e0a\u306e\u4f8b\u3067\u306f)5\u6b21\u5143\u4ee5\u4e0a\u304b\u3089\u306f\u4f8b\u5916\u5410\u304f\u306e\u3067\u30a2\u30a6\u30c8\u3068\u304b\uff0c\u3084\u3063\u3071\u308a\u4f7f\u3044\u306b\u304f\u3055\u306e\u5143\u306b\u306a\u308b\u3093\u3058\u3083\u306a\u3044\u3067\u3057\u3087\u3046\u304b\uff0e
\n\u300c\u4f7f\u3044\u56de\u3057\u306f\u60aa\u3057\u304d\u7fd2\u6163\u300d\u3066\u3053\u3068\u3067const\u6559\u306e\u5e03\u6559\u306b\u3064\u306a\u3052\u308b\u3053\u3068\u3082\u3067\u304d\u306a\u304f\u306f\u306a\u3044\u3067\u3059\u306d\uff0e<\/p>\n
\u3053\u308ctypedef\u3068\u304b\u95a2\u6570\u6e21\u3057\u3068\u304b\u3067\u3061\u3087\u3063\u3068\u8907\u96d1\u306b\u306a\u308b\u3068\uff0c\u305a\u3044\u3076\u3093\u3068\u308f\u304b\u308a\u306b\u304f\u3044\u30d0\u30b0\u306b\u306a\u308a\u307e\u3059\uff0e
\n\u6b21\u306e\u30d7\u30ed\u30b0\u30e9\u30e0\u306f\u3044\u308d\u3093\u306a\u610f\u5473\u3067\u5168\u304f\u7121\u610f\u5473\u3067\u3059\u304c\uff0c\u3061\u3087\u3063\u3068\u30c7\u30d0\u30c3\u30ac\u7121\u3057\u3067\u3069\u3053\u3067\u843d\u3061\u308b\u304b\u5f53\u3066\u3089\u3093\u306a\u3044\u3093\u3058\u3083\u306a\u3044\u3067\u3057\u3087\u3046\u304b\uff0e
\n\u3053\u3053\u3067\u306f\u308f\u3056\u3068\u3089\u3057\u304f\u9593\u9055\u3048\u3066\u308b\u306e\u3067\u3059\u304c\uff0c\u4ee3\u5165\u306b\u4ee3\u5165\u3092\u7d4c\u3066\u305d\u306e\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u306f\u4eca\u4f55\u6b21\u5143\u30d9\u30af\u30c8\u30eb\u306a\u306e\u3063\u3066\u306e\u304c\u66d6\u6627\u306a\u72b6\u614b\u3067\u3084\u3089\u304b\u3057\u3066\u3057\u307e\u3044\u307e\u3059(\u3084\u3089\u304b\u3057\u3066\u3057\u307e\u3044\u307e\u3057\u305f)\uff0e<\/p>\n
\r\n\/\/\u3069\u3053\u3067\u843d\u3061\u308b\u3067\u3057\u3087\u3046\r\n#include <iostream>\r\n#include <boost\/numeric\/ublas\/vector.hpp>\r\n#include <boost\/numeric\/ublas\/io.hpp>\r\nusing namespace boost::numeric::ublas;\r\n\r\ntypedef bounded_vector<double, 3> vec3d;\r\ntypedef bounded_vector<double, 4> vec4d;\r\n\r\nvec4d unit(vec4d v)\r\n{\r\n\tdouble norm = norm_2(v);\r\n\treturn v \/ norm;\r\n}\r\n\r\ndouble angle(vec4d v1, vec4d v2)\r\n{\r\n\tv1 = unit(v1);\r\n\tv2 = unit(v2);\r\n\treturn acos(inner_prod(v1, v2));\r\n}\r\n\r\nint main()\r\n{\r\n\tvec3d v1;\r\n\tv1[0] = 1; v1[1] = 1; v1[2] = 0;\r\n\r\n\tvec4d v2;\r\n\tv2[0] = 0; v2[1] = 1; v2[2] = 1;\r\n\r\n\tstd::cout << angle(v1, v2) << std::endl;\r\n}\r\n<\/pre>\n
\u305d\u3046\uff0c\u30a2\u30db\u306a\u3093\u3067\u3059\u3088\uff0c\u3053\u3093\u306a\u3057\u3087\u30fc\u3082\u306a\u3044\u3053\u3068\u3067\u3064\u307e\u3065\u3044\u305f\u50d5\u304c\u30a2\u30db\u306a\u3093\u3067\u3059\uff0e
\n\u3067\u3082\u306a\u3093\u3068\u306a\u304f\u305d\u308c\u3067\u6e08\u307e\u305b\u305f\u304f\u306a\u304b\u3063\u305f\u2026<\/p>\n
*\u5922\u3092\u898b\u305f\u3044\u306a\u3089Eigen\u3067\u3069\u3046\u305e<\/p>\n
\u3053\u3053\u3067\uff0c\u7dda\u5f62\u4ee3\u6570\u30e9\u30a4\u30d6\u30e9\u30eaEigen\u306e\u7d39\u4ecb\uff08\u3061\u3087
\nuBLAS\u3068\u540c\u3058\u304f\u30d0\u30ea\u30d0\u30ea\u306eC++\u3067\uff0c\u30d8\u30c3\u30c0\u30aa\u30f3\u30ea\u30fc\u306a\u30e9\u30a4\u30d6\u30e9\u30ea\u3067\u3059\uff0e
\n\u884c\u5217\u3084\u30d9\u30af\u30c8\u30eb\u306e\u6027\u8cea(\u6b21\u6570\u3084sparse\u304bdense\u304b\uff0c\u306a\u3069\u3044\u308d\u3044\u308d)\u306b\u4f75\u305b\u305f\u30c7\u30fc\u30bf\u69cb\u9020\u3068\u30a2\u30eb\u30b4\u30ea\u30ba\u30e0\u3092\u6301\u3064\u3053\u3068\uff0c\u4f4e\u6b21\u5143\u30fb\u56fa\u5b9a\u6b21\u5143\u306e\u305f\u3081\u306b\u30d7\u30ea\u30df\u30c6\u30a3\u30d6\u306a\u578b\u306e\u611f\u899a\u3067\u4f7f\u3048\u308b\u30d9\u30af\u30c8\u30eb\u578b\u3092\u3082\u3064\u3053\u3068\u306fuBLAS\u3068\u540c\u3058\u3067\u3059\uff0e
\n\u3055\u3089\u306b
\n-\u771f\u306b\u578b\u30ec\u30d9\u30eb\u3067\u56fa\u5b9a\u9577
\n–\u9055\u3046\u6b21\u6570\u3078\u306e\u4ee3\u5165\u3084\u6f14\u7b97\u306a\u3093\u304b\u306f\u5168\u90e8\u578b\u30ec\u30d9\u30eb\u3067\u30cf\u30b8\u304f
\n-\u884c\u5217\u3082\u56fa\u5b9a\u6b21\u5143\u306b\u3067\u304d\u308b
\n–\u3068\u3044\u3046\u304b\u30d9\u30af\u30c8\u30eb\u81ea\u4f53\u304c1\u884c\u307e\u305f\u306f1\u5217\u306e\u884c\u5217
\n–\u3064\u307e\u308a\u884c\u30d9\u30af\u30c8\u30eb\u3082\u5217\u30d9\u30af\u30c8\u30eb\u3082\u8868\u73fe\u3067\u304d\u308b
\n-\u6b21\u6570\u3092\u578b\u3067\u7ba1\u7406\u3059\u308b\u306e\u3067\u30d5\u30a3\u30fc\u30eb\u30c9\u306b\u4fdd\u6301\u3059\u308b\u5fc5\u8981\u304c\u306a\u304f\uff0cublas::bounded_vector\u3088\u308a\u5e38\u306b4\u30d0\u30a4\u30c8\u5206\u304a\u30c8\u30af
\n-\u7279\u306b\u6b21\u5143\u304c2\u301c4\u306e\u3068\u304d\u306e\u30d9\u30af\u30c8\u30eb\u306e\u5834\u5408\u306e\u6271\u3044\u3084\u3059\u3055\u306b
\n–\u6dfb\u5b57\u3058\u3083\u306a\u304fx,y,z,w\u3068\u3044\u3046\u30e1\u30bd\u30c3\u30c9(\u53c2\u7167\u8fd4\u3057\u306a\u306e\u3067\u4ee3\u5165\u3082\u53ef)\u3067\u30a2\u30af\u30bb\u30b9\u3067\u304d\u308b
\n—\u305f\u3060\u4f8b\u3048\u30702\u6b21\u5143\u306e\u3068\u304dz\u3084w\u306b\u30a2\u30af\u30bb\u30b9\u3059\u308b\u3068\u30b3\u30f3\u30d1\u30a4\u30eb\u30a8\u30e9\u30fc\u3058\u3083\u306a\u304f\u5b9f\u884c\u6642\u30a8\u30e9\u30fc\u306b\u306a\u308b\u3051\u3069\u3053\u308c\u3046\u307e\u304f\u3084\u3063\u3066\u306a\u3093\u3068\u304b\u3067\u304d\u306a\u3044\u3060\u308d\u3046\u304b
\n–\u57fa\u672c\u7684\u306a\u578b\u304ctypedef\u3055\u308c\u3066\u3044\u308b(Vector2d\u3068\u304bMatrix3d\u3068\u304b)
\n–\u30b3\u30f3\u30b9\u30c8\u30e9\u30af\u30bf\u306b\u6b21\u5143\u6570\u3060\u3051\u306e\u5024\u3092\u6e21\u3059\u3053\u3068\u3067\u30d9\u30af\u30c8\u30eb\u306e\u521d\u671f\u5024\u3092\u8a2d\u5b9a\u3067\u304d\u308b<\/p>\n
\u3053\u3093\u306a\u96f0\u56f2\u6c17\uff0e<\/p>\n
\r\n#include <iostream>\r\n#include <Eigen\/Core>\r\n\r\n#include <boost\/format.hpp>\r\n\r\nint main()\r\n{\r\n\tstd::cout << sizeof(Eigen::Matrix<int, 4, 1>) << std::endl;\r\n\tstd::cout << sizeof(Eigen::Matrix<double, 4, 1>) << std::endl;\r\n\tstd::cout << sizeof(Eigen::Matrix<double, 8, 1>) << std::endl;\r\n\tstd::cout << sizeof(Eigen::Matrix<double, 8, 8>) << std::endl;\r\n\r\n\t\/\/Eigen::Vector\r\n\tEigen::Matrix3d m;\r\n\tstd::cout << m << std::endl;\r\n\r\n\tEigen::Vector2d v2d(0, 1);\r\n\t\/\/Eigen::Vector3d v3d = v2d;\t\/\/\u3061\u3083\u3093\u3068\u30b3\u30f3\u30d1\u30a4\u30eb\u30a8\u30e9\u30fc\uff0e(\u30b3\u30d4\u30b3\u30f3\u3058\u3083\u306a\u304f)\u4ee3\u5165\u3067\u3082\u3082\u3061\u308d\u3093\r\n\tEigen::Vector3d v3d(0, 1, 2);\r\n\r\n\t\/\/\u6f14\u7b97\u3082\u578b\u304c\u9055\u3046\u3068\u30a8\u30e9\u30fc\r\n\t\/\/std::cout << (v2d + v3d) << std::endl;\r\n\t\/\/Eigen::Matrix<double, 3, 3>() * Eigen::Matrix<double, 4, 4>();\r\n\r\n\t\/\/\u5358\u306b\u578b\u304c\u9055\u3046\u3068\u30a8\u30e9\u30fc\u306a\u3093\u3058\u3083\u306a\u304f\u3066\uff0c\u7b97\u8853\u7684\u306b\u6709\u52b9\u306a\u578b\u540c\u58eb\u306a\u3089\u554f\u984c\u306a\u3044\r\n\tEigen::Matrix<double, 3, 1>() * Eigen::Matrix<double, 1, 3>();\r\n\r\n\t\/\/Eigen\u3067\u306f\u30d9\u30af\u30c8\u30eb\u306f1\u884c\u307e\u305f\u306f1\u5217\u306e\u884c\u5217\r\n\tEigen::Matrix<double, 3, 1> vec = Eigen::Vector3d(0, 0, 0);\t\/\/\u306a\u306e\u3067\u4ee3\u5165\u3067\u304d\u308b\r\n\t\r\n\r\n\t\/\/4\u6b21\u5143\u307e\u3067\u306f\u30b3\u30f3\u30b9\u30c8\u30e9\u30af\u30bf\u3067\u521d\u671f\u5024\u3092\u6307\u5b9a\u3067\u304d\uff0cx,y,z,w\u3068\u3044\u3046\u30e1\u30f3\u30d0\u3067\u53c2\u7167(\u5de6\u8fba\u306b\u7f6e\u3044\u3066\u4ee3\u5165\u3059\u308b\u3053\u3068\u3082)\u3067\u304d\u308b\r\n\tEigen::Vector2f v2f(0, 1);\r\n\tstd::cout << boost::format("{%1%,%2%}") % v2f.x() % v2f.y() << std::endl;\r\n\r\n\tEigen::Vector3f v3f(0, 1, 2);\r\n\tstd::cout << boost::format("{%1%,%2%,%3%}") % v3f.x() % v3f.y() % v3f.z() << std::endl;\r\n\r\n\tEigen::Vector4f v4f(0, 1, 2, 3);\r\n\tstd::cout << boost::format("{%1%,%2%,%3%,%4%}") % v4f.x() % v4f.y() % v4f.z() % v4f.w() << std::endl;\r\n}\r\n<\/pre>\n
\r\n 16\r\n 32\r\n 64\r\n 512\r\n -7.30135e-42 -1.27833e-41 -1.37662e-41\r\n -1.27833e-41 -7.30002e-42 3.64203e-314\r\n 9.38725e-323 -6.62786e-42 -1.27992e-41\r\n {0,1}\r\n {0,1,2}\r\n {0,1,2,3}\r\n<\/pre>\n
*\u307e\u3068\u3081<\/p>\n
\u3042\u308c\uff0c\u7d50\u5c40\u4f55\u304c\u8a00\u3044\u305f\u304b\u3063\u305f\u3093\u3060\u308d\u3046\uff0e
\n\u3042\uff0c\u601d\u3044\u51fa\u3057\u305f\uff0e<\/p>\n
-uBLAS\u306b\u306f\u4f4e\u6b21\u5143\u3067\u6271\u3044\u3084\u3059\u3044bounded_vector\u304c\u3042\u308b\u3088\uff01
\n-\u3067\u3082\u6b21\u6570\u5468\u308a\u3067\u3061\u3087\u3063\u3068\u30cf\u30de\u308a\u3084\u3059\u3044\u304b\u3082\u3060\u304b\u3089\u6c17\u3092\u3064\u3051\u3066\u306d\uff01
\n-Eigen\u3082\u4fbf\u5229\u3060\u3088\uff01
\n-\u4fbf\u5229\u3060\u304b\u3089\u3063\u3066\uff0c\u6b21\u6570\u3092\u6570\u767e\u3068\u304b\u306b\u3057\u3066\u4e71\u7528\u3059\u308b\u3053\u3068\u306a\u3044\u3088\u3046\u306b\u306d\uff01<\/p>\n
\u3055\u3066\uff0c\u9045\u523b\u306b\u3050\u3060\u3050\u3060\u306a\u5185\u5bb9\u3068\u5fb9\u5e95\u7684\u306b\u30cf\u30fc\u30c9\u30eb\u3092\u4e0b\u3052\u3066\u3057\u307e\u3063\u305f\u306e\u3067\uff0c\u304a\u6b21\u306e[https:\/\/twitter.com\/#!\/yutopp:title=@yutopp]\u3055\u3093\uff0c\u3088\u308d\u3057\u304f\u304a\u9858\u3044\u3057\u307e\u3059\uff01\uff01<\/p>\n","protected":false},"excerpt":{"rendered":"
\u3053\u306e\u8a18\u4e8b\u306f\uff0c[http:\/\/partake.in\/events\/597a0fc3-0e3a-47a3-8fc3-4f32ad846a3d:title=Boost Advent Calendar 2011]\u306e14\u65e5\u76ee\u306e\u8a18\u4e8b […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[24,14,12],"class_list":["post-681","post","type-post","status-publish","format-standard","hentry","category-1","tag-adventcalendar","tag-boost","tag-c"],"_links":{"self":[{"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=\/wp\/v2\/posts\/681","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=681"}],"version-history":[{"count":0,"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=\/wp\/v2\/posts\/681\/revisions"}],"wp:attachment":[{"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=681"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=681"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/daily.belltail.jp\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=681"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}