Dipole broadcasting

.github/environment_ci.yml

@@ -22,6 +22,7 @@ dependencies:
   # testing
   - pytest
   - pytest-cov
+  - sympy
 
   # Building
   - pip

.github/workflows/test_with_conda.yml

@@ -68,7 +68,7 @@ jobs:
 
     - name: "Upload coverage to Codecov"
       if: matrix.python-version == '3.11'
-      uses: codecov/codecov-action@v3
+      uses: codecov/codecov-action@v4
       with:
         token: ${{ secrets.CODECOV_TOKEN }}
         fail_ci_if_error: true

geoana/earthquake/oksar.py

@@ -1657,7 +1657,7 @@ def example():
     #     data=data
     # )
 
-    data = np.fromstring(dinar_file.content, np.float32)
+    data = np.frombuffer(dinar_file.content, np.float32)
     title = 'Dinar, Turkey'
     location = [706216.0606, 4269238.9999]
     location_UTM_zone = 35

geoana/em/base.py

@@ -58,6 +58,26 @@
 
         self._sigma = value
 
+    @property
+    def rho(self):
+        """Electrical resistivity in Ohm m
+
+        Returns
+        -------
+        float
+            Electrical resistivity in Ohm m
+        """
+        return 1/self.sigma
+
+    @rho.setter
+    def rho(self, value):
+        try:
+            value = float(value)
+        except:
+            raise TypeError(f"resistivity must be a number, got {type(value)}")
+
+        self.sigma = 1/value
+
     @property
     def mu(self):
         """Magnetic permeability in H/m
@@ -177,6 +197,8 @@
                 var = np.r_[0., 1., 0.]
             elif var.upper() == 'Z':
                 var = np.r_[0., 0., 1.]
+            else:
+                raise ValueError("Orientation must be one of {'X','Y','Z'}")
         else:
             try:
                 var = np.asarray(var, dtype=float)
@@ -187,7 +209,7 @@
             var = np.squeeze(var)
             if var.shape != (3, ):
                 raise ValueError(
-                    f"orientation must be array_like with shape (3,), got {len(var)}"
+                    f"orientation must be array_like with shape (3,), got {var.shape}"
                 )
 
             # Normalize the orientation

geoana/em/fdem/base.py

@@ -203,13 +203,11 @@ def frequency(self, value):
             value = np.asarray(value, dtype=float)
         except:
             raise TypeError(f"frequencies are not a valid type")
-        value = np.atleast_1d(value)
-
         # Enforce positivity and dimensions
         if (value < 0.).any():
             raise ValueError("All frequencies must be greater than 0")
         if value.ndim > 1:
-            raise TypeError(f"frequencies must be ('*') array")
+            raise TypeError(f"frequencies must have at most 1 dimension.")
 
         self._frequency = value
 

geoana/em/fdem/halfspace.py

@@ -34,7 +34,7 @@ def magnetic_field(self, xy, field="secondary"):
 
         Returns
         -------
-        (n_freq, ..., 3) numpy.array of complex
+        (n_freq, ..., 3) numpy.ndarray of complex
             Magnetic field at all frequencies for the gridded
             locations provided. Output array is squeezed when n_freq and/or
             n_loc = 1.

geoana/em/fdem/layered.py

@@ -300,7 +300,7 @@ def magnetic_field(self, xyz, field="secondary"):
 
         Returns
         -------
-        (n_freq, n_loc, 3) numpy.array of complex
+        (n_freq, n_loc, 3) numpy.ndarray of complex
             Magnetic field at all frequencies for the gridded
             locations provided. Output array is squeezed when n_freq and/or
             n_loc = 1.