9 Nov 2012 vicious   » (Master)

Numerical range

I was fiddling with numerical range of two by two matrices so I modified my root testing python program to do this. The numerical range of A is the set of all values

\frac{v^* A v}{v^* v}

for all nonzero vectors v. This set is a compact set (it can be seen as the image of the mapping $latex v^* A v of vectors on the unit sphere v^* v = 1 which is a compact set). It’s convex which is harder to show. For two by two it is an elliptic disc (could be degenerate).

See the result here, it plugs in random vectors and shows the result. Here’s an example plot for the matrix A=\begin{bmatrix} 1 & i \\ 1 & -1 \end{bmatrix}.

Numerical range of 2x2 matrix

The code is really inefficient and eats up all your cpu. There’s no effort to optimize this.


Syndicated 2012-11-09 00:06:40 from The Spectre of Math

Latest blog entries     Older blog entries

New Advogato Features

New HTML Parser: The long-awaited libxml2 based HTML parser code is live. It needs further work but already handles most markup better than the original parser.

Keep up with the latest Advogato features by reading the Advogato status blog.

If you're a C programmer with some spare time, take a look at the mod_virgule project page and help us with one of the tasks on the ToDo list!