# Tag Info

Accepted

### The point in a plane that minimizes the sum of distances to points outside it

If $x_1$ and $x_2$ are on different sides of $V$, you want $x^*$ to be where the straight line from $x_1$ to $x_2$ intersects $V$. The reason that this is true should be rather obvious, but can be ...

1 vote

### Efficiently compute rows of matrix which are not positive combinations of other rows

$\newcommand{\R}{\mathrm{I\!R}}\newcommand{\cone}{\mathrm{cone}}$I will give an answer regarding the columns of $A$ (I find it more convenient). If you are interested in the rows of $A$ you can apply ...
1 vote
Accepted

### How do Krylov methods performance wise compare to direct methods when no preconditioning is done?

Poorly. The fundamental problem is that Krylov subspace methods have low arithmetic intensity and are incompatible with modern computer architecture In contrast, modern implementations of the QR and ...
1 vote

### Cholesky-like decomposition that works on singular matrices?

An approach using $QR$ decomposition: the following results in a decomposition of the form $A = MM^T$ with $M$ of full rank of size $n \times r$ (with $r$ equal to the rank of $A$), but $M$ is not ...
1 vote

### Computing $(A^TA)A^T$ versus $A^T(AA^T)$ on a Computer with Limited Memory where $A$ is $1000000 \times 2$

The product of an $m\times n$ and an $n \times k$ matrix is an $m \times k$ matrix, and each entry requires $n$ multiplications to compute, for a total of $mnk$ multiplications (assuming the naive ...
1 vote

### Gram matrix of QR decomposition

The projection $\hat x$ of $x$ onto the image of $A$ is $$\hat x = A(A^T A)^{-1} A^T x. \tag{1}$$ so $y = x - 2 A(A^T A)^{-1} A^T x$. This is a well known formula in the topic "linear ...
1 vote
Accepted

### Solution of Diffusion equation with exponential integrator

There are several options: Find a Similarity transform $S$ that diagonalizes $A$, and then $e^{A dt}=S^{-1} e^{Ddt}S$ with $D$ diagonal. Realize that in any case you are dealing with an discrete time ...
1 vote
Accepted

### Arnoldi iteration and minimal poly

As I note in my comment, your approach is incorrect. The key is to note that Arnoldi iteration runs to completion with the initial vector $b$ if and only if the Krylov subspace generated by $b$ is ...

Only top scored, non community-wiki answers of a minimum length are eligible