@@ -2643,12 +2643,12 @@ def run_adjoint_F(self, nd):
2643
2643
preset = 'layers-isotropic' , dtype = np .float64 )
2644
2644
# Run forward operator
2645
2645
rec , u , _ = solver .forward ()
2646
- # printf(norm(u))
2647
2646
assert np .isclose (norm (u ) / Eu , 1.0 )
2648
2647
assert np .isclose (norm (rec ) / Erec , 1.0 )
2649
2648
2650
2649
# Run adjoint operator
2651
2650
srca , v , _ = solver .adjoint (rec = rec )
2651
+
2652
2652
assert np .isclose (norm (v ) / Ev , 1.0 )
2653
2653
assert np .isclose (norm (srca ) / Esrca , 1.0 )
2654
2654
@@ -2657,8 +2657,9 @@ def run_adjoint_F(self, nd):
2657
2657
term2 = norm (rec )** 2
2658
2658
assert np .isclose ((term1 - term2 )/ term1 , 0. , rtol = 1.e-10 )
2659
2659
2660
- @pytest .mark .parametrize ('nd' , [3 ])
2661
- @pytest .mark .parallel (mode = [(4 , 'basic2' )])
2660
+ @pytest .mark .parametrize ('nd' , [1 , 2 , 3 ])
2661
+ @pytest .mark .parallel (mode = [(4 , 'basic' ), (4 , 'basic2' ), (4 , 'diag' ),
2662
+ (4 , 'overlap' ), (4 , 'overlap2' ), (4 , 'full' )])
2662
2663
def test_adjoint_F (self , nd ):
2663
2664
self .run_adjoint_F (nd )
2664
2665
0 commit comments