Skip to content

Commit 0850077

Browse files
author
cpprefjp-autoupdate
committed
update automatically
1 parent c3be1d9 commit 0850077

File tree

8 files changed

+721
-423
lines changed

8 files changed

+721
-423
lines changed

‎reference/linalg.html

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -176,8 +176,8 @@
176176

177177
<p class="text-right"><small>
178178
最終更新日時(UTC):
179-
<span itemprop="datePublished" content="2024-07-05T08:57:02">
180-
2024年07月05日 08時57分02秒
179+
<span itemprop="datePublished" content="2024-07-06T05:57:34">
180+
2024年07月06日 05時57分34秒
181181
</span>
182182
<br/>
183183
<span itemprop="author" itemscope itemtype="http://schema.org/Person">
@@ -458,7 +458,7 @@ <h2>BLAS 2アルゴリズム</h2>
458458
<td>C++26</td>
459459
</tr>
460460
<tr>
461-
<td><code>hermitian_matrix_rank_2_update</code></td>
461+
<td><code><a href="linalg/hermitian_matrix_rank_2_update.html">hermitian_matrix_rank_2_update</a></code></td>
462462
<td>xHER2: ハミルトニアン行列のRank-2更新 (function template)</td>
463463
<td>C++26</td>
464464
</tr>

‎reference/linalg/hermitian_matrix_rank_1_update.html

Lines changed: 5 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -84,7 +84,7 @@
8484

8585
</head>
8686
<body>
87-
<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;a2a66a0d88359edc033efe7cd3c98f8f6e38cae0&#34;, &#34;source&#34;: &#34;#include &lt;array&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 Matrix&gt;\nvoid print_mat(const Matrix&amp; A) {\n for(int i = 0; i &lt; A.extent(0); ++i) {\n for(int j = 0; j &lt; i; ++j) {\n std::cout &lt;&lt; A[j, i] &lt;&lt; &#39; &#39;;\n }\n for(int j = i; j &lt; A.extent(1) - 1; ++j) {\n std::cout &lt;&lt; A[i, j] &lt;&lt; &#39; &#39;;\n }\n std::cout &lt;&lt; A[i, A.extent(1) - 1] &lt;&lt; &#39;\\n&#39;;\n }\n}\n\ntemplate &lt;class Vector&gt;\nvoid init_vec(Vector&amp; v) {\n for (int i = 0; i &lt; v.extent(0); ++i) {\n v[i] = std::complex&lt;double&gt;(0, i);\n }\n}\n\ntemplate &lt;class Matrix&gt;\nvoid init_mat(Matrix&amp; A) {\n for(int i = 0; i &lt; A.extent(0); ++i) {\n A[i,j] = 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\nint main()\n{\n constexpr size_t N = 4;\n\n std::vector&lt;std::complex&lt;double&gt;&gt; A_vec(N * N);\n std::vector&lt;std::complex&lt;double&gt;&gt; x_vec(N);\n std::array&lt;std::complex&lt;double&gt;, N&gt; y_vec;\n\n std::mdspan&lt;\n std::complex&lt;double&gt;,\n std::extents&lt;size_t, N, N&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(), N);\n std::mdspan y(y_vec.data(), N);\n\n init_mat(A);\n init_vec(x);\n init_vec(y);\n\n // (1)\n std::cout &lt;&lt; \&#34;(1)\\n\&#34;;\n std::linalg::matrix_rank_1_update(\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (2)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(2)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n std::execution::par,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (3)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(3)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n -1.0,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (4)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(4)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n std::execution::par,\n -1.0,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n return 0;\n}\n&#34;}], &#34;page_id&#34;: [&#34;reference&#34;, &#34;linalg&#34;, &#34;hermitian_matrix_rank_1_update&#34;]}">
87+
<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;1792939fb57232efc6310058b9ce3580f17d6fa9&#34;, &#34;source&#34;: &#34;#include &lt;array&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 Matrix&gt;\nvoid print_mat(const Matrix&amp; A) {\n for(int i = 0; i &lt; A.extent(0); ++i) {\n for(int j = 0; j &lt; i; ++j) {\n std::cout &lt;&lt; A[j, i] &lt;&lt; &#39; &#39;;\n }\n for(int j = i; j &lt; A.extent(1) - 1; ++j) {\n std::cout &lt;&lt; A[i, j] &lt;&lt; &#39; &#39;;\n }\n std::cout &lt;&lt; A[i, A.extent(1) - 1] &lt;&lt; &#39;\\n&#39;;\n }\n}\n\ntemplate &lt;class Vector&gt;\nvoid init_vec(Vector&amp; v) {\n for (int i = 0; i &lt; v.extent(0); ++i) {\n v[i] = std::complex&lt;double&gt;(0, i);\n }\n}\n\ntemplate &lt;class Matrix&gt;\nvoid init_mat(Matrix&amp; A) {\n for(int i = 0; i &lt; A.extent(0); ++i) {\n A[i,j] = 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\nint main()\n{\n constexpr size_t N = 4;\n\n std::vector&lt;std::complex&lt;double&gt;&gt; A_vec(N * N);\n std::vector&lt;std::complex&lt;double&gt;&gt; x_vec(N);\n\n std::mdspan&lt;\n std::complex&lt;double&gt;,\n std::extents&lt;size_t, N, N&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(), N);\n\n init_mat(A);\n init_vec(x);\n\n // (1)\n std::cout &lt;&lt; \&#34;(1)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (2)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(2)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n std::execution::par,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (3)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(3)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n -1.0,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n // (4)\n init_mat(A);\n std::cout &lt;&lt; \&#34;(4)\\n\&#34;;\n std::linalg::hermitian_matrix_rank_1_update(\n std::execution::par,\n -1.0,\n x,\n A,\n std::linalg::upper_triangle);\n print_mat(A);\n\n return 0;\n}\n&#34;}], &#34;page_id&#34;: [&#34;reference&#34;, &#34;linalg&#34;, &#34;hermitian_matrix_rank_1_update&#34;]}">
8888
<nav class="navbar navbar-default" role="navigation">
8989
<div class="container-fluid">
9090
<div class="navbar-header">
@@ -204,8 +204,8 @@
204204

