Skip to content

Commit

Permalink
zef_nse_sigma: fix indentations with Matlab's smart indent tool
Browse files Browse the repository at this point in the history 
- Use Command/Control + A and Command/Control + I in a code file to run the
  indentation tool.

- This should be done before every commit.
  • Loading branch information
@SeSodesa
SeSodesa committed May 2, 2023
1 parent 46d70c5 commit 2241fdd
Showing 1 changed file with 56 additions and 56 deletions.
112 changes: 56 additions & 56 deletions m/forward_simulation/nse/zef_nse_sigma.m
View file
Expand Up @@ -12,63 +12,63 @@
[~,det] = zef_volume_barycentric(v_nodes,v_tetra);
volume = abs(det)/6;
sigma_aux = sigma_in(active_compartment_ind,1);

if isequal(nse_field.conductivity_model,1)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
I_aux = find(abs(bf_interp(I))> 1e-15);
I = I(I_aux);
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
sigma_out(I_aux,1) = nse_field.blood_conductivity;
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
bf_interp(I(I_aux)) = 1-1e-15;
I_aux = find(abs(bf_interp(I)) <= 1);
I = I(I_aux);
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
m = nse_field.conductivity_exponent;
s = log(1-bf_interp(I).^m)./log(1-bf_interp(I));
sigma_out(active_compartment_ind(I),1) = (1-bf_interp(I)).^s.*background_conductivity + nse_field.blood_conductivity .* bf_interp(I).^m;
end
if isequal(nse_field.conductivity_model,1)

elseif isequal(nse_field.conductivity_model,2)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
%sigma_out(active_compartment_ind(I),1) = nse_field.blood_conductivity.*(1 - 3.*(1-bf_interp(I)).*(nse_field.blood_conductivity-background_conductivity)./(3.*nse_field.blood_conductivity - bf_interp(I).*(nse_field.blood_conductivity-background_conductivity)));
sigma_out(active_compartment_ind(I),1) = background_conductivity.*(1-bf_interp(I)) + nse_field.blood_conductivity.*(bf_interp(I));
end

elseif isequal(nse_field.conductivity_model,3)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
I_aux = find(abs(bf_interp(I))> 1e-15);
I = I(I_aux);
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
sigma_out(I_aux,1) = nse_field.blood_conductivity;
for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
I_aux = find(abs(bf_interp(I))> 1e-15);
I = I(I_aux);
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
bf_interp(I(I_aux)) = 1-1e-15;
I_aux = find(abs(bf_interp(I)) <= 1);
I = I(I_aux);
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
sigma_out(active_compartment_ind(I),1) = background_conductivity.*(1 + (3.*bf_interp(I).*(nse_field.blood_conductivity-background_conductivity))./(3.*background_conductivity + (1 - bf_interp(I)).*(nse_field.blood_conductivity-background_conductivity)));
%sigma_out(active_compartment_ind(I),1) = 1./(background_conductivity./(1-bf_interp(I)) + nse_field.blood_conductivity./bf_interp(I));
end

elseif isequal(nse_field.conductivity_model,4)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
sigma_out(active_compartment_ind(I),1) = (nse_field.blood_conductivity + (background_conductivity - nse_field.blood_conductivity).*(1-(2.*bf_interp(I)./3)))./(1 + (bf_interp(I)./3).*(background_conductivity./nse_field.blood_conductivity-1));
sigma_out(I_aux,1) = nse_field.blood_conductivity;
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
bf_interp(I(I_aux)) = 1-1e-15;
I_aux = find(abs(bf_interp(I)) <= 1);
I = I(I_aux);
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
m = nse_field.conductivity_exponent;
s = log(1-bf_interp(I).^m)./log(1-bf_interp(I));
sigma_out(active_compartment_ind(I),1) = (1-bf_interp(I)).^s.*background_conductivity + nse_field.blood_conductivity .* bf_interp(I).^m;
end

elseif isequal(nse_field.conductivity_model,2)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
%sigma_out(active_compartment_ind(I),1) = nse_field.blood_conductivity.*(1 - 3.*(1-bf_interp(I)).*(nse_field.blood_conductivity-background_conductivity)./(3.*nse_field.blood_conductivity - bf_interp(I).*(nse_field.blood_conductivity-background_conductivity)));
sigma_out(active_compartment_ind(I),1) = background_conductivity.*(1-bf_interp(I)) + nse_field.blood_conductivity.*(bf_interp(I));
end

elseif isequal(nse_field.conductivity_model,3)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
I_aux = find(abs(bf_interp(I))> 1e-15);
I = I(I_aux);
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
sigma_out(I_aux,1) = nse_field.blood_conductivity;
I_aux = find(abs(bf_interp(I)) > 1-1e-15);
bf_interp(I(I_aux)) = 1-1e-15;
I_aux = find(abs(bf_interp(I)) <= 1);
I = I(I_aux);
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
sigma_out(active_compartment_ind(I),1) = background_conductivity.*(1 + (3.*bf_interp(I).*(nse_field.blood_conductivity-background_conductivity))./(3.*background_conductivity + (1 - bf_interp(I)).*(nse_field.blood_conductivity-background_conductivity)));
%sigma_out(active_compartment_ind(I),1) = 1./(background_conductivity./(1-bf_interp(I)) + nse_field.blood_conductivity./bf_interp(I));
end

end

end

elseif isequal(nse_field.conductivity_model,4)

for i = 1 : length(nse_field.capillary_domain_ind)
I = find(domain_labels_aux == nse_field.capillary_domain_ind(i));
volume_aux = sum(volume(I));
background_conductivity = sum(sigma_aux(I).*volume(I))./volume_aux;
sigma_out(active_compartment_ind(I),1) = (nse_field.blood_conductivity + (background_conductivity - nse_field.blood_conductivity).*(1-(2.*bf_interp(I)./3)))./(1 + (bf_interp(I)./3).*(background_conductivity./nse_field.blood_conductivity-1));
end

end

end

0 comments on commit 2241fdd

Please sign in to comment.