A[:,1] * 3 isa Vector</code></pre><pre class="documenter-example-output"><code class="nohighlight hljs ansi">true</code></pre><pre><code class="language-julia hljs">opA[:,1] * 3 isa LinearOperator</code></pre><pre class="documenter-example-output"><code class="nohighlight hljs ansi">true</code></pre><pre><code class="language-julia hljs">opA[:,1] * [3] isa Vector</code></pre><pre class="documenter-example-output"><code class="nohighlight hljs ansi">true</code></pre><p>However, the following returns an error</p><pre><code class="language-julia hljs">A[:,1] * [3]</code></pre><p>This is also true for <code>A[i,:]</code>, which would return a vector and for the scalar <code>A[i,j]</code>. Similarly, <code>opA[1,1]</code> is an operator of size (1,1):"</p><pre><code class="language-julia hljs">(opA[1,1] * [3])[1] - A[1,1] * 3</code></pre><pre class="documenter-example-output"><code class="nohighlight hljs ansi">0.0</code></pre><p>In the same spirit, the operator <code>Matrix</code> always returns a matrix.</p><pre><code class="language-julia hljs">Matrix(opA[:,1])</code></pre><pre class="documenter-example-output"><code class="nohighlight hljs ansi">5×1 Matrix{Float64}:
0 commit comments