Time scales in the pseudospectral simulation of the heat equation
This is a discussion of a pseudospectral simulation of the stochastically forced heat equation and an analysis of the time scales needed for the results to converge.
Numerical Simulations of Fractional Gaussian Fields
I describe the basic properties of fractional Gaussian fields and show some numerical simulations.
Benchmarking Python FFTs
These are some benchmarks on Fast Fourier Transform in Python using different tools: Numpy, Scipy, PyFFTW and PyTorch. The arrays are one-dimensional with complex entries and up to 2^26 elements.
Equipartition Theorem and Analytical Results for a Quartic Potential
In this post I will discuss the application of the equipartition theorem to a non-quadratic Hamiltonian.