This shows you the differences between two versions of the page.
Both sides previous revision Previous revision | Next revision Both sides next revision | ||
docs:utils:poincare [2009/02/23 12:56] gibson |
docs:utils:poincare [2009/02/23 13:00] gibson |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== poincare ====== | ====== poincare ====== | ||
- | compute crossings of a Poincare section in a pre-computed trajectory | + | Compute crossings of a Poincare section in a pre-computed trajectory. The Poincare condition is %%(u(t) - ueqb, e) == 0%%, i.e. the displacement from a given point %%ueqb%% is orthogonal to a direction %%e%%. The program checks for changes in sign in the inner product in the stored trajectory (with large %%dT%%, and when it finds one, goes back and reintegrates with small %%dt%%, and then does quadratic interpolation when the fine-scale trajectory recrosses the section. |
===== usage ===== | ===== usage ===== | ||