Joseph C. Slater
Joseph C. Slater
https://www.reddit.com/r/Python/comments/azbdvr/python_decorators_made_simple_part_1/
def duff_osc_ss(x, params): omega = params['omega'] t = params['cur_time'] xdot = np.array([[x[1]],[-x[0]-.1*x[0]**3-.1*x[1]+1*sin(omega*t)]]) return xdot t, x, e, amps, phases = ms.hb_freq(duff_osc_ss, num_variables = 2, omega = 0.1, eqform='first_order', num_harmonics=5, num_time_steps=51,...
Form of derivative function should be state space by default. Defaulting fails.
https://hgomersall.github.io/pyFFTW/pyfftw/builders/builders.html FFTW is much faster, but license restricted. The code should allow the user to optionally call fftw, but not rely on it. Further, speedups can be obtained with SciPy...
Need to raise exception and create working demo
Calculate Jacobian and use eigen solution for guess? Initial guess is a challenge. Not feasible to be very accurate for hard nonlinearities.
acronym "Science, Math, Technology, Engineering" doesn't discover STEM. (Excited to use this for real)
As of this date, the module does not handle (is not aware) of sparse matrices. It would be prudent for all system matrices from FEA to be considered sparse. This...
Correct negative sign in code to match equation in manuscript
http://www.scipy-lectures.org/advanced/optimizing/index.html#example-of-the-svd