Skip to main content
uhoh's user avatar
uhoh's user avatar
uhoh's user avatar
uhoh
  • Member for 8 years, 8 months
  • Last seen more than a week ago
  • Taipei, Earth
49 votes
2 answers
13k views

What does "symplectic" mean in reference to numerical integrators, and does SciPy's odeint use them?

8 votes
1 answer
472 views

What algorithm does (or did?) Excel use for Bessel functions that is discontinuous at x=8?

4 votes
3 answers
863 views

Convolute a gaussian kernel with a large array of off-grid centroids without looping? (how to make "A Thousand (Gaussian) Points of Light" )

4 votes
1 answer
577 views

Going to try to move some of my scipy/numpy calculation to a new GPU, how to avoid disappointing results?

4 votes
0 answers
178 views

How amenable is this 2D Frenkel–Kontorova-like energy minimization problem in Python to the use of a modest PC + GPU? (Heavy reliance on indexing)

3 votes
0 answers
122 views

How to take the Fourier transform of a Fibonacci chain in a Python script?

3 votes
1 answer
335 views

What problems does softmax() solve and when should I think of using it - in simple terms

2 votes
1 answer
122 views

Need help to fully understand SciPy's odeint's reported step sizes, eval times, # of funct calls & total proc. time (re. question in Astronomy SE)

2 votes
1 answer
240 views

Another way to evaluate the gravitational force from a uniform cube?

2 votes
1 answer
180 views

Interpolating the gradient of a cylindrically symmetric potential field expected obey the Laplace equation, especially near/across r=0

1 vote
1 answer
421 views

2D DFT for lower frequencies only; is there something significantly faster than numpy.fft.fft2 (throwing away high frequencies)?

1 vote
1 answer
179 views

Simple, easy to install and use Python FEM solver (and example) for 2D cylindrical Laplace equation

1 vote
1 answer
344 views

How to generate the convolution of f(x, y) with a parametric function g(t), x(t), y(t) in Python? (Something better than this brute-force sum)

1 vote
1 answer
75 views

Why does this implementation for Eisenstein integer pairs of Euclid's method for finding greatest common denominators get stuck for this one point?