Have some output from a CFD code and need to calculate u/s manifolds, fixed point locations etc. The model is finite difference so the best I can do is interpolate the data to give the right hand side. Would this be supported in PyDSTool? Thanks
If you would like to refer to this comment somewhere else in this project, copy and paste the following link:
Absolutely, although the accuracy of manifold calcs will naturally
depend on the mesh size of your interpolated data. Any "external
input" variable that you use in the RHS is interpolated. e.g. see
pydstool/tests/generator/test_events.py or
pydstool/examples/test_hybrid_extinputs.py, which use the
InterpolateTable generator structure. There's a tutorial for the
saddle manifold calcs too.
Have some output from a CFD code and need to calculate u/s manifolds, fixed
point locations etc. The model is finite difference so the best I can do is
interpolate the data to give the right hand side. Would this be supported in
PyDSTool? Thanks
I see, I wasn't thinking about that. The system I described is set up
for temporal, as you suggest. Probably the easiest route is to add
your own, custom interpolation functions for your data and embed them
either in Python or in C, depending on what solver you're using (you
won't be able to mix python functions in a C-based RHS, or vice
versa). The examples directory shows you both in different contexts.
For instance, read here: https://sourceforge.net/p/pydstool/discussion/472291/thread/2e6181b6
Have some output from a CFD code and need to calculate u/s manifolds, fixed point locations etc. The model is finite difference so the best I can do is interpolate the data to give the right hand side. Would this be supported in PyDSTool? Thanks
Absolutely, although the accuracy of manifold calcs will naturally
depend on the mesh size of your interpolated data. Any "external
input" variable that you use in the RHS is interpolated. e.g. see
pydstool/tests/generator/test_events.py or
pydstool/examples/test_hybrid_extinputs.py, which use the
InterpolateTable generator structure. There's a tutorial for the
saddle manifold calcs too.
On Tue, Feb 14, 2017 at 3:22 PM, James tlr004@users.sf.net wrote:
Spatial interpolation as well as temporal are covered?
I see, I wasn't thinking about that. The system I described is set up
for temporal, as you suggest. Probably the easiest route is to add
your own, custom interpolation functions for your data and embed them
either in Python or in C, depending on what solver you're using (you
won't be able to mix python functions in a C-based RHS, or vice
versa). The examples directory shows you both in different contexts.
For instance, read here:
https://sourceforge.net/p/pydstool/discussion/472291/thread/2e6181b6
On Thu, Feb 16, 2017 at 11:04 AM, James tlr004@users.sf.net wrote: