This tag concerns computational problems central to mathematical and scientific computating. The scope includes algorithms, numerical analysis, optimization, and linear algebra, computational topology, computational geometry, symbolic methods, and inverse problems.

learn more… | top users | synonyms

4
votes
0answers
66 views

On a unique(?) binomial property of $3003$

Given the triangular number, $$T_k = \frac{k(k+1)}{2}$$ and remembering that, $$\binom{n}{m}=\binom{n}{n-m}$$ Excluding $a_0=1$, we then have the six-fold (at least) equalities, $$\begin{aligned} ...
-2
votes
0answers
37 views

What is elnekiti's triangle? (edited) [on hold]

Elementary ceĺular automata shows amazing complex systems such as pascal's triangle is similar to " wolfram rule 90 " , so i looked over youtube searching for extra content and i found this video Here ...
0
votes
0answers
29 views

How to calculate $det(X^T X)$ efficiently, update one column of X each time

$X_{1} = (A, b)$, where $X_{1}$ is a $n\times p$ matrix, $A$ is a $n\times (p-1)$ and $b$ is $n\times1$. First calculate $\det(X_{1}^T X_{1})$, then update $b$ with $c$, st. $X_{2} = (A, c)$ and ...
2
votes
3answers
55 views

What is the value of this function?

Consider the three-variable function defined at the following way for all natural numbers $n,x,y$ : $f(0,x,y) = x+y $ $f(n,x,0) = x$ $f(n,x,y) = f(n-1, $ $ $ $f(n,x,y-1) , $ $ $ $f(n,x,y-1)+y ) $. ...
1
vote
2answers
35 views

(Fast way to) Get a combination given its position in (reverse-)lexicographic order

