show(devnull, obj)

Author: chmerdon

test/runtests.jl

@@ -6,7 +6,6 @@ using ExampleJuggler
 using SparseArrays
 using Aqua
 
-
 @testset "Aqua.jl" begin
     Aqua.test_all(
         ExtendableFEMBase;

test/test_febasis.jl

@@ -23,7 +23,7 @@ end
 
 function test_vertex_values(EG, FEType, order)
     qf = VertexRule(EG, order; T = Rational{Int})
-    @show qf
+    show(devnull, qf)
     basis = get_basis(ON_CELLS, FEType, EG)
     ndofs = get_ndofs(ON_CELLS, FEType, EG)
     basis_vals = zeros(Rational{Int}, ndofs, 1)

test/test_fematrix_and_vector.jl

@@ -7,26 +7,26 @@ function run_fematrix_tests()
         println("===========================")
         xgrid = simplexgrid(0:0.1:1, 0:0.1:1)
         FES1 = FESpace{H1Pk{1, 1, 1}}(xgrid)
-        @show FES1
+        show(devnull, FES1)
         FES2 = FESpace{H1Pk{1, 1, 2}}(xgrid)
-        @show FES2
+        show(devnull, FES2)
         A = FEMatrix(FES1, FES2)
         @test size(A.entries) == (FES1.ndofs, FES2.ndofs)
         @test size(A[1, 1]) == (FES1.ndofs, FES2.ndofs)
-        @show A
+        show(devnull, A)
 
         B = FEMatrix([FES1, FES2])
         @test length(B) == 4
         @test size(B.entries) == (FES1.ndofs + FES2.ndofs, FES1.ndofs + FES2.ndofs)
         @test size(B[1, 2]) == (FES1.ndofs, FES2.ndofs)
-        @show B
+        show(devnull, B)
 
         C = FEMatrix([FES2, FES2], [FES1, FES1])
         @test length(C) == 4
         @test size(C.entries) == (2 * FES2.ndofs, 2 * FES1.ndofs)
         @test size(C[1, 2]) == (FES2.ndofs, FES1.ndofs)
         C.entries.cscmatrix = sprand(2 * FES2.ndofs, 2 * FES1.ndofs, 0.5)
-        @show C
+        show(devnull, C)
 
         b = FEVector([FES1, FES2])
         b.entries .= rand(FES1.ndofs + FES2.ndofs)

test/test_interpolators.jl

@@ -73,7 +73,7 @@ function run_interpolator_tests()
             QP = QPInfos(xgrid)
             qf = VertexRule(EG, order)
             FEB = FEEvaluator(FES, Identity, qf)
-            @show FEB
+            show(devnull, FEB)
             for cell::Int in cells
                 update_trafo!(L2G, cell)
                 update_basis!(FEB, cell)
@@ -106,7 +106,7 @@ function run_interpolator_tests()
         # interpolate
         Solution = FEVector(FES)
         interpolate!(Solution[1], u; bonus_quadorder = order)
-        @show Solution
+        show(devnull, Solution)
 
         # compute error
         error = compute_error(Solution[1], u, order)

test/test_operators.jl

@@ -32,7 +32,7 @@ function test_derivatives2D()
     ## define P2-Courant finite element space
     FEType = H1P2{2, 2}
     FES = FESpace{FEType}(xgrid)
-    @show FES
+    show(devnull, FES)
 
     ## get midpoint quadrature rule for constants
     qf = QuadratureRule{Float64, Triangle2D}(0)
@@ -114,7 +114,7 @@ function test_derivatives3D()
     ## define P2-Courant finite element space
     FEType = H1P2{3, 3}
     FES = FESpace{FEType}(xgrid)
-    @show FES
+    show(devnull, FES)
 
     ## get midpoint quadrature rule for constants
     qf = QuadratureRule{Float64, Tetrahedron3D}(0)

test/test_pointevaluator.jl

@@ -17,7 +17,7 @@ function test_pointevaluation2D()
     Iu = FEVector(FES)
     interpolate!(Iu[1], (result, qpinfo) -> (result[1] = qpinfo.x[1] + qpinfo.x[2]))
     PE = PointEvaluator([(1, Identity)])
-    @show PE
+    show(devnull, PE)
     initialize!(PE, Iu)
     CF = CellFinder(xgrid)
     eval = zeros(Float64, 1)
@@ -39,7 +39,7 @@ function test_pointevaluation3D()
     Iu = FEVector(FES)
     interpolate!(Iu[1], (result, qpinfo) -> (result[1] = qpinfo.x[1] + qpinfo.x[2] + qpinfo.x[3]))
     PE = PointEvaluator([(1, Identity)])
-    @show PE
+    show(devnull, PE)
     initialize!(PE, Iu)
     CF = CellFinder(xgrid)
     eval = zeros(Float64, 1)

test/test_quadrature.jl

@@ -17,7 +17,7 @@ function run_quadrature_tests()
         for order in 1:maxorder1D
             integrand, exactvalue = exact_function(Val(1), order)
             qf = QuadratureRule{Float64, Edge1D}(order)
-            @show qf
+            show(devnull, qf)
             quadvalue = integrate(xgrid, ON_CELLS, integrand, length(exactvalue); force_quadrature_rule = qf)
             println("EG = Edge1D | order = $order ($(qf.name), $(length(qf.w)) points) | error = $(quadvalue - exactvalue)")
             @test isapprox(quadvalue, exactvalue)
@@ -31,7 +31,7 @@ function run_quadrature_tests()
             for order in 1:maxorder2D[j]
                 integrand, exactvalue = exact_function(Val(2), order)
                 qf = QuadratureRule{Float64, EG}(order)
-                @show qf
+                show(devnull, qf)
                 quadvalue = integrate(xgrid, ON_CELLS, integrand, length(exactvalue); force_quadrature_rule = qf)
                 println("EG = $EG | order = $order ($(qf.name), $(length(qf.w)) points) | error = $(quadvalue - exactvalue)")
                 @test isapprox(quadvalue, exactvalue)
@@ -46,7 +46,7 @@ function run_quadrature_tests()
             for order in 1:maxorder3D[j]
                 integrand, exactvalue = exact_function(Val(3), order)
                 qf = QuadratureRule{Float64, EG}(order)
-                @show qf
+                show(devnull, qf)
                 quadvalue = integrate(xgrid, ON_CELLS, integrand, length(exactvalue); force_quadrature_rule = qf)
                 println("EG = $EG | order = $order ($(qf.name), $(length(qf.w)) points) | error = $(quadvalue - exactvalue)")
                 @test isapprox(quadvalue, exactvalue)

test/test_segmentintegrator.jl

@@ -19,7 +19,7 @@ function test_segmentintegrator_nokernel()
     ## init segment integrator
     SI = SegmentIntegrator(Edge1D, [(1, Identity)])
     initialize!(SI, uh)
-    @show SI
+    show(devnull, SI)
 
     ## integrate along line [1/4,1/4] to [3/4,1/4] in first triangle
     ## exact integral should be [3//32,0]
@@ -66,9 +66,7 @@ function test_segmentintegrator_withkernel()
     ## init segment integrator
     SI = SegmentIntegrator(Edge1D, multiply_r!, [(1, Identity)]; bonus_quadorder = 1)
     initialize!(SI, uh)
-    @show SI
-
-    @show xgrid[Coordinates], xgrid[CellNodes]
+    show(devnull, SI)
 
     L2G = L2GTransformer(Triangle2D, xgrid, ON_CELLS)
     update_trafo!(L2G, 1)