All Questions
1,732 questions
3
votes
2
answers
276
views
Practical symmetric equivalent to QR factorization updates
As we know, the QR-factorization $Q\cdot R=A$ of any real symmetric $n \times n$ matrix $A$ with full rank is unconditionally numerically stable. Further, when A is rank-1-updated, the factorization ...
- 51
5
votes
2
answers
352
views
Locating the maximum point $x_n$ of $f_n(x):=e^{-1/x}\Bigl(1+\frac{1}{n^2 x^n} \Bigr)$ in $(0,1)$
I am trying to observe the behavior of $x_n \in (0,1)$ defined such that the function
\begin{equation}
f_n(x):=e^{-1/x}\Bigl(1+\frac{1}{n^2 x^n} \Bigr)
\end{equation}
attains its maximum inside the ...
- 3,477
1
vote
0
answers
96
views
On optimizing a multivariate quadratic function subject to certain conditions
The problem is to maximize $f(x_1,x_2,\cdots,x_n)=\sum\limits_{i=1}^{n}\Big(x_i-k_i\Big)^2$ for $n\ge 3$ subject to the conditions (1) $\sum\limits_{i=1}^{n}x_i=\sum\limits_{i=1}^{n}k_i\le n(n-1)$ ...
- 61
4
votes
0
answers
714
views
What is the computational complexity of Arnoldi algorithm for diagonalization?
What is the space and time computational complexity of finding $k$ eigenvalues of an $N\times N$ matrix using the iterative Arnoldi algorithm?
I know that exact diagonalization scales like $O(N^3)$, ...
- 141
1
vote
0
answers
335
views
Closed-form solution of a particular linear program
(Note: I asked a similar question at math.stackexchange but the present one is more precise.)
I have a linear program of the form:
$$\text{minimize} \space\space x_1 \space\space \text{subject to:}$$
$...
- 990
2
votes
1
answer
84
views
Pressure integrated by parts in finite element method
Most FEM texts or tutorials apply FEMs on diffusion equations where the 2nd spatial derivative is integrated by parts during weak formulation. For convection diffusion equations, there is a first ...
- 159
6
votes
2
answers
539
views
Optimal polynomial approximation of rational function $\frac{1}{1-x}$
I've been working on the following polynomial approximation problem. I want to find the optimal Chebyshev approximation of the rational function $\frac{1}{1-x}$ on the real interval $x\in[-\rho, \rho]$...
- 63
1
vote
1
answer
131
views
Numerical methods for systems of trilinear polynomials
I have some large system of particular non-linear polynomial equations:
each equation mentions at most three variables
no variable appears with a degree larger than 1.
I'm not an expert in this area ...
- 121
1
vote
1
answer
169
views
Best projection on non-convex discrete set with two constraints
I want to compute the projection of a vector $\left( x\right) _{1\leq
i,j\leq n}\in \lbrack 0,1]^{n\times n}$ on the following discrete set
$$
S=\left\{ x\in \{0,1\}^{n\times n}:x_{i,j}+x_{j,i}\leq 1;\...
- 47
2
votes
1
answer
61
views
Counting the number of pair of d-uplets with upper bounded distance
Consider two d-uplets $u = (u_1,...,u_d)$ and $v = (v_1, ..., v_d)$ both living in $\mathbb{N}^d$ with $d$ a positive integer. They both verify $$(*) \sum_{i=1}^d u_i = \sum_{i=1}^d v_i = k$$ with $k$ ...
- 55
0
votes
1
answer
147
views
Is there a redundant constraint in linear programming? [closed]
From wikipedia:
But... Why do we need the $x\ge 0$ part? We can instead do $-x\le 0$, and thus saving a line in the definition (which is not a big deal but nevertheless nice).
(In order to do that, ...
- 3
1
vote
1
answer
155
views
Finding minimax approximation of a permutation equivariant polynomial
Is there any known method to approximate a given permutation-equivariant smooth function $f: \mathbb{R}^{n} \to \mathbb{R}^{n}$ as multivariable polynomial function $p: \mathbb{R}^{n} \to \mathbb{R}^{...
- 2,215
0
votes
0
answers
272
views
Finding the eigenvectors of a submatrix
Let $A=(a_{kl})$ be a matrix in $M_n(\mathbb{R})$ when $n$ is even. Let $B=(b_{kl})$ be the symmetric $2n$ by $2n$ matrix whose entries are given by,
$b_{k,l}=a_{kl}$ if $1\leq k,l\leq n$.
$b_{n+k,l}=...
- 4,058
1
vote
0
answers
47
views
Nash Equilibria change linearly in (some) game parameters. Already known / follows from a more general result?
EDIT: The key thing that I am wondering about is the linearity of the P2 strategy, not the constancy of P1. (The latter is straightforward.)
Question: Is the following result already known? Or is it a ...
0
votes
1
answer
83
views
Combining Dantzig-Wolfe and Benders decomposition
I'm now solving an LP that has a few coupling rows (as in Dantzig-Wolfe decomposition) and a few coupling columns (as in Benders decomposition) simultaneously; other rows and columns are block-angular....
- 3
4
votes
0
answers
99
views
Interpolation on Sobolev space on $[0, 1]^d$ over finite meshes
Let $\Omega = [0, 1]^d$ and suppose that $f \colon \Omega \to \mathbb{R}$ lies in order $m > d/2$ Sobolev space; i.e.,
$$
\|f\|_{H^m(\Omega)}^2 = \sum_{|\alpha| \leq m} \|D^\alpha f\|_{L^2(\Omega)}^...
- 380
1
vote
0
answers
56
views
Error bounds for Sobolev space norm approximation on a finite grid
Suppose that $f : [0, 1] \to \mathbb{R}$ is an element of the order-$k$ Sobolev space,
\begin{multline}
f^{(k - 1)}~\text{is absolutely continuous},\quad \|f\|_{W^k}^2 := \int_0^1 f^{(k)}(x)^2 \, dx &...
- 380
0
votes
1
answer
36
views
Benefit of adding a trivial constraint to ILPs
let ILP be an integer linear program with constraints-matrix $\boldsymbol{\mathrm{M}}\in\mathbb{Z}^{m\times n}$ and cost vector $\boldsymbol{\mathrm{c}}\in\mathbb{Z}^n$,
${\boldsymbol{\mathrm{x}}^*}\...
- 13.2k
0
votes
0
answers
94
views
Boolean operation on n dimensional polyhedron
A polyhedron in $R^n$ is defined by a set of half-planes: $P = \{x \in R^n \mid Ax - b \le 0\}$.
Given a set of polyhedra in $R^n$, $ P_1, P_2, \dotsc, P_k$, is there an algorithm/implementation that ...
- 21
0
votes
1
answer
93
views
How quickly can this IQP or its MILP relaxation be solved
Let $A\in\{0,1\}^{(n,n)}$ be a $n$ by $n$ boolean matrix (in particular think of an adjacency matrix of a graph), and consider the following optimization problem:
$$\begin{align*}&&\max_{P\in\{...
- 71
2
votes
0
answers
88
views
Kolmogorov $\epsilon$-entropy, $n$-width, and $\epsilon$-capacity and applications
What is the relationship between Kolmogorov $\epsilon$-entropy, Kolmogorov $n$-width, and Kolmogorov $\epsilon$-capacity of a set $M$ in a metric space $X$? (The $\epsilon$-capacity here is the ...
- 131
1
vote
0
answers
57
views
Uniform bound on the measure of $\Omega_\delta = \Omega \cap \delta\mathbb Z^d$ if $\Omega$ is an open bounded set with Lipschitz boundary
Let $\Omega \subset \mathbb R^d$ be an open bounded set with Lipschitz boundary. Let us consider $\Omega_\delta = \Omega \cap \delta\mathbb Z^d$ for $\delta >0$. I want to say that the measure of $\...
- 131
1
vote
0
answers
240
views
How to numerically compute the operator norm of an operator acting on a matrix algebra?
Let $M_n(C)$ denote the $n\times n$ matrices with complex entries acting on the Hilbert space $C^n$. As norm on $M_n(C)$ we take the operator norm, i.e. the largest eigenvalue of its absolute value. ...
- 11
8
votes
2
answers
1k
views
Mathematics of sustainable development and energy sobriety in the classroom
Faculty members are encouraged to highlight the connection between the courses we teach and climate change, and raise awareness of the issue in our lectures, across subjects in my university. I am ...
0
votes
1
answer
539
views
Method for (binary) optimization under constraints
I would like to know if there is a method to solve the Problem.
Problem:
Maximize the following function: $$f(p_{1,i},p_{2,i},\dotsc,p_{m,i})=\sum_{i=1}^{n}\begin{bmatrix}p_{1,i} & p_{2,i} & \...
- 3
3
votes
0
answers
182
views
Rate of uniform approximation by piecewise constant functions
Definitions and Notation:
Fix a positive constant $M>0$ with positive integers $m,n$ and the standard orthonormal basis $e_1,\dots,e_n$ of $\mathbb{R}^n$.
For every positive integer $N$, define the ...
- 5,405
3
votes
1
answer
162
views
Approximation in Bochner spaces
Is there any result like the Bramble-Hilbert lemma for Bochner spaces?
More specifically: let $H$ be a (e.g.) Hilbert space, $I\subset \mathbb R$ a bounded interval, and $L \in \mathcal L(H^k(I;H), Y)...
- 235
3
votes
0
answers
147
views
Chebyshev-like polynomials [closed]
In some approximation problems I'm working on, the errors turned out to be polynomials of various degrees whose graphs on the interval $[-1,1]$ look like this:
As you can see, these things look a bit ...
- 649
1
vote
0
answers
59
views
How do I incorporate Ito's lemma into the solution for a finite-horizon stochastic cake-eating problem?
I'm interested in finite-horizon, continuous-time cake-eating problems in which the agent has a time-horizon $W$ over which to eat the cake, and then chooses an optimal consumption path $\{h_t\}_0^W$, ...
- 11
0
votes
0
answers
63
views
Direct (first-order ?) algorithm to minimize $u(x) := \|x-a\|_C + r\|x\|_p$
Fix $a \in \mathbb R^n$, $r \ge 0$, $p \in \{1,2\}$, and a positive-definite matrix $C$ of order $n$. Define $u:\mathbb R^n \to \mathbb R$ by $u(x) := \|x-a\|_C + r\|x\|_p$, where $\|z\|_C := \sqrt{z^\...
- 6,853
5
votes
1
answer
284
views
Unbounded solution but bounded Euler discretization
Is there an ordinary differential equation in $\mathbb{R}^d$ induced by a gradient vector field with unbounded solutions, for which the difference equations obtained by using the forward Euler method ...
- 311
0
votes
1
answer
143
views
$\mathrm{ILP}$-formulation for Minimum Maximal Matching (MMM) Problem
Despite some online searching I couldn't find examples of dedicated Integer Linear Programs ($\mathrm{ILP}$s) for determining smallest matchings, that are not contained in a larger one.
It seems that ...
- 13.2k
0
votes
0
answers
54
views
Numerically expanding a function in a rational-power "basis"
I have some scientific code which interfaces with a library which accepts real functions specified as any number of additive terms with exponential powers. For instance, it is capable of accepting ...
- 173
11
votes
3
answers
726
views
Can computers find zeros of order $2$?
We assume we are given an entire function $f: \mathbb C \to \mathbb C$ with $f(0)=1$ and $f'(0)=0$ and $f$ is real on the real axis.
We assume (as a fact about $f$, that we want to demonstrate ...
- 124
1
vote
1
answer
181
views
Linear programming with "nice" matrices
Consider the following linear programming problem
\begin{array}{ll}
\text{minimize} & \mathrm 1^{\top} \mathrm x\\
\text{subject to} & v\le \mathrm A \mathrm x \le u\\
& \mathrm x \geq ...
- 650
1
vote
1
answer
703
views
Calculating the eigenvalues of the Laplacian numerically
I am trying to find the eigenvalues of the Laplacian operator, or in other words, solve the Helmholtz equation
$\nabla^2f=\lambda f$
on a compact 2D domain (comes from a quantum mechanics particle-in-...
- 3,738
1
vote
0
answers
59
views
Functional approximation with derivatives
I am trying to solve a functional approximation problem.
Consider a set of measurements of a d-dimensional state $\mathrm x \in \mathbb{R}^d$, together with velocities $\dot{\mathrm x}$ and ...
0
votes
0
answers
108
views
Solving a nonlinear equation maybe with Lambert W function
Can you please help me solve the following nonlinear equation?
\begin{equation}
\boldsymbol{z} \odot\left(\boldsymbol{\Gamma}^{\top} \boldsymbol{y}\right)=(\beta)^{\frac{1}{m-1}}\left(\frac{m-1}{...
- 11
3
votes
0
answers
69
views
how to efficiently find level sets (using a modification of a root-finding algorithm)?
I'm trying to find a set of points $\{ a_i | f(a_i) = c_i \}_{i=1}^k$ where $f$ and $\{ c_i \}_{i=1}^k$ are given in sorted order. All $c_i > 0$, $f$ is continuous and monotonically increasing, $f(...
- 131
6
votes
1
answer
913
views
Resultant of linear combinations of Chebyshev polynomials of the second kind
The Chebyshev polynomial $U_n(x)$ of the second kind is characterized by
$$
U_n(\cos\theta)=\frac{\sin(n+1)\theta}{\sin(\theta)}.
$$
It seems that
$$\operatorname*{Res}_x \left( U_n(x)+tU_{n-1}(x),\...
- 437
3
votes
1
answer
175
views
Is there a classical textbook/reference on numerical discretization schemes?
I found that it is relatively easy to find a book that discusses Euler discretization or Runge-Kutta discretization, but I am not aware of one that is well-known and/or common knowledge (i.e., field-...
- 479
2
votes
1
answer
177
views
Obtain 3D function from 2D slices [closed]
I am given a graph of a motor's torque vs RPM values at two different current draws, 730A and 300A (see graph). I need to obtain the 3D function to find current as a function of torque and RPM by ...
- 31
2
votes
1
answer
877
views
Interpreting mincost flow dual variables
Consider the task of finding flow of size $b$ with minimum possible cost.
It may be formulated as linear programming in a following way:
$$\boxed{\begin{gather}
\min\limits_{f_{ij} \in \mathbb R} &...
- 1,171
3
votes
0
answers
122
views
Preconditioners for $Ax=y$ that rely on hierarchical statistical modeling
Solving $Ax=y$ exactly can be done as:
fit a linear autoregressive model by treating rows of $A$ as data
apply this model to $A^T y$
Imperfect predictive model corresponds to an approximate inverse ...
- 2,442
1
vote
1
answer
59
views
Does norm of discrepancy decrease monotonously in CGLS/CGNR
I am the author of the package for tomographic reconstruction https://github.com/kulvait/KCT_cbct I have implemented CGLS/CGNR , algorithm which applies conjugate gradients on normal equation
$$
A^\...
- 151
0
votes
1
answer
64
views
Round Robin volleyball Tournament [closed]
Consider a set of N teams (N even number) that must make a
Round Robin Tournament. To each pair i; j, i ≠ j, of teams there is associated level
of interest si,j ∈ {1;2;3} of the match between them (1 =...
0
votes
1
answer
123
views
Explicit expression of Padé–Hermite approximant of type I
It is well known that the Padé approximants $(P,Q)$ of an analytic function in the neighborhood of $0$ can be expressed as a quotient of Hankel determinants built on the coefficients of the function $...
- 3,996
1
vote
0
answers
114
views
How is the quasipotential in Freidlin-Wentzell theory of large deviations affected by $C^1$ transformations?
I have a 3D differential equation I'm interested in studying the potential landscape using quasipotentials, described in depth in this paper. I need to calculate the potential landscape several times, ...
- 73
0
votes
0
answers
114
views
Degeneracies in linear combination of tensor product of Pauli matrices
Let $P_i \in \{I,X,Y,Z\}^{\otimes n} $, that is $P_i = \bigotimes_{i =1 }^n \sigma_i$ with $\sigma_i \in \{I,X,Y,Z\}$, where
$$
I = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \hspace{1cm} X =...
- 131
1
vote
0
answers
76
views
Newton-Raphson and bisection for solving for a system of nonlinear equations?
We know that for a single equation root-finding, we can use the Newton's method, or a combination of Newton with bisection to guarantee convergence. Can we use Newton+bisection for a system of ...
- 11