I am trying to solve the time-harmonic elastic wave problem using the Galerkin method.
But the results are not correct.
I am attaching a jupyter-notebook with the case that I tried. Could you tell me what I am doing wrong?
In addition, I would like to include a pointwise source in this problem, could you also help me with this?
Finally, how can I use the complex conjugate to measure L2-norm error?
gfu = GridFunction(fes) u_exact = CoefficientFunction(r_source) u = CoefficientFunction(gfu) #Error analysis print ("Displacement L2-error:", sqrt (Integrate ( (u-u_exact)*(u-u_exact), mesh)))
Thank you in advance for your time.