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