
%landscape AB: 2 basins connected by a saddle point
%base eigenvectors at acute angles======
clear f;
x0=[-10; 10];
y0=[-10; 10];
R1=[4 ; 4];
R2=[10 ; 10];
w=[1;1];
theta = pi*[0.15; 0.55];
gamma=[1;1];
parms=[x0 y0 R1 R2 theta gamma w];

landscape.AB=parms;

%landscape AB2 another version of the
% accessible basin landscape
% these are closer than AB
clear f;
x0=[-5; 5];
y0=[-5; 5];
R1=[5 ; 5];
R2=[9 ; 9];
w=[1;1];
theta = pi*[0.15; 0.55];
gamma=[1;1];
parms=[x0 y0 R1 R2 theta gamma w];

landscape.AB2=parms;
%=========================================
% %landscape SB:  single basin 
theta=0.2*pi
%theta=0;
R1=3; R2=5;
x0=0; y0=0; w=1;
R=[cos(theta) sin(theta); -sin(theta) cos(theta)];
S=[R1^2 0 ; 0 R2^2];
G=R'*S^-1*R;
gamma=1;  %really flat 0.25, normal 0.75
parms=[];
parms=[x0 y0 R1 R2 theta gamma w];
landscape.SB=parms;
%========================================
% %landscape SB2:  single basin 
% a skinny version so that the 
% principal axes are visible
theta=0.2*pi
%theta=0;
R1=1; R2=5;
x0=0; y0=0; w=1;
R=[cos(theta) sin(theta); -sin(theta) cos(theta)];
S=[R1^2 0 ; 0 R2^2];
G=R'*S^-1*R;
gamma=1;  %really flat 0.25, normal 0.75
parms=[];
parms=[x0 y0 R1 R2 theta gamma w];
landscape.SB2=parms;
%=========================================
%landscape DB:  disjoint basins 
% base eigenvectors at acute angles======
clear f;
x0=[-12; 8];
y0=[-12; 12];
R1=[2 ; 5];
R2=[5 ; 10];
w=[1;4];
theta = pi*[0.1; 0.35];
gamma=[4; 7];
parms=[x0 y0 R1 R2 theta gamma w];
landscape.DB=parms;
%=========================================

%$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
%landscape MB: Multiple Basins
%complicated landscape
% landscape parameters======
clear f;
x0=[-12;   5;   3;    15;  -3];
y0=[-15;  12;  -10;    5;  2];
R1=[  6;   5;   5;    10;  2];
R2=[  14;  10;   3;    8;  3];
w= [  3;   4;   2;    3;  1];
%w= [  1;   1;   1;    1;  1];
theta = pi*...
[  0.15; .25;   0; -0.1; 0.3];
gamma=...
1*[     2;   3;   1;    2;   1;];
parms=[x0 y0 R1 R2 theta gamma w];
landscape.MB=parms;
%=========================================

%-------------------------------------------
%landscape OB: 2 orthogonal basins
%orthogonal to axes
clear parms
clear f;
x0=[-10; 10];
y0=[-10; 10];
R1=[4 ; 12];
R2=[12 ; 4];
theta=[0;0];
w=[1;1];   %different weights cause interesting offset!theta = pi*[0.0; 0.0];
gamma=[9;9];
parms=[x0 y0 R1 R2 theta gamma w];
landscape.OB=parms
%=========================================
