Website build

Author: colinjcotter

Computationallinearalgebracourse.pdf

@@ -193,11 +193,11 @@ <h2><span class="section-number">2.2. </span>QR factorisation by classical Gram-
 <p>The <a class="reference internal" href="cla_utils.html#cla_utils.exercises2.GS_classical" title="cla_utils.exercises2.GS_classical"><code class="xref py py-func docutils literal notranslate"><span class="pre">cla_utils.exercises2.GS_classical()</span></code></a> function has been
 left unimplemented. It should implement the classical Gram-Schmidt
 algorithm above, using Numpy slice notation so that only one Python
-for loop is used. The function should work “in place” by making a
-copy of <span class="math notranslate nohighlight">\(A\)</span> and then changing the values in the copy, without
-introducing additional intermediate arrays (you will need to create
-a new array to store <span class="math notranslate nohighlight">\(R\)</span>). The test script <code class="docutils literal notranslate"><span class="pre">test_exercises2.py</span></code>
-in the <code class="docutils literal notranslate"><span class="pre">test</span></code> directory will test this function.</p>
+for loop is used. The function should work “in place” by changing
+the values in <span class="math notranslate nohighlight">\(A\)</span>, without introducing additional intermediate
+arrays (you will need to create a new array to store <span class="math notranslate nohighlight">\(R\)</span>). The test
+script <code class="docutils literal notranslate"><span class="pre">test_exercises2.py</span></code> in the <code class="docutils literal notranslate"><span class="pre">test</span></code> directory will test
+this function.</p>
 </div></div></section>
 <section id="projector-interpretation-of-gram-schmidt">
 <h2><span class="section-number">2.3. </span>Projector interpretation of Gram-Schmidt<a class="headerlink" href="#projector-interpretation-of-gram-schmidt" title="Permalink to this headline">¶</a></h2>
@@ -307,9 +307,8 @@ <h2><span class="section-number">2.4. </span>Modified Gram-Schmidt<a class="head
 algorithm above, using Numpy slice notation where possible.
 What is the minimal number of Python
 for loops possible?</p>
-<p>The function should work “in place” by making a
-copy of <span class="math notranslate nohighlight">\(A\)</span> and then changing those values, without introducing
-additional intermediate arrays. The test script
+<p>The function should work “in place” by changing the values of <span class="math notranslate nohighlight">\(A\)</span>,
+without introducing additional intermediate arrays. The test script
 <code class="docutils literal notranslate"><span class="pre">test_exercises2.py</span></code> in the <code class="docutils literal notranslate"><span class="pre">test</span></code> directory will test this
 function.</p>
 </div></div><div class="proof proof-type-exercise" id="id6">
@@ -682,7 +681,7 @@ <h2><span class="section-number">2.6. </span>Householder triangulation<a class="
 done so, you will need to modified
 <a class="reference internal" href="cla_utils.html#cla_utils.exercises3.householder" title="cla_utils.exercises3.householder"><code class="xref py py-func docutils literal notranslate"><span class="pre">cla_utils.exercises3.householder()</span></code></a> to use the <code class="docutils literal notranslate"><span class="pre">kmax</span></code>
 argument. You may make use of the built-in triangular solve
