Questions tagged [numerical-linear-algebra]
Questions on the various algorithms used in linear algebra computations (matrix computations).
3,669 questions
0 votes
0 answers
31 views
Efficient sparse matrix vector multiplication at low 30-90% unstructured sparsity. Relevant use cases?
I have been working on an algorithm for efficient matrix–vector multiplication on GPUs in the regime of moderate sparsity (30%–90%), with no structural assumptions on the sparsity pattern. The ...
1 vote
2 answers
112 views
Inversion of an ill-conditioned triangular matrix
Let $L$ be a lower triangular matrix where all elements on and below the diagonal are positive. The matrix is also highly ill-conditioned. I want to solve $Lx=1$. Additionally, I know that all ...
0 votes
0 answers
72 views
Ratio of eigenvector amplitudes
I am studying a generalized eigenvalue problem which can be partitioned as $$ \left[ {\begin{array}{*{20}{c}} {{A_{11}}}&{{A_{12}}}\\ {{A_{21}}}&{{A_{22}}} \end{array}} \right]\left[ {\begin{...
3 votes
1 answer
86 views
Factorizing a matrix as a product of band matrices
Let $A$ be $n\times n$ matrix. Can $A$ be written as $A=B_1\cdots B_k$ where $B_i$ are band matrices with constant bandwidth, and $k=O(n)$? Is it always possible? Is there an efficient algorithm for ...
0 votes
0 answers
42 views
Stability of specific floating point operation algorithm
Problem: We want to determine whether the following algorithm is stable or not. Data is $x_{1},x_{2} \in \mathbb{C}$, Solution is $x_{1}(x_{2}+1)$, computed as $\text{fl}(x_{1}) \otimes (\text{fl}(x_{...
0 votes
0 answers
24 views
Inequality of Perturbation of Linear System [duplicate]
Problem: Let $A \in \mathbb{R}^{m \times n}$ be a nonsingular matrix. Consider the system $Ax = b$. If $b$ is perturbed by $\delta b$, suppose that the solution of the original system is perturbed by $...
2 votes
0 answers
69 views
Using eigenvalue to estimate the accuracy of the solution of a system
I'm studying for an exam and I came into the following question regarding the use of eigenvalues to estimate a solution. I'm having trubles understating what the request is, can somebody give me a ...
0 votes
0 answers
20 views
$LU$ decomposition and Leading principal submatrix [duplicate]
Problem Setup: Let $A \in \mathbb{C}^{m \times m}$ be nonsingular. For each $k$ with $1 \leq k \leq m$, the upper-left $k \times k$ block $A_{1:k,1:k}$ is nonsingular. We want to show that $A$ has an ...
2 votes
0 answers
85 views
How to use Exercise 2.1 to solve Exercise 2.5(a)? (Numerical Linear Algebra by Lloyd N. Trefethen and David Bau, III.)
I am reading Numerical Linear Algebra by Lloyd N. Trefethen and David Bau, III. On p.16 in this book: Exercise 2.5. Let $S\in\mathbb{C}^{m\times m}$ be skew-hermitian, i.e., $S^*=-S$. (a) Show by ...
0 votes
2 answers
60 views
Condition Number of $A$
Let $A = \begin{bmatrix} 1.0000 & 2.0000\\ 1.0001 & 2.0000 \end{bmatrix}$. Suppose we wish to find $Ax = b$, where $b = (3.0000,3.0001)^T$. Instead of $x$, we obtain $x' = ...
0 votes
0 answers
82 views
Geometric Effect of Householder Reflector
We have the Householder reflector $F = \begin{bmatrix} -c & s\\ s &c \end{bmatrix}$, where $c$ stands for $\cos{\theta}$, $s$ stands for $\sin{\theta}$ for some $\theta$. We want to know its ...
0 votes
0 answers
51 views
Nozero diagonal entries imply full rank
Let $A$ be an $m \times n$ matrix with $m \geq n$. And let $A = \hat{Q}\hat{R}$ be the reduced $QR$ factorization. Suppose that all the diagonal matrix entries $\hat{R}$ are nonzero. We want to show ...
1 vote
0 answers
41 views
Bandwidth of a matrix and reducing fill-in during LU factorization
The bandwidth of a matrix $M=(m_{ij})$ is defined as the maximum value of $|i-j|$ such that $m_{ij}$ is nonzero. During LU factorization for solving a linear system, reducing the bandwidth of a matrix ...
1 vote
1 answer
82 views
I need help with solving a double integral exp of a quadratic variable
I have a double integral, which I want to solve: $\mathcal{I}_{au}=\int_0^{t'}e^{(a-c)t''}\int_0^{t''}e^{(-b+c)t'''}\text{exp}\big(-\frac{\lambda^2}{2}\left( S_2(t''')^2+S_1(t')^2+S_1(t'')^2+2S_{12}t'...
0 votes
1 answer
176 views
How do you compute the solution to least-squares problem if neither $A^TA$ nor $AA^T$ nor $A$ are invertible?
For a least-squares problem find $x$ such that $\|Ax - b\|_2, A \in \mathbb{R}^{m \times n}$ is minimized, the solution of is captured by the pseudo-inverse, $$x = A^\dagger b$$ There exists three ...