cpprefjp
diff --git a/‎reference/linalg/hermitian_matrix_vector_product.html
Lines changed: 18 additions & 17 deletions b/‎reference/linalg/hermitian_matrix_vector_product.html
Lines changed: 18 additions & 17 deletions
@@ -84,7 +84,7 @@
 
     </head>
     <body>
-<header data-kunai-mdinfo="{&#34;meta&#34;: {&#34;header&#34;: [&#34;linalg&#34;], &#34;id-type&#34;: [&#34;function template&#34;], &#34;namespace&#34;: [&#34;std::linalg&#34;], &#34;cpp&#34;: [&#34;cpp26&#34;]}, &#34;sources&#34;: [{&#34;id&#34;: &#34;57af623807ea7b4e9b34bde27a1f2ccd9add2517&#34;, &#34;source&#34;: &#34;#include &lt;array&gt;\n#include &lt;complex&gt;\n#include &lt;iostream&gt;\n#include &lt;linalg&gt;\n#include &lt;mdspan&gt;\n#include &lt;vector&gt;\n\ntemplate &lt;class Vector&gt;\nvoid print(const Vector&amp; v, const std::string&amp; name) {\n  for (int i = 0; i &lt; v.extent(0); ++i) {\n    std::cout &lt;&lt; name &lt;&lt; \&#34;[\&#34; &lt;&lt; i &lt;&lt; \&#34;]\&#34; &lt;&lt; \&#34; = \&#34; &lt;&lt; v[i] &lt;&lt; &#39;\\n&#39;;\n  }\n}\n\nint main()\n{\n  constexpr size_t N = 4;\n  constexpr size_t M = 4;\n\n  std::vector&lt;std::complex&lt;double&gt;&gt; A_vec(N*M);\n  std::vector&lt;double&gt; x_vec(M);\n  std::array&lt;double, N&gt; y_vec, z_vec;\n\n  std::mdspan&lt;\n    std::complex&lt;double&gt;,\n    std::extents&lt;size_t, N, M&gt;,\n    std::linalg::layout_blas_packed&lt;\n      std::linalg::upper_triangle_t,\n      std::linalg::row_major_t&gt;\n  &gt; A(A_vec.data());\n  std::mdspan x(x_vec.data(), M);\n  std::mdspan y(y_vec.data(), N);\n  std::mdspan z(z_vec.data(), N);\n\n  for(int i = 0; i &lt; A.extent(0); ++i) {\n    for(int j = i; j &lt; A.extent(1); ++j) {\n      A[i,j] = std::complex&lt;double&gt;(i, j);\n    }\n  }\n\n  for(int j = 0; j &lt; x.extent(0); ++j) {\n    x[j] = j;\n  }\n\n  // (1)\n  std::cout &lt;&lt; \&#34;(1)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(A, std::linalg::upper_triangle, x, y);\n  print(y, \&#34;y\&#34;);\n\n  // (2)\n  std::cout &lt;&lt; \&#34;(2)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(std::execution::par, A, std::linalg::upper_triangle, x, y);\n  print(y, \&#34;y\&#34;);\n\n  for(int i = 0; i &lt; y.extent(0); ++i) {\n    y[i] = -i;\n  }\n\n  // (3)\n  std::cout &lt;&lt; \&#34;(3)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(A, std::linalg::upper_triangle, x, y, z);\n  print(z, \&#34;z\&#34;);\n\n  // (4)\n  std::cout &lt;&lt; \&#34;(4)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(std::execution::par, A, std::linalg::upper_triangle, x, y, z);\n  print(z, \&#34;z\&#34;);\n\n  return 0;\n}\n&#34;}], &#34;page_id&#34;: [&#34;reference&#34;, &#34;linalg&#34;, &#34;hermitian_matrix_vector_product&#34;]}">
+<header data-kunai-mdinfo="{&#34;meta&#34;: {&#34;header&#34;: [&#34;linalg&#34;], &#34;id-type&#34;: [&#34;function template&#34;], &#34;namespace&#34;: [&#34;std::linalg&#34;], &#34;cpp&#34;: [&#34;cpp26&#34;]}, &#34;sources&#34;: [{&#34;id&#34;: &#34;5d6dc4179e7e08da018b9a61c9f163db70984b26&#34;, &#34;source&#34;: &#34;#include &lt;array&gt;\n#include &lt;complex&gt;\n#include &lt;iostream&gt;\n#include &lt;linalg&gt;\n#include &lt;mdspan&gt;\n#include &lt;vector&gt;\n\ntemplate &lt;class Vector&gt;\nvoid print(const Vector&amp; v, const std::string&amp; name) {\n  for (int i = 0; i &lt; v.extent(0); ++i) {\n    std::cout &lt;&lt; name &lt;&lt; \&#34;[\&#34; &lt;&lt; i &lt;&lt; \&#34;]\&#34; &lt;&lt; \&#34; = \&#34; &lt;&lt; v[i] &lt;&lt; &#39;\\n&#39;;\n  }\n}\n\nint main()\n{\n  constexpr size_t N = 4;\n  constexpr size_t M = 4;\n\n  std::vector&lt;std::complex&lt;double&gt;&gt; A_vec(N*M);\n  std::vector&lt;double&gt; x_vec(M);\n  std::array&lt;double, N&gt; y_vec, z_vec;\n\n  std::mdspan&lt;\n    std::complex&lt;double&gt;,\n    std::extents&lt;size_t, N, M&gt;,\n    std::linalg::layout_blas_packed&lt;\n      std::linalg::upper_triangle_t,\n      std::linalg::row_major_t&gt;\n  &gt; A(A_vec.data());\n  std::mdspan x(x_vec.data(), M);\n  std::mdspan y(y_vec.data(), N);\n  std::mdspan z(z_vec.data(), N);\n\n  for(int i = 0; i &lt; A.extent(0); ++i) {\n    A[i,i] = std::complex&lt;double&gt;(i, 0);\n    for(int j = i + 1; j &lt; A.extent(1); ++j) {\n      A[i,j] = std::complex&lt;double&gt;(i, j);\n    }\n  }\n\n  for(int j = 0; j &lt; x.extent(0); ++j) {\n    x[j] = j;\n  }\n\n  // (1)\n  std::cout &lt;&lt; \&#34;(1)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(A, std::linalg::upper_triangle, x, y);\n  print(y, \&#34;y\&#34;);\n\n  // (2)\n  std::cout &lt;&lt; \&#34;(2)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(std::execution::par, A, std::linalg::upper_triangle, x, y);\n  print(y, \&#34;y\&#34;);\n\n  for(int i = 0; i &lt; y.extent(0); ++i) {\n    y[i] = -i;\n  }\n\n  // (3)\n  std::cout &lt;&lt; \&#34;(3)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(A, std::linalg::upper_triangle, x, y, z);\n  print(z, \&#34;z\&#34;);\n\n  // (4)\n  std::cout &lt;&lt; \&#34;(4)\\n\&#34;;\n  std::linalg::hermitian_matrix_vector_product(std::execution::par, A, std::linalg::upper_triangle, x, y, z);\n  print(z, \&#34;z\&#34;);\n\n  return 0;\n}\n&#34;}], &#34;page_id&#34;: [&#34;reference&#34;, &#34;linalg&#34;, &#34;hermitian_matrix_vector_product&#34;]}">
   <nav class="navbar navbar-default" role="navigation">
     <div class="container-fluid">
       <div class="navbar-header">
