Numerical algorithms for problems in analysis and algebra, scientific computation

**1**

vote

**1**answer

109 views

### What is a one-parameter Newton's method?

The Newton's method that I know is defined as follows:
$$x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)}$$
However, I've recently encountered a paper that talks about a one-parameter family of Newton's method ...

**0**

votes

**1**answer

118 views

### Generating random variables from the Cantor Distribution [closed]

I am looking for a method (exact, if possible, but at least asymptotically correct) for generating random variates from a Cantor Distribution? It seems like its abstract definition prevents this. In ...

**0**

votes

**0**answers

26 views

### Column Subset Selection implementations

Are there readily available implementations of algorithms for the CSSP - Column Subset Selection Problem?

**0**

votes

**0**answers

39 views

### Comparing Calculation Error in Divergent Numerical Methods

I'm not an expert in numerical methods, but I'm doing a simulation based on non-linear differential equations (General Relativity), there solutions has singularities, thus at some points numerical ...

**0**

votes

**1**answer

82 views

### What is exponentially fitted osculating straight line?

While reading an article about iterative methods for solving nonlinear equations I can't understand what is exponentially fitted osculating straight line. Could someone please briefly explain this ...

**3**

votes

**2**answers

270 views

### Analytic Solution to SDEs

Are there any example of SDEs with constant diffusion terms, other than the Ornstein Uhlenbeck process, which have exact solutions? I'm thinking of something of the form:
\begin{equation}
dX_t = ...

**-2**

votes

**1**answer

286 views

### BDF2 and TR-BDF2: what is better? [closed]

What method of numerical solving ODEs is better? BDF2 or TR-BDF2?
Namely, what advantages has TR-BDF2 over BDF2? Is TR-BDF2 more accurate? Does it damp "ringing" better?
The BDF2 method requires the ...

**0**

votes

**0**answers

28 views

### non-coherent estimation problem

