Contents
Construct solution vector x0 from block variables
x0 = zeros(sum(ntv), 1);
idx = 1;
for j = 1:Nblock
for i = 1:NVAR(j)
variable_name = list_var{j}{i};
len = nt(j);
x0(idx:idx+len-1) = reshape(eval(sprintf('%s%d', variable_name, j)), len, 1);
idx = idx + len;
end
end
Construct full solution vector x0full
x0full = zeros(ntA * NVA, 1);
idx = 1;
for j = 1:Nblock
for i = 1:NVAR(j)
variable_name = list_var{j}{i};
x0full(idx:idx+ntA-1) = reshape(eval(sprintf('%s%dfull', variable_name, j)), ntA, 1);
idx = idx + ntA;
end
end