function [ relative_errors ] = calculate_relative_errors( results_EDFM_matrix, results_pEDFM_case1, results_pEDFM_case2, r_EDFMs)
%CALCULATE_RELATIVE_ERRORS 此处显示有关此函数的摘要
%   此处显示详细说明
% % % EDFM_result = results_EDFM_matrix(:,1);
% % % pEDFM_result_case1 = results_pEDFM_case1(:,2);
% % % modified_EDFM_result = results_pEDFM_case2(:,3);
% % % new_EDFM_result = results_EDFMs_matrix(:,4);
% % % LGR_result = results_LGR(:,1);

EDFMs_cell_mid_point_coordinates = r_EDFMs.cell_mid_coordinates;
% LGR_cell_mid_point_coordinates = r_LGR.cell_mid_coordinates;

num_of_cells = size(EDFMs_cell_mid_point_coordinates,1);
% the_reference_solution_for_comparison = zeros(num_of_cells,1);

% % % for i = 1:num_of_cells
% % %     [a,b]=min((LGR_cell_mid_point_coordinates(:,1)- EDFMs_cell_mid_point_coordinates(i,1)).^2+(LGR_cell_mid_point_coordinates(:,2)- EDFMs_cell_mid_point_coordinates(i,2)).^2);
% % %     the_reference_solution_for_comparison(i,1) = LGR_result(b(1),1);
% % % end
% 相对误差计算
fracture_perm = [ 100, 700, 4000, 20000, 100000];
relative_errors = zeros(2,5);
for j = 1:5
relative_errors(1,j) = norm(results_EDFM_matrix(1:num_of_cells,j)-results_pEDFM_case2(1:num_of_cells,j))/norm(results_pEDFM_case2(1:num_of_cells,j));
relative_errors(2,j) = norm(results_pEDFM_case1(1:num_of_cells,j)-results_pEDFM_case2(1:num_of_cells,j))/norm(results_pEDFM_case2(1:num_of_cells,j));

% relative_errors(1,2) = norm(pEDFM_result(1:num_of_cells,1)-the_reference_solution_for_comparison(:,1))/norm(the_reference_solution_for_comparison(:,1));
% relative_errors(1,3) = norm(modified_EDFM_result(1:num_of_cells,1)-the_reference_solution_for_comparison(:,1))/norm(the_reference_solution_for_comparison(:,1));
% relative_errors(1,4) = norm(new_EDFM_result(1:num_of_cells,1)-the_reference_solution_for_comparison(:,1))/norm(the_reference_solution_for_comparison(:,1));
end
% relative_errors(1,5) = norm(LGR_result(:,1)-the_reference_solution_for_comparison(:,1))/norm(the_reference_solution_for_comparison(:,1));


end