Hi,
I've a dynamical systsem (Aircraft) with 8 state variables (U1,U2,U3..U8) (and eight first order ODEs ...F(1),F(2),..F(8)) and three parameters (PAR1, PAR2,PAR3). I used one of those as continuation parameter (PAR1) and the other two (PAR2 and PAR3) as free parameters, while i performed continuation and bifurcation analysis, i could only see a partial set of states (6)
Hope the figures below explain my question better. In the PyPlaut and bifurcation diagram, i could only see six states(U(1)..U(6)) displayed. But when i enter 'print start(1)' as shown in the attached picture, i see U(1)...U(7), U(8). I have 8 states in my model. I'm not able to understand where is it storing the data about how to access the same. That data is not there in fort.7 and fort.8.
Thanks,
Kumara
take the help of MATLAB or contact author of spectral methods Using MATLAB.
@Biswanath Rath: I'm sorry i could not understand your answer. I'm using AUTO-07P software and i don't understand how Matlab is related to this. Thanks.
Once you extract the AUTO files and before building using ./configure go to the path ..../auto/07p/src and open the file io.f90, in the file in the subroutine WRLINE you can find the lines nearly at the beginning of subroutine (line 1052 in auto07p-0.8 & around line 1095 in auto07p-0.9)
IF(N1.GT.7)THEN
N1=7
N2=0
ELSEIF(N1+N2.GT.7)THEN
N2=7-N1
ENDIF
Here the no.7 limits the size of your output printing (L2-NORM + 6 States (i guess)) increase to 10 i guess you can have 9th order system. in general increase 7 to n+1, n being the dimension of your state vector and then install using ./configure and make.
The attached file is an execution of auto in Windows using MinGw of an 8th order system with matrix being [BR, PT, TY, LAB, P(1), L2-NORM, U(1), U(2), U(3), U(4), U(5), U(6), U(7), U(8)]
if you are looking for windows instillation check
https://people.unica.it/massimilianogrosso/auto-on-windows/
Hope this is helpful.......
Also, During Periodic Computations you will get NAN or INF or some very small no printed for all the states with index greater than 7, This is also because of the limitations in no of variables allowed to be printed.
To overcome this you have to open ..../auto/07p/src/bvp.f90 and go to Subroutine STPLBV, and increase the size of the allocated array UMX to a required number. later in the same subroutine you will find the line
IF(N2.GT.7)N2=7
at line no. 1235, Increase the no. 7 to required to compute the max of orbits for all the states.
Also, The Parameter Over specification using the Subroutine PVLS (explained in demo pvl) can be used to print additional parameters, the computations added in PVLS will not be used in Continuation. It is barely to display additional parameters.
It is explained under Parameter Over specification topic in the manual