fullwave25 1.2.6.dev9__tar.gz → 1.2.6.dev10__tar.gz

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/.bumpversion.toml

@@ -1,5 +1,5 @@
 [tool.bumpversion]
-current_version = "1.2.6-dev9"
+current_version = "1.2.6-dev10"
 parse = """(?x)
     (?P<major>0|[1-9]\\d*)\\.
     (?P<minor>0|[1-9]\\d*)\\.

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: fullwave25
-Version: 1.2.6.dev9
+Version: 1.2.6.dev10
 Summary: Fullwave 2.5: Ultrasound wave propagation simulation with heterogeneous power law attenuation modelling capabilities
 Author: Masashi Sode, Gianmarco Pinton
 Maintainer: Masashi Sode

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/fullwave/__init__.py

@@ -60,7 +60,7 @@ try:
     __version__ = version("fullwave")
 except PackageNotFoundError:
     # Update via bump-my-version, not manually
-    __version__ = "1.2.6-dev9"
+    __version__ = "1.2.6-dev10"
 
 VERSION = __version__  # for convenience
 logger.info("Fullwave version: %s", __version__)

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/fullwave/utils/pulse.py

@@ -61,6 +61,7 @@ def gaussian_modulated_sinusoidal_signal(
     dt = duration / nt
     t_offset = ncycles / f0 + delay_sec
 
+    phase_correction = 0.0
     if i_layer is not None:
         if dt_for_layer_delay is None:
             error_msg = "dt_for_layer_delay must be provided if i_layer is provided"
@@ -68,7 +69,17 @@ def gaussian_modulated_sinusoidal_signal(
         if cfl_for_layer_delay is None:
             error_msg = "cfl_for_layer_delay must be provided if i_layer is provided"
             raise ValueError(error_msg)
-        t_offset += (dt_for_layer_delay / cfl_for_layer_delay) * i_layer
+        layer_delay_sec = (dt_for_layer_delay / cfl_for_layer_delay) * i_layer
+        t_offset += layer_delay_sec
+        # Correct carrier phase for fractional-sample delays.
+        # A continuous time shift of a fractional number of samples changes
+        # the carrier phase at each discrete sample point, causing sign
+        # inversion and wrong peak positions. Remove the fractional phase
+        # offset so the carrier stays on the same discrete phase grid as
+        # the base (i_layer=0) signal.
+        delay_samples = layer_delay_sec / dt
+        fractional_samples = delay_samples - round(delay_samples)
+        phase_correction = 2.0 * np.pi * f0 * fractional_samples * dt
 
     t = np.arange(nt, dtype=dtype) * dt - t_offset
 
@@ -87,5 +98,5 @@ def gaussian_modulated_sinusoidal_signal(
     else:
         env = np.exp(-(a_sq**drop_off))
 
-    # Compute final signal
-    return env * np.sin(w_t) * p0
+    # Compute final signal with phase-corrected carrier
+    return env * np.sin(w_t + phase_correction) * p0

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "fullwave25"
-version = "1.2.6-dev9" # Update via bump-my-version, not manually
+version = "1.2.6-dev10" # Update via bump-my-version, not manually
 description = "Fullwave 2.5: Ultrasound wave propagation simulation with heterogeneous power law attenuation modelling capabilities"
 readme = "README.md"
 requires-python = ">=3.10"

{fullwave25-1.2.6.dev9 → fullwave25-1.2.6.dev10}/uv.lock

@@ -735,7 +735,7 @@ wheels = [
 
 [[package]]
 name = "fullwave25"
-version = "1.2.6.dev9"
+version = "1.2.6.dev10"
 source = { editable = "." }
 dependencies = [
     { name = "matplotlib" },