Flip wavelet domain for the Lp norm constraint.

TOMOFAST · Oct 16, 2023 · 44dae12 · 44dae12
1 parent e8054dc
commit 44dae12
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/inversion/joint_inverse_problem.F90 b/src/inversion/joint_inverse_problem.F90
@@ -190,7 +190,8 @@ subroutine joint_inversion_initialize(this, par, nnz_sensit, myrank)
   ! Determine if we perform inversion using wavelet domain for the model update vector.
   this%WAVELET_DOMAIN = .true.
   if (this%add_cross_grad .or. this%add_clustering .or. &
-      this%add_damping_gradient(1) .or. this%add_damping_gradient(2)) then
+      this%add_damping_gradient(1) .or. this%add_damping_gradient(2) .or. &
+      par%norm_power /= 2.d0) then
     this%WAVELET_DOMAIN = .false.
   endif