你已经派生过 gprMax
镜像自地址
https://gitee.com/sunhf/gprMax.git
已同步 2025-08-07 04:56:51 +08:00
81 行
2.9 KiB
Python
81 行
2.9 KiB
Python
# Copyright (C) 2015-2023: The University of Edinburgh, United Kingdom
|
|
# Authors: Craig Warren, Antonis Giannopoulos, and John Hartley
|
|
#
|
|
# This file is part of gprMax.
|
|
#
|
|
# gprMax is free software: you can redistribute it and/or modify
|
|
# it under the terms of the GNU General Public License as published by
|
|
# the Free Software Foundation, either version 3 of the License, or
|
|
# (at your option) any later version.
|
|
#
|
|
# gprMax is distributed in the hope that it will be useful,
|
|
# but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
# GNU General Public License for more details.
|
|
#
|
|
# You should have received a copy of the GNU General Public License
|
|
# along with gprMax. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
import itertools
|
|
import logging
|
|
from operator import add
|
|
from pathlib import Path
|
|
|
|
import matplotlib.gridspec as gridspec
|
|
import matplotlib.pyplot as plt
|
|
import numpy as np
|
|
|
|
from testing.diff_output_files import diff_output_files
|
|
|
|
logger = logging.getLogger(__name__)
|
|
|
|
# Create/setup plot figure
|
|
# Plot colours from http://tools.medialab.sciences-po.fr/iwanthue/index.php
|
|
colorIDs = ["#79c72e", "#5774ff", "#ff7c2c", "#4b4e80", "#d7004e", "#007545"]
|
|
colors = itertools.cycle(colorIDs)
|
|
lines = itertools.cycle(("--", ":", "-.", "-"))
|
|
markers = ["o", "d", "^", "s", "*"]
|
|
|
|
fn = Path(__file__)
|
|
basename = "pml_basic"
|
|
PMLIDs = ["off", "x0", "y0", "z0", "xmax", "ymax", "zmax"]
|
|
maxerrors = []
|
|
|
|
for x, PMLID in enumerate(PMLIDs):
|
|
file1 = fn.parent.joinpath(basename + str(x + 1) + "_CPU.h5")
|
|
file2 = fn.parent.joinpath(basename + str(x + 1) + "_GPU.h5")
|
|
time, datadiffs = diff_output_files(file1, file2)
|
|
|
|
# Print maximum error value
|
|
start = 150
|
|
maxerrors.append(f": {np.amax(np.amax(datadiffs[start::, :])):.1f} [dB]")
|
|
print(f"{PMLID}: Max. error {maxerrors[x]}")
|
|
|
|
# Plot diffs
|
|
gs = gridspec.GridSpec(3, 2, hspace=0.3, wspace=0.3)
|
|
fig, ax = plt.subplots(figsize=(20, 10), facecolor="w", edgecolor="w")
|
|
ax.remove()
|
|
fig.suptitle(f"{PMLID}")
|
|
|
|
outputs = ["Ex", "Ey", "Ez", "Hx", "Hy", "Hz"]
|
|
for i, output in enumerate(outputs):
|
|
if i < 3:
|
|
ax = plt.subplot(gs[i, 0])
|
|
else:
|
|
ax = plt.subplot(gs[i - 3, 1])
|
|
ax.plot(time[start::], datadiffs[start::, i], color=next(colors), lw=2)
|
|
ax.set_xticks(np.arange(0, 1800, step=200))
|
|
ax.set_xlim([0, 1600])
|
|
ax.set_yticks(np.arange(-400, 80, step=40))
|
|
ax.set_ylim([-400, 40])
|
|
ax.set_axisbelow(True)
|
|
ax.grid(color=(0.75, 0.75, 0.75), linestyle="dashed")
|
|
ax.set_xlabel("Time [iterations]")
|
|
ax.set_ylabel(f"{output} error [dB]")
|
|
|
|
# Save a PDF/PNG of the figure
|
|
fig.savefig(basename + "_diffs_" + PMLID + ".pdf", dpi=None, format="pdf", bbox_inches="tight", pad_inches=0.1)
|
|
# fig.savefig(basename + "_diffs_" + PMLID + ".png", dpi=150, format='png', bbox_inches='tight', pad_inches=0.1)
|
|
|
|
plt.show()
|