Skip to main content

Questions tagged [implementation-details]

For questions regarding the underlying implementation of built-in functions.

Filter by
Sorted by
Tagged with
1 vote
1 answer
101 views

On `ContinuedFraction`

The code p = Integrate[E^-x^2, {x, 0, 1}] p // ContinuedFraction produces Does this mean that (i) Mathematica knows that a terminating or periodic continued ...
Iosif Pinelis's user avatar
5 votes
1 answer
132 views

Understanding NDSolve's implementation of the method of lines

In order to better understand the MethodOfLines flag, I tried comparing the results of NDSolve for a driven heat equation with a ...
phonon's user avatar
  • 105
1 vote
1 answer
131 views

A request on the kernel of the function ”SpheroidalEigenvalue“ in Mathemitica 11

I want to understand the algorithm of the function "SpheroidalEigenvlaue" in Mathematica 11, especially the code of the following command. Then I will try to use this algorithm to reproduce ...
amon xu's user avatar
  • 41
0 votes
0 answers
98 views

How is the analytical continuation for the HurwitzZeta function implemented?

Following up on this question, I am trying to understand the implementation details of the HurwitzZeta[x,y] function in Mathematica, particularly when the first ...
stefan_chem's user avatar
3 votes
1 answer
151 views

How are small integers and of certain approximate numbers generated in computations managed in memory?

In the documentation tutorial/SomeNotesOnInternalImplementation, it says: Most raw objects such as strings and numbers are allocated separately; however, unique ...
rnotlnglgq's user avatar
  • 3,740
5 votes
2 answers
281 views

What do the arguments of a generated BSplineFunction mean?

For example, if we evaluate this: BSplineFunction@{{0,100},{200,50},{200,0}} we'll get ...
rnotlnglgq's user avatar
  • 3,740
12 votes
3 answers
811 views

FEM doesn't output exactly the same solution for exactly the same code?

Consider this toy example: ...
xzczd's user avatar
  • 66.8k
13 votes
2 answers
332 views

How is the zoom tooltip in the "Coordinates Tool" implemented?

I notice that the built-in "Coordinates Tool" has a very efficient zoom tooltip displaying an enlarged portion of the image along with coordinates and row/column indices: How to get the ...
Alexey Popkov's user avatar
4 votes
2 answers
181 views

Interior point method for unconstrained problems

I was recently working on applying the FindMinimum function to $f(x)=(x+1)(x-1)$, and chose the interior point method as the 'method', which yielded $x=0$. I would like to better understand how the ...
manofwar's user avatar
24 votes
2 answers
2k views

What happens if you don't Reap what you've Sown?

As I understand it, when using Sow[expr] you throw the expr on some private stack which you can ...
Gert's user avatar
  • 1,530
4 votes
1 answer
218 views

Inverse Laplace Algorithm used in Mathematica [closed]

I have a general question. I want to know what algorithm (The name of this numerical method) is already used to calculate a numerical Inverse Laplace in Mathematica. I used the function ...
Ali AlCapone's user avatar
2 votes
0 answers
90 views

What is LatticeReduce minimizing?

According to the manual page on LatticeReduce, The product of the norms of lattice vectors will decrease. The claim seems to be true for small lattices, but, I cannot confirm this for higher ...
A. Kato's user avatar
  • 406
5 votes
2 answers
96 views

Which model behind SequencePredictorFunction?

Is the SequencePredictorFunction based on Hidden Markov Model? If so, how does the model estimate the number of hidden states?
andrea simonetti's user avatar
3 votes
2 answers
413 views

Solving hyperbolic equation with parallelization in python by elucidating Mathematica algorithm [closed]

So! I have what seems like a rather basic hyperbolic equation to be solved written in Mathematica, with an integer input ‘n’, for which I am seeking the First {x,y} Integer solution to satisfy the ...
Travis Arlen McCracken's user avatar
2 votes
1 answer
83 views

Does SessionSubmit[ScheduledTask[expr,spec]] create a separate thread?

I don't quite understand how SessionSubmit[ScheduledTask[expr,spec]] works. We can create a background task, which can go without interfering with the main loop. If ...
Kirill Vasin's user avatar
  • 1,245
7 votes
1 answer
175 views

Are derivatives from interpolating functions created by NDSolve using derivatives of the spline or numerics from solving DE

For example given the code: ...
akozi's user avatar
  • 873
2 votes
0 answers
112 views

Theory behind CountRoots

According to this page, one can compute the number of zeros of a meromorphic function over an interval in Mathematica. I am working with univariate functions of x, which are meromorphic elementary ...
James Cheung's user avatar
4 votes
1 answer
107 views

