small addition for SymTracelessTensorValue

Antoinemarteau · Aug 8, 2024 · 98fbe44 · 98fbe44
1 parent 1543f7e
commit 98fbe44
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/src/Fields/AutoDiff.jl b/src/Fields/AutoDiff.jl
@@ -94,7 +94,7 @@ function divergence(f::Function,x::Point,fx::TensorValue{3,3})
    )
 end
 
-function divergence(f::Function,x::Point,fx::SymTensorValue{2})
+function divergence(f::Function,x::Point,fx::AbstractSymTensorValue{2})
   g(x) = SVector(f(x).data)
   a = ForwardDiff.jacobian(g,get_array(x))
   VectorValue(
@@ -103,7 +103,7 @@ function divergence(f::Function,x::Point,fx::SymTensorValue{2})
   )
 end
 
-function divergence(f::Function,x::Point,fx::SymTensorValue{3})
+function divergence(f::Function,x::Point,fx::AbstractSymTensorValue{3})
   g(x) = SVector(f(x).data)
   a = ForwardDiff.jacobian(g,get_array(x))
   VectorValue(

diff --git a/test/VisualizationTests/VtkTests.jl b/test/VisualizationTests/VtkTests.jl
@@ -96,6 +96,7 @@ writevtk(trian,f,nsubcells=10, cellfields=[
   "v2"=>x->VectorValue(1,2),
   "v"=>x->VectorValue(1,2,3),
   "s"=>x->SymTensorValue(1.0,2.0,3.0),
+  "q"=>x->SymTracelessTensorValue(1.0,2.0),
   "c"=>x->SymFourthOrderTensorValue(1,2,3, 1,2,3, 1,2,3),
   "t"=>x->TensorValue(1,2,3,4),])