@@ -204,8 +204,8 @@
 
     <p class="text-right"><small>
       最終更新日時(UTC):
-      <span itemprop="datePublished" content="2024-06-26T02:45:26">
-        2024年06月26日 02時45分26秒
+      <span itemprop="datePublished" content="2024-06-26T03:41:50">
+        2024年06月26日 03時41分50秒
       </span>
       <br/>
       <span itemprop="author" itemscope itemtype="http://schema.org/Person">
@@ -316,7 +316,7 @@ <h2>備考</h2>
 </ul>
 <h2>例</h2>
 <p><strong>[注意] 処理系にあるコンパイラで確認していないため、間違っているかもしれません。</strong></p>
-<p><div class="yata" id="57af623807ea7b4e9b34bde27a1f2ccd9add2517"><div class="codehilite"><pre><span></span><code><span class="cp">#include</span> <span class="cpf"><a href="../array.html">&lt;array&gt;</a></span><span class="cp"></span>
+<p><div class="yata" id="5d6dc4179e7e08da018b9a61c9f163db70984b26"><div class="codehilite"><pre><span></span><code><span class="cp">#include</span> <span class="cpf"><a href="../array.html">&lt;array&gt;</a></span><span class="cp"></span>
 <span class="cp">#include</span> <span class="cpf"><a href="../complex.html">&lt;complex&gt;</a></span><span class="cp"></span>
 <span class="cp">#include</span> <span class="cpf"><a href="../iostream.html">&lt;iostream&gt;</a></span><span class="cp"></span>
 <span class="cp">#include</span> <span class="cpf"><a href="../linalg.html">&lt;linalg&gt;</a></span><span class="cp"></span>
