PyNEC 1.7.4__tar.gz → 1.8.1__tar.gz

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

@@ -1,16 +1,15 @@
 Metadata-Version: 2.4
 Name: PyNEC
-Version: 1.7.4
+Version: 1.8.1
 Summary: Python Antenna Simulation Module (nec2++) object-oriented interface
 Author-email: Tim Molteno <tim@physics.otago.ac.nz>
-License: GPLv3
+License-Expression: GPL-3.0-only
 Project-URL: Homepage, http://github.com/tmolteno/python-necpp
 Project-URL: Repository, http://github.com/tmolteno/python-necpp
 Keywords: nec2,nec2++,antenna,electromagnetism,radio
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Topic :: Scientific/Engineering
 Classifier: Topic :: Communications :: Ham Radio
-Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
 Classifier: Programming Language :: Python :: 3
 Classifier: Intended Audience :: Science/Research
 Requires-Python: >=3.9

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

@@ -1,16 +1,15 @@
 Metadata-Version: 2.4
 Name: PyNEC
-Version: 1.7.4
+Version: 1.8.1
 Summary: Python Antenna Simulation Module (nec2++) object-oriented interface
 Author-email: Tim Molteno <tim@physics.otago.ac.nz>
-License: GPLv3
+License-Expression: GPL-3.0-only
 Project-URL: Homepage, http://github.com/tmolteno/python-necpp
 Project-URL: Repository, http://github.com/tmolteno/python-necpp
 Keywords: nec2,nec2++,antenna,electromagnetism,radio
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Topic :: Scientific/Engineering
 Classifier: Topic :: Communications :: Ham Radio
-Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
 Classifier: Programming Language :: Python :: 3
 Classifier: Intended Audience :: Science/Research
 Requires-Python: >=3.9

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

@@ -11,6 +11,7 @@ PyNEC.egg-info/dependency_links.txt
 PyNEC.egg-info/requires.txt
 PyNEC.egg-info/top_level.txt
 example/antenna_util.py
+example/basic_usage.py
 example/context_clean.py
 example/dipole.py
 example/impedance_plot.py
@@ -36,7 +37,6 @@ necpp_src/src/PowerBudget.h
 necpp_src/src/RadiationInput.h
 necpp_src/src/XGetopt.cpp
 necpp_src/src/XGetopt.h
-necpp_src/src/atlas_check.cpp
 necpp_src/src/c_evlcom.cpp
 necpp_src/src/c_evlcom.h
 necpp_src/src/c_geometry.cpp
@@ -56,6 +56,7 @@ necpp_src/src/matrix_algebra.h
 necpp_src/src/misc.cpp
 necpp_src/src/misc.h
 necpp_src/src/nec2cpp.h
+necpp_src/src/nec_card_parser.h
 necpp_src/src/nec_context.cpp
 necpp_src/src/nec_context.h
 necpp_src/src/nec_debug.h

{pynec-1.7.4 → pynec-1.8.1}/PyNEC.py

@@ -198,8 +198,8 @@ class c_geometry(object):
     def arc(self, tag_id, segment_count, rada, ang1, ang2, rad):
         return _PyNEC.c_geometry_arc(self, tag_id, segment_count, rada, ang1, ang2, rad)
 
-    def helix(self, s, hl, a1, b1, a2, b2, rad, segment_count, tag_id):
-        return _PyNEC.c_geometry_helix(self, s, hl, a1, b1, a2, b2, rad, segment_count, tag_id)
+    def helix(self, tag_id, segment_count, s, hl, a1, b1, a2, b2, rad):
+        return _PyNEC.c_geometry_helix(self, tag_id, segment_count, s, hl, a1, b1, a2, b2, rad)
 
     def scale(self, xw1):
         return _PyNEC.c_geometry_scale(self, xw1)
@@ -213,6 +213,9 @@ class c_geometry(object):
     def sp_card(self, ns, in_x1, in_y1, in_z1, in_x2, in_y2, in_z2):
         return _PyNEC.c_geometry_sp_card(self, ns, in_x1, in_y1, in_z1, in_x2, in_y2, in_z2)
 
+    def sc_card(self, i2, x3, y3, z3, x4, y4, z4):
+        return _PyNEC.c_geometry_sc_card(self, i2, x3, y3, z3, x4, y4, z4)
+
     def gx_card(self, card_int_1, card_int_2):
         return _PyNEC.c_geometry_gx_card(self, card_int_1, card_int_2)
 

{pynec-1.7.4 → pynec-1.8.1}/PyNEC_wrap.cxx

@@ -3620,13 +3620,13 @@ SWIGINTERN void c_geometry_move(c_geometry *self,nec_float rox_deg,nec_float roy
       nec_float rox_rad = degrees_to_rad(rox_deg);
       nec_float roy_rad = degrees_to_rad(roy_deg);
       nec_float roz_rad = degrees_to_rad(roz_deg);
-      
+
       return self->move( rox_rad, roy_rad, roz_rad, xs, ys, zs, its, nrpt, itgi );
     }
 SWIGINTERN void c_geometry_arbitrary_shaped_patch(c_geometry *self,nec_float ax1,nec_float ay1,nec_float az1,nec_float ax2_deg,nec_float ay2_deg,nec_float az2){
       nec_float ax2_rad = degrees_to_rad(ax2_deg);
       nec_float ay2_rad = degrees_to_rad(ay2_deg);
-      
+
       return self->patch( 0, 1, ax1, ay1, az1, ax2_rad, ay2_rad, az2, 0, 0, 0, 0, 0, 0 );
     }
 SWIGINTERN void c_geometry_rectangular_patch(c_geometry *self,nec_float ax1,nec_float ay1,nec_float az1,nec_float ax2,nec_float ay2,nec_float az2,nec_float ax3,nec_float ay3,nec_float az3){
@@ -6786,20 +6786,20 @@ fail:
 SWIGINTERN PyObject *_wrap_c_geometry_helix(PyObject *self, PyObject *args) {
   PyObject *resultobj = 0;
   c_geometry *arg1 = (c_geometry *) 0 ;
-  nec_float arg2 ;
-  nec_float arg3 ;
+  int arg2 ;
+  int arg3 ;
   nec_float arg4 ;
   nec_float arg5 ;
   nec_float arg6 ;
   nec_float arg7 ;
   nec_float arg8 ;
-  int arg9 ;
-  int arg10 ;
+  nec_float arg9 ;
+  nec_float arg10 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  double val2 ;
+  int val2 ;
   int ecode2 = 0 ;
-  double val3 ;
+  int val3 ;
   int ecode3 = 0 ;
   double val4 ;
   int ecode4 = 0 ;
@@ -6811,9 +6811,9 @@ SWIGINTERN PyObject *_wrap_c_geometry_helix(PyObject *self, PyObject *args) {
   int ecode7 = 0 ;
   double val8 ;
   int ecode8 = 0 ;
-  int val9 ;
+  double val9 ;
   int ecode9 = 0 ;
-  int val10 ;
+  double val10 ;
   int ecode10 = 0 ;
   PyObject *swig_obj[10] ;
   
@@ -6824,16 +6824,16 @@ SWIGINTERN PyObject *_wrap_c_geometry_helix(PyObject *self, PyObject *args) {
     SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "c_geometry_helix" "', argument " "1"" of type '" "c_geometry *""'"); 
   }
   arg1 = reinterpret_cast< c_geometry * >(argp1);
-  ecode2 = SWIG_AsVal_double(swig_obj[1], &val2);
+  ecode2 = SWIG_AsVal_int(swig_obj[1], &val2);
   if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "c_geometry_helix" "', argument " "2"" of type '" "nec_float""'");
+    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "c_geometry_helix" "', argument " "2"" of type '" "int""'");
   } 
-  arg2 = static_cast< nec_float >(val2);
-  ecode3 = SWIG_AsVal_double(swig_obj[2], &val3);
+  arg2 = static_cast< int >(val2);
+  ecode3 = SWIG_AsVal_int(swig_obj[2], &val3);
   if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), "in method '" "c_geometry_helix" "', argument " "3"" of type '" "nec_float""'");
+    SWIG_exception_fail(SWIG_ArgError(ecode3), "in method '" "c_geometry_helix" "', argument " "3"" of type '" "int""'");
   } 
-  arg3 = static_cast< nec_float >(val3);
+  arg3 = static_cast< int >(val3);
   ecode4 = SWIG_AsVal_double(swig_obj[3], &val4);
   if (!SWIG_IsOK(ecode4)) {
     SWIG_exception_fail(SWIG_ArgError(ecode4), "in method '" "c_geometry_helix" "', argument " "4"" of type '" "nec_float""'");
@@ -6859,16 +6859,16 @@ SWIGINTERN PyObject *_wrap_c_geometry_helix(PyObject *self, PyObject *args) {
     SWIG_exception_fail(SWIG_ArgError(ecode8), "in method '" "c_geometry_helix" "', argument " "8"" of type '" "nec_float""'");
   } 
   arg8 = static_cast< nec_float >(val8);
-  ecode9 = SWIG_AsVal_int(swig_obj[8], &val9);
+  ecode9 = SWIG_AsVal_double(swig_obj[8], &val9);
   if (!SWIG_IsOK(ecode9)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode9), "in method '" "c_geometry_helix" "', argument " "9"" of type '" "int""'");
+    SWIG_exception_fail(SWIG_ArgError(ecode9), "in method '" "c_geometry_helix" "', argument " "9"" of type '" "nec_float""'");
   } 
-  arg9 = static_cast< int >(val9);
-  ecode10 = SWIG_AsVal_int(swig_obj[9], &val10);
+  arg9 = static_cast< nec_float >(val9);
+  ecode10 = SWIG_AsVal_double(swig_obj[9], &val10);
   if (!SWIG_IsOK(ecode10)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode10), "in method '" "c_geometry_helix" "', argument " "10"" of type '" "int""'");
+    SWIG_exception_fail(SWIG_ArgError(ecode10), "in method '" "c_geometry_helix" "', argument " "10"" of type '" "nec_float""'");
   } 
-  arg10 = static_cast< int >(val10);
+  arg10 = static_cast< nec_float >(val10);
   {
     try {
       (arg1)->helix(arg2,arg3,arg4,arg5,arg6,arg7,arg8,arg9,arg10);
@@ -7142,6 +7142,97 @@ fail:
 }
 
 
+SWIGINTERN PyObject *_wrap_c_geometry_sc_card(PyObject *self, PyObject *args) {
+  PyObject *resultobj = 0;
+  c_geometry *arg1 = (c_geometry *) 0 ;
+  int arg2 ;
+  nec_float arg3 ;
+  nec_float arg4 ;
+  nec_float arg5 ;
+  nec_float arg6 ;
+  nec_float arg7 ;
+  nec_float arg8 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  int val2 ;
+  int ecode2 = 0 ;
+  double val3 ;
+  int ecode3 = 0 ;
+  double val4 ;
+  int ecode4 = 0 ;
+  double val5 ;
+  int ecode5 = 0 ;
+  double val6 ;
+  int ecode6 = 0 ;
+  double val7 ;
+  int ecode7 = 0 ;
+  double val8 ;
+  int ecode8 = 0 ;
+  PyObject *swig_obj[8] ;
+  
+  (void)self;
+  if (!SWIG_Python_UnpackTuple(args, "c_geometry_sc_card", 8, 8, swig_obj)) SWIG_fail;
+  res1 = SWIG_ConvertPtr(swig_obj[0], &argp1,SWIGTYPE_p_c_geometry, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "c_geometry_sc_card" "', argument " "1"" of type '" "c_geometry *""'"); 
+  }
+  arg1 = reinterpret_cast< c_geometry * >(argp1);
+  ecode2 = SWIG_AsVal_int(swig_obj[1], &val2);
+  if (!SWIG_IsOK(ecode2)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "c_geometry_sc_card" "', argument " "2"" of type '" "int""'");
+  } 
+  arg2 = static_cast< int >(val2);
+  ecode3 = SWIG_AsVal_double(swig_obj[2], &val3);
+  if (!SWIG_IsOK(ecode3)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode3), "in method '" "c_geometry_sc_card" "', argument " "3"" of type '" "nec_float""'");
+  } 
+  arg3 = static_cast< nec_float >(val3);
+  ecode4 = SWIG_AsVal_double(swig_obj[3], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), "in method '" "c_geometry_sc_card" "', argument " "4"" of type '" "nec_float""'");
+  } 
+  arg4 = static_cast< nec_float >(val4);
+  ecode5 = SWIG_AsVal_double(swig_obj[4], &val5);
+  if (!SWIG_IsOK(ecode5)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode5), "in method '" "c_geometry_sc_card" "', argument " "5"" of type '" "nec_float""'");
+  } 
+  arg5 = static_cast< nec_float >(val5);
+  ecode6 = SWIG_AsVal_double(swig_obj[5], &val6);
+  if (!SWIG_IsOK(ecode6)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode6), "in method '" "c_geometry_sc_card" "', argument " "6"" of type '" "nec_float""'");
+  } 
+  arg6 = static_cast< nec_float >(val6);
+  ecode7 = SWIG_AsVal_double(swig_obj[6], &val7);
+  if (!SWIG_IsOK(ecode7)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode7), "in method '" "c_geometry_sc_card" "', argument " "7"" of type '" "nec_float""'");
+  } 
+  arg7 = static_cast< nec_float >(val7);
+  ecode8 = SWIG_AsVal_double(swig_obj[7], &val8);
+  if (!SWIG_IsOK(ecode8)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode8), "in method '" "c_geometry_sc_card" "', argument " "8"" of type '" "nec_float""'");
+  } 
+  arg8 = static_cast< nec_float >(val8);
+  {
+    try {
+      (arg1)->sc_card(arg2,arg3,arg4,arg5,arg6,arg7,arg8);
+    }
+    catch (nec_exception* nex)  {
+      SWIG_exception(SWIG_RuntimeError,nex->get_message().c_str());
+    }
+    catch (const char* message) {
+      SWIG_exception(SWIG_RuntimeError,message);
+    }
+    catch (...) {
+      SWIG_exception(SWIG_RuntimeError,"Unknown exception");
+    }       
+  }
+  resultobj = SWIG_Py_Void();
+  return resultobj;
+fail:
+  return NULL;
+}
+
+
 SWIGINTERN PyObject *_wrap_c_geometry_gx_card(PyObject *self, PyObject *args) {
   PyObject *resultobj = 0;
   c_geometry *arg1 = (c_geometry *) 0 ;
@@ -12210,6 +12301,7 @@ static PyMethodDef SwigMethods[] = {
 	 { "c_geometry_reflect", _wrap_c_geometry_reflect, METH_VARARGS, NULL},
 	 { "c_geometry_generate_cylindrical_structure", _wrap_c_geometry_generate_cylindrical_structure, METH_VARARGS, NULL},
 	 { "c_geometry_sp_card", _wrap_c_geometry_sp_card, METH_VARARGS, NULL},
+	 { "c_geometry_sc_card", _wrap_c_geometry_sc_card, METH_VARARGS, NULL},
 	 { "c_geometry_gx_card", _wrap_c_geometry_gx_card, METH_VARARGS, NULL},
 	 { "c_geometry_move", _wrap_c_geometry_move, METH_VARARGS, NULL},
 	 { "c_geometry_arbitrary_shaped_patch", _wrap_c_geometry_arbitrary_shaped_patch, METH_VARARGS, NULL},

pynec-1.8.1/example/basic_usage.py

@@ -0,0 +1,135 @@
+#!/usr/bin/env python3
+"""
+Basic PyNEC usage example — demonstrates:
+
+  1. Creating a dipole antenna with geo.wire()  (Issue #22: wire args)
+  2. Getting feedpoint impedance                 (Issue #22: feedpoint impedance)
+  3. Getting structure currents                  (Issue #22: segment currents)
+  4. Using sc_card with multiple_patch           (Issue #28: SC card)
+
+Run:  python basic_usage.py
+"""
+
+import numpy as np
+from PyNEC import *
+
+# ── 1. Wire arguments ──────────────────────────────────────────────────
+#
+# geo.wire(tag_id, segment_count,
+#          xw1, yw1, zw1,   # start point (meters)
+#          xw2, yw2, zw2,   # end point   (meters)
+#          rad,             # wire radius (meters)
+#          rdel=1.0,        # segment length ratio (1.0 = uniform)
+#          rrad=1.0)        # radius taper ratio   (1.0 = no taper)
+#
+# Rule of thumb: use 12-20 segments per wavelength.
+
+context = nec_context()
+geo = context.get_geometry()
+
+freq_mhz = 300.0  # 300 MHz
+wavelength = 300.0 / freq_mhz
+length = wavelength / 2  # half-wave dipole
+radius = length / 1000
+
+geo.wire(
+    tag_id=1,
+    segment_count=21,
+    xw1=0,
+    yw1=-length / 2,
+    zw1=0,  # bottom end
+    xw2=0,
+    yw2=+length / 2,
+    zw2=0,  # top end
+    rad=radius,
+)
+
+context.geometry_complete(0)
+context.fr_card(0, 1, freq_mhz, 0)
+
+# Excite segment 11 (center)
+context.ex_card(0, 1, 11, 0, 1.0, 0.0, 0.0, 0.0, 0.0, 0.0)
+
+# Request radiation pattern (triggers calculation)
+context.rp_card(0, 90, 1, 0, 5, 0, 0, 0.0, 90.0, 1.0, 1.0, 0.0, 0.0)
+
+
+# ── 2. Feedpoint impedance ─────────────────────────────────────────────
+#
+# get_input_parameters(index) returns a nec_antenna_input object.
+# Available methods:
+#   .get_impedance()   -> [complex, ...]   impedance at each source
+#   .get_current()     -> [complex, ...]   current at each source
+#   .get_voltage()     -> [complex, ...]   voltage at each source
+#   .get_power()       -> [float, ...]     power at each source
+#   .get_frequency()   -> float            frequency in Hz
+#   .get_tag()         -> [int, ...]       segment tag numbers
+#   .get_segment()     -> [int, ...]       segment numbers
+
+ipt = context.get_input_parameters(0)
+Z = ipt.get_impedance()[0]
+V = ipt.get_voltage()[0]
+I = ipt.get_current()[0]
+P = ipt.get_power()[0]
+
+print("── Feedpoint ──")
+print(f"  Frequency: {ipt.get_frequency():.1f} MHz")
+print(f"  Impedance: {Z.real:.1f} + {Z.imag:+.1f}j Ω")
+print(
+    f"  VSWR (50Ω): {abs((1 + abs((Z - 50) / (Z + 50))) / (1 - abs((Z - 50) / (Z + 50)))):.3f}"
+)
+print()
+
+
+# ── 3. Structure currents ──────────────────────────────────────────────
+#
+# get_structure_currents(index) returns a nec_structure_currents object.
+# Available methods:
+#   .get_current()                  -> [complex, ...]   currents (A)
+#   .get_current_segment_number()   -> [int, ...]       segment numbers
+#   .get_current_segment_tag()      -> [int, ...]       tag numbers
+#   .get_current_segment_center_x() -> [float, ...]     x coords (m)
+#   .get_current_segment_center_y() -> [float, ...]     y coords (m)
+#   .get_current_segment_center_z() -> [float, ...]     z coords (m)
+#   .get_current_segment_length()   -> [float, ...]     segment lengths (m)
+#   .get_n()                        -> int              total wire segments
+#   .get_m()                        -> int              total patches
+#
+# Also for charge densities: get_q_density(), get_q_density_segment_*()
+# And for patches: get_patch_number(), get_patch_center_*(), get_patch_e_*()
+
+sc = context.get_structure_currents(0)
+currents = sc.get_current()
+tags = sc.get_current_segment_tag()
+
+print("── Segment currents (first 5) ──")
+for i in range(min(5, len(currents))):
+    c = currents[i]
+    mag = abs(c)
+    phase = np.angle(c, deg=True)
+    print(f"  seg {i:2d}  tag={tags[i]:2d}  |I|={mag:.4f} A  ∠={phase:+.1f}°")
+print()
+
+
+# ── 4. SC card with multiple_patch ─────────────────────────────────────
+#
+# multiple_patch() wraps the SM card.  Immediately after it, you can add
+# SC continuation cards via geo.sc_card() to define a row of patches.
+# This example is adapted from necpp/tests/test_multiple_sc_cards.py.
+#
+#   geo.multiple_patch(nx, ny,
+#       ax1, ay1, az1,   # corner 1
+#       ax2, ay2, az2,   # corner 2
+#       ax3, ay3, az3)   # corner 3
+#
+#   geo.sc_card(i2,
+#       x3, y3, z3,      # corner 3 of next patch
+#       x4, y4, z4)      # corner 4 of next patch
+
+print("── SC card example ──")
+ctx2 = nec_context()
+g2 = ctx2.get_geometry()
+
+g2.sp_card(3, 0.019, -0.001424, 0.078830, 0.019, 0.001424, 0.078830)
+g2.sc_card(3, 0.019, 0.001424, 0.076180, 0.019, -0.001424, 0.076180)
+print("  SP + SC card added successfully.")