fixes for discontinuous function plots

Project.toml

@@ -1,7 +1,7 @@
 name = "PlutoVista"
 uuid = "646e1f28-b900-46d7-9d87-d554eb38a413"
 authors = ["Juergen Fuhrmann <juergen.fuhrmann@wias-berlin.de>"]
-version = "0.8.17"
+version = "0.8.18"
 
 [deps]
 ColorSchemes = "35d6a980-a343-548e-a6ea-1d62b119f2f4"
@@ -16,7 +16,7 @@ UUIDs = "cf7118a7-6976-5b1a-9a39-7adc72f591a4"
 ColorSchemes = "3"
 Colors = "0.12,1"
 DocStringExtensions = "0.8,0.9"
-GridVisualizeTools = "0.2.1,0.3"
+GridVisualizeTools = "0.3"
 HypertextLiteral = "0.9"
 Pluto = "0.19"
 julia = "1.6"

src/plutovtkplot.jl

@@ -283,9 +283,11 @@ contour!(p::PlutoVTKPlot,X,Y,f; kwargs...)=tricontour!(p,triang(X,Y)...,vec(f);k
 $(SIGNATURES)
 
 Plot isosurfaces given by `levels` and contour maps on planes given by the `*planes` parameters
-for piecewise linear function on  tetrahedral mesh. 
+for piecewise linear function on  tetrahedral mesh.
+
+`pts,tets,func` can be arrays of grid data, allowing to plot discontinouos functions
 """
-function tetcontour!(p::PlutoVTKPlot, pts, tets,func; kwargs...)
+function tetcontour!(p::PlutoVTKPlot, pts, tets,func; parentpts=pts, kwargs...)
     reset!(p)
 
     default_args=(colormap=:viridis,
@@ -326,12 +328,12 @@ function tetcontour!(p::PlutoVTKPlot, pts, tets,func; kwargs...)
     if facecolormap==nothing
         facecolormap=bregion_cmap(nbregions)
     end
-        
+
     @views for idim=1:3
-        xyzmin[idim]=minimum(pts[idim,:])
-        xyzmax[idim]=maximum(pts[idim,:])
+        xyzmin[idim]=minimum(parentpts[idim,:])
+        xyzmax[idim]=maximum(parentpts[idim,:])
     end 
-
+    
     xplanes=args[:xplanes] 
     yplanes=args[:yplanes] 
     zplanes=args[:zplanes]  
@@ -400,7 +402,7 @@ function tetcontour!(p::PlutoVTKPlot, pts, tets,func; kwargs...)
 
     if args[:outlinealpha]>0 && faces!=nothing
         parameter!(p,"outline",1)
-        outline!(p,pts,faces,facemarkers,facecolormap,nbregions,xyzmin,xyzmax;alpha=args[:outlinealpha])
+        outline!(p,parentpts,faces,facemarkers,facecolormap,nbregions,xyzmin,xyzmax;alpha=args[:outlinealpha])
     else
         parameter!(p,"outline",0)
     end