Hello dear all
I have solved a transitional problem on a square domain. However, I would like to do the hilbert transform of the solution on this same domain. Only, I did not find a fenics library for that. I have tried the scipy.signal library, but I get no result. Please, has anyone done a similar maneuver?
thanks in advance?
from dolfin import *
import numpy as np
from scipy.signal import hilbert2, hilbert
# after solving have the solution
I4f1 = inner(f01, E_out1 * f01)
I4f11 = np.asarray(I4f1)
#%%%%%% phase singularities %%%%% hilbert transfort %%%%%%#
im_I1 = ufl.imag(hilbert2(I4f11))
r_I1 = ufl.real(hilbert2(I4f11))
PH1 = project(ufl.atan(im_I1/r_I1),Mh1)
AM1 = project ((im_I1**2+r_I1**2)**0.5,Mh1)
Thanks