Contents
close all; clear all; global s0 s1 u0 u1 Psi par xi um1 um2 sig;
first some CPs at R=26
sd0='FSS'; sp0='pt17'; sd1='p1'; sp1='pt38'; flip=0;
fn=setfnflip(sd0,sp0,sd1,sp1,flip); opt=[]; opt=ocstanopt(opt);
opt.rhoi=1; opt.t1=100; opt.start=1; opt.tv=[]; opt.nti=10; opt.start=1;
call iscnat and plot
alvin=[0.1 0.25 0.5 0.75 1]; opt.Nmax=200; opt.msw=1; opt.Stats_step='on';
[alv,vv,sol,udat,tlv,tv,uv]=iscnat(alvin,[],[],opt,fn); v=[15,30]; vegsolplot(sol,v,4,fn);
save('17to38','sol','fn','alv','vv');
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
26.0000 0.0100 0.1000 1.0000 1.1000 0.3000
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
26.0000 0.0100 0.1000 1.0000 1.1000 0.3000
getting Psi, c-evals at i=160, mu=-0.196209, mu(i+1)=-0.196209,
c-evals at i=162, mu=-0.187169, mu(i+1)=-0.187169,
c-evals at i=166, mu=0.157165, mu(i+1)=0.157165,
c-evals at i=168, mu=0.166198, mu(i+1)=0.166198,
done in 0.245191s, n/2=164, d=0, suggested T=350.32
...
p1/p11 to FSS
sd0='FSS'; sp0='pt17'; sd1='p1'; sp1='pt11'; flip=1;
fn=setfnflip(sd0,sp0,sd1,sp1,flip);
[alv,vv,sol,udat,tlv,tv,uv]=iscnat(alvin,[],[],opt,fn);vegsolplot(sol,v,4,fn);
save('11to17','sol','fn','alv','vv');
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
26.0000 0.0100 0.1000 1.0000 1.1000 0.3000
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
26.0000 0.0100 0.1000 1.0000 1.1000 0.3000
getting Psi, c-evals at i=41, mu=-1923.48, mu(i+1)=-1923.48,
done in 0.215624s, n/2=164, d=0, suggested T=169.475
===> Non linear iteration = 1, Linear iteration = 1
N= 10, Maximum relative error = 0.00496996, order = 2
al=0.1, flag=0
...
R=20:
sd0='FSS'; sp0='pt28'; sd1='p1'; sp1='pt49'; flip=0;
fn=setfnflip(sd0,sp0,sd1,sp1,flip);
[alv,vv,sol,udat,tlv,tv,uv]=iscnat(alvin,[],[],opt,fn); vegsolplot(sol,v,15,fn);
save('28to49','sol','fn','alv','vv');
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
20.0000 0.0100 0.1000 1.0000 1.1000 0.3000
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
20.0000 0.0100 0.1000 1.0000 1.1000 0.3000
getting Psi, c-evals at i=160, mu=-0.205547, mu(i+1)=-0.205547,
c-evals at i=162, mu=-0.189794, mu(i+1)=-0.189794,
c-evals at i=166, mu=0.159789, mu(i+1)=0.159789,
c-evals at i=168, mu=0.175532, mu(i+1)=0.175532,
done in 0.209442s, n/2=164, d=0, suggested T=469.408
...
R=10:
sd0='FSS'; sp0='pt44'; sd1='p1'; sp1='pt65'; flip=0;
fn=setfnflip(sd0,sp0,sd1,sp1,flip); opt.msw=0;
[alv,vv,sol,udat,tlv,tv,uv]=iscnat(alvin,[],[],opt,fn); vegsolplot(sol,v,10,fn);
save('44to65','sol','fn','alv','vv');
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
10.0000 0.0100 0.1000 1.0000 1.1000 0.3000
ans =
Columns 1 through 7
0.0300 0.0010 0.5000 0.0300 0.0050 0.9000 0.0010
Columns 8 through 13
10.0000 0.0100 0.1000 1.0000 1.1000 0.3000
getting Psi, c-evals at i=161, mu=-0.184986, mu(i+1)=-0.184986,
c-evals at i=167, mu=0.154975, mu(i+1)=0.154975,
done in 0.218367s, n/2=164, d=0, suggested T=13995.9
===> Non linear iteration = 1, Linear iteration = 1
N= 10, Maximum relative error = 1.06602, order = 2
...