-
Notifications
You must be signed in to change notification settings - Fork 0
/
update_matrices.m
21 lines (19 loc) · 1.63 KB
/
update_matrices.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
function [ matrices ] = matrices_linearized_t_dep(t, matrices,params )
%MATRICES_LINEARIZED_T_DEP Summary of this function goes here
% Detailed explanation goes here
[matrices.adim.ephase.dl.Dmod,matrices.adim.ephase.dl.Emod,matrices.adim.ephase.dl.Umod,matrices.adim.ephase.dl.Vmod] = ...
(matrices_bdrys_wrt_internal(...
matrices.adim.ephase.D1_noBC(1,:),...
matrices.adim.ephase.D1_noBC(end,:),...
current_profile_adim(t,params),...
current_profile_adim(t,params),...
[0,0]...
));
matrices.adim.ephase.dl.D2_BC = ...
(matrices.adim.ephase.D2_noBC(2:end-1,2:end-1)...
+ matrices.adim.ephase.D2_noBC(2:end-1,1) * matrices.adim.ephase.dl.Dmod...
+ matrices.adim.ephase.D2_noBC(2:end-1,end) * matrices.adim.ephase.dl.Emod);
matrices.adim.ephase.dl.exogeneous_bdry = ...
(matrices.adim.ephase.D2_noBC(2:end-1,1) * matrices.adim.ephase.dl.Umod...
+ matrices.adim.ephase.D2_noBC(2:end-1,end) * matrices.adim.ephase.dl.Vmod);
end