-algorithm <a class="reference external" href="http://scipy.github.io/devdocs/reference/generated/scipy.linalg.solve_triangular.html#scipy.linalg.solve_triangular" title="(in SciPy v1.8.0.dev0+1767.0785d6f)"><code class="xref py py-func docutils literal notranslate"><span class="pre">scipy.linalg.solve_triangular()</span></code></a> (we shall consider
+algorithm <a class="reference external" href="http://scipy.github.io/devdocs/reference/generated/scipy.linalg.solve_triangular.html#scipy.linalg.solve_triangular" title="(in SciPy v1.8.0.dev0+1859.14341ea)"><code class="xref py py-func docutils literal notranslate"><span class="pre">scipy.linalg.solve_triangular()</span></code></a> (we shall consider
 triangular matrix algorithms briefly later). The test script
 <code class="docutils literal notranslate"><span class="pre">test_exercises3.py</span></code> in the <code class="docutils literal notranslate"><span class="pre">test</span></code> directory will also test this
 function.</p>
@@ -794,7 +793,7 @@ <h2><span class="section-number">2.7. </span>Application: Least squares problems
 appropriate augmented matrix <span class="math notranslate nohighlight">\(\hat{A}\)</span>, calling
 <a class="reference internal" href="cla_utils.html#cla_utils.exercises3.householder" title="cla_utils.exercises3.householder"><code class="xref py py-func docutils literal notranslate"><span class="pre">cla_utils.exercises3.householder()</span></code></a> and extracting appropriate
 subarrays using slice notation, before using
-<a class="reference external" href="http://scipy.github.io/devdocs/reference/generated/scipy.linalg.solve_triangular.html#scipy.linalg.solve_triangular" title="(in SciPy v1.8.0.dev0+1767.0785d6f)"><code class="xref py py-func docutils literal notranslate"><span class="pre">scipy.linalg.solve_triangular()</span></code></a> to solve the resulting upper triangular
+<a class="reference external" href="http://scipy.github.io/devdocs/reference/generated/scipy.linalg.solve_triangular.html#scipy.linalg.solve_triangular" title="(in SciPy v1.8.0.dev0+1859.14341ea)"><code class="xref py py-func docutils literal notranslate"><span class="pre">scipy.linalg.solve_triangular()</span></code></a> to solve the resulting upper triangular
 system, before returning the solution <span class="math notranslate nohighlight">\(x\)</span>. The test script
 <code class="docutils literal notranslate"><span class="pre">test_exercises3.py</span></code> in the <code class="docutils literal notranslate"><span class="pre">test</span></code> directory will also test this
 function.</p>

L2_QR_factorisation.html

@@ -6,15 +6,14 @@
     <meta charset="utf-8" />
     <meta name="viewport" content="width=device-width, initial-scale=1.0" />
     <title>cla_utils.exercises1 &#8212; Computational linear algebra course 2020.0 documentation</title>
-    <link rel="stylesheet" href="../../_static/pygments.css" type="text/css" />
-    <link rel="stylesheet" href="../../_static/fenics.css" type="text/css" />
+    <link rel="stylesheet" type="text/css" href="../../_static/pygments.css" />
+    <link rel="stylesheet" type="text/css" href="../../_static/fenics.css" />
     <link rel="stylesheet" type="text/css" href="../../_static/proof.css" />
-    <script id="documentation_options" data-url_root="../../" src="../../_static/documentation_options.js"></script>
+    <script data-url_root="../../" id="documentation_options" src="../../_static/documentation_options.js"></script>
     <script src="../../_static/jquery.js"></script>
     <script src="../../_static/underscore.js"></script>
     <script src="../../_static/doctools.js"></script>
     <script src="../../_static/proof.js"></script>
-    <script async="async" src="https://cdn.jsdelivr.net/npm/mathjax@2/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
     <link rel="index" title="Index" href="../../genindex.html" />
     <link rel="search" title="Search" href="../../search.html" />
 <!--[if lte IE 6]>
@@ -133,7 +132,7 @@ <h1>Source code for cla_utils.exercises1</h1><div class="highlight"><pre>
 
 <div class="viewcode-block" id="rank2"><a class="viewcode-back" href="../../cla_utils.html#cla_utils.exercises1.rank2">[docs]</a><span class="k">def</span> <span class="nf">rank2</span><span class="p">(</span><span class="n">u1</span><span class="p">,</span> <span class="n">u2</span><span class="p">,</span> <span class="n">v1</span><span class="p">,</span> <span class="n">v2</span><span class="p">):</span>
     <span class="sd">&quot;&quot;&quot;</span>
-<span class="sd">    Return the rank2 matrix A = u1*v2^* + u2*v2^*.</span>
+<span class="sd">    Return the rank2 matrix A = u1*v1^* + u2*v2^*.</span>
 
 <span class="sd">    :param u1: m-dimensional numpy array</span>
 <span class="sd">    :param u2: m-dimensional numpy array</span>
@@ -186,7 +185,7 @@ <h1>Source code for cla_utils.exercises1</h1><div class="highlight"><pre>
     </div>
     <div class="footer" role="contentinfo">
         &#169; Copyright 2020, Colin J. Cotter.
-      Created using <a href="https://www.sphinx-doc.org/">Sphinx</a> 3.4.3.
+      Created using <a href="https://www.sphinx-doc.org/">Sphinx</a> 4.2.0.
     </div>
   </body>
 </html>

_modules/cla_utils/exercises1.html

@@ -129,7 +129,7 @@ <h1>Source code for cla_utils.exercises2</h1><div class="highlight"><pre>
 
     <span class="k">raise</span> <span class="ne">NotImplementedError</span>
 
-    <span class="k">return</span> <span class="n">Q</span><span class="p">,</span> <span class="n">R</span></div>
+    <span class="k">return</span> <span class="n">R</span></div>
 
 <div class="viewcode-block" id="GS_modified"><a class="viewcode-back" href="../../cla_utils.html#cla_utils.exercises2.GS_modified">[docs]</a><span class="k">def</span> <span class="nf">GS_modified</span><span class="p">(</span><span class="n">A</span><span class="p">):</span>
     <span class="sd">&quot;&quot;&quot;</span>

_modules/cla_utils/exercises2.html

@@ -161,11 +161,11 @@ for code blocks in pseudo-code in these notes.)
    The :func:`cla_utils.exercises2.GS_classical` function has been
    left unimplemented. It should implement the classical Gram-Schmidt
    algorithm above, using Numpy slice notation so that only one Python
-   for loop is used. The function should work "in place" by making a
-   copy of `A` and then changing the values in the copy, without
-   introducing additional intermediate arrays (you will need to create
-   a new array to store `R`). The test script ``test_exercises2.py``
-   in the ``test`` directory will test this function.
+   for loop is used. The function should work "in place" by changing
+   the values in `A`, without introducing additional intermediate
+   arrays (you will need to create a new array to store `R`). The test
+   script ``test_exercises2.py`` in the ``test`` directory will test
+   this function.
 
 Projector interpretation of Gram-Schmidt
 ----------------------------------------
@@ -308,9 +308,8 @@ in `A` with the `v` s and eventually the `q` s.
    What is the minimal number of Python
    for loops possible?
 
-   The function should work "in place" by making a
-   copy of `A` and then changing those values, without introducing
-   additional intermediate arrays. The test script
+   The function should work "in place" by changing the values of `A`,
+   without introducing additional intermediate arrays. The test script
    ``test_exercises2.py`` in the ``test`` directory will test this
    function.
 

_sources/L2_QR_factorisation.rst.txt

@@ -7,3 +7,11 @@ Errata for 2021/22
 Here we'll list all the updates to the exercises and coursework Classroom
 repositories for Autumn 2021. An up to date version is in the
 `public repository <https://github.com/comp-lin-alg/comp-lin-alg-course>`_.
+
+1. The docstring for :code:`cla_utils.rank2` was updated to show the
+   correct indices. Now reads: "Return the rank2 matrix A = u_1v_1^* +
+   u_2*v_2^*." See the code on the course website for the up-to-date
+   version.
+
+2. The skeleton code for exercises2.GS_classical was amended as the
+   function should just return R, not Q and R.

_sources/errata.rst.txt

@@ -88,23 +88,9 @@ build the :doc:`cla_utils` web documentation.
 As you do the exercises, **commit your code** to your repository. This
 will build up your computational exercise solution sets. You should
 commit code early and often - small commits are easier to understand
-and debug than large ones. **Never** commit back to the ``master``
-branch of your fork, that should always remain a clean copy of the
-main repository.
-
-Pull requests for feedback
---------------------------
-
-#. Click on the ``New pull request`` button at the top of your
-   repository page on GitHub.
-#. Make sure **left** dropdown box ("base") is set to ``master``.
-#. Make sure **right** dropdown box ("compare") is set to ``implementation``.
-#. Type a suitable title in the title box. For example 
-   ``Request for feedback 30/1/19``.
-#. If you have any comments you would like to pass on to the lecturer
-   (for example questions about how you should have done a particular
-   exercise) then type these in the ``Description`` box.
-#. Click ``Create pull request``.
+and debug than large ones. **Do not commit to the feedback branch**
+This branch is just there so that we can provide feedback on your
+changes to the main branch.
 
 
 Testing your work

_sources/exercises.rst.txt

