% Bree Aldridge, 2/15/2006
% script to integrate over the phase-space grid and save the end trajectory
% locations

%phase-space grid to launch trajectories where the first two numbers
%indicate the span of the grid in each direction while the last bumber is
%the numgber of point taken along that direction.  
c8var=linspace(100,350000,9);   %casp-8
c8avar=linspace(100,100000,9);  %casp-8*
c3var=linspace(100,350000,9);   %casp-3
c3avar=linspace(100,100000,9);  %casp-3*
c8c3avar=linspace(100,25000,3); %casp-8:casp-3*
c8ac3var=c8c3avar;              %casp-8*:casp-3
xiapvar=linspace(100,100000,9); %XIAP
c3axvar=c8c3avar;               %casp-3*:XIAP

%time span - ending at 6 hours
t=linspace(0,6*60*60,500);

%solve odes across the grid.
%programming note: The new Matlab versions (6.5 and up) have greatly
%improved for-loop performance.  We have implemented a vectorized version
%of this script, which ended up slowing down the run-time (presumibly
%because of memory constraints).  These for loops are easily
%parallelizable if increased grid resolution is required.  
IC=[c8var(1) c8avar(1) c3var(1) c3avar(1) c8c3avar(1) c8ac3var(1) xiapvar(1) c3axvar(1)]; %initial conditions
for a=1:length(c8var)
    IC(1)=c8var(a);
    for b=1:length(c8avar)
        IC(2)=c8avar(b);
        for c=1:length(c3var)
            IC(3)=c3var(c);
            for d=1:length(c3avar)
                IC(4)=c3avar(d);
                for e=1:length(c8c3avar)
                    IC(5)=c8c3avar(e);
                    for f=1:length(c8ac3var)
                        IC(6)=c8ac3var(f);
                        for g=1:length(xiapvar)
                            IC(7)=xiapvar(g);
                            for h=1:length(c3axvar)
                                IC(8)=c3axvar(h);                                
                                [T,Y]=ode15s(@casp3modelodes,t,IC,[]);
                                stor(a,b,c,d,e,f,g,h,:)=Y(end,:); %store the end locations of the trajectories
                            end
                        end
                    end
                end
            end
        end
    end
end

%saves the .mat (matlab workspace) including the grid and ending trajectory
%locations
save integration_results;