This question is the inverse of the Fast way to get a position of combination (without repetitions). Given all $\left(\!\!\binom{n}{k}\!\!\right)$ combinations without repetitions (in either ...
0
votes
0answers
48 views

Need to solve for an angle, do I need to use numerical methods??

I need to run a simulation on MATLAB where I need to solve two equations for two unknowns. The equations look something like this, $$x_{comp} = G\cdot \cos^2 b \cdot \cos(a+b).$$ $$z_{comp} = ...
1
vote
1answer
57 views

Expected Utility Method and a Repeated Game Solution

I am trying to replicate Bruce B. de Mesquita's (BDM) results on political game theory for prediction. Based on where actors stand on issues, their capabilities, salience, BDM's method attempts to ...
2
votes
3answers
101 views

Fast way to get a position of combination (without repetitions)

This question has an inverse: (Fast way to) Get a combination given its position in (reverse-)lexicographic order What would be the most efficient way to ...
0
votes
1answer
46 views

Constructive Induction to derive and prove the formula for a geometric sequence

$\sum_{i=1}^nr^i = a \cdot (b^n) + c$ Base Case: n=1, this holds Inductive Hypothesis: Assume for $n = k$, $k \ge 1$ that $\sum_{i=1}^k r^i = a * (b^k) + c$. Inductive Step: Prove for $n = k+1$ ...
0
votes
0answers
25 views

In bipartite directed graph, how to efficiently find nontrivial edge values such that for each node, sum of edge values is zero?

Suppose I have a set of sources U and a set of sinks V and a set of directed edges E going from elements of U to elements of V. I want to find a vector z of size |E| (each value of z is assigned to ...
7
votes
1answer
83 views

Are the unit partial quotients of $\pi, \log(2), \zeta(3) $ and other constants $all$ governed by $H=0.415\dots$?

Khinchin showed that given the simple continued fraction of a real number, $$r = a_0+\cfrac{1}{a_1+\cfrac{1}{a_2+\cfrac{1} {\ddots}}}$$ then it is almost always true that the partial quotients $a_i$ ...
0
votes
1answer
28 views

Seeking the Recommendation on Complexity Theory books

S.E advisers, I am a rising college junior in US with a major in mathematics and an aspiring applied mathematician in the fields of theoretical computing. I just recently got a research project on ...
-1
votes
0answers
158 views

need mathematical expression for the below [on hold]

Want to convert below algorithm into a mathematical model:- General points 1. Let there be a Connected Directed Graph. G = (V, E) V vertices or nodes E edges. This graph can be seen as a network ...
0
votes
0answers
14 views

How can I measure the dispersion of a discrete time system?

I have a discrete time system: x[n+1]=A*x[n]+B[n], with B[n] as an excitation signal. How can I measure the dispersion of this system? I guess this is related to some property of A, which depicts the ...
2
votes
2answers
85 views

Verify input is the sum of other numbers

I have a relationship: 4000k + 2500j + 400g = n, k >= 0, 0 <= j <= k, 0 <= g <= j I have to, given n, verify ...
2
votes
1answer
70 views

How do computers generate primes so quickly?

From what I understand, when a computer encrypts a file using an encryption standard like RSA, one of the steps is to generate two large primes, and multiply them together. I have created RSA keys on ...
6
votes
5answers
259 views

Computing as many digits as possible of $\sqrt{2}$ with a pen and a paper in 5 minutes

You have to compute as many digits as possible of $\sqrt{2}$ with a pen and a paper (an eraser if you're lucky...) in 5 minutes. What will you do? What is your justification for doing it? The ...
0
votes
0answers
23 views

Unknown variable in formula - binomial coefficient? [duplicate]

I'm currently researching a filter and don't quite understand one of the equations used there, since it contains a variable I don't know how to calculate: $$(1)\,\,a^{m, k}_s = \frac{c^{k, ...
0
votes
0answers
18 views

hyperoperation sequence with non-integer values of n

This probably has a very simple answer of some sort, but I'm not a mathematician. For the hyperoperation sequence: $$G(n,a,b)$$ ...there are obvious defined values for positive integer values of $n$ ...
0
votes
0answers
29 views

deconvolution of exp($x^2$)

I would like to know whether we can get the function of type exp($x^2$) by convoluting any functions. That is which function convolution gives exp($x^2$). Thanks in advance
1
vote
1answer
25 views

computing a fixed interest rate

I've been struggling for hours now with understanding a Topcoder problem, Autoloan , but i cannot grasp the way of computing it from a mathematical point of view. The excerpt goes as follows: A ...
0
votes
0answers
19 views

Analytical expressions for the orthogonalization of a specific set of vectors

I would like to know whether analytical or closed-form expressions could be obtained for the orthogonalization of a set of vectors in the following setting. Let $x_t$ be a vector indexed as a time ...
0
votes
1answer
17 views

Calculate height for rows of rectangles within a given width

I have an array of rectangles, all of the same height, but with different widths. Imagine they are on a single line with a uniform gap between them as shown below... |XXX| |X| |XXXXX| |XXX| ...
1
vote
2answers
61 views

Can car traffic be managed by mathematical formula? [closed]

How car traffic is managed? Is it managed by mathematical algorithm? Or by human(operator)? If it's by operator, can it be managed mathematically? Or is it by physics? By what theories/formula? ...
2
votes
1answer
24 views

Approximate recursively defined error in fixed point iteration

Problem: With an initial guess of $x_0$, the fixed point iteration is given by $$x_{k+1} = e^{-x_k}, \mbox{ for } k=0,1,2,...$$ If $x^*$ is the exact solution, then the approximation error is $$e_n = ...
1
vote
2answers
76 views

Is there any sort algorithm quicker than Quicksort given a random array of integers?

How can we proof (mathematically) that any complexity of sorting algorithm that sorts a random array of integers is no better than $O(n\log n)$?
1
vote
0answers
17 views

Is there a way to delineate the parameter of highest influence in a system of differential equations?

So I have a system of nonlinear ordinary differential equations dependent on parameters. These equations can traditionally be solved numerically with robust methods and the solution is well defined. ...
0
votes
0answers
6 views

State-of-Art library or a method for paralell matrix inversion?

Do you have reference to a computer library or a paper regarding a state-of-art method to obtain inverse of a matrix in parallel? Thanks, Mojo.
8
votes
0answers
74 views

Generalizing Bellard's “exotic” formula for $\pi$ to $m=11$

Bellard's "exotic" pi formula has the form, $$a\pi+b = \sum_{n=1}^\infty \dfrac{P(n)}{{\displaystyle \binom{mn}{2n}2^{n-1}}}$$ where $a,b,m$ are integers and he uses $m=7$. However, it seems there ...
0
votes
0answers
16 views

FFT differential equations

Given a generical differential equation what is the procedure to solve it using fft command. Can anyone explain me how to do it? For example: $$\frac{d^2y}{dt^2}+10\cdot \frac{d\:y}{dt}=-5\cdot ...
0
votes
0answers
9 views

How to solve the Helmholtz equation in a triangular region?

Suppose we take the Dirichlet boundary condition, namely the function must vanish on the boundary of the triangle. How about a general n-polygon?
-1
votes
1answer
28 views

What is transpose multiplier and forward multiplier?

For linear system X = A*s, we define the forward and transpose multiplies Af and At as follows: Af = @(s) A*s; At = @(s) A'*s; I want to know what is forward ...
1
vote
0answers
36 views

Simplifying the Generalized Eigenvalue Problem

Let $\Sigma_1$, $\Sigma_2$ be symmetric positive-definite real $n\times n$ matrices. We want to solve the generalized eigenvalue problem $$ \Sigma_1V=\Lambda\Sigma_2V, $$ where $\Lambda$ is the ...
0
votes
0answers
18 views

Using Orthogonal Collocation to solve Coupled Ordinary Differential Equations

I am trying to solve six first order coupled ODE's, two of these are associated with a heat balance of a catalyst pellet, and four are mass balances. I have been trying to solve these equations using ...
3
votes
1answer
53 views

Prime numbers distribution theorem

I'm trying to understand Gauss' theorem: $$ \frac{\pi(x) }{x/\ln x} \to 1 $$ for large $x$. I've taken the list of first 1000 prime numbers from Utah university site, saved them to file ...
1
vote
0answers
21 views

periodic boundary conditions and the FEM

I am trying to set up the mass matrix for a 1D system which I want to solve using finite elements. So the mass matrix is defined as $$ M = \int{NN^T}dL, $$ where $N$ is the finite element linear ...
0
votes
0answers
5 views

What is the behavior of the spatial median in high-dimensional spaces?

I am a photographer who is investigating a technique known as image stacking, in which multiple images of the same subject are combined to reduce noise (by CLT). Commonly used techniques are mean and ...
1
vote
1answer
24 views

How to make Poisson voronoi diagram

I am facing a problem as follows : I want to make poisson voronoi diagram & for this I have to appropriately choose some generating points. Sources on internet(WIKI) refers that these points have ...
1
vote
0answers
20 views

Calculate pairwise cosine distance only returning the lower triangular matrix

I have a matrix, where each row is a feature vector. I would like to find out the pairwise cosine distance between all of these feature vectors. The cosine value between all rows in a matrix could be ...
0
votes
2answers
21 views

Special properties in the direct solving of sparse symmetric linear systems

In the area of computational solving of large sparse linear systems, some solvers specialize only on symmetric sparse matrices, be it positive definite or indefinite as compared to general ...
1
vote
1answer
25 views

Computational complexity of the algorithm

Make an analysis of the computational complexity of the algorithm below, where it is given by the number of elementary operations that the algorithm performs (assignment is not considered). Where ...
1
vote
1answer
48 views

inverse of $AQ^{-1}A'$

Suppose that $A$ is a $m\times n$ full row rank sparse matrix, and $Q$ is an $n\times n$ symmetric positive definite sparse matrix with $m<n$. Besides, $m$ is about $10^5$, and $n$ is about $10^6$. ...
6
votes
2answers
86 views

How to determine whether a point is inside a closed region or not?

Take the following parametric equation of an implicit curve as an example: $$ \left\{\quad \begin{array}{rl} x=& 9 \sin 2 t+5 \sin 3 t \\ y=& 9 \cos 2 t-5 \cos 3 t \\ \end{array} \right. $$ ...
1
vote
0answers
34 views

Cholesky decomposition and rotation matrix inverse

I implemented three methods for inversion of a matrix, all are classic. I wanted to test for the most generalized method, while taking efficiency into account. For Cholesky decomposition, which is ...
4
votes
1answer
39 views

Compute finite series

The problem is to count the sum of the finite series $$\sum_{k=0}^{k_0} \frac{a_k}{b_k}$$ I need to count this series in binary with some precision, that would output $n$ correct binary digits after ...
0
votes
0answers
12 views

Literature Reference for transformations through vector spaces

I am trying to understand the transformations through vector spaces: Problem 1. Let's say we have orthonormal basis $B=\{v_1, v_2, \ldots, v_n\}$ spanning the vector space $V$ and basis $B_1=\{u_1, ...
0
votes
0answers
36 views

Algorithm for finding zero of an odd function from n-sphere -> R^n

There is a well-known Borsuk-Ulam theorem stating that each continuous mapping $f : S^n \rightarrow \mathbb{R}^n$ that is odd in sence of $f(v) = -f(-v)$ for each $v \in S^n$ (where $-v$ denotes the ...
0
votes
1answer
24 views

Analysis of iterative optimization methods using lyapunov analysis

In analysis of iterative methods, is it possible that we have to use two time-lagged version of the time-varying system to analyze its convergence? (that is, we construct the evolution of x^k, ...
1
vote
0answers
21 views

Why steepest descent gives a wrong direction search?

I have to minimize the function $ƒ(x_1,x_2)=(x_1-1)^2+x_2^3-x_1x_2$. The initial point is $[1,1]^T$. The gradient of this function is $∇ƒ(x_1,x_2)=[2(x_1-1)-x_2,3x_2^2-x1]$. This gradient evaluated ...
8
votes
1answer
91 views

How to find this number, which is probably a very big prime or a product of big primes?

Let $\mathcal{N}(n)$ be the next prime greater than $n$. Which is the smallest natural number $n>0\;$ such that: $\mathcal N(2\cdot 3\cdot 5\cdot 7\cdot 11\cdot n)−2\cdot 3\cdot 5\cdot 7\cdot ...