Skip to content

Commit 2f985fa

Browse files
waltsimsclaude
andcommitted
Revert unnecessary nabla1 conditional — not worth the complexity
Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
1 parent 7a4ce97 commit 2f985fa

1 file changed

Lines changed: 3 additions & 6 deletions

File tree

kwave/solvers/kspace_solver.py

Lines changed: 3 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -427,12 +427,9 @@ def _setup_physics_operators(self):
427427
self._absorption = lambda div_u: 0 if no_absorption else -2 * alpha_np * self.c0 * self.rho0 * div_u
428428
self._dispersion = lambda rho: 0
429429
else: # Power-law with fractional Laplacian
430-
if no_absorption:
431-
self._absorption = lambda div_u: 0
432-
else:
433-
tau = -2 * alpha_np * self.c0 ** (alpha_power - 1)
434-
nabla1 = self._fractional_laplacian(alpha_power - 2)
435-
self._absorption = lambda div_u: tau * self._diff(self.rho0 * div_u, nabla1)
430+
tau = -2 * alpha_np * self.c0 ** (alpha_power - 1)
431+
nabla1 = self._fractional_laplacian(alpha_power - 2)
432+
self._absorption = (lambda div_u: 0) if no_absorption else (lambda div_u: tau * self._diff(self.rho0 * div_u, nabla1))
436433

437434
if no_dispersion:
438435
self._dispersion = lambda rho: 0

0 commit comments

Comments
 (0)