PyNEC 1.8.1__tar.gz → 2.2.1__tar.gz

{pynec-1.8.1 → pynec-2.2.1}/MANIFEST.in

@@ -2,7 +2,7 @@ include LICENCE.txt
 
 include necpp_src/src/*.h
 
-include necpp_src/config.h
+include necpp_src/src/config.h.in
 
 include example/*.py
 

{pynec-1.8.1 → pynec-2.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PyNEC
-Version: 1.8.1
+Version: 2.2.1
 Summary: Python Antenna Simulation Module (nec2++) object-oriented interface
 Author-email: Tim Molteno <tim@physics.otago.ac.nz>
 License-Expression: GPL-3.0-only
@@ -16,8 +16,12 @@ Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENCE.txt
 Requires-Dist: numpy
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Requires-Dist: pytest-cov; extra == "test"
 Dynamic: license-file
 
+<!-- Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz) -->
 # Python NEC2++ Module
 
 This module wraps the C++ API for antenna simulation of nec2++. It is easier to work with, and more powerful than the C-style API wrapper. Works with Python 2.7 and 3+.

{pynec-1.8.1 → pynec-2.2.1}/PyNEC.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PyNEC
-Version: 1.8.1
+Version: 2.2.1
 Summary: Python Antenna Simulation Module (nec2++) object-oriented interface
 Author-email: Tim Molteno <tim@physics.otago.ac.nz>
 License-Expression: GPL-3.0-only
@@ -16,8 +16,12 @@ Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENCE.txt
 Requires-Dist: numpy
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Requires-Dist: pytest-cov; extra == "test"
 Dynamic: license-file
 
+<!-- Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz) -->
 # Python NEC2++ Module
 
 This module wraps the C++ API for antenna simulation of nec2++. It is easier to work with, and more powerful than the C-style API wrapper. Works with Python 2.7 and 3+.

{pynec-1.8.1 → pynec-2.2.1}/PyNEC.egg-info/SOURCES.txt

@@ -29,7 +29,6 @@ example/test_rp2.py
 example/test_se.py
 example/test_structure_currents.py
 example/test_surface_patch_currents.py
-necpp_src/config.h
 necpp_src/src/AntennaInput.h
 necpp_src/src/BaseInput.h
 necpp_src/src/CurrentInput.h
@@ -46,6 +45,7 @@ necpp_src/src/c_ggrid.h
 necpp_src/src/c_plot_card.cpp
 necpp_src/src/c_plot_card.h
 necpp_src/src/common.h
+necpp_src/src/config.h.in
 necpp_src/src/electromag.cpp
 necpp_src/src/electromag.h
 necpp_src/src/libNEC.cpp
@@ -75,5 +75,13 @@ necpp_src/src/nec_structure_currents.h
 necpp_src/src/nec_wire.h
 necpp_src/src/safe_array.h
 necpp_src/src/typesafe_stdint.h
+tests/test_charge_densities.py
 tests/test_examples.py
-tests/test_get_gain.py
+tests/test_get_gain.py
+tests/test_input_parameters.py
+tests/test_near_field.py
+tests/test_norm_rx_pattern.py
+tests/test_patches.py
+tests/test_radiation_pattern.py
+tests/test_structure_currents.py
+tests/test_structure_excitation.py

pynec-2.2.1/PyNEC.egg-info/requires.txt

@@ -0,0 +1,5 @@
+numpy
+
+[test]
+pytest
+pytest-cov

{pynec-1.8.1 → pynec-2.2.1}/README.md

@@ -1,3 +1,4 @@
+<!-- Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz) -->
 # Python NEC2++ Module
 
 This module wraps the C++ API for antenna simulation of nec2++. It is easier to work with, and more powerful than the C-style API wrapper. Works with Python 2.7 and 3+.

{pynec-1.8.1 → pynec-2.2.1}/example/antenna_util.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #
 # Some antenna utility functions
 #

{pynec-1.8.1 → pynec-2.2.1}/example/basic_usage.py

@@ -1,4 +1,5 @@
 #!/usr/bin/env python3
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 """
 Basic PyNEC usage example — demonstrates:
 

{pynec-1.8.1 → pynec-2.2.1}/example/context_clean.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 # Note: explicit zeroes are blanks. All other values should be specified symbolically.
 
 # Currently these contain only the subset of cards that I needed

{pynec-1.8.1 → pynec-2.2.1}/example/dipole.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 import numpy as np
 import scipy.optimize
 import pylab as plt

{pynec-1.8.1 → pynec-2.2.1}/example/impedance_plot.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #
 #  Plot of reflection coefficient vs antenna length for a fixed base height.
 #

{pynec-1.8.1 → pynec-2.2.1}/example/logperiodic_opt.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 import numpy as np
 import scipy.optimize
 import matplotlib.pyplot as plt

{pynec-1.8.1 → pynec-2.2.1}/example/monopole.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #
 #  Simple vertical monopole antenna simulation using python-necpp
 #  pip install PyNEC

{pynec-1.8.1 → pynec-2.2.1}/example/monopole_realistic_ground_plane.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 import numpy as np
 import scipy.optimize
 import matplotlib.pyplot as plt

{pynec-1.8.1 → pynec-2.2.1}/example/optimized.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #
 #  Automatically tune antenna
 #

{pynec-1.8.1 → pynec-2.2.1}/example/radiation_pattern.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #
 #  Simple vertical monopole antenna simulation using python-necpp
 #  pip install necpp

{pynec-1.8.1 → pynec-2.2.1}/example/test_ai.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_charge_densities.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_ne_nh.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_nrp.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_rp.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 
 
 from PyNEC import *

{pynec-1.8.1 → pynec-2.2.1}/example/test_rp2.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 from PyNEC import *
 
 #creation of a nec context

{pynec-1.8.1 → pynec-2.2.1}/example/test_se.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_structure_currents.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/example/test_surface_patch_currents.py

@@ -1,3 +1,4 @@
+# Copyright (c) 2008-2026 Tim Molteno (tim@elec.ac.nz)
 #instructions to compile :
 #
 #swig -c++ -python PyNEC.i 

{pynec-1.8.1 → pynec-2.2.1}/necpp_src/src/c_evlcom.cpp

@@ -186,7 +186,7 @@ void c_evlcom::gshank( nec_complex start, nec_complex dela, complex_array& sum,
 				if (den < denm)
 					den=denm;
 				a1=q1[i][j]-a1;
-				amg=fabs(real(a1)+fabs(imag(a1)));
+				amg=fabs(real(a1))+fabs(imag(a1));
 				if (amg > den)
 				{
 					brk = true;
@@ -563,54 +563,52 @@ void c_evlcom::evlua( nec_complex *erv, nec_complex *ezv,
 	gshank(cp1,delta,ans,6,sum,0,bk,bk);
 	rmis=m_rho*(real(m_ck1)-m_ck2);
 	
-	bool jump = false;
-	if ( (rmis >= 2.*m_ck2) && (m_rho >= 1.e-10) )
-	{
-		if (m_zph >= 1.e-10)
-		{
-			bk=nec_complex(-m_zph,m_rho)*(m_ck1-cp3);
-			rmis = -real(bk)/fabs(imag(bk));
-			if (rmis > 4.*m_rho/m_zph)
-				jump = true;
-		}
-		
-		if ( false == jump )
-		{
-			/* integrate up between branch cuts, then to + infinity */
-			cp1=m_ck1- nec_complex(0.1,+0.2);
-			cp2=cp1+.2;
-			bk=nec_complex(0.,del);
-			gshank(cp1,bk,sum,6,ans,0,bk,bk);
-			m_contour_a=cp1;
-			m_contour_b=cp2;
-			rom1(6,ans,1);
-			for (int i = 0; i < 6; i++ )
-				ans[i] -= sum[i];
-			
-			gshank(cp3,bk,sum,6,ans,0,bk,bk);
-			gshank(cp2,delta2,ans,6,sum,0,bk,bk);
-		}
-		
-		jump = true;	
-	} /* if ( (rmis >= 2.*m_ck2) || (m_rho >= 1.e-10) ) */
-	else
+	bool jump;
+	if ( (rmis < 2.*m_ck2) || (m_rho < 1.e-10) )
+		jump = true;
+	else if ( m_zph < 1.e-10 )
 		jump = false;
