cell_values.jl
会调用 interpolations.jl
中的 value(...)
函数,value(...)
函数是用来计算参考单元在积分点处的值:
for (qp, ξ) in enumerate(quad_rule.points)
for i in 1:n_func_basefuncs
dNdξ[i, qp], N[i, qp] = gradient(ξ -> value(func_interpol, i, ξ), ξ, :all)
end
for i in 1:n_geom_basefuncs
dMdξ[i, qp], M[i, qp] = gradient(ξ -> value(geom_interpol, i, ξ), ξ, :all)
end
end