{"id":2971,"date":"2017-07-13T17:42:54","date_gmt":"2017-07-13T08:42:54","guid":{"rendered":"https:\/\/na-inet.jp\/weblog2\/?p=2971"},"modified":"2017-07-13T17:45:54","modified_gmt":"2017-07-13T08:45:54","slug":"%e7%a0%b4%e7%b6%bb%e3%81%99%e3%82%8b%e5%8c%ba%e9%96%93%e6%bc%94%e7%ae%97%e3%81%ae%e4%ba%8b%e4%be%8b","status":"publish","type":"post","link":"https:\/\/na-inet.jp\/weblog2\/2017\/07\/13\/%e7%a0%b4%e7%b6%bb%e3%81%99%e3%82%8b%e5%8c%ba%e9%96%93%e6%bc%94%e7%ae%97%e3%81%ae%e4%ba%8b%e4%be%8b\/","title":{"rendered":"\u7834\u7dbb\u3059\u308b\u533a\u9593\u6f14\u7b97\u306e\u4e8b\u4f8b"},"content":{"rendered":"<p><br \/>\n\u3000\u7cbe\u5ea6\u4fdd\u8a3c\u8a08\u7b97\uff08\u3068\u3044\u3046\u8a00\u8449\u306f\u3042\u307e\u308a\u597d\u304d\u3067\u306f\u306a\u3044\u306e\u3060\u3051\u308c\u3069\uff09\u306e\u4e8b\u4f8b\u3082\u7d39\u4ecb\u3057\u3088\u3046\u304b\u3068\uff0c<a href=\"http:\/\/arblib.org\/\" target=\"_blank\">ARB<\/a>\uff08\u534a\u5f84\uff09\u3068<a href=\"https:\/\/gforge.inria.fr\/frs\/?group_id=157\" target=\"_blank\">MPFI<\/a>\uff08\u533a\u9593\uff09\u3092\u4f7f\u3063\u3066\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u5199\u50cf\\(x_{i+1} := 4x_i(1-x_i)\\)\u306e\u8a08\u7b97\u3092\u3055\u305b\u3066\u307f\u305f\u3068\u3053\u308d\uff0c\u6841\u6570\u304c\u8db3\u308a\u306a\u3044\u3068\u7834\u7dbb\u3059\u308b\u3053\u3068\u304c\u5206\u304b\u3063\u305f\uff08\u4eca\u66f4\uff09\u3002\u3064\u30fc\u3053\u3068\u3067\uff0c128bits\u8a08\u7b97\u306e\u4e8b\u4f8b\uff08\u4e0a\uff09\u3068256bits\u306e\u4e8b\u4f8b\uff08\u4e0b\uff09\u3092\u3070\u3002<\/p>\n<p><a href=\"https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843.jpg\"><img loading=\"lazy\" decoding=\"async\" src=\"\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843-212x300.jpg\" alt=\"\" width=\"212\" height=\"300\" class=\"alignnone size-medium wp-image-2974\" srcset=\"https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843-212x300.jpg 212w, https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843-106x150.jpg 106w, https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843-768x1087.jpg 768w, https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843-723x1024.jpg 723w, https:\/\/na-inet.jp\/wordpress\/wp-content\/uploads\/2017\/07\/2017-07-13_170843.jpg 837w\" sizes=\"auto, (max-width: 212px) 100vw, 212px\" \/><\/a><\/p>\n<p>\u3000\u6709\u52b9\u6841\u6570\u304c\u306a\u304f\u306a\u3063\u305f\u9014\u7aef\u306b\u533a\u9593\u3084\u534a\u5f84\u304c\u6025\u6fc0\u306b\u30c7\u30ab\u304f\u306a\u3063\u3066\u7834\u7dbb\u3057\u3066\u3044\u308b\u3053\u3068\u304c\u4e00\u76ee\u3067\u308f\u304b\u308b\u306e\u306f\u3044\u3044\u3068\u3057\u3066\uff0c\u7cbe\u5ea6\u306f\u3068\u3082\u304b\u304f\u8a08\u7b97\u306f\u6700\u5f8c\u307e\u3067\u5b9f\u884c\u3057\u305f\u3044\u6642\u306b\u306f\u4f7f\u3048\u3093\u3002\u3068\u308a\u3042\u3048\u305a\uff0c\u3053\u3046\u3044\u3046\u6841\u843d\u3061\u304c\u9177\u3044\u4e8b\u4f8b\u306b\u306f\u533a\u9593\u6f14\u7b97\u3060\u308d\u3046\u3068\u305d\u3046\u3067\u306a\u304b\u308d\u3046\u3068\u6841\u6570\u5897\u3084\u3057\u3066\u3084\u3093\u306a\u3055\u3044\u3068\u3044\u3046\u7d50\u8ad6\u306b\u306a\u308b\u306a\u3002<\/p>\n<p>\u3000\u3042\u3068\uff0c\u5b9f\u4fc2\u65702\u6b21\u65b9\u7a0b\u5f0f\u306e\u4e8b\u4f8b\u3082\u4f5c\u3063\u3066\u307f\u305f\u304c\uff0c\u5224\u5225\u5f0f\u3067\u5b9f\u6570\u89e3\u3068\u8907\u7d20\u6570\u89e3\u306e\u533a\u5225\u3092\u3059\u308bif\u6587\u3092\u3069\u3046\u3057\u3088\u3046\u3068\u304b\uff0c\u5b9f\u7528\u306b\u5f53\u305f\u3063\u3066\u306f\u3084\u306f\u308a\u30e1\u30f3\u30c9\u30af\u30b5\u30a4\u3053\u3068\u3092\u8003\u3048\u306a\u3044\u3068\u3044\u3051\u306a\u3044\u306e\u3067\u305d\u306e\u8fba\u3092\u3069\u3046\u304b\u3057\u3066\u5de5\u592b\u3059\u308b\u3068\u7814\u7a76\u306b\u306a\u308b\u8a33\u3067\uff0c\u3053\u3063\u304b\u3089\u5148\u306f\u7cbe\u5ea6\u4fdd\u8a3c\u306e\u65b9\u3005\u306b\u304a\u4efb\u305b\u3057\u3088\u3046\u3063\u3068\u3002\u307e\u3041\u300c\u6700\u521d\u304b\u3089\u8907\u7d20\u4fc2\u6570\u3067\u3084\u308c\u300d\u3068\u3044\u3046\u7d50\u8ad6\u306b\u306a\u308b\u3093\u3060\u308d\u3046\u3051\u3069\uff0c\u305d\u308c\u306fARB\u306e\u6f14\u7fd2\u554f\u984c\u306b\u3057\u3066\u304a\u3051\u3070\u3044\u3044\u304b\u3002<\/p>\n<p>\u3000\u3064\u30fc\u3053\u3068\u3067\uff0c\u4f7f\u7528\u3057\u305f\u30d7\u30ed\u30b0\u30e9\u30e0(logistic_arb.c, logisitc_mpfi.c)\u3092\u6700\u5f8c\u306b\u63b2\u8f09\u3057\u3066\u304a\u304f\u3002\u4f7f\u7528\u6642\u306b\u306fARB, MPFI\u306e\u4ed6\uff0cMPFR\/GMP\u304c\u5fc5\u8981\u3002<\/p>\n<pre class=\"lang:c decode:true \" title=\"logisitc_arb.c\" >\/\/ logistic \u5199\u50cf\r\n\/\/ ARB\u7248\r\n#include &lt;stdio.h&gt;\r\n#include \"arb.h\"\r\n\r\nint main()\r\n{\r\n\tint i;\r\n\tslong prec;\r\n\tarb_t x[102];\r\n\r\n\t\/\/ \u521d\u671f\u5316\r\n\tfor(i = 0; i &lt; 102; i++)\r\n\t\tarb_init(x[i]);\r\n\r\n\tprintf(\"prec(bits) = \"); scanf(\"%ld\", &amp;prec);\r\n\r\n\t\/\/ \u521d\u671f\u5024\r\n\t\/\/x[0] = 0.7501;\r\n\tarb_set_str(x[0], \"0.7501\", prec);\r\n\r\n\tfor(i = 0; i &lt;= 90; i++)\r\n\t{\r\n\t\tif((i % 10) == 0)\r\n\t\t{\r\n\t\t\tprintf(\"%5d, \", i);\r\n\t\t\tarb_printd(x[i], 17);\r\n\t\t\tprintf(\"\\n\");\r\n\t\t}\r\n\r\n\t\t\/\/x[i + 1] = 4 * x[i] * (1 - x[i]);\r\n\t\tarb_sub_ui(x[i + 1], x[i], 1UL, prec);\r\n\t\tarb_neg(x[i + 1], x[i + 1]);\r\n\t\tarb_mul(x[i + 1], x[i + 1], x[i], prec);\r\n\t\tarb_mul_ui(x[i + 1], x[i + 1], 4UL, prec);\r\n\r\n\t}\r\n\r\n\t\/\/ \u6d88\u53bb\r\n\tfor(i = 0; i &lt; 102; i++)\r\n\t\tarb_clear(x[i]);\r\n\r\n\treturn 0;\r\n}\r\n<\/pre>\n<pre class=\"lang:c decode:true \" title=\"logistic_mpfi.c\" >\/\/ logistic \u5199\u50cf\r\n\/\/ MPFR &amp; MPFI\u7248\r\n#include &lt;stdio.h&gt;\r\n#include \"mpfr.h\"\r\n#include \"mpfi.h\"\r\n\r\nint main()\r\n{\r\n\tint i;\r\n\tunsigned long prec;\r\n\tmpfi_t x[102];\r\n\tmpfr_t relerr;\r\n\r\n\tprintf(\"prec(bits) = \"); scanf(\"%ld\", &amp;prec);\r\n\tmpfr_set_default_prec(prec);\r\n\r\n\t\/\/ \u521d\u671f\u5316\r\n\tmpfr_init(relerr);\r\n\tfor(i = 0; i &lt; 102; i++)\r\n\t\tmpfi_init(x[i]);\r\n\r\n\t\/\/ \u521d\u671f\u5024\r\n\t\/\/x[0] = 0.7501;\r\n\tmpfi_set_str(x[0], \"0.7501\", 10);\r\n\r\n\tfor(i = 0; i &lt;= 90; i++)\r\n\t{\r\n\t\tif((i % 10) == 0)\r\n\t\t{\r\n\t\t\tprintf(\"%5d, \", i);\r\n\t\t\tmpfi_out_str(stdout, 10, 17, x[i]);\r\n\t\t\tmpfi_diam(relerr, x[i]);\r\n\t\t\tmpfr_printf(\"%10.3RNe\\n\", relerr);\r\n\t\t}\r\n\r\n\t\t\/\/x[i + 1] = 4 * x[i] * (1 - x[i]);\r\n\t\tmpfi_ui_sub(x[i + 1], 1UL, x[i]);\r\n\t\tmpfi_mul(x[i + 1], x[i + 1], x[i]);\r\n\t\tmpfi_mul_ui(x[i + 1], x[i + 1], 4UL);\r\n\r\n\t}\r\n\r\n\t\/\/ \u6d88\u53bb\r\n\tmpfr_clear(relerr);\r\n\tfor(i = 0; i &lt; 102; i++)\r\n\t\tmpfi_clear(x[i]);\r\n\r\n\treturn 0;\r\n}\r\n<\/pre>\n","protected":false},"excerpt":{"rendered":"<p>\u3000\u7cbe\u5ea6\u4fdd\u8a3c\u8a08\u7b97\uff08\u3068\u3044\u3046\u8a00\u8449\u306f\u3042\u307e\u308a\u597d\u304d\u3067\u306f\u306a\u3044\u306e\u3060\u3051\u308c\u3069\uff09\u306e\u4e8b\u4f8b\u3082\u7d39\u4ecb\u3057\u3088\u3046\u304b\u3068\uff0cARB\uff08\u534a\u5f84\uff09\u3068MPFI\uff08\u533a\u9593\uff09\u3092\u4f7f\u3063\u3066\u30ed\u30b8\u30b9\u30c6\u30a3\u30c3\u30af\u5199\u50cf\u306e\u8a08\u7b97\u3092\u3055\u305b\u3066\u307f\u305f\u3068\u3053\u308d\uff0c\u6841\u6570\u304c\u8db3\u308a\u306a\u3044\u3068\u7834\u7dbb\u3059\u308b\u3053\u3068\u304c\u5206\u304b\u3063\u305f\uff08\u4eca\u66f4\uff09\u3002\u3064\u30fc\u3053 &hellip; <\/p>\n<p class=\"link-more\"><a href=\"https:\/\/na-inet.jp\/weblog2\/2017\/07\/13\/%e7%a0%b4%e7%b6%bb%e3%81%99%e3%82%8b%e5%8c%ba%e9%96%93%e6%bc%94%e7%ae%97%e3%81%ae%e4%ba%8b%e4%be%8b\/\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;\u7834\u7dbb\u3059\u308b\u533a\u9593\u6f14\u7b97\u306e\u4e8b\u4f8b&#8221;<\/span><\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_hide_featured_image":false,"_jetpack_featured_image_duplicate_dismissed":0,"_jetpack_feature_clip_id":0,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_publicize_message":"","jetpack_publicize_feature_enabled":true,"jetpack_social_post_already_shared":false,"jetpack_social_options":{"image_generator_settings":{"template":"highway","default_image_id":0,"font":"","enabled":false},"version":2},"jetpack_post_was_ever_published":false},"categories":[12,7],"tags":[],"class_list":["post-2971","post","type-post","status-publish","format-standard","hentry","category-mpna","category-7"],"jetpack_publicize_connections":[],"jetpack_sharing_enabled":true,"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/posts\/2971","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/comments?post=2971"}],"version-history":[{"count":5,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/posts\/2971\/revisions"}],"predecessor-version":[{"id":2977,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/posts\/2971\/revisions\/2977"}],"wp:attachment":[{"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/media?parent=2971"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/categories?post=2971"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/na-inet.jp\/weblog2\/wp-json\/wp\/v2\/tags?post=2971"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}