View Source Integrator.SampleEqns (Integrator v0.1.3)

Functions to be used in testing

Summary

Functions

The Euler equations of a rigid body without external forces. This is a standard test problem proposed by Krogh for solvers intended for nonstiff problems [see below]. Based on "rigidode.m" from Matlab/Octave. The analytical solutions are Jacobian elliptic functions.

Simulates a point mass or particle falling through pass affected by gravity. Used for comparisons with the Matlab/Octave ballode.m routine

Functions

Link to this function

euler_equations(arg1, x)

View Source
@spec euler_equations(Nx.t(), Nx.t()) :: Nx.t()

The Euler equations of a rigid body without external forces. This is a standard test problem proposed by Krogh for solvers intended for nonstiff problems [see below]. Based on "rigidode.m" from Matlab/Octave. The analytical solutions are Jacobian elliptic functions.

See rigidode.m in Matlab.

Shampine, L. F., and M. K. Gordon, Computer Solution of Ordinary Differential Equations, W.H. Freeman & Co., 1975

Link to this function

falling_particle(arg1, x)

View Source
@spec falling_particle(Nx.t(), Nx.t()) :: Nx.t()

Simulates a point mass or particle falling through pass affected by gravity. Used for comparisons with the Matlab/Octave ballode.m routine

@spec van_der_pol_fn(Nx.t(), Nx.t()) :: Nx.t()

From octave b

fvdp = @(t,x) [x(2); (1 - x(1)^2) * x(2) - x(1)];