Add shape_symmetric_gradient for PointValues (#1325)

Defines generic `shape_symmetric_gradient` for `AbstractValues`

Author: henrij22

CHANGELOG.md

@@ -6,6 +6,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
 ## [Next] - xxxx-xx-xx
+ - Make `shape_symmetric_gradient` work for `PointValues` ([#1325])
 
 ## [v1.4.0] - 2026-04-20
 

src/FEValues/CellValues.jl

@@ -100,7 +100,6 @@ shape_hessian_type(cv::CellValues) = shape_hessian_type(cv.fun_values)
 @propagate_inbounds shape_value(cv::CellValues, q_point::Int, i::Int) = shape_value(cv.fun_values, q_point, i)
 @propagate_inbounds shape_gradient(cv::CellValues, q_point::Int, i::Int) = shape_gradient(cv.fun_values, q_point, i)
 @propagate_inbounds shape_hessian(cv::CellValues, q_point::Int, i::Int) = shape_hessian(cv.fun_values, q_point, i)
-@propagate_inbounds shape_symmetric_gradient(cv::CellValues, q_point::Int, i::Int) = shape_symmetric_gradient(cv.fun_values, q_point, i)
 
 # Access quadrature rule values
 get_quadrature_rule(cv::CellValues) = cv.qr

src/FEValues/FacetValues.jl

@@ -93,7 +93,6 @@ get_fun_values(fv::FacetValues) = @inbounds fv.fun_values[getcurrentfacet(fv)]
 @propagate_inbounds shape_value(fv::FacetValues, q_point::Int, i::Int) = shape_value(get_fun_values(fv), q_point, i)
 @propagate_inbounds shape_gradient(fv::FacetValues, q_point::Int, i::Int) = shape_gradient(get_fun_values(fv), q_point, i)
 @propagate_inbounds shape_hessian(fv::FacetValues, q_point::Int, i::Int) = shape_hessian(get_fun_values(fv), q_point, i)
-@propagate_inbounds shape_symmetric_gradient(fv::FacetValues, q_point::Int, i::Int) = shape_symmetric_gradient(get_fun_values(fv), q_point, i)
 
 """
     getcurrentfacet(fv::FacetValues)

src/FEValues/FunctionValues.jl

@@ -109,7 +109,6 @@ getnquadpoints(funvals::FunctionValues) = size(funvals.Nx, 2)
 @propagate_inbounds shape_value(funvals::FunctionValues, q_point::Int, base_func::Int) = funvals.Nx[base_func, q_point]
 @propagate_inbounds shape_gradient(funvals::FunctionValues, q_point::Int, base_func::Int) = funvals.dNdx[base_func, q_point]
 @propagate_inbounds shape_hessian(funvals::FunctionValues{2}, q_point::Int, base_func::Int) = funvals.d2Ndx2[base_func, q_point]
-@propagate_inbounds shape_symmetric_gradient(funvals::FunctionValues, q_point::Int, base_func::Int) = symmetric(shape_gradient(funvals, q_point, base_func))
 
 function_interpolation(funvals::FunctionValues) = funvals.ip
 function_difforder(::FunctionValues{DiffOrder}) where {DiffOrder} = DiffOrder

src/FEValues/common_values.jl

@@ -135,7 +135,7 @@ shape_gradient(fe_v::AbstractValues, q_point::Int, base_function::Int)
 Return the symmetric gradient of shape function `base_function` evaluated in
 quadrature point `q_point`.
 """
-function shape_symmetric_gradient end
+@propagate_inbounds shape_symmetric_gradient(cv::AbstractValues, q_point::Int, base_func::Int) = symmetric(shape_gradient(cv, q_point, base_func))
 
 """
     shape_divergence(fe_v::AbstractValues, q_point::Int, base_function::Int)

test/test_pointevaluation.jl

@@ -513,8 +513,19 @@ end
     @test function_value(pvv, uv) ≈ function_value(cvv, 1, uv)
     @test function_gradient(pvv, uv) ≈ function_gradient(cvv, 1, uv)
     @test function_symmetric_gradient(pvv, uv) ≈ function_symmetric_gradient(cvv, 1, uv)
+
+    @test shape_gradient(pvv, 1, 1) ≈ shape_gradient(cvv, 1, 1)
+    @test shape_divergence(pvv, 1, 1) ≈ shape_divergence(cvv, 1, 1)
+    @test shape_curl(pvv, 1, 1) ≈ shape_curl(cvv, 1, 1)
+    @test shape_symmetric_gradient(pvv, 1, 1) ≈ shape_symmetric_gradient(cvv, 1, 1)
+
     reinit!(pvv, x, ξ₂)
     @test function_value(pvv, uv) ≈ function_value(cvv, 2, uv)
     @test function_gradient(pvv, uv) ≈ function_gradient(cvv, 2, uv)
     @test function_symmetric_gradient(pvv, uv) ≈ function_symmetric_gradient(cvv, 2, uv)
+
+    @test shape_gradient(pvv, 1, 1) ≈ shape_gradient(cvv, 2, 1)
+    @test shape_divergence(pvv, 1, 1) ≈ shape_divergence(cvv, 2, 1)
+    @test shape_curl(pvv, 1, 1) ≈ shape_curl(cvv, 2, 1)
+    @test shape_symmetric_gradient(pvv, 1, 1) ≈ shape_symmetric_gradient(cvv, 2, 1)
 end