matrix exponential solution with abritrary precision floats

ukaea / decay-rate-ode-solver

Solver for a hyper-stiff set of ODEs

GNU Lesser General Public License v3.0

1 stars 1 forks source link

matrix exponential solution with abritrary precision floats #6

Open jwscook opened 5 years ago

jwscook commented 5 years ago

See ArbNumerics.jl and https://gist.github.com/JeffreySarnoff/b01361f8eada74974c5f4e25eae071ef

jwscook commented 5 years ago

This seems to be working really nicely, but takes a long time to solve. It essentially is solving the equations analytically. I trust a matrix exponential answer via ArbNumerics to be the most reliable.

jwscook commented 5 years ago

in python e.g.

from mpmath import mp
mp.exp(mp.randmatrix(3, 3))