...

Author: xtalax

test/robin.jl

@@ -21,10 +21,6 @@ for i in 1:5
     @test Q_L[1,:] ≈ [-1 / (1-al[i]*dx_l[i]/bl[i]); zeros(5i-1)]
     @test Q_L[5i+2,:] ≈ [zeros(5i-1); 1 / (1+ar[i]*dx_r[i]/br[i])]
 
-    #Test the test
-    @test (u[1]*(bl[i]/dx_l[i]))/(al[i]-bl[i]/dx_l[i]) ≈ -u[1] / (1-al[i]*dx_l[i]/bl[i])
-    @test (u[end]*(br[i]/dx_r[i]))/(ar[i]+br[i]/dx_r[i]) ≈ u[end] / (1+ar[i]*dx_r[i]/br[i])
-
     #Check that Q_b is computed correctly
     @test Q_b ≈ [cl[i]/(al[i]-bl[i]/dx_l[i]); zeros(5i); cr[i]/(ar[i]+br[i]/dx_r[i])]
 
@@ -35,8 +31,10 @@ for i in 1:5
     Qextended = Q*u
     CorrectQextended = [(cl[i]+(bl[i]/dx_l[i])*u[1])/(al[i]-bl[i]/dx_l[i]); u; (cr[i]+ (br[i]/dx_r[i])*u[5i])/(ar[i]+br[i]/dx_r[i])]
 
+
     @test length(Qextended) ≈ 5i+2
-    #@test Qextended ≈ CorrectQextended
+
+    #test the test    
     @test (u[1]*(bl[i]/dx_l[i]))/(al[i]-bl[i]/dx_l[i]) ≈ -u[1] / (1-al[i]*dx_l[i]/bl[i])
     @test (u[end]*(br[i]/dx_r[i]))/(ar[i]+br[i]/dx_r[i]) ≈ u[end] / (1+ar[i]*dx_r[i]/br[i])