-	
+	else
+	{
+		bk=nec_complex(-m_zph,m_rho)*(m_ck1-cp3);
+		rmis = -real(bk)/fabs(imag(bk));
+		if (rmis > 4.*m_rho/m_zph)
+			jump = true;
+		else
+			jump = false;
+	}
+
 	if ( false == jump )
+	{
+		/* integrate up between branch cuts, then to + infinity */
+		cp1=m_ck1- nec_complex(0.1,+0.2);
+		cp2=cp1+.2;
+		bk=nec_complex(0.,del);
+		gshank(cp1,bk,sum,6,ans,0,bk,bk);
+		m_contour_a=cp1;
+		m_contour_b=cp2;
+		rom1(6,ans,1);
+		for (int i = 0; i < 6; i++ )
+			ans[i] -= sum[i];
+
+		gshank(cp3,bk,sum,6,ans,0,bk,bk);
+		gshank(cp2,delta2,ans,6,sum,0,bk,bk);
+	}
+	else
 	{
 		/* integrate below branch points, then to + infinity */
 		for (int i = 0; i < 6; i++ )
 			sum[i] = -ans[i];
-		
+
 		rmis=real(m_ck1)*1.01;
 		if ( (m_ck2+1.) > rmis )
 			rmis=m_ck2+1.;
-		
+
 		bk=nec_complex(rmis,.99*imag(m_ck1));
 		delta=bk-cp3;
 		delta *= del/abs(delta);
 		gshank(cp3,delta,ans,6,sum,1,bk,delta2);
-	} /* if ( false == jump ) */
+	}
 	
 	ans[5] *= m_ck1;
 	