205205
<p class="text-right"><small>
206206
最終更新日時(UTC):
207-
<span itemprop="datePublished" content="2024-07-04T09:27:11">
208-
2024年07月04日 09時27分11秒
207+
<span itemprop="datePublished" content="2024-07-06T05:58:05">
208+
2024年07月06日 05時58分05秒
209209
</span>
210210
<br/>
211211
<span itemprop="author" itemscope itemtype="http://schema.org/Person">
@@ -306,7 +306,7 @@ <h2>備考</h2>
306306
<p>(3), (4)は$A \leftarrow A - xx^T$を行うために用意された。</p>
307307
<h2>例</h2>
308308
<p><strong>[注意] 処理系にあるコンパイラで確認していないため、間違っているかもしれません。</strong></p>
309-
<p><div class="yata" id="a2a66a0d88359edc033efe7cd3c98f8f6e38cae0"><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>
309+
<p><div class="yata" id="1792939fb57232efc6310058b9ce3580f17d6fa9"><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>
310310
<span class="cp">#include</span> <span class="cpf"><a href="../iostream.html">&lt;iostream&gt;</a></span><span class="cp"></span>
311311
<span class="cp">#include</span> <span class="cpf"><a href="../linalg.html">&lt;linalg&gt;</a></span><span class="cp"></span>
312312
<span class="cp">#include</span> <span class="cpf"><a href="../mdspan.html">&lt;mdspan&gt;</a></span><span class="cp"></span>
@@ -348,7 +348,6 @@ <h2>例</h2>
348348

349349
<span class="n"><a href="../vector/vector.html">std::vector</a></span><span class="o">&lt;</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;&gt;</span> <span class="n">A_vec</span><span class="p">(</span><span class="n">N</span> <span class="o">*</span> <span class="n">N</span><span class="p">);</span>
350350
<span class="n"><a href="../vector/vector.html">std::vector</a></span><span class="o">&lt;</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;&gt;</span> <span class="n">x_vec</span><span class="p">(</span><span class="n">N</span><span class="p">);</span>
351-
<span class="n"><a href="../array/array.html">std::array</a></span><span class="o">&lt;</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">N</span><span class="o">&gt;</span> <span class="n">y_vec</span><span class="p">;</span>
352351

353352
<span class="n"><a href="../mdspan/mdspan.html">std::mdspan</a></span><span class="o">&lt;</span>
354353
<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>
@@ -358,15 +357,13 @@ <h2>例</h2>
358357
<span class="n">std</span><span class="o">::</span><span class="n">linalg</span><span class="o">::</span><span class="n">row_major_t</span><span class="o">&gt;</span>
359358
<span class="o">&gt;</span> <span class="n">A</span><span class="p">(</span><span class="n">A_vec</span><span class="p">.</span><span class="n">data</span><span class="p">());</span>
360359
<span class="n"><a href="../mdspan/mdspan.html">std::mdspan</a></span> <span class="nf">x</span><span class="p">(</span><span class="n">x_vec</span><span class="p">.</span><span class="n">data</span><span class="p">(),</span> <span class="n">N</span><span class="p">);</span>
361-
<span class="n"><a href="../mdspan/mdspan.html">std::mdspan</a></span> <span class="nf">y</span><span class="p">(</span><span class="n">y_vec</span><span class="p">.</span><span class="n">data</span><span class="p">(),</span> <span class="n">N</span><span class="p">);</span>
362360

363361
<span class="n">init_mat</span><span class="p">(</span><span class="n">A</span><span class="p">);</span>
364362
<span class="n">init_vec</span><span class="p">(</span><span class="n">x</span><span class="p">);</span>
365-
<span class="n">init_vec</span><span class="p">(</span><span class="n">y</span><span class="p">);</span>
366363

367364
<span class="c1">// (1)</span>
368365
<span class="n"><a href="../iostream/cout.html">std::cout</a></span> <span class="o">&lt;&lt;</span> <span class="s">"(1)</span><span class="se">\n</span><span class="s">"</span><span class="p">;</span>
369-
<span class="n">std</span><span class="o">::</span><span class="n">linalg</span><span class="o">::</span><span class="n">matrix_rank_1_update</span><span class="p">(</span>
366+
<span class="n">std</span><span class="o">::</span><span class="n">linalg</span><span class="o">::</span><span class="n">hermitian_matrix_rank_1_update</span><span class="p">(</span>
370367
<span class="n">x</span><span class="p">,</span>
371368
<span class="n">A</span><span class="p">,</span>
372369
<span class="n">std</span><span class="o">::</span><span class="n">linalg</span><span class="o">::</span><span class="n">upper_triangle</span><span class="p">);</span>

0 commit comments

Comments
 (0)