Skip to content
Snippets Groups Projects
Select Git revision
  • master default protected
  • revert-c23056ed
2 results
Blame Permalink
plot_conductive_y.m 1.32 KiB 
close all;
clear all;
clc;
hold on;

% equi-dimensional MFD
mfd = csvread('mfd/mfd_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(mfd(:, 6), mfd(:, 1), 'Color', 'k', 'LineWidth', 1);

% Box-DFM
dfm = csvread('boxdfm/boxdfm_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(dfm(:, 14), dfm(:, 4), 'Color', 'g', 'LineWidth', 1);

% CC-DFM
ccdfm = csvread('ccdfm/ccdfm_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(ccdfm(:, 4), ccdfm(:, 1), 'Color', 'c', 'LineWidth', 1);

% EDFM
edfm = csvread('edfm/edfm_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(edfm(:, 4), edfm(:, 1), 'Color', 'r', 'LineWidth', 1);

% Mortar-DFM
mortar = csvread('mortardfm/mortardfm_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(mortar(:, 4), mortar(:, 1), 'Color', [1 0.7 0], 'LineWidth', 1);

% P-XFEM
p_xfem = csvread('pxfem/pxfem_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(p_xfem(:, 3), p_xfem(:, 1), 'Color', 'b', 'LineWidth', 1);

% D-XFEM
d_xfem_fine = csvread('dxfem/dxfem_geiger_conductive_y_eq_0.7.csv', 1, 0);
plot(d_xfem_fine(:, 3), d_xfem_fine(:, 1), 'Color', 'm', 'LineWidth', 1);

% format the plot
ylim([1 1.51]);
xlabel('arc length', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
ylabel('pressure', 'FontSize', 12, 'FontWeight', 'bold', 'Color', 'k');
set(gca, 'FontSize', 12);
legend('reference', 'Box-DFM', 'CC-DFM', 'EDFM', 'mortar-DFM', 'P-XFEM', 'D-XFEM');