Skip to content
Snippets Groups Projects
Commit f7aae77f authored by Marco Ratto's avatar Marco Ratto
Browse files

removed useless elements + reduce loops for efficiency

(cherry picked from commit f2271264)
parent 3ddafb16
No related branches found
No related tags found
No related merge requests found
...@@ -181,13 +181,13 @@ d2P1 = zeros(ns,ns,k,k); ...@@ -181,13 +181,13 @@ d2P1 = zeros(ns,ns,k,k);
for ii = 1:k for ii = 1:k
dAi = dA(:,:,ii); dAi = dA(:,:,ii);
dFi = dP0(Z,Z,ii); dFi = dP0(Z,Z,ii);
d2Omi = d2Om(:,:,ii); % d2Omi = d2Om(:,:,ii);
diFi = -iF*dFi*iF; diFi = -iF*dFi*iF;
dKi = dK0(:,:,ii); dKi = dK0(:,:,ii);
for jj = 1:k for jj = 1:ii,
dAj = dA(:,:,jj); dAj = dA(:,:,jj);
dFj = dP0(Z,Z,jj); dFj = dP0(Z,Z,jj);
d2Omj = d2Om(:,:,jj); % d2Omj = d2Om(:,:,jj);
dFj = dP0(Z,Z,jj); dFj = dP0(Z,Z,jj);
diFj = -iF*dFj*iF; diFj = -iF*dFj*iF;
dKj = dK0(:,:,jj); dKj = dK0(:,:,jj);
...@@ -218,6 +218,9 @@ for ii = 1:k ...@@ -218,6 +218,9 @@ for ii = 1:k
d2K(:,:,ii,jj) = d2Kij; %#ok<NASGU> d2K(:,:,ii,jj) = d2Kij; %#ok<NASGU>
d2P1(:,:,ii,jj) = d2AtmpA + d2Omij; %#ok<*NASGU> d2P1(:,:,ii,jj) = d2AtmpA + d2Omij; %#ok<*NASGU>
d2S(:,:,ii,jj) = d2Fij; d2S(:,:,ii,jj) = d2Fij;
d2K(:,:,jj,ii) = d2Kij; %#ok<NASGU>
d2P1(:,:,jj,ii) = d2AtmpA + d2Omij; %#ok<*NASGU>
d2S(:,:,jj,ii) = d2Fij;
% d2iS(:,:,ii,jj) = d2iF; % d2iS(:,:,ii,jj) = d2iF;
end end
end end
...@@ -244,13 +247,13 @@ d2P1 = zeros(ns,ns,k,k); ...@@ -244,13 +247,13 @@ d2P1 = zeros(ns,ns,k,k);
for ii = 1:k for ii = 1:k
dAi = dA(:,:,ii); dAi = dA(:,:,ii);
dFi = Z*dP0(:,:,ii)*Z; dFi = Z*dP0(:,:,ii)*Z;
d2Omi = d2Om(:,:,ii); % d2Omi = d2Om(:,:,ii);
diFi = -iF*dFi*iF; diFi = -iF*dFi*iF;
dKi = dK0(:,:,ii); dKi = dK0(:,:,ii);
for jj = 1:k for jj = 1:ii,
dAj = dA(:,:,jj); dAj = dA(:,:,jj);
dFj = Z*dP0(:,:,jj)*Z; dFj = Z*dP0(:,:,jj)*Z;
d2Omj = d2Om(:,:,jj); % d2Omj = d2Om(:,:,jj);
dFj = Z*dP0(:,:,jj)*Z; dFj = Z*dP0(:,:,jj)*Z;
diFj = -iF*dFj*iF; diFj = -iF*dFj*iF;
dKj = dK0(:,:,jj); dKj = dK0(:,:,jj);
...@@ -282,6 +285,9 @@ for ii = 1:k ...@@ -282,6 +285,9 @@ for ii = 1:k
d2P1(:,:,ii,jj) = d2AtmpA + d2Omij; %#ok<*NASGU> d2P1(:,:,ii,jj) = d2AtmpA + d2Omij; %#ok<*NASGU>
d2S(:,:,ii,jj) = d2Fij; d2S(:,:,ii,jj) = d2Fij;
% d2iS(:,:,ii,jj) = d2iF; % d2iS(:,:,ii,jj) = d2iF;
d2K(:,:,jj,ii) = d2Kij; %#ok<NASGU>
d2P1(:,:,jj,ii) = d2AtmpA + d2Omij; %#ok<*NASGU>
d2S(:,:,jj,ii) = d2Fij;
end end
end end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment