io fix

Author: TomMunoz3772

electroacPy/circuitSolver/solver.py

@@ -145,7 +145,7 @@ def countNodesAndSources(self):
                 nbsource += comp.vsource  # add number of sources
                 # print("nbsource", nbsource)
                 # print("comp.np", comp.np)
-                if comp.vsource == 2:
+                if comp.vsource == 2: # check number of sources, if two = CCVS
                     self.source_id[comp.np] = nbsource-1
                     self.source_id[comp.np1] = nbsource
                 else:

electroacPy/io.py

@@ -290,33 +290,48 @@ def storePressureMeshResults(acoustic_study):
     Nfft = len(study.frequency)
     nRad = study.Ns
     nCoeff = len(study.p_mesh[0, 0].coefficients)
-    nCoeffV = len(study.u_mesh[0, 0].coefficients)
+    
+    # get max u_mesh len (because acoustic radiators can be of different size, i.e. mesh elements)
+    radSurf_elements = []
+    for i in range(np.shape(study.u_mesh)[1]):
+        tmp_length = len(study.u_mesh[0, i].coefficients)
+        radSurf_elements.append(tmp_length)
+    nCoeffV = np.max(radSurf_elements)
 
     # store pressure
     pressureMesh = np.zeros([Nfft, nRad, nCoeff], dtype=complex)
     velocityMesh = np.zeros([Nfft, nRad, nCoeffV], dtype=complex)
     for freq in range(Nfft):
         for rad in range(nRad):
+            tmp_length = radSurf_elements[rad]
             pressureMesh[freq, rad, :] = study.p_mesh[freq, rad].coefficients
-            velocityMesh[freq, rad, :] = study.u_mesh[freq, rad].coefficients
+            velocityMesh[freq, rad, :tmp_length] = study.u_mesh[freq, rad].coefficients
     return pressureMesh, velocityMesh
 
 def storePressureMeshResults_ADM(acoustic_study):
     study = acoustic_study
     Nfft = len(study.frequency)
     nRad = study.Ns
     nCoeff = len(study.p_mesh[0, 0].coefficients)
-    nCoeffV = len(study.u_mesh[0, 0].coefficients)
-    Yn_c        = np.zeros((Nfft, nCoeff), dtype=complex)
-
+    # nCoeffV = len(study.u_mesh[0, 0].coefficients)
+    Yn_c   = np.zeros((Nfft, nCoeff), dtype=complex)
+    
+    # get max u_mesh len (because acoustic radiators can be of different size, i.e. mesh elements)
+    radSurf_elements = []
+    for i in range(np.shape(study.u_mesh)[1]):
+        tmp_length = len(study.u_mesh[0, i].coefficients)
+        radSurf_elements.append(tmp_length)
+    nCoeffV = np.max(radSurf_elements)
+    
     # store pressure
     pressureMesh = np.zeros([Nfft, nRad, nCoeff], dtype=complex)
     velocityMesh = np.zeros([Nfft, nRad, nCoeffV], dtype=complex)
     for freq in range(Nfft):
         Yn_c[freq, :] = study.Yn_c[freq]
         for rad in range(nRad):
+            tmp_length = radSurf_elements[rad]
             pressureMesh[freq, rad, :] = study.p_mesh[freq, rad].coefficients
-            velocityMesh[freq, rad, :] = study.u_mesh[freq, rad].coefficients
+            velocityMesh[freq, rad, :tmp_length] = study.u_mesh[freq, rad].coefficients
     return pressureMesh, velocityMesh, Yn_c
 
 # LOAD PRESSURE

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "electroacPy"
-version = "2025.05.12"
+version = "2025.07.01"
 authors = [
   { name="Tom Munoz", email="munozt@outlook.fr" },
 ]