{"id":3038,"date":"2020-04-25T13:41:00","date_gmt":"2020-04-25T04:41:00","guid":{"rendered":"https:\/\/www.ktech.biz\/jp\/?p=3038"},"modified":"2020-07-30T17:07:54","modified_gmt":"2020-07-30T08:07:54","slug":"c","status":"publish","type":"post","link":"https:\/\/www.ktech.biz\/jp\/manual-api\/c\/","title":{"rendered":"XLPack API: C\/C++ \u4f7f\u7528\u4f8b"},"content":{"rendered":"

C\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(1), Visual C)<\/h3>\n

Visual C (Visual Studio (Windows\u7248)) \u3092\u4f7f\u7528\u3057\u305f\u4f8b\u3067\u3059. \u884c\u5217\u3092\u5217\u512a\u5148(VBA\u3084Fortran\u3068\u540c\u3058, \u901a\u5e38\u306eC\/C++\u306e\u683c\u7d0d\u9806\u756a\u3068\u306f\u884c\u3068\u5217\u304c\u53cd\u5bfe)\u3067\u683c\u7d0d\u3057\u3066\u3044\u308b\u3053\u3068\u306b\u6ce8\u610f\u3057\u3066\u304f\u3060\u3055\u3044.
\na\u3068b\u306b1\u6b21\u5143\u914d\u5217\u3092\u4f7f\u3044\u307e\u3059(\u30ea\u30d5\u30a1\u30ec\u30f3\u30b9\u30de\u30cb\u30e5\u30a2\u30eb\u306f\u4e0b\u306eClang\u306e\u4f8b\u306e\u3088\u3046\u306b\u53ef\u5909\u95772\u6b21\u5143\u914d\u5217\u3067\u8868\u8a18\u3057\u3066\u3044\u307e\u3059\u304cVisual C\u306e\u5834\u5408\u306b\u306f\u8aad\u307f\u66ff\u3048\u3066\u304f\u3060\u3055\u3044).<\/p>\n

\n
#include <stdio.h>\n#include \"cnumlib.h\"\n\nint main(void)\n{\n    double a[3][3] = {\n        { 0.2, -0.32, -0.8 },\n        { -0.11, 0.81, -0.92 },\n        { -0.93, 0.37, -0.29 } };\n    double b[] = { -0.3727, 0.4319, -1.4247 };\n    int n = 3, nrhs = 1, lda = 3, ldb = 3, ipiv[3], info;\n\n    dgesv(n, nrhs, lda, (double *)a, ipiv, ldb, b, &info);\n    printf(\"x = %f  %f  %f\\n\", b[0], b[1], b[2]);\n    printf(\"info = %d\\n\", info);\n    return 0;\n}<\/code><\/pre>\n<\/div>\n
\u5b9f\u884c\u7d50\u679c<\/p>\n
\n
x = 0.860000  0.640000  0.510000\ninfo = 0<\/code><\/pre>\n<\/div>\n
<\/code><\/p>\n
<\/p>\n
C\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(1), Clang)<\/h3>\n
Clang \u3067\u306fC99\u306e\u53ef\u5909\u9577\u914d\u5217\u304c\u4f7f\u3048\u308b\u306e\u3067a\u3068b\u306b\u306f2\u6b21\u5143\u914d\u5217\u3092\u4f7f\u3044\u307e\u3059.<\/p>\n
\n
#include <stdio.h>\n#include \"cnumlib.h\"\n\nint main(void)\n{\n    double a[3][3] = {\n        { 0.2, -0.32, -0.8 },\n        { -0.11, 0.81, -0.92 },\n        { -0.93, 0.37, -0.29 } };\n    double b[] = { -0.3727, 0.4319, -1.4247 };\n    int n = 3, nrhs = 1, lda = 3, ldb = 3, ipiv[3], info;\n\n    dgesv(n, nrhs, lda, a, ipiv, ldb, (double (*)[ldb])b, &info);\n    printf(\"x = %f  %f  %f\\n\", b[0], b[1], b[2]);\n    printf(\"info = %d\\n\", info);\n    return 0;\n}<\/code><\/pre>\n<\/div>\n
C\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(1), LAPACKE \u884c\u512a\u5148)<\/h3>\n
LAPACKE\u3067\u306f\u7b2c1\u30d1\u30e9\u30e1\u30fc\u30bf\u3067\u884c\u512a\u5148\u306e\u914d\u5217\u3092\u4f7f\u3046\u3053\u3068\u3092\u6307\u5b9a\u3067\u304d\u307e\u3059. \u3053\u308c\u306b\u3088\u308a\u884c\u5217\u306e\u8a18\u8ff0\u3092\u6570\u5b66\u8868\u8a18\u3068\u540c\u3058\u306b\u3059\u308b\u3053\u3068\u304c\u3067\u304d\u307e\u3059. \u5185\u90e8\u51e6\u7406\u3068\u3057\u3066\u306f, \u8a08\u7b97\u958b\u59cb\u524d\u306b\u914d\u5217\u306e\u8ee2\u7f6e\u3092\u884c\u3044, \u8a08\u7b97\u7d42\u4e86\u5f8c\u306b\u518d\u3073\u8ee2\u7f6e\u3092\u3057\u307e\u3059. ldb = 1 \u3068\u306a\u3063\u3066\u3044\u308b\u3053\u3068\u306b\u6ce8\u610f\u3057\u3066\u304f\u3060\u3055\u3044. \u3053\u308c\u306f, a\u3060\u3051\u3067\u306a\u304fb\u3082\u8ee2\u7f6e\u51e6\u7406\u3092\u884c\u3046\u4ed5\u69d8\u306b\u306a\u3063\u3066\u304a\u308a, b\u3068\u3057\u3066\u306f\u6a2a\u30d9\u30af\u30c8\u30eb\u306e\u5165\u529b\u3092\u671f\u5f85\u3055\u308c\u3066\u3044\u308b\u304b\u3089\u3067\u3059.<\/p>\n
\n
#include <stdio.h>\n#include \"lapacke.h\"\n\nint main(void)\n{\n    double a[3][3] = {\n        { 0.2, -0.11, -0.93 },\n        { -0.32, 0.81, 0.37 },\n        { -0.8, -0.92, -0.29 } };\n    double b[] = { -0.3727, 0.4319, -1.4247 };\n    int n = 3, nrhs = 1, lda = 3, ldb = 1, ipiv[3], info;\n\n    info = LAPACKE_dgesv(LAPACK_ROW_MAJOR, n, nrhs, (double *)a, lda, ipiv, b, ldb);\n    printf(\"x = %f  %f  %f\\n\", b[0], b[1], b[2]);\n    printf(\"info = %d\\n\", info);\n    return 0;\n}<\/code><\/pre>\n<\/div>\n
C++\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(1))<\/h3>\n
\u30d8\u30c3\u30c0\u30fc\u30d5\u30a1\u30a4\u30eb\u306f cnumlib (XLPack Lite\u3067\u306f cnumlib_lite) \u3092\u4f7f\u3046\u3053\u3068\u304c\u3067\u304d\u307e\u3059(info\u5f15\u6570\u306b\u6ce8\u610f). \u884c\u5217\u306fClang\u306e\u5834\u5408\u3067\u3082\u4e00\u6b21\u5143\u914d\u5217\u306b\u5165\u308c\u307e\u3059.<\/p>\n
\n
#include <iostream>\n#include \"cnumlib\"\nusing namespace std;\n\nint main(void)\n{\n    double a[3][3] = {\n        { 0.2, -0.32, -0.8 },\n        { -0.11, 0.81, -0.92 },\n        { -0.93, 0.37, -0.29 } };\n    double b[] = { -0.3727, 0.4319, -1.4247 };\n    int n = 3, nrhs = 1, lda = 3, ldb = 3, ipiv[3], info;\n\n    dgesv(n, nrhs, lda, (double *)a, ipiv, ldb, b, info);\n    cout << \"x = \" << b[0] << \", \" << b[1] << \", \" << b[2] << endl;\n    cout << \"info = \" << info << endl;\n    return 0;\n}<\/code><\/pre>\n<\/div>\n
C\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(2))<\/h3>\n
\u4f8b\u984c(2)\u3067\u306ff(x)\u306e\u7a4d\u5206\u3092qk15\u3092\u4f7f\u3063\u3066\u6c42\u3081\u307e\u3059. qk15\u306ff(x)\u3092\u5b9a\u7fa9\u3059\u308b\u5916\u90e8\u95a2\u6570\u3092\u5fc5\u8981\u3068\u3057\u307e\u3059. \u5916\u90e8\u95a2\u6570\u306fC\u3067\u4f5c\u6210\u3059\u308b\u3053\u3068\u304c\u3067\u304d, qk15\u306f\u5fc5\u8981\u306a\u3068\u304d\u306b\u305d\u308c\u3092\u547c\u3073\u51fa\u3057\u307e\u3059.<\/p>\n
\n
#include <stdio.h>\n#include \"cnumlib.h\"\n\ndouble f(double x)\n{\n    return 1\/(1 + x*x);\n}\n\nvoid test_qk15()\n{\n    double a, b, result, abserr, resabs, resasc;\n\n    a = 0; b = 4;\n    qk15(f, a, b, &result, &abserr, &resabs, &resasc);\n    printf(\"result = %g, abserr = %g\\n\", result, abserr);\n}<\/code><\/pre>\n<\/div>\n
\u5b9f\u884c\u7d50\u679c<\/p>\n
\n
result = 1.32582, abserr = 0.00148272<\/code><\/pre>\n<\/div>\n
<\/code><\/code><\/p>\n","protected":false},"excerpt":{"rendered":"
C\u30d7\u30ed\u30b0\u30e9\u30e0\u4f8b (\u4f8b\u984c(1), Visual C) Visual C (Visual Studio (Windows\u7248)) \u3092\u4f7f\u7528\u3057\u305f\u4f8b\u3067\u3059. \u884c\u5217\u3092\u5217\u512a\u5148(VBA\u3084Fortran\u3068\u540c\u3058, \u901a\u5e38\u306eC\/C++\u306e\u683c\u7d0d\u9806\u756a\u3068\u306f […]<\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"footnotes":""},"categories":[16],"tags":[],"class_list":["post-3038","post","type-post","status-publish","format-standard","hentry","category-manual-api"],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/posts\/3038","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/comments?post=3038"}],"version-history":[{"count":2,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/posts\/3038\/revisions"}],"predecessor-version":[{"id":3093,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/posts\/3038\/revisions\/3093"}],"wp:attachment":[{"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/media?parent=3038"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/categories?post=3038"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.ktech.biz\/jp\/wp-json\/wp\/v2\/tags?post=3038"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}