1
vote
1answer
26 views

Solving Ax = b where A is composed of diagonal blocks

I would like to solve the equation $Ax=b$ where $x\in\mathbb{R}^n$ and $A$ is of the form: $$A= \begin{bmatrix} D_1 & D_2 &D_3 \\ D_2 & D_4 & D_5 \\ D_3 & D_5 & D_6 ...
0
votes
1answer
25 views

Proving the equality: $D^{-1}+D^{-1}CPBD^{-1}=(D-CA^{-1}B)^{-1}$ where $P=(A-BD^{-1}C)^{-1}$

I am trying to prove the following equality that I need to use as an intermediate step to solve one of my problems. The equality is the following: $D^{-1}+D^{-1}CPBD^{-1}=(D-CA^{-1}B)^{-1}$, where ...
2
votes
0answers
44 views

linear algebra formulation help

I asked this on mathoverflow as well and apologies for cross-posting. I am trying to compute this so-called bending energy matrix. The bending energy of a thin plate in 3D is given by: $$ BE = ...
1
vote
1answer
246 views

Jordan Canonical form 2x2 matrix

Compute the Jordan Canonical form of A = $\begin{bmatrix}i & 1\\1 & -1\end{bmatrix}$. My (feeble) attempt: After I compute the characteristic polynomial, which gives me $x^2=0$, the ...
0
votes
0answers
116 views

Row & Column Removal and Rank Reduction

I have a problem involving a n x n square, real matrix $K$ which is initially full rank and is not positive definite. In each iteration of my program, I have to remove a row and the corresponding ...
4
votes
1answer
167 views

Eigenvector of a sparse structured matrix corresponding to the eigenvalue 1

I have a matrix with the following sparsity pattern: $M = \begin{bmatrix} \ast &\ast &0 &0 &0 &0 &0 &0\\ 0 & 0 &\ast &\ast &0 &0 &0 &0 \\ 0 ...
1
vote
1answer
59 views

Help regarding a weird Matrix

Hi I have a matrix of the following form arising by discretization of a system of PDEs. I am working to get the invertibility of the Matrix. Can some one help me or at least give me some reference on ...
5
votes
2answers
937 views

Are there any decompositions of a symmetric matrix that allow for the inversion of any submatrix?

I am given a $J \times J$ symmetric matrix $\Sigma$. For a subset of $\{1, ..., J\}$, $v$, let $\Sigma_{v}$ denote the associated square submatrix. I am in need of an efficient scheme for inverting ...
0
votes
1answer
87 views

Iterative Method for a special kind of Sparse Matrix

I've the following problem.I've a sparse square Matrix $\bf M$. I can write $\bf M$ as: $${\mathbf M} = \begin{bmatrix}\mathbf A_{11} & \dots & \mathbf A_{1n} \\ \vdots & \ddots & ...
5
votes
0answers
76 views

Is there a way to exploit the fact that the covariance matrix has a blocked structure to more easily compute the multivariate normal density?

I'm trying to minimize the (negative) multivariate normal log likelihood (dropping constants): $$ \log |\boldsymbol\Sigma|\,+(\mathbf{x}-\boldsymbol\mu)^{\rm ...