@@ -351,7 +351,8 @@ <h2>例</h2>
   <span class="n"><a href="../mdspan/mdspan.html">std::mdspan</a></span> <span class="nf">z</span><span class="p">(</span><span class="n">z_vec</span><span class="p">.</span><span class="n">data</span><span class="p">(),</span> <span class="n">N</span><span class="p">);</span>
 
   <span class="k">for</span><span class="p">(</span><span class="kt">int</span> <span class="n">i</span> <span class="o">=</span> <span class="mi">0</span><span class="p">;</span> <span class="n">i</span> <span class="o">&lt;</span> <span class="n">A</span><span class="p">.</span><span class="n">extent</span><span class="p">(</span><span class="mi">0</span><span class="p">);</span> <span class="o">++</span><span class="n">i</span><span class="p">)</span> <span class="p">{</span>
-    <span class="k">for</span><span class="p">(</span><span class="kt">int</span> <span class="n">j</span> <span class="o">=</span> <span class="n">i</span><span class="p">;</span> <span class="n">j</span> <span class="o">&lt;</span> <span class="n">A</span><span class="p">.</span><span class="n">extent</span><span class="p">(</span><span class="mi">1</span><span class="p">);</span> <span class="o">++</span><span class="n">j</span><span class="p">)</span> <span class="p">{</span>
+    <span class="n">A</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">i</span><span class="p">]</span> <span class="o">=</span> <span class="n">std</span><span class="o">::</span><span class="n">complex</span><span class="o">&lt;</span><span class="kt">double</span><span class="o">&gt;</span><span class="p">(</span><span class="n">i</span><span class="p">,</span> <span class="mi">0</span><span class="p">);</span>
+    <span class="k">for</span><span class="p">(</span><span class="kt">int</span> <span class="n">j</span> <span class="o">=</span> <span class="n">i</span> <span class="o">+</span> <span class="mi">1</span><span class="p">;</span> <span class="n">j</span> <span class="o">&lt;</span> <span class="n">A</span><span class="p">.</span><span class="n">extent</span><span class="p">(</span><span class="mi">1</span><span class="p">);</span> <span class="o">++</span><span class="n">j</span><span class="p">)</span> <span class="p">{</span>
       <span class="n">A</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">std</span><span class="o">::</span><span class="n">complex</span><span class="o">&lt;</span><span class="kt">double</span><span class="o">&gt;</span><span class="p">(</span><span class="n">i</span><span class="p">,</span> <span class="n">j</span><span class="p">);</span>
     <span class="p">}</span>
   <span class="p">}</span>
@@ -391,24 +392,24 @@ <h2>例</h2>
 <h3>出力</h3>
 <p><pre><code>(1)
 y[0] = (0,14)
-y[1] = (6,14)
-y[2] = (11,11)
-y[3] = (14,0)
+y[1] = (6,13)
+y[2] = (11,7)
+y[3] = (14,-9)
 (2)
 y[0] = (0,14)
-y[1] = (6,14)
-y[2] = (11,11)
-y[3] = (14,0)
+y[1] = (6,13)
+y[2] = (11,7)
+y[3] = (14,-9)
 (3)
 z[0] = (0,14)
-z[1] = (5,14)
-z[2] = (9,11)
-z[3] = (11,0)
+z[1] = (5,13)
+z[2] = (9,7)
+z[3] = (11,-9)
 (4)
 z[0] = (0,14)
-z[1] = (5,14)
-z[2] = (9,11)
-z[3] = (11,0)
+z[1] = (5,13)
+z[2] = (9,7)
+z[3] = (11,-9)
 </code></pre></p>
 <h2>バージョン</h2>
 <h3>言語</h3>