Does RandomArrayLayer implement the reparameterization trick?

I want to backpropagate through random operations (e.g. learn the mean and variance of a Gaussian random variable). I assumed that ...
Ian's user avatar
  • 101
6 votes
1 answer
452 views

Is K is a function? [duplicate]

Some of the capitalized single letters are functions. Starting a fresh kernel and typing in the 24 capital letters, the following letters turn black since they are special functions, ...
Kai's user avatar
  • 2,099
6 votes
1 answer
189 views

How does FindRoot decide if a solution has converged?

I am solving a 1D non-linear differential equation using the finite element method with NDSolve. From the documentation I understand that the equation is discretized and then solved with FindRoot, I ...
Eddi's user avatar
  • 147
12 votes
1 answer
772 views

Mathematica internal source code for a numerical method

I read this (Six Reasons Why the Wolfram Language Is (Like) Open Source), which inspired me to re-investigate Mathematica's implementation of some numerical methods. I would like to see the code that ...
user avatar
1 vote
0 answers
148 views

How is MatrixPower computed?

How does Mathematica compute MatrixPower[m,n]? The documentation states it is equivalent to Dot @@ ConstantArray[m, {n}], but ...
GregH's user avatar
  • 1,919
6 votes
1 answer
263 views

Why does NMinimize using NelderMead without PostProcessing always stop at the nearest decade of iterations, e.g. 100, 110, 120

Why does NMinimize using NelderMead without PostProcessing always stop at the nearest decade of iterations, e.g. 100, 110, 120... Perhaps Catch Throw can prevent this e.g., see this answer, but that ...
Carl's user avatar
  • 717
4 votes
2 answers
289 views

