fixes show method crash for FEVectorBlock + additional show functions (#45)

Author: chmerdon

CHANGELOG.md

@@ -1,8 +1,9 @@
 # CHANGES
 
-## v1.2.1 July 09, 2025
+## v1.3.0 July 09, 2025
   - some bugfixes related to new template parameter from 1.2.0
   - @show of FEVectorBlock does not crash anymore
+  - added show functions for FEBasisEvaluator, PointEvaluator, SegmentEvaluator
 
 ## v1.2.0 July 07, 2025
   - major documentation and docstring overhaul

Project.toml

@@ -1,7 +1,7 @@
 name = "ExtendableFEMBase"
 uuid = "12fb9182-3d4c-4424-8fd1-727a0899810c"
 authors = ["Christian Merdon <[email protected]>", "Patrick Jaap <[email protected]>"]
-version = "1.2.1"
+version = "1.3.0"
 
 [deps]
 DiffResults = "163ba53b-c6d8-5494-b064-1a9d43ac40c5"

src/feevaluator.jl

@@ -25,6 +25,17 @@ struct SingleFEEvaluator{T <: Real, TvG <: Real, TiG <: Integer, operator, FETyp
     compressiontargets::Array{Int, 1}            # some operators allow for compressed storage (e.g. SymmetricGradient)
 end
 
+function Base.show(io::IO, FEB::SingleFEEvaluator)
+    println(io, "SingleFEEvaluator")
+    println(io, "-----------------")
+    println(io, "Operator:   ", typeof(FEB).parameters[4])
+    println(io, "FESpace:    ", typeof(FEB.FE))
+    println(io, "Element:    ", typeof(FEB.FE).parameters[3])
+    println(io, "Geometry:   ", typeof(FEB.L2G).parameters[3])
+    println(io, "Quadrature: ", length(FEB.xref), " points")
+    return
+end
+
 
 Base.getindex(FEB::FEEvaluator, c, dof, qp) = FEB.cvals[c, dof, qp]
 

src/fevector.jl

@@ -55,9 +55,13 @@ $(TYPEDSIGNATURES)
 Custom `show` function for `FEVectorBlock` that prints some information and the view of that block.
 """
 function Base.show(io::IO, ::MIME"text/plain", FEB::FEVectorBlock)
-    @printf(io, "block %s [%d:%d] = ", FEB.name, FEB.offset + 1, FEB.last_index)
     return show(io, view(FEB))
 end
+function Base.show(io::IO, FEB::FEVectorBlock)
+    @printf(io, "block %s [%d:%d] = ", FEB.name, FEB.offset + 1, FEB.last_index)
+    show(io, view(FEB))
+    return nothing
+end
 
 """
 $(TYPEDEF)

src/point_evaluator.jl

@@ -13,6 +13,16 @@ mutable struct PointEvaluator{Tv <: Real, UT, KFT <: Function}
     parameters::Dict{Symbol, Any}
 end
 
+function Base.show(io::IO, PE::PointEvaluator)
+    println(io, "PointEvaluator")
+    println(io, "--------------")
+    println(io, "Unknowns: ", PE.u_args)
+    println(io, "Operators: ", PE.ops_args)
+    println(io, "Kernel function: ", typeof(PE.kernel))
+    println(io, "Parameters: ", PE.parameters)
+    return
+end
+
 default_peval_kwargs() = Dict{Symbol, Tuple{Any, String}}(
     :name => ("PointEvaluator", "name for operator used in printouts"),
     :resultdim => (0, "dimension of result field (default = length of operators)"),

src/reconstructionoperators.jl

@@ -91,3 +91,13 @@ function update_basis!(FEBE::FEReconstEvaluator)
     end
     return nothing
 end
+
+function Base.show(io::IO, FEB::FEReconstEvaluator)
+    println(io, "FEReconstEvaluator (reconstruction operator)")
+    println(io, "-------------------------------------------")
+    println(io, "Reconstruction FE space: ", typeof(FEB.FEB.FE))
+    println(io, "Reconstruction operator: ", typeof(FEB.FEB).parameters[4])
+    println(io, "Underlying SingleFEEvaluator:")
+    show(io, FEB.FEB)
+    return
+end

src/segment_integrator.jl

@@ -11,6 +11,17 @@ mutable struct SegmentIntegrator{Tv <: Real, UT, KFT <: Function, EG}
     parameters::Dict{Symbol, Any}
 end
 
+function Base.show(io::IO, SI::SegmentIntegrator)
+    println(io, "SegmentIntegrator")
+    println(io, "-----------------")
+    println(io, "Domain geometry: ", segment_geometry(SI))
+    println(io, "Unknowns: ", SI.u_args)
+    println(io, "Operators: ", SI.ops_args)
+    println(io, "Kernel function: ", typeof(SI.kernel))
+    println(io, "Parameters: ", SI.parameters)
+    return
+end
+
 segment_geometry(::SegmentIntegrator{Tv, UT, KFT, EG}) where {Tv, UT, KFT, EG} = EG
 
 

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,6 +23,7 @@ end
 
 function test_vertex_values(EG, FEType, order)
     qf = VertexRule(EG, order; T = Rational{Int})
+    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

@@ -2,28 +2,31 @@ function run_fematrix_tests()
 
     @testset "FEMatrixVector" begin
         println("\n")
-        println("=======================")
-        println("Testing FEMatrix&VEctor")
-        println("=======================")
+        println("===========================")
+        println("Testing FEMatrix & FEVector")
+        println("===========================")
         xgrid = simplexgrid(0:0.1:1, 0:0.1:1)
         FES1 = FESpace{H1Pk{1, 1, 1}}(xgrid)
+        show(devnull, FES1)
         FES2 = FESpace{H1Pk{1, 1, 2}}(xgrid)
+        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(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(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)