| [sfac,sdfac,dfac,dqfac,kprint]=check2(sfac,sdfac,dfac,dqfac,kprint); |
function [sfac,sdfac,dfac,dqfac,kprint]=check2(sfac,sdfac,dfac,dqfac,kprint);
persistent ca cdotc cdotu csize1 csize2 ct10x ct10y ct6 ct7 ct8 cx cx1 cy cy1 da db dc dpar dparam ds dsize dsize1 dsize2 dt10x dt10y dt19x dt19xa dt19xb dt19xc dt19xd dt19y dt19ya dt19yb dt19yc dt19yd dt2 dt7 dt8 dt9x dt9y dtx dty dx dx1 dx2 dy dy1 dy2 firstCall i incxs incys j ki kn kni kpar ksize lens lenx leny mx my ns qc sa sb sc sparam ss ssize ssize1 ssize2 ssize3 st7b stx sty sx sy ; if isempty(firstCall),firstCall=1;end;
if isempty(i), i=0; end;
global combla_2; if isempty(combla_2), combla_2=0; end;
global combla_4; if isempty(combla_4), combla_4=0; end;
global combla_5; if isempty(combla_5), combla_5=0; end;
if isempty(j), j=0; end;
if isempty(ki), ki=0; end;
if isempty(kn), kn=0; end;
if isempty(kni), kni=0; end;
if isempty(kpar), kpar=0; end;
if isempty(ksize), ksize=0; end;
if isempty(lenx), lenx=0; end;
if isempty(leny), leny=0; end;
global combla_6; if isempty(combla_6), combla_6=0; end;
if isempty(mx), mx=0; end;
if isempty(my), my=0; end;
global combla_3; if isempty(combla_3), combla_3=0; end;
global combla_1; if isempty(combla_1), combla_1=0; end;
if isempty(sa), sa=0; end;
if isempty(sb), sb=0; end;
if isempty(sc), sc=0; end;
if isempty(ss), ss=0; end;
% common :: ;
global combla_7; if isempty(combla_7), combla_7=false; end;
%% common /combla/ nprint , icase , n , incx , incy , mode , pass;
%% common /combla/ combla_1 , combla_2 , combla_3 , combla_4 , combla_5 , combla_6 , combla_7;
if isempty(incxs), incxs=zeros(1,4); end;
if isempty(incys), incys=zeros(1,4); end;
if isempty(lens), lens=zeros(4,2); end;
if isempty(ns), ns=zeros(1,4); end;
if isempty(sx), sx=zeros(1,7); end;
if isempty(sy), sy=zeros(1,7); end;
if isempty(stx), stx=zeros(1,7); end;
if isempty(sty), sty=zeros(1,7); end;
if isempty(ssize1), ssize1=zeros(1,4); end;
if isempty(ssize2), ssize2=zeros(14,2); end;
if isempty(ssize), ssize=zeros(1,7); end;
if isempty(qc), qc=zeros(1,30); end;
if isempty(sparam), sparam=zeros(1,5); end;
if isempty(st7b), st7b=zeros(4,4); end;
if isempty(ssize3), ssize3=zeros(1,4); end;
if isempty(dx), dx=zeros(1,7); end;
if isempty(da), da=0; end;
if isempty(dx1), dx1=zeros(1,7); end;
if isempty(dy1), dy1=zeros(1,7); end;
if isempty(dy), dy=zeros(1,7); end;
if isempty(dt7), dt7=zeros(4,4); end;
if isempty(dt8), dt8=zeros(7,4,4); end;
if isempty(dx2), dx2=zeros(1,7); end;
if isempty(dy2), dy2=zeros(1,7); end;
if isempty(dt2), dt2=zeros(4,4,2); end;
if isempty(dparam), dparam=zeros(1,5); end;
if isempty(dpar), dpar=zeros(5,4); end;
if isempty(dt10x), dt10x=zeros(7,4,4); end;
if isempty(dt10y), dt10y=zeros(7,4,4); end;
if isempty(db), db=0; end;
if isempty(dsize1), dsize1=zeros(1,4); end;
if isempty(dsize2), dsize2=zeros(7,2); end;
if isempty(dsize), dsize=zeros(1,7); end;
if isempty(dc), dc=0; end;
if isempty(ds), ds=0; end;
if isempty(dt9x), dt9x=zeros(7,4,4); end;
if isempty(dt9y), dt9y=zeros(7,4,4); end;
if isempty(dtx), dtx=zeros(1,7); end;
if isempty(dty), dty=zeros(1,7); end;
if isempty(dt19x), dt19x=zeros(7,4,16); end;
if isempty(dt19xa), dt19xa=zeros(7,4,4); end;
if isempty(dt19xb), dt19xb=zeros(7,4,4); end;
if isempty(dt19xc), dt19xc=zeros(7,4,4); end;
if isempty(dt19xd), dt19xd=zeros(7,4,4); end;
if isempty(dt19y), dt19y=zeros(7,4,16); end;
if isempty(dt19ya), dt19ya=zeros(7,4,4); end;
if isempty(dt19yb), dt19yb=zeros(7,4,4); end;
if isempty(dt19yc), dt19yc=zeros(7,4,4); end;
if isempty(dt19yd), dt19yd=zeros(7,4,4); end;
%!!EQUIVALENCE (dt19x(1,1,1),dt19xa(1,1,1))
%!!EQUIVALENCE (dt19x(1,1,5),dt19xb(1,1,1))
%!!EQUIVALENCE (dt19x(1,1,9),dt19xc(1,1,1))
%!!EQUIVALENCE (dt19x(1,1,13),dt19xd(1,1,1))
%!!EQUIVALENCE (dt19y(1,1,1),dt19ya(1,1,1))
%!!EQUIVALENCE (dt19y(1,1,5),dt19yb(1,1,1))
%!!EQUIVALENCE (dt19y(1,1,9),dt19yc(1,1,1))
%!!EQUIVALENCE (dt19y(1,1,13),dt19yd(1,1,1))
if isempty(cx), cx=zeros(1,7); end;
if isempty(ca), ca=0; end;
if isempty(cx1), cx1=zeros(1,7); end;
if isempty(cy1), cy1=zeros(1,7); end;
if isempty(cy), cy=zeros(1,7); end;
if isempty(ct6), ct6=zeros(4,4); end;
if isempty(ct7), ct7=zeros(4,4); end;
if isempty(ct8), ct8=zeros(7,4,4); end;
if isempty(csize1), csize1=zeros(1,4); end;
if isempty(csize2), csize2=zeros(7,2); end;
if isempty(ct10x), ct10x=zeros(7,4,4); end;
if isempty(ct10y), ct10y=zeros(7,4,4); end;
if isempty(cdotc), cdotc=0; end;
if isempty(cdotu), cdotu=0; end;
if firstCall, sa =[.3]; end;
if firstCall, da =[.3d0]; end;
if firstCall, ca =[complex(.4,-.7)]; end;
if firstCall, db =[.25d0]; end;
if firstCall, sb=[.1]; end;
if firstCall, incxs=[1,2,-2,-1]; end;
if firstCall, incys=[1,-2,1,-2]; end;
if firstCall, lens([1:4],1)=[1,1,2,4]; end;
if firstCall, lens([1:4],2)=[1,1,3,7]; end;
if firstCall, ns=[0,1,2,4]; end;
if firstCall, sc =[.8]; end;
if firstCall, ss =[.6]; end;
if firstCall, dc =[.8d0]; end;
if firstCall, ds=[.6d0]; end;
if firstCall, dx1=[.6d0,.1d0,-.5d0,.8d0,.9d0,-.3d0,-.4d0]; end;
if firstCall, dy1=[.5d0,-.9d0,.3d0,.7d0,-.6d0,.2d0,.8d0]; end;
if firstCall, dx2=[1.0d0,.01d0,.02d0,1.0d0,.06d0,2.0d0,1.0d0]; end;
if firstCall, dy2=[1.0d0,.04d0,-.03d0,-1.0d0,.05d0,3.0d0,-1.0d0]; end;
if firstCall, cx1=[complex(.7,-.8),complex(-.4,-.7),complex(-.1,-.9),complex(.2,-.8),complex(-.9,-.4),complex(.1,.4),complex(-.6,.6)]; end;
if firstCall, cy1=[complex(.6,-.6),complex(-.9,.5),complex(.7,-.6),complex(.1,-.5),complex(-.1,-.2),complex(-.5,-.3),complex(.8,-.7)]; end;
if firstCall, dt2=[0.25d0,1.25d0,1.2504d0,0.2498d0,0.25d0,1.25d0,0.24d0,0.2492d0,0.25d0,1.25d0,0.31d0,0.2518d0,0.25d0,1.25d0,1.2497d0,0.2507d0,0.0d0,2.0d0,2.0008d0,-.0004d0,0.0d0,2.0d0,-.02d0,-.0016d0,0.0d0,2.0d0,.12d0,.0036d0,0.0d0,2.0d0,1.9994d0,.0014d0]; end;
if firstCall, dt7=[0.0d0,.30d0,.21d0,.62d0,0.0d0,.30d0,-.07d0,.85d0,0.0d0,.30d0,-.79d0,-.74d0,0.0d0,.30d0,.33d0,1.27d0]; end;
if firstCall, st7b=[.1,.4,.31,.72,.1,.4,.03,.95,.1,.4,-.69,-.64,.1,.4,.43,1.37]; end;
if firstCall, ct7=[complex(0.,0.),complex(-.06,-.90),complex(.65,-.47),complex(-.34,-1.22),complex(0.,0.),complex(-.06,-.90),complex(-.59,-1.46),complex(-1.04,-.04),complex(0.,0.),complex(-.06,-.90),complex(-.83,.59),complex(.07,-.37),complex(0.,0.),complex(-.06,-.90),complex(-.76,-1.15),complex(-1.33,-1.82)]; end;
if firstCall, ct6=[complex(0.,0.),complex(.90,0.06),complex(.91,-.77),complex(1.80,-.10),complex(0.,0.),complex(.90,0.06),complex(1.45,.74),complex(.20,.90),complex(0.,0.),complex(.90,0.06),complex(-.55,.23),complex(.83,-.39),complex(0.,0.),complex(.90,0.06),complex(1.04,0.79),complex(1.95,1.22)]; end;
if firstCall, dt8=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,-.87d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,-.87d0,.15d0,.94d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.35d0,-.9d0,.48d0,0.0d0,0.0d0,0.0d0,0.0d0,.38d0,-.9d0,.57d0,.7d0,-.75d0,.2d0,.98d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.35d0,-.72d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.38d0,-.63d0,.15d0,.88d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,-.9d0,.33d0,0.0d0,0.0d0,0.0d0,0.0d0,.68d0,-.9d0,.33d0,.7d0,-.75d0,.2d0,1.04d0]; end;
if firstCall, ct8=[complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(-1.55,.5),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(-1.55,.5),complex(.03,-.89),complex(-.38,-.96),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.07,-.89),complex(-.9,.5),complex(.42,-1.41),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.78,.06),complex(-.9,.5),complex(.06,-.13),complex(.1,-.5),complex(-.77,-.49),complex(-.5,-.3),complex(.52,-1.51),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.07,-.89),complex(-1.18,-.31),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.78,.06),complex(-1.54,.97),complex(.03,-.89),complex(-.18,-1.31),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(-.9,.5),complex(.05,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.32,-1.41),complex(-.9,.5),complex(.05,-.6),complex(.1,-.5),complex(-.77,-.49),complex(-.5,-.3),complex(.32,-1.16)]; end;
if firstCall, dt9x=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,-.46d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,-.46d0,-.22d0,1.06d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.66d0,.1d0,-.1d0,0.0d0,0.0d0,0.0d0,0.0d0,.96d0,.1d0,-.76d0,.8d0,.90d0,-.3d0,-.02d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.06d0,.1d0,-.1d0,0.0d0,0.0d0,0.0d0,0.0d0,.90d0,.1d0,-.22d0,.8d0,.18d0,-.3d0,-.02d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,.26d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.78d0,.26d0,-.76d0,1.12d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt9y=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,-.78d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,-.78d0,.54d0,.08d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,-.9d0,-.12d0,0.0d0,0.0d0,0.0d0,0.0d0,.64d0,-.9d0,-.30d0,.7d0,-.18d0,.2d0,.28d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,-1.08d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.64d0,-1.26d0,.54d0,.20d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,-.9d0,.18d0,0.0d0,0.0d0,0.0d0,0.0d0,.04d0,-.9d0,.18d0,.7d0,-.18d0,.2d0,.16d0]; end;
if firstCall, dt10x=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,.7d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.3d0,.1d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,.8d0,.1d0,-.6d0,.8d0,.3d0,-.3d0,.5d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,.1d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,.1d0,.3d0,.8d0,-.9d0,-.3d0,.5d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,.3d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,.3d0,-.6d0,.8d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt10y=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,.8d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.5d0,-.9d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,-.4d0,-.9d0,.9d0,.7d0,-.5d0,.2d0,.6d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.5d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.4d0,.9d0,-.5d0,.6d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,-.9d0,.1d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,-.9d0,.1d0,.7d0,-.5d0,.2d0,.8d0]; end;
if firstCall, ct10x=[complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(-.9,.5),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(-.9,.5),complex(.7,-.6),complex(.1,-.5),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.6),complex(-.4,-.7),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.8,-.7),complex(-.4,-.7),complex(-.1,-.2),complex(.2,-.8),complex(.7,-.6),complex(.1,.4),complex(.6,-.6),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.9,.5),complex(-.4,-.7),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.1,-.5),complex(-.4,-.7),complex(.7,-.6),complex(.2,-.8),complex(-.9,.5),complex(.1,.4),complex(.6,-.6),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(.7,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(.7,-.6),complex(-.1,-.2),complex(.8,-.7),complex(0.,0.),complex(0.,0.),complex(0.,0.)]; end;
if firstCall, ct10y=[complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(-.4,-.7),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(-.4,-.7),complex(-.1,-.9),complex(.2,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.1,-.9),complex(-.9,.5),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.6,.6),complex(-.9,.5),complex(-.9,-.4),complex(.1,-.5),complex(-.1,-.9),complex(-.5,-.3),complex(.7,-.8),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.1,-.9),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(-.6,.6),complex(-.9,-.4),complex(-.1,-.9),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.6,-.6),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(-.9,.5),complex(-.4,-.7),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(.7,-.8),complex(-.9,.5),complex(-.4,-.7),complex(.1,-.5),complex(-.1,-.9),complex(-.5,-.3),complex(.2,-.8)]; end;
if firstCall, dt19xa=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,3.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,2.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,-.4d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,.8d0,0.0d0,0.0d0,0.0d0,-.8d0,3.8d0,-2.2d0,-1.2d0,0.0d0,0.0d0,0.0d0,-.9d0,2.8d0,-1.4d0,-1.3d0,0.0d0,0.0d0,0.0d0,3.5d0,-.4d0,-2.2d0,4.7d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt19xb=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.1d0,-3.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.3d0,.1d0,-2.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.3d0,.1d0,-2.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,.8d0,.9d0,-.3d0,-.4d0,-2.0d0,.1d0,1.4d0,.8d0,.6d0,-.3d0,-2.8d0,-1.8d0,.1d0,1.3d0,.8d0,0.0d0,-.3d0,-1.9d0,3.8d0,.1d0,-3.1d0,.8d0,4.8d0,-.3d0,-1.5d0]; end;
if firstCall, dt19xc=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,0.0d0,0.0d0,0.0d0,0.0d0,4.8d0,.1d0,-3.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.3d0,.1d0,-2.0d0,0.0d0,0.0d0,0.0d0,0.0d0,2.1d0,.1d0,-2.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,.8d0,.9d0,-.3d0,-.4d0,-1.6d0,.1d0,-2.2d0,.8d0,5.4d0,-.3d0,-2.8d0,-1.5d0,.1d0,-1.4d0,.8d0,3.6d0,-.3d0,-1.9d0,3.7d0,.1d0,-2.2d0,.8d0,3.6d0,-.3d0,-1.5d0]; end;
if firstCall, dt19xd=[.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.8d0,-1.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.9d0,-.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,3.5d0,.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.6d0,.1d0,-.5d0,.8d0,0.0d0,0.0d0,0.0d0,-.8d0,-1.0d0,1.4d0,-1.6d0,0.0d0,0.0d0,0.0d0,-.9d0,-.8d0,1.3d0,-1.6d0,0.0d0,0.0d0,0.0d0,3.5d0,.8d0,-3.1d0,4.8d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt19ya=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,-4.8d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,-.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,3.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,.7d0,0.0d0,0.0d0,0.0d0,.7d0,-4.8d0,3.0d0,1.1d0,0.0d0,0.0d0,0.0d0,1.7d0,-.7d0,-.7d0,2.3d0,0.0d0,0.0d0,0.0d0,-2.6d0,3.5d0,-.7d0,-3.6d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt19yb=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,0.0d0,0.0d0,0.0d0,0.0d0,4.0d0,-.9d0,-.3d0,0.0d0,0.0d0,0.0d0,0.0d0,-.5d0,-.9d0,1.5d0,0.0d0,0.0d0,0.0d0,0.0d0,-1.5d0,-.9d0,-1.8d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,.7d0,-.6d0,.2d0,.8d0,3.7d0,-.9d0,-1.2d0,.7d0,-1.5d0,.2d0,2.2d0,-.3d0,-.9d0,2.1d0,.7d0,-1.6d0,.2d0,2.0d0,-1.6d0,-.9d0,-2.1d0,.7d0,2.9d0,.2d0,-3.8d0]; end;
if firstCall, dt19yc=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,4.0d0,-6.3d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-.5d0,.3d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-1.5d0,3.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,.7d0,0.0d0,0.0d0,0.0d0,3.7d0,-7.2d0,3.0d0,1.7d0,0.0d0,0.0d0,0.0d0,-.3d0,.9d0,-.7d0,1.9d0,0.0d0,0.0d0,0.0d0,-1.6d0,2.7d0,-.7d0,-3.4d0,0.0d0,0.0d0,0.0d0]; end;
if firstCall, dt19yd=[.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,0.0d0,0.0d0,0.0d0,0.0d0,.7d0,-.9d0,1.2d0,0.0d0,0.0d0,0.0d0,0.0d0,1.7d0,-.9d0,.5d0,0.0d0,0.0d0,0.0d0,0.0d0,-2.6d0,-.9d0,-1.3d0,0.0d0,0.0d0,0.0d0,0.0d0,.5d0,-.9d0,.3d0,.7d0,-.6d0,.2d0,.8d0,.7d0,-.9d0,1.2d0,.7d0,-1.5d0,.2d0,1.6d0,1.7d0,-.9d0,.5d0,.7d0,-1.6d0,.2d0,2.4d0,-2.6d0,-.9d0,-1.3d0,.7d0,2.9d0,.2d0,-4.0d0]; end;
if firstCall, ssize1=[0.,.3,1.6,3.2]; end;
if firstCall, dsize1=[0.0d0,.3d0,1.6d0,3.2d0]; end;
if firstCall, ssize3=[.1,.4,1.7,3.3]; end;
if firstCall, csize1=[complex(0.,0.),complex(.9,.9),complex(1.63,1.73),complex(2.90,2.78)]; end;
if firstCall, ssize2=[0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17,1.17]; end;
if firstCall, dsize2=[0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,0.0d0,1.17d0,1.17d0,1.17d0,1.17d0,1.17d0,1.17d0,1.17d0]; end;
if firstCall, csize2=[complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(0.,0.),complex(1.54,1.54),complex(1.54,1.54),complex(1.54,1.54),complex(1.54,1.54),complex(1.54,1.54),complex(1.54,1.54),complex(1.54,1.54)]; end;
if firstCall, dpar=[-2.0d0,0.0d0,0.0d0,0.0d0,0.0d0,-1.0d0,2.0d0,-3.0d0,-4.0d0,5.0d0,0.0d0,0.0d0,2.0d0,-3.0d0,0.0d0,1.0d0,5.0d0,2.0d0,0.0d0,-4.0d0]; end;csize2=reshape(csize2,[7,2]);ct10x=reshape(ct10x,[7,4,4]);ct10y=reshape(ct10y,[7,4,4]);ct6=reshape(ct6,[4,4]);ct7=reshape(ct7,[4,4]);ct8=reshape(ct8,[7,4,4]);dpar=reshape(dpar,[5,4]);dsize2=reshape(dsize2,[7,2]);dt10x=reshape(dt10x,[7,4,4]);dt10y=reshape(dt10y,[7,4,4]);dt19xa=reshape(dt19xa,[7,4,4]);dt19xb=reshape(dt19xb,[7,4,4]);dt19xc=reshape(dt19xc,[7,4,4]);dt19xd=reshape(dt19xd,[7,4,4]);dt19ya=reshape(dt19ya,[7,4,4]);dt19yb=reshape(dt19yb,[7,4,4]);dt19yc=reshape(dt19yc,[7,4,4]);dt19yd=reshape(dt19yd,[7,4,4]);dt2=reshape(dt2,[4,4,2]);dt7=reshape(dt7,[4,4]);dt8=reshape(dt8,[7,4,4]);dt9x=reshape(dt9x,[7,4,4]);dt9y=reshape(dt9y,[7,4,4]);lens=reshape(lens,[4,2]);ssize2=reshape(ssize2,[14,2]);st7b=reshape(st7b,[4,4]);
firstCall=0;
%Barrowes - commented out calls to stest with complex variables
dt19x(:,:,[1:4])=dt19xa;
dt19x(:,:,[5:8])=dt19xb;
dt19x(:,:,[9:12])=dt19xc;
dt19x(:,:,[13:16])=dt19xd;
dt19y(:,:,[1:4])=dt19ya;
dt19y(:,:,[5:8])=dt19yb;
dt19y(:,:,[9:12])=dt19yc;
dt19y(:,:,[13:16])=dt19yd;
for ki = 1 : 4;
combla_4 = fix(incxs(ki));
combla_5 = fix(incys(ki));
mx = fix(abs(combla_4));
my = fix(abs(combla_5));
for kn = 1 : 4;
combla_3 = fix(ns(kn));
ksize = fix(min(2,kn));
lenx = fix(lens(kn,mx));
leny = fix(lens(kn,my));
for i = 1 : 7;
sx(i) = dx1(i);
sy(i) = dy1(i);
dx(i) = dx1(i);
dy(i) = dy1(i);
cx(i) = cx1(i);
cy(i) = cy1(i);
end; i = fix(7+1);
if( combla_2==2 )
[dumvar1,dumvar2,dumvar3,ssize1(sub2ind(size(ssize1),max(kn,1)):end),sfac,kprint]=stest(1,real(dsdot(combla_3,sx,combla_4,sy,combla_5)),real(dt7(kn,ki)),ssize1(sub2ind(size(ssize1),max(kn,1)):end),sfac,kprint);
elseif( combla_2==3 ) ;
[dumvar1,dumvar2,st7b(sub2ind(size(st7b),kn,ki):end),ssize3(sub2ind(size(ssize3),max(kn,1)):end),sfac,kprint]=stest(1,sdsdot(combla_3,sb,sx,combla_4,sy,combla_5),st7b(sub2ind(size(st7b),kn,ki):end),ssize3(sub2ind(size(ssize3),max(kn,1)):end),sfac,kprint);
elseif( combla_2==4 ) ;
[dumvar1,dumvar2,dt7(sub2ind(size(dt7),kn,ki):end),dsize1(sub2ind(size(dsize1),max(kn,1)):end),dfac,kprint]=dtest(1,ddot(combla_3,dx,combla_4,dy,combla_5),dt7(sub2ind(size(dt7),kn,ki):end),dsize1(sub2ind(size(dsize1),max(kn,1)):end),dfac,kprint);
elseif( combla_2==5 ) ;
combla_6 = 1;
[dumvar1,dumvar2,dumvar3,dumvar3,dqfac,kprint]=dtest(1,dqdoti(combla_3,db,qc,dx2,combla_4,dy2,combla_5),dt2(sub2ind(size(dt2),kn,ki,1):end),dt2(sub2ind(size(dt2),kn,ki,1):end),dqfac,kprint); dt2(sub2ind(size(dt2),kn,ki,1):end)=dumvar3;
elseif( combla_2==6 ) ;
combla_6 = 1;
[dumvar1,dumvar2,dumvar3,dumvar3,dqfac,kprint]=dtest(1,dqdoti(combla_3,db,qc,dx2,combla_4,dy2,combla_5),dt2(sub2ind(size(dt2),kn,ki,1):end),dt2(sub2ind(size(dt2),kn,ki,1):end),dqfac,kprint); dt2(sub2ind(size(dt2),kn,ki,1):end)=dumvar3;
combla_6 = 2;
[dumvar1,dumvar2,dumvar3,dumvar3,dqfac,kprint]=dtest(1,dqdota(combla_3,-db,qc,dx2,combla_4,dy2,combla_5),dt2(sub2ind(size(dt2),kn,ki,2):end),dt2(sub2ind(size(dt2),kn,ki,2):end),dqfac,kprint); dt2(sub2ind(size(dt2),kn,ki,2):end)=dumvar3;
elseif( combla_2==7 ) ;
%Barrowes - commented out calls to stest with complex variables
%!! CALL STEST(2,CDOTC(N,cx,INCx,cy,INCy),ct6(kn,ki),csize1(kn),%!!Sfac,Kprint)
elseif( combla_2==8 ) ;
%Barrowes - commented out calls to stest with complex variables
%!! CALL STEST(2,CDOTU(N,cx,INCx,cy,INCy),ct7(kn,ki),csize1(kn),%!!Sfac,Kprint)
elseif( combla_2==9 ) ;
[combla_3,sa,sx,combla_4,sy,combla_5]=saxpy(combla_3,sa,sx,combla_4,sy,combla_5);
for j = 1 : leny;
sty(j) = dt8(j,kn,ki);
end; j = fix(leny+1);
[leny,sy,sty,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint]=stest(leny,sy,sty,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint);
elseif( combla_2==10 ) ;
[combla_3,da,dx,combla_4,dy,combla_5]=daxpy(combla_3,da,dx,combla_4,dy,combla_5);
[leny,dy,dt8(sub2ind(size(dt8),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint]=dtest(leny,dy,dt8(sub2ind(size(dt8),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint);
elseif( combla_2==11 ) ;
[combla_3,ca,cx,combla_4,cy,combla_5]=caxpy(combla_3,ca,cx,combla_4,cy,combla_5);
%Barrowes - commented out calls to stest with complex variables
%!! CALL STEST(2*leny,cy,ct8(1,kn,ki),csize2(1,ksize),Sfac,Kprint)
elseif( combla_2==12 || combla_2==13 || combla_2==16 ||combla_2==17 ) ;
error(['stop encountered in original fortran code ',char(10),';']);
elseif( combla_2==14 ) ;
for i = 1 : 7;
sx(i) = dx1(i);
sy(i) = dy1(i);
stx(i) = dt9x(i,kn,ki);
sty(i) = dt9y(i,kn,ki);
end; i = fix(7+1);
[combla_3,sx,combla_4,sy,combla_5,sc,ss]=srot(combla_3,sx,combla_4,sy,combla_5,sc,ss);
[lenx,sx,stx,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint]=stest(lenx,sx,stx,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint);
[leny,sy,sty,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint]=stest(leny,sy,sty,ssize2(sub2ind(size(ssize2),1,ksize):end),sfac,kprint);
elseif( combla_2==15 ) ;
for i = 1 : 7;
dx(i) = dx1(i);
dy(i) = dy1(i);
end; i = fix(7+1);
[combla_3,dx,combla_4,dy,combla_5,dc,ds]=drot(combla_3,dx,combla_4,dy,combla_5,dc,ds);
[lenx,dx,dt9x(sub2ind(size(dt9x),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint]=dtest(lenx,dx,dt9x(sub2ind(size(dt9x),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint);
[leny,dy,dt9y(sub2ind(size(dt9y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint]=dtest(leny,dy,dt9y(sub2ind(size(dt9y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,ksize):end),dfac,kprint);
elseif( combla_2==18 ) ;
kni = fix(kn + 4.*(ki-1));
for kpar = 1 : 4;
for i = 1 : 7;
sx(i) = dx1(i);
sy(i) = dy1(i);
stx(i) = dt19x(i,kpar,kni);
sty(i) = dt19y(i,kpar,kni);
end; i = fix(7+1);
for i = 1 : 5;
sparam(i) = dpar(i,kpar);
end; i = fix(5+1);
combla_6 = fix(fix(sparam(1)));
for i = 1 : lenx;
ssize(i) = stx(i);
end; i = fix(lenx+1);
if((kpar==2) &&(kni==7) )
ssize(1) = 2.4e0;
end;
if((kpar==3) &&(kni==8) )
ssize(5) = 1.8e0;
end;
[combla_3,sx,combla_4,sy,combla_5,sparam]=srotm(combla_3,sx,combla_4,sy,combla_5,sparam);
[lenx,sx,stx,ssize,sfac,kprint]=stest(lenx,sx,stx,ssize,sfac,kprint);
sty_orig=sty; [leny,sy,sty,dumvar4,sfac,kprint]=stest(leny,sy,sty,sty,sfac,kprint); sty(dumvar4~=sty_orig)=dumvar4(dumvar4~=sty_orig);
end; kpar = fix(4+1);
elseif( combla_2==19 ) ;
kni = fix(kn + 4.*(ki-1));
for kpar = 1 : 4;
for i = 1 : 7;
dx(i) = dx1(i);
dy(i) = dy1(i);
dtx(i) = dt19x(i,kpar,kni);
dty(i) = dt19y(i,kpar,kni);
end; i = fix(7+1);
for i = 1 : 5;
dparam(i) = dpar(i,kpar);
end; i = fix(5+1);
combla_6 = fix(fix(dparam(1)));
for i = 1 : lenx;
dsize(i) = dtx(i);
end; i = fix(lenx+1);
if((kpar==2) &&(kni==7) )
dsize(1) = 2.4d0;
end;
if((kpar==3) &&(kni==8) )
dsize(5) = 1.8d0;
end;
[combla_3,dx,combla_4,dy,combla_5,dparam]=drotm(combla_3,dx,combla_4,dy,combla_5,dparam);
[lenx,dx,dtx,dsize,dfac,kprint]=dtest(lenx,dx,dtx,dsize,dfac,kprint);
dty_orig=dty; [leny,dy,dty,dumvar4,dfac,kprint]=dtest(leny,dy,dty,dty,dfac,kprint); dty(dumvar4~=dty_orig)=dumvar4(dumvar4~=dty_orig);
end; kpar = fix(4+1);
elseif( combla_2==20 ) ;
for i = 1 : 7;
sty(i) = dt10y(i,kn,ki);
end; i = fix(7+1);
[combla_3,sx,combla_4,sy,combla_5]=scopy(combla_3,sx,combla_4,sy,combla_5);
[leny,sy,sty,ssize2(sub2ind(size(ssize2),1,1):end),dumvar5,kprint]=stest(leny,sy,sty,ssize2(sub2ind(size(ssize2),1,1):end),1.,kprint);
elseif( combla_2==21 ) ;
[combla_3,dx,combla_4,dy,combla_5]=dcopy(combla_3,dx,combla_4,dy,combla_5);
[leny,dy,dt10y(sub2ind(size(dt10y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),dumvar5,kprint]=dtest(leny,dy,dt10y(sub2ind(size(dt10y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),1.0d0,kprint);
elseif( combla_2==22 ) ;
[combla_3,cx,combla_4,cy,combla_5]=ccopy(combla_3,cx,combla_4,cy,combla_5);
%Barrowes - commented out calls to stest with complex variables
%!! CALL STEST(2*leny,cy,ct10y(1,kn,ki),ssize2(1,1),1.,Kprint)
elseif( combla_2==23 ) ;
[combla_3,sx,combla_4,sy,combla_5]=sswap(combla_3,sx,combla_4,sy,combla_5);
for i = 1 : 7;
stx(i) = dt10x(i,kn,ki);
sty(i) = dt10y(i,kn,ki);
end; i = fix(7+1);
[lenx,sx,stx,ssize2(sub2ind(size(ssize2),1,1):end),dumvar5,kprint]=stest(lenx,sx,stx,ssize2(sub2ind(size(ssize2),1,1):end),1.,kprint);
[leny,sy,sty,ssize2(sub2ind(size(ssize2),1,1):end),dumvar5,kprint]=stest(leny,sy,sty,ssize2(sub2ind(size(ssize2),1,1):end),1.,kprint);
elseif( combla_2==24 ) ;
[combla_3,dx,combla_4,dy,combla_5]=dswap(combla_3,dx,combla_4,dy,combla_5);
[lenx,dx,dt10x(sub2ind(size(dt10x),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),dumvar5,kprint]=dtest(lenx,dx,dt10x(sub2ind(size(dt10x),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),1.0d0,kprint);
[leny,dy,dt10y(sub2ind(size(dt10y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),dumvar5,kprint]=dtest(leny,dy,dt10y(sub2ind(size(dt10y),1,kn,ki):end),dsize2(sub2ind(size(dsize2),1,1):end),1.0d0,kprint);
elseif( combla_2==25 ) ;
[combla_3,cx,combla_4,cy,combla_5]=cswap(combla_3,cx,combla_4,cy,combla_5);
%Barrowes - commented out calls to stest with complex variables
%!! CALL STEST(2*lenx,cx,ct10x(1,kn,ki),ssize2(1,1),1.,Kprint)
%!! CALL STEST(2*leny,cy,ct10y(1,kn,ki),ssize2(1,1),1.,Kprint)
else;
[dumvar1,dumvar2,dumvar3,ssize1(sub2ind(size(ssize1),max(kn,1)):end),sfac,kprint]=stest(1,sdot(combla_3,sx,combla_4,sy,combla_5),real(dt7(kn,ki)),ssize1(sub2ind(size(ssize1),max(kn,1)):end),sfac,kprint);
end;
end; kn = fix(4+1);
end; ki = fix(4+1);
return;
end %subroutine check2
|
|
