# Tagged Questions

Numerical algorithms for problems in analysis and algebra, scientific computation

375 views

### sign-flipping inverse

Consider this matrix: $Z=\begin{bmatrix}23.9 & -7 & -17 \\\\ -7 & 23.9 & -17 \\\\ -17 & -17 & 33.9 \end{bmatrix}$ Its inverse is entrywise negative (you can check...) and ...
337 views

### Fast inversion of a special kind of matrices - approximations are ok

Suppose I have a stochastic matrix $M$ (with thousands or millions of stochastic column vectors), which I split into two matrices: $D$ containing only the diagonal entries of $M$, and $R$ containing ...
432 views

### Update the inverse of sum of two symmetric matrices

There are two invertible symmetric matrices A and B, of which B is a block diagonal. A and B have the same dimensions. I need to iteratively calculate the inverse of M = s * A + B, where s is a ...
13 views

### Numerical evaluation of orthogonal polynomials [on hold]

I've written some Matlab procedures that evaluate orthogonal polynomials, and as a sanity check I was trying to ensure that their dot product would be zero. But, while I'm fairly sure there's not ...
37 views

### What's the advantage of majorization-minimization (MM) algorithm [on hold]

The majorization-minimization (MM) algorithm is a framework for convex and nonconvex optimization. When applied to nonconvex optimization, the MM algorithm solves a sequence of convex problems to ...
1k views

### How to solve a fifth degree polynomIal

Charles Hermite have created a method using elliptic functions to solve fifth degree polynomial, to get around the theory of Galois. Can someone explain me it and give a simple example? Tank you.
2k views

### What is… A Grossone?

Y. Sergeyev developed a positional system for representing infinite numbers using a basic unit called a "grossone", as well as what he calls an "infinity computer". The mathematical value of this ...
373 views

### Discrete Wavelet Transform and L2 Basis

Using the mother wavlet $phi$ one obtains an orthonormal basis $\phi_{j,k}(x):=2^{j/2}\,\phi(2^j\,x-k)$of L^2 (on the unit interval say). Given a function $f$ on can calculate the coefficients using ...
6k views

### There must be a good introductory numerical analysis course out there!

Background As a numerical analyst, I've frequently taught the 'Introductory Numerical Analysis' class. Such courses are found in many major universities; the audience typically consists of reluctant ...
89 views

### Evaluate a Function to Full Machine Precision [closed]

If we want to evaluate $$f(x)=\frac{e^x-1-x}{x^2}$$ then we have to observe its large relative error as $x\to 0$. My question is that how can we find a method so that we can compute $f(x)$ to full ...
246 views

### Why Householder reflection is better than Givens rotation in dense linear algebra?

It’s obvious that Givens rotation works better with sparse matrices. But I don’t know why Householder reflection is better for dense matrices. Does it require less computations? Or it’s numerically ...
207 views

### Numerical solution of SDEs with colored noise

I am trying to numerically solve an SDE with both white and colored noise that models a non-linear circuit: $$dX_t = f(X_t) dt + \sigma_w dW + \sigma_c dC$$ where $W$ is a standard Brownian motion ...
579 views

### Numerical solution to Fisher-Kolmogorov equation

Can you use the Crank-Nicolson method to get a numerical approximation to the fisher-kolmogorov equation? If not what would be the easiest way to model the equation using matlab? Thanks and sorry its ...
25 views

### state-of-art numerical contour (complex) integration method when contour is square and available values are evenly spaced

What is current state-of-art for numerical contour integration method (for $f(z)$ with $z$ being complex number and $f$ complex-valued) when contour is square on complex plane, and one only has ...
22 views

### numerical differentiation of sum of one-dimensional sinusoids with angular frequency close to Nyquist one

Suppose that $f(t) = \sum_i C_i e^{i\omega_i t}$, and $f$ is sampled at certain sampling angular frequency $\omega_s$. All $\omega_i$s are very close to $\omega_s/2$, and thus standard finite ...
49 views

### Injectivity of vector functions: Numerical Verification

Problem Setup Let $f:A\rightarrow B$, be a continuous function, $A\subset\Re^{n}$,$B\subset\Re^{m}$, $m\geq n$ and $A, B$ compact. The function $f(\cdot)$ can only be evaluated numerically. ...
103 views

### Reference Request: Variational Problem

I want to solve approximately the following variational problem: Given a function $c:[-1,1]^2\rightarrow [0,1]$, constants $p_1...p_n\in \mathbb{R}^+$, $\alpha_1...\alpha_n\in \mathbb{R}$, and ...
111 views

21k views

### Why were matrix determinants once such a big deal?

I have been told that the study of matrix determinants once comprised the bulk of linear algebra. Today, few textbooks spend more than a few pages to define it and use it to compute a matrix inverse. ...
401 views

### “Fractally self-similar” numbers

This is another question about visualization of Ford circles, the previous one being Confusion with practically implementing rational approximations. Here is an output of zooming into Ford circles at ...
204 views

### Confusion with practically implementing rational approximations

Writing a program visualizing Ford circles I've encountered a seemingly purely programmatic puzzle but then gradually realized there are some mathematical aspects of it which I don't understand. Let ...
190 views

### Existence of nonlinear equation

How can we prove that equation (1) has solutions for every $p$. I mean, is there an analytic method that can be used to show that there exist solutions for every $p$ for this nonlinear equation: ...
4k views

### Meaning of $\Subset$ notation

The symbol $\Subset$ (occurring in places where $\subseteq$ could occur syntactically) comes up frequently in a paper I'm reading. The paper lives at the intersection of a few areas of math, and I ...
21 views

### Monotonicity per dimension of multivariate scattered data

For my thesis, I am working on interpolation using the RBF method (Radial Basis Functions). Before interpolating, I want some a priori insight into the data, for example check in which dimensions it ...
68 views

### Difference stencils approximating Laplacian

Let $\Delta$ be the Laplace operator on the interval $[0,1]\subset \mathbb{R}$. Divide $[0,1]$ into small intervals of size $h$ to get an equidistant grid. One can approminate $-\Delta$ on this grid ...
125 views

### The condition number of a scaled Vandermonde matrix

Let $V(x_1,..,x_n)$ be the Vandermonde matrix induced by $x_1,..,x_n$, and let $\tilde{V} := V(\frac{x_1}{h},...,\frac{x_n}{h})$. My intuition says that the condition number should be invariant under ...
130 views

### Iterative Root Finding

Consider a function $f(x)=g(x,h(x))$, which we know has a unique root $x^*$. The functions $f$, $g$ and $h$ are all continuous in $x$ and behave nicely. Iteratively solving $g(x_{i+1}, h(x_i))=0$ with ...
99 views

When the term Gaussian Quadrature appears in most Literatures, does it actually refer to Gauss-Legendre Quadrature. In other words, do they implicitly admit that they use the Legendre orthogonal ...
Consider Galois fields $\mathbb{F}_{2^n}$ and $\mathbb{F}_{2^k}$, where $n=km$, and $\mathbb{F}_{2^k}$ is a ground field of $\mathbb{F}_{2^n}$. I’d appreciate pointers to papers or suggestions on: ...