From 457bc833d0409061b8100349d16e8a5a1bcb9757 Mon Sep 17 00:00:00 2001 From: Craig Warren Date: Tue, 11 Oct 2016 13:25:36 +0100 Subject: [PATCH] Added check when processing materials to set averagable to False if conductivity is infinite, i.e. pec --- gprMax/materials.py | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/gprMax/materials.py b/gprMax/materials.py index 51bb157c..5361c7a8 100644 --- a/gprMax/materials.py +++ b/gprMax/materials.py @@ -157,6 +157,10 @@ def process_materials(G): materialsdata = [['\nID', '\nName', '\nType', '\neps_r', 'sigma\n[S/m]', '\nDelta eps_r', 'tau\n[s]', 'omega\n[Hz]', 'delta\n[Hz]', 'gamma\n[Hz]', '\nmu_r', 'sigma*\n[S/m]', 'Dielectric\nsmoothable']] for material in G.materials: + # Set material averaging to False if infinite conductivity, i.e. pec + if material.se == float('inf'): + material.averagable = False + # Calculate update coefficients for material material.calculate_update_coeffsE(G) material.calculate_update_coeffsH(G)