Hi dear all, i want to use Dirac delta function in my domain not just at special point like x0. below is a example of this
function MatrixOf(fem_params) x0 = VectorValue(0,300) # Position of the field to be optimized δ = 1 return assemble_matrix(fem_params.U, fem_params.V) do u, v ∫((x->(1/(2*π)*exp(-norm(x - x0)^2 / 2 / δ^2))) * (∇(u) ⋅ ∇(v)) )fem_params.dΩ_c end end
and to be clear it finds the matrix around x0. however i want to assemble it for each node. i have derived the node coordinate by below line
however i cant make a loop for gridap to go through this vector and find o for each and then assemble it for all.
thanks for any suggestion.