RandomReal vs Random[Real

Let's say I need random numbers with increasing precision to simulate a particular function, so I try: ...
Syed's user avatar
  • 55.4k
2 votes
1 answer
93 views

My construction of additive compound matrices doesn't work sometimes; why?

Using the multiplicative compound matrix --see definition at Are compound matrices implemented in mathematica? one may define also additive compound matrices ...
florin's user avatar
  • 1,798
3 votes
1 answer
181 views

Construction of an additive compound matrix

I want to construct an additive compound matrix which has the following form as output : ...
Rim ADENANE's user avatar
1 vote
0 answers
67 views

Can LinearOptimization (ex. LinearProgramming) support other methods (Ellipsoid method)?

Can LinearOptimization (ex. LinearProgramming) support methods other than "InteriorPoint", "Simplex", "RevisedSimplex"? Specifically, I need the Ellipsoid Method. Are ...
Vepir's user avatar
  • 622
1 vote
0 answers
138 views

Spelunking yields <<kernel function>> [closed]

This is in the spirit of the old but excellent thread: What is the most convenient way to read definitions of in-memory symbols when we don't have the source files? (Spelunking tools) Using ...
qreus's user avatar
  • 343
2 votes
1 answer
209 views

How FunctionConvexity is implemented in Wolfram Mathematica?

I am wondering how FunctionConvexity was implemented in Wolfram Mathematica. I am trying to prove the convexity of the function $A(v)$: $$ A(v) = -\dfrac{k-1}{\...
Oleh's user avatar
  • 81
2 votes
1 answer
207 views

How is this asymptotic expansion of an integral calculated?

I am strongly impressed by this example from New in 13 as = AsymptoticIntegrate[ (t^10 + 3) Exp[I λ (t^5 + t + 1)], {t, -2, 2}, {λ, Infinity, 2} ] <...
user64494's user avatar
  • 26.6k
7 votes
1 answer
221 views

Exactly recreating image NetEncoder output

I need to run an exported NetChain in ONNX, surprisingly the model exports and runs successfully, but there's a discrepancy in the image loading stage (the NetEncoder'ing that must be written manually)...
user5601's user avatar
  • 3,625
1 vote
0 answers
104 views

How does Mathematica compute PseudoInverse? [closed]

What is the algorithm that is used in computing the PseudoInverse in Mathematica?
User101's user avatar
  • 613
1 vote
0 answers
156 views

How does Mathematica numerically evaluate Gaussian Hypergeometric function $_2F_1(a,b,c,z)$ where $a,b,c \in \mathbb{C}$ [closed]

Gaussian hypergeometric function $_2F_1(a,b,c,z)$ usually takes $a,b,c \in \mathbb{R}$ and $z \in \mathbb{C}$. But using Barnes integral representation, one can analytically continue the $_2F_1(a,b,c,...
Mockingbird's user avatar
4 votes
1 answer
72 views

`NetNestOperator` with custom input / output ports

I'd like to use NetNestOperator on a NetGraph that has custom input / output port names, such as ...
Fidel I. Schaposnik's user avatar
1 vote
0 answers
48 views

How is the Riemann-Siegel Z function evaluated in Mathematica?

How is the Riemann-Siegel Z function evaluated in Mathematica? Is it calculated using the first so many terms of its power series or by other means?
Gabi23's user avatar
  • 295
5 votes
2 answers
788 views

How does mathematica numerically evaluate the Meijer G-function?

The Meijer G-function is defined as a contour integral in the complex plane. Mathematica is able to numerically evaluate such a function quickly and accurately. How does she do that?
LeastSquare's user avatar
5 votes
0 answers
63 views

What are the valid options for the "ParametricCaching" option in ParametricNDSolve?

I am using ParametricNDSolve as part of the calculation of an objective function for an optimization, so I am trying to strike a balance between memory usage and ...
Michael Seifert's user avatar
2 votes
2 answers
139 views

How does Plot choose the points in a figure?

Taking a simple example, ...
Mark_Phys's user avatar
  • 491
7 votes
1 answer
733 views

Is anyone aware of Mathematica implementation of OpenAI's GPT-3 API?

Is anyone aware of Mathematica's implementation of OpenAI's GPT-3 API?
blackace's user avatar
  • 415
19 votes
2 answers
1k views

How do the functions Reap and Sow actually work behind the scenes?

I have read that using Reap and Sow to generate a list of unknown length is much more efficient than using AppendTo over and over. I believe AppendTo works by creating a new list that is one element ...
Cat Bisque's user avatar
1 vote
1 answer
77 views

AsymptoticOutputTracker to Matlab-Simulink

What is the structural scheme for AsymptoticOutputTracker in Mathematica? I need this in order to transfer the feedback signals received with the help of ...
dtn's user avatar
  • 2,404
11 votes
3 answers
315 views

Why does ReliefPlot not have a version for functions?

TL, DR: Use the LightingAngle option to DensityPlot to achieve the same shadows effect as ...
Gravifer's user avatar
  • 872
6 votes
1 answer
186 views

GroebnerBasis internals and runtime dependence on variable list ordering

I have a problem I have reduced to asking for a Gröbner basis. For some reason, Mathematica is able to solve this in minutes, while other programs more dedicated to these types of calculations run for ...
RedCurry's user avatar
9 votes
1 answer
213 views

Automatic Method for FindMinimum with constraints

What is the default method used by FindMinimum with constraints? The documentation says: "Currently, the only method available for constrained optimization is ...
VarunM's user avatar
  • 101
4 votes
0 answers
249 views

NMinimize differential evolution: how do search points and initial points really work?

I am puzzled by the settings of a specific method for NMinimize called DifferentialEvolution, given for example here. The ...
Kiro's user avatar
  • 1,521
3 votes
0 answers
83 views

Algorithm used by Mathematica for evaluating partial sums

Today, while using mathematica, I entered the command Sum[1/Factorial[n], {n, 0, x}] and found that: $$\sum_{x\geq n\geq0}\frac{1}{n!}=\frac{e\Gamma(x+1,1)}{\Gamma(...
user avatar
3 votes
0 answers
195 views

Is there any way to change the box-formatting behavior thoroughly?

Background The box-formmating behavior is affected by the downvalues of the symbol MakeBoxes. But the formatting behavior is not exactly ...
asd1dsa's user avatar
  • 573
4 votes
1 answer
409 views

Implementation of quantum Uhlmann Fidelity in Mathematica

I want to implement the quantum Uhlmann fidelity $$F(\rho, \sigma) := \mathrm{tr} \left[\sqrt{\sqrt{\rho} \sigma \sqrt{\rho}} \right]^2.$$ in Mathematica as a measure of "closeness" between ...
MrDerDart's user avatar
  • 201
6 votes
1 answer
145 views

Why Does Subsets[...,{n}] not Output a Packed List, Even Though it Doesn't Unpack?

Assume list is packed. I expect Subsets[] is a structural operation because it depends on the number of elements, not on what ...
Just Some Old Man's user avatar
2 votes
0 answers
52 views

How algebraic (or more general) expression is represented (memory level, data structure) in Mathematica?

This is not a question about the usage of Mathematica, it is more about implementation. Mathematica is a very huge system, it can provide a series of very advanced operation against a Mathematica ...
Audra Jacot's user avatar

1
2 3 4 5
7