{"id":131,"date":"2014-07-02T11:04:55","date_gmt":"2014-07-02T02:04:55","guid":{"rendered":"http:\/\/cs-tklab.na-inet.jp\/?p=131"},"modified":"2015-05-21T16:17:09","modified_gmt":"2015-05-21T07:17:09","slug":"mathjax%e3%81%a8syntax-highlighter%e3%81%ae%e3%83%86%e3%82%b9%e3%83%88","status":"publish","type":"post","link":"https:\/\/kd-tklab.na-inet.jp\/?p=131","title":{"rendered":"MathJax\u3068Syntax Highlighter\u306e\u30c6\u30b9\u30c8"},"content":{"rendered":"

\u3000\u3068\u308a\u3042\u3048\u305a\uff0cMathJax\u306e\u30c6\u30b9\u30c8\u3002\u30c7\u30d5\u30a9\u30eb\u30c8\u3092Display mode\u306b\u3057\u3066
\n\\[ y(s) = f(s) + \\int^b_a g(s, t) y(t) dt \\]
\n\u6570\u5f0f\u304c\u51fa\u308b\u304b\u3057\u3089\uff1f\u3000\u3064\u3044\u3067\u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a\u3082\"y'(t)\u30fb\u30fb\u30fb\u307e\u3041\u3053\u3093\u306a\u3082\u3093\u304b\u3002<\/p>\n

\u3000\u3067\uff0c\u30b3\u30fc\u30c9\u4f8b\u3082\u4e00\u3064\u3002$n$\u6b21\u6b63\u65b9\u884c\u5217\"A\in\mathbb{R}^{n\times\u3068\"n\"\u6b21\u5143\u30d9\u30af\u30c8\u30eb\"\mathbf{x},\u304b\u3089\u4f5c\u3089\u308c\u308b\u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f
\n\\[ A\\mathbf{x} = \\mathbf{b} \\]
\n\u3092\u672a\u77e5\u30d9\u30af\u30c8\u30eb\"\mathbf{x}\"\u306b\u3064\u3044\u3066\u89e3\u304fLAPACK\u30d7\u30ed\u30b0\u30e9\u30e0\u3002DGESV\u3092\u4f7f\u3046\u6700\u3082\u7c21\u5358\u306a\u3082\u306e\u3002\u89e3\u306e\u8868\u793a\u3092\u30b3\u30e1\u30f3\u30c8\u306b\u3057\u3066\u3042\u308b\u306e\u306f\u3067\u304b\u3044\u6b21\u5143\u6570\u306e\u30d9\u30f3\u30c1\u30de\u30fc\u30af\u7528\u306b\u4f7f\u3063\u305f\u304b\u3089\u3067\u3059\u3002\u9069\u5b9c\u306f\u305a\u3057\u3066\u304a\u4f7f\u3044\u4e0b\u3055\u3044\u307e\u305b\u3002<\/p>\n

\/*************************************************\/\r\n\/* LAPACK Tutorial                               *\/\r\n\/* Sample Program of DGESV                       *\/\r\n\/*                        with Intel Math Kernel *\/\r\n\/* Last Update: 2013-01-21 (Mon) T.Kouya         *\/\r\n\/\/ icc -I\/opt\/intel\/include -I\/opt\/intel\/mkl\/include -L\/opt\/intel\/mkl\/lib\/intel64 -L\/opt\/intel\/lib\/intel64 blas1_mkl.c -Wl,--start-group \/opt\/intel\/mkl\/lib\/intel64\/libmkl_intel_lp64.a \/opt\/intel\/mkl\/lib\/intel64\/libmkl_sequential.a \/opt\/intel\/mkl\/lib\/intel64\/libmkl_core.a -Wl,--end-group -lpthread\r\n\/*************************************************\/\r\n#include <stdio.h>\r\n#include <stdlib.h>\r\n#include "mkl.h" \/\/ for Intel Math Kernel Library\r\n#include "mkl_cblas.h"\r\n#include "mkl_lapacke.h"\r\n\r\n#define min(i, j) (((i) < (j)) ? (i) : (j))\r\n#define max(i, j) (((i) > (j)) ? (i) : (j))\r\n\r\nint main()\r\n{\r\n\tMKL_INT i, j, dim; \/\/ dimension of vectors\r\n\tMKL_INT *pivot, info;\r\n\tdouble *ma, *vx, *vb;\r\n\tdouble alpha, beta;\r\n\tchar str_mkl_version[1024];\r\n\tdouble running_time;\r\n\r\n\t\/\/ print MKL version\r\n\tMKL_Get_Version_String(str_mkl_version, 1024);\r\n\tprintf("%s\\n", str_mkl_version);\r\n\r\n\t\/\/ input dimension\r\n\tprintf("Dim = "); scanf("%d", &dim);\r\n\r\n\tif(dim <= 0)\r\n\t{\r\n\t\tprintf("Illigal dimension! (dim = %d)\\n", dim);\r\n\t\treturn EXIT_FAILURE;\r\n\t}\r\n\r\n\t\/\/ Initialize\r\n\tpivot = (MKL_INT *)calloc(sizeof(MKL_INT), dim);\r\n\tma = (double *)calloc(sizeof(double), dim * dim);\r\n\tvx = (double *)calloc(sizeof(double), dim);\r\n\tvb = (double *)calloc(sizeof(double), dim);\r\n\r\n\t\/\/ input va and vb\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tfor(j = 0; j < dim; j++)\r\n\t\t{\r\n\t\t\t\/\/ A = hirbert matrix\r\n\t\t\t\/\/*(ma + i * dim + j) = 1.0 \/ (double)(i + j + 1);\r\n\t\t\t\/\/ A = frank matrix\r\n\t\t\t*(ma + i * dim + j) = (double)(dim - max(i, j));\r\n\t\t}\r\n\t\t\/\/ x  = answer\r\n\t\t*(vx + i) = (double) (i + 1);\r\n\t}\r\n\r\n\t\/\/ b := A * x\r\n\talpha = 1.0;\r\n\tbeta = 0.0;\r\n\tcblas_dgemv(CblasRowMajor, CblasNoTrans, dim, dim, alpha, ma, dim, vx, 1, beta, vb, 1);\r\n\r\n\t\/\/ print\r\n\/*\tprintf("A = \\n");\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tprintf("%3d: ", i);\r\n\t\tfor(j = 0; j < dim; j++)\r\n\t\t\tprintf("%10f ", *(ma + i * dim + j));\r\n\t\tprintf("\\n");\r\n\r\n\t}\r\n\r\n\t\/\/ print\r\n\tprintf("x = \\n");\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tprintf("%3d: ", i);\r\n\t\tprintf("%10f ", *(vx + i));\r\n\t\tprintf("\\n");\r\n\r\n\t}\r\n\r\n\t\/\/ print\r\n\tprintf("b = \\n");\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tprintf("%3d: ", i);\r\n\t\tprintf("%10f ", *(vb + i));\r\n\t\tprintf("\\n");\r\n\t}\r\n*\/\r\n\r\n\t\/\/ solve A * X = C -> C := X\r\n\trunning_time = dsecnd(); \/\/ start \r\n\tinfo = LAPACKE_dgesv(LAPACK_ROW_MAJOR, dim, 1, ma, dim, pivot, vb, 1);\r\n\trunning_time = dsecnd() - running_time; \/\/ end\r\n\r\n\t\/\/ error occurs if info > 0\r\n\tif(info > 0)\r\n\t{\r\n\t\tprintf("pivot ma[%d][%d] is zero!\\n", info, info);\r\n\t\treturn EXIT_FAILURE;\r\n\t}\r\n\r\n\tprintf("Dimension: %d, Running Time(s): %g\\n", dim, running_time);\r\n\r\n\t\/\/ print\r\n\/*\tprintf("calculated x = \\n");\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tprintf("%3d: ", i);\r\n\t\tprintf("%10f ", *(vb + i));\r\n\t\tprintf("\\n");\r\n\r\n\t}\r\n\r\n\t\/\/ diff\r\n\tprintf("x - calculated x = \\n");\r\n\tfor(i = 0; i < dim; i++)\r\n\t{\r\n\t\tprintf("%3d: ", i);\r\n\t\tprintf("%10f ", fabs(*(vx + i) - *(vb + i)));\r\n\t\tprintf("\\n");\r\n\t}\r\n*\/\r\n\r\n\t\/\/ free\r\n\tfree(pivot);\r\n\tfree(ma);\r\n\tfree(vb);\r\n\tfree(vx);\r\n\r\n\treturn EXIT_SUCCESS;\r\n}<\/pre>\n","protected":false},"excerpt":{"rendered":"
\u3000\u3068\u308a\u3042\u3048\u305a\uff0cMathJax\u306e\u30c6\u30b9\u30c8\u3002\u30c7\u30d5\u30a9\u30eb\u30c8\u3092Display mode\u306b\u3057\u3066 \\[ y(s) = f(s) + \\int^b_a g(s, t) y(t) dt \\] \u6570\u5f0f\u304c\u51fa\u308b\u304b\u3057\u3089\uff1f\u3000\u3064\u3044\u3067\u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a\u3082\u30fb\u30fb … <\/p>\n
“MathJax\u3068Syntax Highlighter\u306e\u30c6\u30b9\u30c8” \u306e<\/span>\u7d9a\u304d\u3092\u8aad\u3080<\/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_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":[2],"tags":[],"class_list":["post-131","post","type-post","status-publish","format-standard","hentry","category-2"],"jetpack_publicize_connections":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"jetpack_shortlink":"https:\/\/wp.me\/p6Oxer-27","jetpack-related-posts":[{"id":2493,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=2493","url_meta":{"origin":131,"position":0},"title":"\u9ad8\u6027\u80fd\u8a08\u7b97\u30d7\u30ed\u30b0\u30e9\u30df\u30f3\u30b0","author":"T.Kouya","date":"2023\u5e747\u670811\u65e5","format":false,"excerpt":"\u5e78\u8c37\u667a\u7d00 \u7b2c\u4e00\u90e8\uff1aPython\u306b\u3088\u308b\u65b9\u7a0b\u5f0f\u6c42\u89e3\u30d7\u30ed\u30b0\u30e9\u30df\u30f3\u30b0 \u300cPython\u6570\u5024\u8a08\u7b97\u30d7\u30ed\u30b0\u30e9\u30df\u30f3\u30b0\u300d\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":1769,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=1769","url_meta":{"origin":131,"position":1},"title":"PC\u30b1\u30fc\u30b9\u30d5\u30a1\u30f3\u306e\u4ea4\u63db","author":"T.Kouya","date":"2018\u5e744\u670818\u65e5","format":false,"excerpt":"\u3000\u672c\u7814\u7a76\u5ba4\u306e\u30c7\u30b9\u30af\u30c8\u30c3\u30d7\u30de\u30b7\u30f3\u306f\u5927\u6982\u53e4\u307c\u3051\u305fPC\u30b1\u30fc\u30b9\u306b\u53ce\u307e\u3063\u3066\u3044\u307e\u3059\u3002\u65e2\u306b\uff13\u4e16\u4ee3\u4ee5\u4e0a\u66f4\u65b0\u3092\u91cd\u306d\u3066\u3044\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"","src":"https:\/\/i0.wp.com\/cs-tklab.na-inet.jp\/blog\/wp-content\/uploads\/2018\/04\/IMG_4041-300x241.jpg?resize=350%2C200&ssl=1","width":350,"height":200},"classes":[]},{"id":328,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=328","url_meta":{"origin":131,"position":2},"title":"\u6a21\u64ec\u8b1b\u7fa9\u300c\u30b3\u30f3\u30d4\u30e5\u30fc\u30bf\u306f\u3069\u3046\u3084\u3063\u3066\u8a08\u7b97\u3057\u3066\u3044\u308b\u306e\u304b\uff1f\u300d\uff08\u6539\u984c\u300c\u6570\u5b66\u3068\u8a08\u7b97\u3000\uff5e\u7406\u8ad6\u3068\u5b9f\u8df5\uff5e\u300d\uff09","author":"T.Kouya","date":"2016\u5e749\u67089\u65e5","format":false,"excerpt":"\u6982\u8981 \u3000\u6570\u5b66\uff08\uff1d\u7406\u8ad6\uff09\u3068\u8a08\u7b97\uff08\uff1d\u5b9f\u8df5\uff09\u306f\u5168\u304f\u6027\u8cea\u306e\u7570\u306a\u308b\u3082\u306e\u3067\uff0c\u8a08\u7b97\u306e\u52b9\u7387\u6027\u3084\u6b63\u78ba\u6027\u3092\u8ffd\u53ca\u3059\u308b\u3068\u6570\u5b66\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"basis_computing1","src":"https:\/\/i0.wp.com\/cs-tklab.na-inet.jp\/blog\/wp-content\/uploads\/2015\/05\/basis_computing1-300x195.png?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":586,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=586","url_meta":{"origin":131,"position":3},"title":"[\u6a21\u64ec\u8b1b\u7fa92] \u5927\u898f\u6a21\u8a08\u7b97\u3068CPU\uff06GPU","author":"T.Kouya","date":"2015\u5e749\u670829\u65e5","format":false,"excerpt":"\u3000\u3053\u308c\u306f\u6a21\u64ec\u8b1b\u7fa9\u300c\u6570\u5b66\u3068\u8a08\u7b97\u3000\uff5e\u7406\u8ad6\u3068\u5b9f\u8df5\uff5e\u300d\u306e\u7d9a\u304d\u306b\u306a\u308a\u307e\u3059\u3002\u307e\u305a\u305d\u3061\u3089\u3092\u8aad\u3093\u3067\u9802\u3044\u3066\u304b\u3089\uff0c\u3053\u3061\u3089\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"lu","src":"https:\/\/i0.wp.com\/cs-tklab.na-inet.jp\/blog\/wp-content\/uploads\/2015\/09\/lu-98x300.png?resize=350%2C200","width":350,"height":200},"classes":[]},{"id":2146,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=2146","url_meta":{"origin":131,"position":4},"title":"\u300c\u9023\u7acb\u4e00\u6b21\u65b9\u7a0b\u5f0f\u3092\u30ac\u30a6\u30b9\u306e\u6d88\u53bb\u6cd5\u3067\u89e3\u3051\u308b\u30a2\u30d7\u30ea\u300d\u30b3\u30e1\u30f3\u30c8","author":"T.Kouya","date":"2019\u5e746\u67087\u65e5","format":false,"excerpt":"\u3000\u4eca\u65e5\u306e\u300c\u5fdc\u7528\u7dda\u5f62\u4ee3\u6570\u300d\u3067\u30b3\u30e1\u30f3\u30c8\u3092\u6c42\u3081\u3089\u308c\u305f\u306e\u3067\uff0c\u3061\u3087\u3063\u3068\u89e6\u3063\u3066\u307f\u307e\u3057\u305f\u3002URL\u306f\u516c\u958b\u3057\u3066\u3044\u3044\u306e\u304b\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"","src":"","width":0,"height":0},"classes":[]},{"id":266,"url":"https:\/\/kd-tklab.na-inet.jp\/?p=266","url_meta":{"origin":131,"position":5},"title":"\u5927\u5b66\u898b\u5b66\u4f1a\u300c\u9ad8\u6027\u80fd\u8a08\u7b97\u7814\u7a76\u5ba4\u300d","author":"T.Kouya","date":"2017\u5e743\u670822\u65e5","format":false,"excerpt":"\u9ad8\u6027\u80fd\u8a08\u7b97\u3068Web\u30a2\u30d7\u30ea\u30b1\u30fc\u30b7\u30e7\u30f3 \u9759\u5ca1\u7406\u5de5\u79d1\u5927\u5b66\u3000\u30b3\u30f3\u30d4\u30e5\u30fc\u30bf\u30b7\u30b9\u30c6\u30e0\u5b66\u79d1\u3000\u5e78\u8c37\u3053\u3046\u3084\u3000\u667a\u7d00\u3068\u3082\u306e\u2026","rel":"","context":"\u30a8\u30c3\u30bb\u30a4","block_context":{"text":"\u30a8\u30c3\u30bb\u30a4","link":"https:\/\/kd-tklab.na-inet.jp\/?cat=7"},"img":{"alt_text":"try_mpfr_pi_3321","src":"https:\/\/i0.wp.com\/cs-tklab.na-inet.jp\/blog\/wp-content\/uploads\/2015\/05\/try_mpfr_pi_3321-300x129.png?resize=350%2C200","width":350,"height":200},"classes":[]}],"amp_enabled":true,"_links":{"self":[{"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/posts\/131","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=131"}],"version-history":[{"count":18,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/posts\/131\/revisions"}],"predecessor-version":[{"id":170,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=\/wp\/v2\/posts\/131\/revisions\/170"}],"wp:attachment":[{"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=131"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=131"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/kd-tklab.na-inet.jp\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=131"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}