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

2
votes
1answer
56 views

The “computability” of fundamental physical constants

I would like to ask if any of the fundamental physical quantities like the speed of light or plancks constant (all measured according to a common standard of of units) can be classified as computable ...
0
votes
1answer
87 views

Discretize an ellipsoid given its semi-major axes and orientation

An ellipsoid centered at the origin can be defined by the solutions to $$ \mathbf{x}^\text{T} A \mathbf{x} = 1 $$ where $A$ is symmetric and positive-definite. The eigenvectors of $A$ define the ...
2
votes
1answer
46 views

Generating Eulerian digraphs/isographs

I would like to be able to quickly generate (all) non-isomorphic isographs (that is, digraphs where each node has the same indegree and outdegree - also called "balanced networks" in the distributed ...
0
votes
2answers
79 views

logarithm and exponent computation performance

Using glibc on a x86 processor, which takes more CPU time? $a\ log\ b$ or $b^a$? For which values of $a$ is one faster than the other? Optional: Does the base used matter? See also: What algorithm ...
0
votes
1answer
277 views

Calculating eigenvectors and eigenvalues of a 2x2 complex matrix

I've previously asked elsewhere, http://stackoverflow.com/questions/21118820/non-trivial-eigenvectors-of-a-22-matrix-in-code, how to calculate the eigenvectors and eigenvalues of a 2x2 matrix in a ...
1
vote
1answer
55 views

Is it possible to reduce a lambda expression to it's smallest equivalent form?

In the Untyped Lambda Calculus, is it possible to reduce any arbitrary expression to it's smallest equivalent form? (defined as an expression with the smallest number of lambda terms) If so, is there ...
2
votes
2answers
199 views

Galois group command for Magma online calculator?

I need to test if a family of 7th deg and 13 deg equations are solvable. I'm new to Magma, so my apologies, but what would I type in, http://magma.maths.usyd.edu.au/calc/ to determine the Galois ...
3
votes
0answers
70 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
51 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 ...
5
votes
3answers
288 views

FFT with powers of 3

Classic Fast Fourier Transfrom (FFT) works fine, when $n$ is power of 2. How to generalize FFT procedure when $n$ is power of 3? Is it possible to easily modify the algorithm and preserve its ...
1
vote
2answers
36 views

Subtlety about the definition of B-splines

I came across the following definition for the zero'th order B-spline $$b_0(x) = \left\{ \begin{array}{lr} 0 & |x|>1/2\\ 1 & |x|<1/2\\ 1/2& |x|=1/2. \end{array} ...
2
votes
0answers
62 views

How calculators compute. [duplicate]

I would like to teach a class on the "magic" behind the calculator, so I would like to generate a list of "algorithms" for how a calculator computes the things we want it to. I will get the ball ...
3
votes
2answers
167 views

Calling GAP's StructureDescription from SAGE

Given a sage group object $G$, I want to obtain its structure description using GAP, as follows: sage: gap.StructureDescription(G) The command works fine in ...
2
votes
3answers
69 views

Satisfying equality between logarithmic expressions

Apologies in advance for any misused terminology, or if this is the wrong place for the question (I think it's okay though). I am given a group of logarithmic expressions such as: $- (a \log(a) + ...
1
vote
0answers
54 views

What exactly is 'computer mathematics'?

I'm looking at some potential things to study next semester and I see a full B.sc. degree called 'Computer mathematics'. It says it's a hybrid between computer-science and mathematics. Does anyone ...
0
votes
1answer
65 views

Which strings belong to the regular set represented by the regular expression (1∗01∗0)?

I know the string should be like 1…101…10, but not sure how to describe it. Can anyone help me?
3
votes
1answer
338 views

Software for numerical solution of a non-linear ODE system?

I have been given a nonlinear system of ODEs which has arisen out of a colleague's engineering research: $$\begin{array}{rcl} \dot{x}_0&=&x_1\\ ...
1
vote
0answers
194 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
46 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$ ...
1
vote
3answers
49 views

Linear algebra computable

I'm working in a linear algebra software, but I have a problem. If I try to reduce the matrix to reduced echelon form, many resources will be consumed. That is, if I perform many basic operations ...
1
vote
2answers
348 views

Minimum number of iterations in Newton's method to find a square root

I am writing an algorithm that evaluates the square root of a positive real number $y$. To do this I am using the Newton-Raphton method to approximate the roots to $f(x)=x^2-y$. The $n^{th}$ iteration ...
0
votes
2answers
149 views

Inverse Percentage.

Sorry for asking this foolish question. Here is the data i have. I purchased the product as $5 and additional fee is 2%. So Here is the total dollor ...
4
votes
0answers
64 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 ...
3
votes
3answers
80 views

A good way to find $a_{50000}$ where $a_n$ is a number in the form of $2^j\cdot 3^k$

Letting $A=\{2^j\cdot 3^k| j,k \ \text{are non-negative integers} \}$, let us define $a_n$ as the $n$-th element of $A$ in ascending order. We can see $$a_1=1, ...
0
votes
1answer
36 views

computing and algebra

I have a question: Let $f(x)=\sqrt{x^2 +1} - 1$. When $x=10^{-3}$ compute $f(x)$ working to 5 sf. Show algebriacally $f(x)=\frac{x^2}{\sqrt{x^2+1}+1}.$ After desperately rearranging I'm just going ...
1
vote
0answers
91 views

Lamé's equation - elliptic PDE?

I have a specific PDE in 3-dimensional space + time (the right one). u(x,t) is the unknown function (values are in R^3) and F(x,t) the right hand side, mu and lambda are positive constants. Now ...
1
vote
0answers
28 views

Nonlinear optimization using parallel input/output

I have a system that accepts a vector and returns a function value. The goal is to change the elements of the vector such that the function value is minimized using a derivative-free solver, eg. using ...
2
votes
0answers
57 views

The Command “TzGoGo” in GAP

I am learning GAP and would like to ask one question about a command called "TzGoGo": If $P$ is a finite presentation of a group $G$, then will the eventual result of the command "TzGoGo(P)" be ...
3
votes
1answer
95 views

GAP Responding Time

I am running GAP 4.6.5 on my six-year-old computer and sometimes it takes like forever for GAP to respond to my simple commands. An easy example is as follows: ...
2
votes
1answer
91 views

Need help for a beginner to floating-point arithmetic

I have this question that I need to complete and I literally have zero idea on what to do. I basically need someone to talk me through it and will appreciate all the help I can get. The question ...
1
vote
1answer
49 views

Survey/encyclopaedia/website of mathematical theorems connected

Is there, or is someone creating a survey/encyclopaedia/website of mathematical theorems which connects theorems together with their assumptions (axioms, other theorems, hypotheses etc.)? I'm thinking ...
0
votes
1answer
131 views

Multi layer perceptron activation function

How can you show that the Fourier series approximation of a function (so $f(x)=\sum\limits_{n=0}^{\infty} (a_n cos(nx) + b_n sin(nx))$ can be approximated to arbitrary precision by a feedforward ...
1
vote
1answer
74 views

Ackermann function and primitive recursiveness

If we define $b_n(m) := a(n,m)$ for all $n$ and $m \in \mathbb{N}$. For which $n$ is the function $b_n$ primitive recursive and for which $n$ it is not a primitive recursive function? Can anyone ...
2
votes
0answers
118 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 ...
4
votes
1answer
138 views

What is the upper bound on the error of a matrix multiplication

When both A and B are n x n upper-triangular matrices, the entries of C = AB are defined as follows: $$ c_{ij} = \begin{cases} \sum _{k=i}^ja_{ik}b_{kj} & 1\leq i\leq j\leq n \\0 & 1\leq j\lt ...
0
votes
1answer
212 views

Simultaneous recursion

I have no idea how to even start proving the following theorem: If $f_0, f_1: \mathbb{N}^r \rightarrow \mathbb{N}$ and $g_0, g_1: \mathbb{N}^{r+3} \rightarrow \mathbb{N}$ are primitive recursive, ...
3
votes
1answer
113 views

Algorithm for Finitely Presented Torsion-Free Nilpotent Groups

I am studying some finitely presented, torsion-free and nilpotent groups $G$ and need to consider the following question: Let $H$ be a subgroup of $G$ and suppose that $H$ is generated by ...
2
votes
3answers
105 views

Websites/Software for Group computation

Anyone knows a website or software that helps to do computations in a group? For example, by inputting generators and relations in the group, can we tell when two particular elements in the group ...
2
votes
1answer
92 views

Is WolframAlpha computing this radical correctly?

Is WolframAlpha computing this radical correctly? $$\sqrt{\frac{1}{1 + {10}^{-375}}}$$ When I double-check again, the inequality: $$\sqrt{\frac{1}{1 + {10}^{-x}}} > 1$$ leads to a ...
3
votes
0answers
137 views

$\pi$, disjunctive numbers, and finite sequences of given length

It is an open problem whether the number $\pi$ is disjunctive in base $10$, i.e., whether every finite sequence appears (at least once) in the base $10$ expansion of $\pi$. Of course, every sequence ...
2
votes
0answers
436 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 ...
5
votes
2answers
1k views

Fast Matlab Code for hypergeometric function $_2F_1$

I am looking for a good numerical algorithm to evaluate the hypergeometric function $_2F_1$ in Matlab (hypergeom in Matlab is very slow). I looked across the ...
0
votes
2answers
68 views

Generating a random number of higher range

I had a discussion with my friend about writing a function (a computer program) which generates 9 values randomly using a random number generator which generates 4 values, i.e I have a PRNG rand(4) ...
3
votes
0answers
51 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 ...
4
votes
1answer
423 views

Wolfram Alpha error?

I was seeing some equations in WA, and i got with http://www.wolframalpha.com/input/?i=%28k%2B1%29%5E2%3E%3D4%28k-1%29%5E2 Let's manually solve the equation $$(k+1)^2\ge4(k-1)^2$$ ...
12
votes
1answer
631 views

Why are there mathematicians that do not use computers?

I was watching a video on Andrew Wiles and his proof of Fermat's Last Theorem and I quite liked the video, especially the complexity of the proof only to prove a simple concept which can be understood ...
2
votes
1answer
149 views

Is there a computer programm or CAS (maybe GAP?) that can calculate with projective (indecomposable) A-modules (A is a finite dimensional k-algebra)?

I have the following question(s): I have an "Algebra-With-One" $R$ as a subalgebra of a full matrix algebra in GAP. Furthermore, I have 5 primitive orthogonal idempotents $e_1,...,e_5$, which sum up ...
6
votes
1answer
198 views

how to calculate this complementary Bessel function?

I am trying to calculate this complementary Bessel function $$\Psi(a,b,\gamma)=\int_0^\infty\Phi({a\over \sqrt{u}}+b\sqrt{u}){u^{\gamma-1}e^{-u}\over \Gamma(\gamma)}du$$ where $\Phi$ is the standard ...
2
votes
0answers
89 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 ...
1
vote
1answer
66 views

Is there a good strategy for computing eigenspace corresponding to $1$ of a matrix with entries of trigonom

For example, say $A= \left ( \begin{matrix} \cos x & -\sin x & 0 \\ \cos y \sin x & \cos x \cos y & -\sin y \\ \sin x \sin y & \sin y \cos x & \cos y \end{matrix} \right)$. ...