From 98fbe44b17699935edf9b8f4978027bbfeb877ec Mon Sep 17 00:00:00 2001 From: Antoine Marteau Date: Thu, 8 Aug 2024 16:22:52 +1000 Subject: [PATCH] small addition for SymTracelessTensorValue --- src/Fields/AutoDiff.jl | 4 ++-- test/VisualizationTests/VtkTests.jl | 1 + 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/src/Fields/AutoDiff.jl b/src/Fields/AutoDiff.jl index 5c6ada102..6bfa9d12d 100644 --- 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 index 055da65eb..af9ad91c4 100644 --- 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),])