search

JuliaDynamics / Attractors.jl

Find attractors (and their basins) of dynamical systems. Perform global continuation. Study global stability (a.k.a. non-local, or resilience). Also tipping points functionality.
MIT License
32 stars 6 forks source link

Irregular grid #88

Closed StefanVaylBX2023 closed 1 year ago

StefanVaylBX2023 commented 1 year ago

Initial implementation, but very strange behavior of mapper for default test of new implementation of regular grid, need to be discussed

Datseris commented 1 year ago

i don't know what "very strange" means, but I hope you will explain on our meeting today

StefanVaylBX2023 commented 1 year ago

So, if you look into the code, for just base case with regular grid, literally nothing was changed in terms of algorithm structure, just basin_nfo now has this parameter grid_nfo which has the same as before .minima, .maxima and .steps. You can now explicitly call mapper.bsn_nfo.grid_nfo.minima/.maxima/.steps and correct values will be outputted, but for some reason calls to mapper(ics) doesn't work properly as it only possible outputs are -1,0,1. But in the same time basin_of_attraction returns 3 collect values for attractors. This is the point where I am a bit stacked and cannot test other features. New addition of basin_cell_index supposed to works the following way: 1) initial check for condition to be inside of grid 2) for each axis subtract the value of condition along this axis to find the closest one and output its index, on paper works just fine, but I cannot test it properly given that mappers seems have some bug.

Datseris commented 1 year ago

? @StefanVaylBX2023 what happened here?