Skip to content
Snippets Groups Projects
Commit a1c7b565 authored by Ivar Stefansson's avatar Ivar Stefansson
Browse files

Delete EvaluateAllFractures.m

parent 0df36c1b
No related branches found
No related tags found
No related merge requests found
% Evaluates the norm in each of the fractures and adds the computed values
% to the input file. This file must contain a cell array with fracxandp for
% each of the fractures, see the Description for details on how each
% fracxandp should be constructed.
clear
load geiger_permeable_reference.mat
load anna_reference_lefttoright.mat
% load geiger_blocking_reference.mat
%load Hydrocoin_reference.mat
%coarsefile = 'coarse file.mat';
%coarsefile = '../hydrocoin/results/edfm_fractures.mat';
coarsefile = '../geiger/results/permeable_edfm.mat';
%coarsefile = '../geiger/results/blocking_edfm.mat';
coarsefile = '/home/elle/Dropbox/2015-fractures-comparison/new_case_4/results_bis/edfm_fractures.mat';
load(coarsefile)
nfrac = length(fracxandp);
Erel_f_vec = zeros(nfrac,1);
E2_f = zeros(nfrac,1);
dP2_f = zeros(nfrac,1);
FRAC = true;
%%
for i = 1:nfrac
xandp = fracxandp{i};
[Erel_f_vec(i),E2_f(i),dP2_f(i)] = LineNormFucntion(xandp,p,t,x,FRAC);
end
Erel_f = sqrt(sum(E2_f)/sum(dP2_f))
save(coarsefile,'Erel_f','E2_f','dP2_f','-append')
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment