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

11
votes
0answers
106 views

On the prime-generating polynomial $m^2+m+234505015943235329417$

In 2009, J. Waldvogel and Peter Leikauf found the remarkable Euler-like polynomial, $$F(m)=m^2+m+234505015943235329417$$ which is prime for $m=0\to20$, but composite for $m=21$. Define, ...
11
votes
0answers
145 views

Product of primes mod n

Let $n$ be an odd composite number. I'm trying to compute $$ f(n)=\prod_{n/2<p<n}p\pmod n $$ where $p$ ranges over the primes in the indicated region. Can this be done (significantly) faster ...
8
votes
0answers
91 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 ...
5
votes
0answers
85 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} ...
5
votes
0answers
104 views

A summation involving the ceiling function

I'm trying to find a better method of calculating the sum $$\sum_{k=1}^N\lceil ak\rceil^2$$ where $a$ is an irrational number. So far, my only idea is to somehow use a best rational approximation. ...
5
votes
0answers
181 views

A property of the subgroups lattices

Let $G$ be a finite group. Consider all the subgroups $H$ such that its subgroups lattice $\mathcal{L}(H)$ is distributive (i.e. the group $H$ is cyclic, by Ore's theorem), and among them, let $\{ ...
4
votes
0answers
69 views

Finding every $n$ such that $n\times$ ('reverse' number of $n)=m^2$ such as $1584\times 4851={2772}^2$

Let $r(n)$ be the 'reverse' number of $n$ in the decimal system. For example, $r(1234)=4321$. Then, here is my question. Question : Can we find every $n(\in\mathbb N)$, which is not a square ...
4
votes
0answers
537 views

Obtain a contradiction

Motivation : The motivation is to show that the equation $x^{2b}.x^{2a} +(3-x^{2b}) x^{a} + (1-s^2)=0 $ has no solutions in integers for any values of $x,b,a,s$ ( choosen as per the constraints ...
4
votes
0answers
82 views

is there a computationally efficient formula for computing the mutual information between two continuous variables?

I need to compute the mutual information between two continuous variables. Below is an equation shown to compute the mutual information between a variable $X$ and $Y$. $I(X;Y) = \int_Y \int_X ...
3
votes
0answers
53 views

Is there a systematic way of “discovering” an algebra from observations of its universe?

I am faced with the following situation: I have a finite set of some $m$ positive integers $Q^m \in \mathbb{N}$ These integers go through a series of $N$ possible black boxes that transform them. ...
3
votes
0answers
118 views

Computing the density of a set of multiples

I have a finite* set $A=a_1<\cdots<a_r$ of positive integers. Define $B$ as the set of positive integer multiples of $A$ and $$ A_1=\frac{1}{a_1} $$ $$ ...
3
votes
0answers
67 views

A Free Boundary Problem

Is there any special way to solve such a problem. Any idea would be appreciated. At least does anybody know which method is useful to solve this problem numerically? Is it even solvable numerically? ...
3
votes
0answers
26 views

Has there been work on computational group theory applications to computing colimits of crosses n-cubes of groups?

I'm trying to compute homotopy groups of a few spaces using crossed n-cubes of groups. I'm able to describe a few colimits in terms of quotients of induced crossed modules and nonabelian tensor ...
3
votes
0answers
411 views

Nerve Theorem: Is the finite union of closed convex sets triangulable?

My Question: Let $A_1, \ldots, A_k \subseteq \mathbb{R}^n$ be closed convex sets. Is the union $\bigcup_{i=1}^k A_i$ triangulable$^1$? If so, why? Background: I'm trying to better understand the ...
3
votes
0answers
83 views

Computing the decomposition of a representation of $S_n$

I have an explicitly defined representation of the symmetric group that I would like to decompose into irreducibles. How to do this most easily? The best approach I have so far is as follows: Find a ...
3
votes
0answers
64 views

Algorithm for primary decomposition of ideals in a power series ring over a field

Let $K$ be a field such that there exists an algorithm for factoring a polynomial over $K$ into the product of irreducible polynomials. For example, the field of rational numbers $\mathbb{Q}$ is such ...
3
votes
0answers
60 views

Supremum over unitary group action

Let $A$ and $B$ are two given Hermitian operators on matrix algebra $M_n(\mathbb{C})$. $A$ is positive semi-definite with unit trace. I want to know the general method for calculating the following ...
3
votes
0answers
402 views

“Green Globs” question

When I was in high school geometry, we had a fun little game on the computer called Green Globs (the website for the software is http://www.greenglobs.net/index.html). A number of targets (globs) are ...
2
votes
0answers
18 views

Given a set of integers $S$, what is the maximum integer that is a product of one or more integers from $S$ not exceeding $X$?

Given a set of integers $S$, which will contain no more than $100$ integers. Now, what would be the fastest approach to find $M$ which is a product of one or more integers from $S$ (and multiple usage ...
2
votes
0answers
45 views

A recursive problem with GAP concerning lists and an iterator loop

I have the following question concerning a list algorithm in GAP: Let $L_1$ be a non-empty list with certain objects as entries. I wrote a program and called it helping_program_1. The Input for ...
2
votes
0answers
21 views

terms of taylor expansions of multiple variables at the origin

By the fundamental theorem of symmetric polynomials, $X_1,X_2,\cdots,X_n$ are polynomials of $ e_1,\cdots,e_n$ and $$ \mathbb{Z}[ e_1,\cdots,e_n]=\mathbb{Z}[X_1,X_2,\cdots,X_n]. $$ We define a ...
2
votes
0answers
54 views

Summation of n terms of a series$ 1+2+8+64+…$

In one of my problem , I got a series as $1$,$2$,$8$,$64$,$1024$...and so on. can we really get a sum expression for that series$???$ If yes, then what is the expression $f$ $?$ or the sum of $n$ ...
2
votes
0answers
73 views

I have a finite permutation group and access to a computer algebra system, how can I recognize the structure of the group?

I have a fixed permutation group $G$, and cannot tell which finite group it really is in a "human readable" way. I also have GAP. Is there a step by step computation to give me the structure of this ...
2
votes
0answers
17 views

How to compute Generalized Group Inverse?

Given a transition matrix $P \in \mathbb{R}^{n \times n}$, i.e. $\sum_j P_{ij} = 1$ and $P_{ij} \geq 0$ for all $i,j$. One can show that there exists a unique group inverse $B$ of $A:= I - P$ which ...
2
votes
0answers
34 views

Good method for finding roots that *usually* fall within an interval?

I've been using Brent's method to find the roots of a monotonic, nonlinear, non-differentiable function. The roots often fall within a known interval, but Brent's method fails if they occasionally ...
2
votes
0answers
46 views

Computer program to simplify formulas

What is the computer program that attempts to simplify sums of binomial coefficients, factorials, etc.? Possibly Zeilberger wrote it, but I'm unsure. If so, possibly it was talked about in his A=B ...
2
votes
0answers
40 views

Quantitatively comparing event trains of different lengths for Poissonness

I have a parameterized, effectively black box process that generates a series of events (simulated action potentials). Different parameter values often lead to different numbers of events. How can I ...
2
votes
0answers
481 views

Two-layer Perceptron for XOR

I'm reading Neural Networks for Pattern Recognition by Christopher M. Bishop. It's for a physics class, but I think the problem is closer to mathematics so I'm asking here instead of PSE. Chapter 4 of ...
2
votes
0answers
55 views

LLL and factoring polynomials in $\Bbb Z[x]$

Given a degree $2k$ reducible polynomial $f(x)=\sum_{i=0}^{2k}a_ix^i\in\Bbb Z[x]$ with $gcd(a_{2k},\dots,a_0)=1$ that is known to be of the form $f_1(x)f_2(x)$ with $deg(f_i(x))=\frac{deg(f(x)}{2}=k$ ...
2
votes
0answers
172 views

Using GAP to compute the abelianization of a subgroup

Let $K$ be the group generated by four elements $x_1,\cdots,x_4$ with relations that each generator commutes with all its conjugates. (An equivalent relation is, any simple commutator with repeated ...
2
votes
0answers
661 views

Solving large, sparse system of linear equations

I have a system of linear equations as follows: $$(A+I)x=B$$ where $I$ is the $n\times n$ identity matrix, $A$ is a $n\times n$ matrix such that the first and last rows are blank, and, for every ...
2
votes
0answers
123 views

Solving a particular system of Diophantine equations in $n$ variables (Frobenius equations)

I have a particular system of linear Diophantine equations in $n$ variables for which I need to find all nonnegative integer solutions. Specifically, they are Frobenius equations, meaning the ...
2
votes
0answers
36 views

Computing a particular finite set of quaternion matrices.

Let $B = \left(\frac{-1,-11}{\mathbb{Q}}\right)$ be a choice of quaternion algebra ramifying at $11$ and consider the maximal order ...
2
votes
0answers
379 views

How to convert a hologram into an image?

Suppose one knows in full detail the phase and intensity of monochromatic light in a plane. This is basically what a hologram records, at least for some section of a plane. By using this as the ...
2
votes
0answers
488 views

Logistic regression algorithm in Casio and Texas Instruments calculators

When using logistic regression on a Casio or Texas Instruments calculator, the output is of the form $$f(x) = \frac{c}{1+ae^{-bx}} $$ The problem I have (when teaching in a class where both types of ...
2
votes
0answers
467 views

Clarification on different types of solutions: analytical, closed form, iterative, algorithmic,

When doing computation, there are different types of solutions: analytical, closed form, iterative, algorithmic, exact, approximate, ... (there are probably more than I just listed and please don't ...
1
vote
0answers
22 views

Trying to learn about kernel PCA but cannot understand some math.

I'm trying to learn about kernel PCA by reading through the paper of it's creators (I assume) "Nonlinear Component Analysis as a Kernel Eigenvalue Problem", Bernhard Schölkopf, Alexander Smola, ...
1
vote
0answers
17 views

Integrating Wishart density

I have several points $\textbf{s} = s_1,...,s_n$ which follow Wishart distribution. In one of my problem, I have to integrate this Wishart pdf over a ball of radius $r$ at origin in $\mathbf{R}^2$ ...
1
vote
0answers
10 views

Efficiently compute the Determinant of a Banded Matrix

So I've got a large (~ 2 million x 2 million) positive semi-definite, banded, square matrix that I need to find the determinant of. What is the correct way to efficiently compute the determinant of a ...
1
vote
0answers
17 views

I need help with this Standard Uniform random variable problem

Click here for problem image I need help on how to do b,c and d. Can I please have an idea of how to start it because i have no idea how to do this.
1
vote
0answers
31 views

for what value of p does the integral converge

For what value of p does $\int_0 ^1 x^p dx$ converge? So I used improper integral value for this, looking at the result $$\lim_{p\to 0+} {1\over p+1}(x^{p+1}-p^{p+1})$$ *** I made a mistake, I got ...
1
vote
0answers
54 views

Find the different values of this integral when all paths of integration are possible

This is the question.. I only know how to do the question from infinity to infinity.. enter image description here Find the different values of this integral when all paths of integration are ...
1
vote
0answers
42 views

Constructing triangulations algorithmically

I am developing a Python package for computations in algebraic topology (namely: cohomology and Massey products on manifolds). Basically all the stuff I'm doing requires an explicit triangulation of ...
1
vote
0answers
20 views

Algorithm detect simple curves using Voronoi diagram or Delaunay triangulation?

I wonder if there is algorithm/method to determine if closed (or even non closed) curve is simple or not, using the mathematics from the field of computational geometry? Especially I wonder if exist ...
1
vote
0answers
37 views

Algorithm for min max of a function

I want to solve a problem of the form $$ \min_{x \in X} \max_{y \in Y} f(x,y) $$ where $X \subset \mathbb{R}^d$ and $Y \subset \mathbb{R}^k$. Typically $d$ is large ($>100$) and $k$ is small ...
1
vote
0answers
19 views

Minimum vertex cover of vertex disjoint odd holes and antiholes

I am interested in knowing whether the minimum vertex cover of a graph that can be written as the union of vertex-disjoint odd holes and odd antiholes can be found exactly, in polynomial time. I could ...
1
vote
0answers
41 views

How can we calculate the tensor product of Lagrange basis polynomials?

Given data points $(x_i,y_i)$, the Lagrange basis polynomials are $$\mathcal l_j(x):=\sum_{i\ne j}\frac{x-x_i}{x_j-x_i}\;.$$ I'm reading a text targeting Smolyak's algorithm. In this text, they use ...
1
vote
0answers
18 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. ...
1
vote
0answers
54 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 ...
1
vote
0answers
42 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 ...