@@ -131,7 +131,7 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 <dl class="py function">
 <dt class="sig sig-object py" id="cla_utils.exercises1.rank2">
 <span class="sig-prename descclassname"><span class="pre">cla_utils.exercises1.</span></span><span class="sig-name descname"><span class="pre">rank2</span></span><span class="sig-paren">(</span><em class="sig-param"><span class="n"><span class="pre">u1</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">u2</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">v1</span></span></em>, <em class="sig-param"><span class="n"><span class="pre">v2</span></span></em><span class="sig-paren">)</span><a class="reference internal" href="_modules/cla_utils/exercises1.html#rank2"><span class="viewcode-link"><span class="pre">[source]</span></span></a><a class="headerlink" href="#cla_utils.exercises1.rank2" title="Permalink to this definition">¶</a></dt>
-<dd><p>Return the rank2 matrix A = u1*v2^* + u2*v2^*.</p>
+<dd><p>Return the rank2 matrix A = u1*v1^* + u2*v2^*.</p>
 <dl class="field-list simple">
 <dt class="field-odd">Parameters</dt>
 <dd class="field-odd"><ul class="simple">

cla_utils.html

@@ -52,6 +52,14 @@ <h1>Errata for 2021/22<a class="headerlink" href="#errata-for-2021-22" title="Pe
 <p>Here we’ll list all the updates to the exercises and coursework Classroom
 repositories for Autumn 2021. An up to date version is in the
 <a class="reference external" href="https://github.com/comp-lin-alg/comp-lin-alg-course">public repository</a>.</p>
+<ol class="arabic simple">
+<li><p>The docstring for <code class="code docutils literal notranslate"><span class="pre">cla_utils.rank2</span></code> was updated to show the
+correct indices. Now reads: “Return the rank2 matrix A = u_1v_1^* +
+u_2*v_2^*.” See the code on the course website for the up-to-date
+version.</p></li>
+<li><p>The skeleton code for exercises2.GS_classical was amended as the
+function should just return R, not Q and R.</p></li>
+</ol>
 </section>
 
 

errata.html

@@ -124,24 +124,9 @@ <h2>How to do the computational exercises<a class="headerlink" href="#how-to-do-
 <p>As you do the exercises, <strong>commit your code</strong> to your repository. This
 will build up your computational exercise solution sets. You should
 commit code early and often - small commits are easier to understand
-and debug than large ones. <strong>Never</strong> commit back to the <code class="docutils literal notranslate"><span class="pre">master</span></code>
-branch of your fork, that should always remain a clean copy of the
-main repository.</p>
-</section>
-<section id="pull-requests-for-feedback">
-<h2>Pull requests for feedback<a class="headerlink" href="#pull-requests-for-feedback" title="Permalink to this headline">¶</a></h2>
-<ol class="arabic simple">
-<li><p>Click on the <code class="docutils literal notranslate"><span class="pre">New</span> <span class="pre">pull</span> <span class="pre">request</span></code> button at the top of your
-repository page on GitHub.</p></li>
-<li><p>Make sure <strong>left</strong> dropdown box (“base”) is set to <code class="docutils literal notranslate"><span class="pre">master</span></code>.</p></li>
-<li><p>Make sure <strong>right</strong> dropdown box (“compare”) is set to <code class="docutils literal notranslate"><span class="pre">implementation</span></code>.</p></li>
-<li><p>Type a suitable title in the title box. For example
-<code class="docutils literal notranslate"><span class="pre">Request</span> <span class="pre">for</span> <span class="pre">feedback</span> <span class="pre">30/1/19</span></code>.</p></li>
-<li><p>If you have any comments you would like to pass on to the lecturer
-(for example questions about how you should have done a particular
-exercise) then type these in the <code class="docutils literal notranslate"><span class="pre">Description</span></code> box.</p></li>
-<li><p>Click <code class="docutils literal notranslate"><span class="pre">Create</span> <span class="pre">pull</span> <span class="pre">request</span></code>.</p></li>
-</ol>
+and debug than large ones. <strong>Do not commit to the feedback branch</strong>
+This branch is just there so that we can provide feedback on your
+changes to the main branch.</p>
 </section>
 <section id="testing-your-work">
 <h2>Testing your work<a class="headerlink" href="#testing-your-work" title="Permalink to this headline">¶</a></h2>