On Mon, Jun 15, 2015 at 8:17 PM, Piotr Held <

[hidden email]> wrote:

Basically this is emulating an event. Instead of terminating the

program you can ask for a handle to a function that will plot/print

the data provided every time in a while (could be another parameter,

something like "ShowEvery" and an integer indicating the number of

iterations to wait before showing a result). This is how it is done in

the odepkg package[1], so you can continue with the calculation but

from time to time update a plot with the current values. How do you

see this?

You could implement this part in an m-file and have a call to the C++

code only to do the calculation.

Some thing like

(assume ShowEvery was given)

while (calculation is not finished)

L = __calculate_lyap__ (..., ShowEvery);

print/plot L;

endwhile

The last argument of calculate_lyap is just indicating how many

iterations should be calculated. This means that at each call all the

necessary data has to be provided. It would be easier if one could

have a reference to the variables being updated, but I do not know if

this is possible (that is updating a plot/print from within the C++

code). I will check and report back.

Also looking at the code...we need to discuss the gram_schmidt in

there. It might be convenient to use qr (but this is just a future

improvement of TISEAN!).

[1]

http://sourceforge.net/p/octave/odepkg/ci/default/tree/inst/integrate_functions/integrate_adaptive.mline 280 and so on