@@ -841,3 +839,45 @@ void hankel( nec_complex z, nec_complex *h0, nec_complex *h0p )
 }
 
 
+
+/* Numerical integration convergence tests used by Sommerfeld integration. */
+
+void test(
+        nec_float f1r, nec_float f2r, nec_float *tr,
+        nec_float f1i, nec_float f2i, nec_float *ti,
+        nec_float dmin )
+{
+    static nec_float _min_val = 1.0e-37;
+    
+    nec_float den = fabs( f2r);
+    nec_float temp_tr = fabs( f2i);
+    
+    if( den < temp_tr)
+        den = temp_tr;
+    if( den < dmin)
+        den = dmin;
+    
+    if( den < _min_val) {
+        *tr = 0.0;
+        *ti = 0.0;
+        return;
+    }
+    
+    *tr= fabs((f1r - f2r)/ den);
+    *ti= fabs((f1i - f2i)/ den); 
+}
+
+nec_float test_simple( nec_float f1r, nec_float f2r, nec_float dmin )
+{
+    static nec_float _min_val = 1.0e-37;
+    
+    nec_float den = fabs(f2r);
+    
+    if( den < dmin)
+        den = dmin;
+    if (den < _min_val)   {
+        return 0.0;
+    }
+    
+    return fabs((f1r - f2r) / den);    
+}

{pynec-1.8.1 → pynec-2.2.1}/necpp_src/src/c_evlcom.h

@@ -61,3 +61,8 @@ private:
 	/*! \brief Flag to select Bessel or Hankel function form (was jh) */
 	bool m_bessel_flag;
 };
+
+/* Numerical integration convergence tests (used by Sommerfeld integration) */
+void test(nec_float f1r, nec_float f2r, nec_float *tr, nec_float f1i,
+  nec_float f2i, nec_float *ti, nec_float dmin);
+nec_float test_simple( nec_float f1r, nec_float f2r, nec_float dmin );