I have the following signals
$$\left[\begin{array}{c} y_{mn} \\ y_{nm}\end{array}\right] =\left[\begin{array}{c} x_{n} \\ x_{m}\end{array}\right]h_{nm} +\left[\begin{array}{c} e_{mn} \\ ...

**7**

votes

**1**answer

204 views

### Approximation theory on the disc

Chebyshev theory provides a very effective method for approximating continuous real valued functions on the unit interval. Is there something similar for continuous real valued functions on the ...

**0**

votes

**0**answers

83 views

### Eigenvalue problem (finite difference operator)

Consider an arbitrary elliptic (perhaps, degenerate) finite difference operator
$$L_{i,j,k}=-\Delta_{i,j,k}+\alpha_{i,j,k}\frac{\partial}{\partial x}_{i,j,k}+\beta_{i,j,k}\frac{\partial}{\partial ...

**0**

votes

**0**answers

47 views

### Shooting method for non-oscillatory solutions

I am having some trouble with using the shooting method :
I am given a system of ode's representing initial value problem, and I know I want one of the unknowns (=u) to vanish at infinity (which I ...

**2**

votes

**1**answer

165 views

### Divergence of the Lagrange interpolation on the Chebyshev nodes

Faber theorem states that for every $\lbrace x_k^{(n)} \rbrace$ there exists a continuous $f$ function such that $\| f - L_n \|_{\infty} \not\rightarrow 0$, where $L_n$ is interpolation polynomial on ...

**0**

votes

**1**answer

130 views

### A kind of Discrete Fourier Transform

Given a $z\in \mathbb{C}^N$, the DFT of $z$ is given for every $k\in [0,N-1]_\mathbb{N}$ by
$$DFT_z(k)=\frac{1}{N} \sum_{j=0}^{N-1} z_j\, \omega^{-k j}$$ where I have denoted by $\omega$ the $N$-th ...

**0**

votes

**0**answers

137 views

### $Ax=b$ in a function space (again)

Let
$X$ be compact Hausdorff topological space,
$C(X)$ denote the algebra of complex-valued continuous functions on $X$,
$b\in \mathbb{C}^m$,
$\mathbf{A}\in C(X)^{m\times n}$,
Let ${\mathbb{C}}^n$ ...

**0**

votes

**0**answers

32 views

### Multiple integral of the resolvent kernel

I am trying to implent Volterra equations using resolvent kernel.To do this, the iterative kernel
$$K_i(x,y) = \int\limits_x^y K_1(y,t)K_{i-1}(t, x)dt. $$
should be calculated. However, it is not ...

**1**

vote

**1**answer

53 views

### Convergent algorithm for dividing a body into two regions of equal volume

Let $\Omega \subset R^3$ be a bounded open region. It is well known that there exists a smooth surface $\Gamma$ with minimum area and constant mean curvature which is orthogonal to $\partial \Omega$ ...

**1**

vote

**1**answer

223 views

### Residues and values of Riemann Zeta function at some points

I need the following computational results for proving something.
Let $1/2 + i\gamma_0$, be the first nontrivial zero of Riemann zeta function, $\zeta(s)$,
i.e. $\gamma_0\sim 14.134...$.
1) what is ...

**0**

votes

**0**answers

75 views

### Rational dynamical system with nonnegative paramaters

let $A$ be a rational system of the form :$\begin{cases} x_{n+1}=\frac{\alpha_{1}}{y_{n}} \\ y_{n+1}=\frac{\alpha_{2}}{z_{n}} \\ ...

**3**

votes

**0**answers

108 views

### error estimate of linear interpolation in high dimension

Consider convex functions $f,g$ on $[0,1]^d$. Let $x_1,\cdots,x_n$ be $n\geq d+1$ fixed point in $[0,1]^d$ that is equally 'distributed' in the sense that
$$c_1\leq ...

**4**

votes

**1**answer

181 views

### Resolvent of a triangular matrix

Suppose $A$ is a triangular matrix. What is the most efficient known algorithm to compute the polynomial (in $x$) matrix $(xI-A)^{-1}$?
Of course, $(xI-A)^{-1}= N(x)/p_A(x)$, where $p_A$ is the ...

**0**

votes

**0**answers

44 views

### Bound on change of function given bound on Hessian

Suppose I have very some smooth function $F(x)$, and let $x_0 = \text{argmin}_x F(x)$. I would like to bound $F(x) - F(x_0)$ from above, in terms of the gradient $\nabla f(x)$ and the Hessian matrix ...

**-2**

votes

**2**answers

65 views

### Systems of ODEs that fulfill a matrix relationship at steady state [closed]

It is well known that for a system of linear ODE $$x'(t) = A(t) \cdot x(t) + b(t)$$
with initial condition $x(t_0) = x_0$, that for a solution at any other time point $t_1$, $x(t_1) = (z_1, \ldots, ...

**4**

votes

**0**answers

65 views

### Estimating polynomial approximation error in high dimension

Question
Let $x \in [-1, 1]^d \subset \mathbb{R}^d$ be a $d$-dimensional variable and assume that -- given $n$ -- I have a way of computing a polynomial $p_n(x)$ of degree $n$ that approximates a ...

**0**

votes

**0**answers

36 views

### Probability of close approach for multivariate normal variables

The following problem comes from a physical model of two groups of particles in three dimensions. I need to know the probability that the two groups of particles approach each other within some ...

**2**

votes

**0**answers

70 views

### What is the computational complexity to compute the integral numerically?

Given $$\int_{\Delta}\frac{P_1(x_1,x_2,\dots,x_n)}{P_2(x_1,x_2,\dots,x_n)}$$ where $P_i$ is polynomial(that is $P_1(x_1,x_2,\dots,x_n), P_2(x_1,x_2,\dots,x_n)$ are polynomial) whose coefficients are ...

**1**

vote

**1**answer

69 views

### Condition Number and CFL Condition in Finite difference Methods

when applying a Finite Difference scheme for an IVP, two factors come to mind when considering stability:
One factor would be the condition number of the approximation operator. The other factor ...

**1**

vote

**0**answers

97 views

### Frozen coefficient method (von Neumann stability analysis)

Earlier it was considered that frozen coefficients method for Neumann stability analysis for finite difference scheme is more heuristic than rigorous. But I have read some information in a book by ...

**1**

vote

**2**answers

156 views

### iterative solution better than analytic solution? [closed]

My supervisor and I were discussing a specific optimisation problem this afternoon.
To be simple: solve for $R$ in the equation $Rx=y$, where $x$, $y$ are made of samples in two difference ...

**8**

votes

**1**answer

400 views

### Algorithm to produce random number with a gamma distribution

I'd like to produce pseudo-random numbers with different distributions for a Monte Carlo simulation.
I've got the poisson distribution working nicely with an algorithm from Knuth. I'm having trouble ...

**6**

votes

**1**answer

163 views

### Computing certain integrals over high-dimensional polyhedra

Let $\delta>0$ be a small real number and consider the $k$-dimensional region consisting of points for which
$$\delta\leq x_1\leq x_2\leq\ldots \leq x_k$$
and
$$x_1+\ldots+x_k\leq 1.$$
I am ...

**1**

vote

**0**answers

66 views

### General reparameterization of a b-spline

Say I have a bspline function (or curve) of order $k_1$, defined over some knot vector
$\mathbf{t} = \{ t_i\}_1^{n_1}$, i.e. $$f(x) = \sum_i a^i B_{i,k_1}(x).$$
Do you know of a process of finding ...

**3**

votes

**1**answer

75 views

### The link and equivalence between variant definition of computation model and computational complexity over reals

To unify the numerical computation and classic computability theory, or to pave a foundation for the numerical computation, mathematicians present variant computation model and computational ...

**6**

votes

**2**answers

219 views

### Conditions for convergence of Euler's method

It is known that a sufficient and necessary condition for
$$\dot y(t) = f(y(t), t), \quad t > 0, \quad y(0) = y_0$$
to have a unique solution is $f$ Lipschitz in $y$ and continuous in $t$. However, ...

**3**

votes

**1**answer

126 views

### Selecting Rays for Simulated Radon Transform

I have the task of determining approximations of a 2D function $f: (x,y)\in \mathbb{R}^2\mapsto\mathbb{R}$ from integrals along lines, i.e. from its Radon transform $R(\phi,\tau)[f(x,y)]$ and, because ...

**3**

votes

**0**answers

132 views

### Numerical Methods for stochastic PDE, from rough paths to backward equations

this question is about some literary references regarding the state of the art in terms of numerical methods for SPDE's. In particular,
Have the numerical implications, if any, of the results in ...

**2**

votes

**1**answer

142 views

### books on very large scale linear optimization

Recently in my material science research, I have encountered problems of very large scale linear optimization. I read the introductory book "Introduction to Linear Optimization (Athena Scientific ...

**0**

votes

**1**answer

75 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 ...

**1**

vote

**1**answer

46 views

### Partial Constraint of Low Rank Matrix

Suppose $X \in \mathbb{R}^{m \times n}$ is a rank $r$ matrix. Let $\Omega$ be a generic subset of $\in \{1, \ldots, m\} \times \{1, \ldots, n\}$ of cardinality $r(m + n - r)$. Denote by $X_{\Omega}$ ...

**4**

votes

**1**answer

164 views

### Special Function, Series Expansion, or Simpler Form of a Certain Infinite Product?

$\prod _{n=1}^{\infty } \left(1+a (c+n)^b\right)$ where a > 0, b < -1, and c >= 0
Is there a special function, series expansion, or other simpler (or maybe just interesting) representation of ...

**0**

votes

**1**answer

99 views

### Could somebody recomends a good book or article about numerical methods for Stochastic Partial Differential Equations

Could somebody recomend a good book or article about numerical methods for Stochastic Partial Differential Equations. I'm looking for a good introductory material thanks.

**0**

votes

**0**answers

41 views

### Numerical method for self-consistency of one-dimensional probability density function

I have an integral equation for self-consistency of one-dimensional probability density function, like this
$$\rho_x(x) = \frac{1}{|a|}\int \int \rho_x\left(\frac{s-b}{a}\right) \rho_P(p) ...

**0**

votes

**1**answer

257 views

### Improving Newton's Inequalities using the Taylor Theorem

Newton's inequalities say that if $f(x) = \sum \binom{n}{k} a_k x^k$ is a polynomial with all real roots then $ a_k^2 > a_{k-1}a_{k+1}$.
The proof this result uses that if $f(x)$ has all real ...

**2**

votes

**1**answer

51 views

### IVP accuracy - scheme accuracy Vs. derivative accuracy?

General Question: If I have an IVP with periodic and continuous initial condition, which rules the accuracy of the scheme - the manner in which we approximate spatial derivative or the acuuracy of the ...

**7**

votes

**1**answer

162 views

### Compute only selected components of an eigenvector

I am wondering whether it is possible to compute portions of the eigenvectors of a given (possibly very big) matrix. More formally, consider the eigenvalue problem $\mathbf{Ax} = \lambda \mathbf{x}$, ...

**7**

votes

**3**answers

428 views

### Multiprecision numerical evaluation of integral: Sage vs. PARI/GP vs. mpmath

I am trying to compute thousands of integrals of the below type, that comes up in a conformal mapping problem, to as many accurate digits as possible (preferably 50+):
$$
\int_{-1}^1\textrm{d}t ...

**7**

votes

**4**answers

481 views

### Numerical integration of legendre polynomials

I hope that numerical questions are also permitted here.
I want to expand a smooth functions $f \in C^{\infty}$in terms of Legendre polynomials. Thus I need to calculate integrals of the form ...

**2**

votes

**1**answer

66 views

### Monotonicity of Trapezoid Approximations

Here's a numerical analysis question which may not be very important, especially in practice, but has been bugging me.
Suppose $f$ is a continuous function on an interval $[a,b]$. Let $T_n(f)$ be ...

**-2**

votes

**1**answer

113 views

### reducing an n-order differential equation to a first order system of equations using either sagemath or sympy [closed]

I want to reduce a n-order ordinary differential equation into a first order system of equations. This is in preparation for numerical analysis. I use both Sympy and Sagemath for Computer Algebra, but ...

**0**

votes

**1**answer

185 views

### Proving that the eigenvalues of a certain matrix product are positive

Let $A$ be an $m \times n$ matrix, and define:
\begin{align*}
U &= {\rm diag} \{ \frac{1}{\beta_j} \}, \beta_j = \sum_{k=1}^m |a_{kj}|, j = 1 \dots n \\
V &= {\rm diag} \{ \frac{1}{\alpha_i} ...

**1**

vote

**1**answer

141 views

### Estimating the volume of a union of balls

Let $\{ B_i \}_{i=1}^n$ be a set of $n$ ball in the unit cube $C$ of dimension $d$.
If I want to estimate
$$
\frac{ \lambda \left( \cup B_i \right) }{\lambda\left( C \right) }, \tag{1}
$$
where ...