Dear Sir/Madam .
I really adore and appreciate your effort for getting back to me.Thank a ton ! !
 Yes I was naively wrong with the variation formulation :
The weak form of the time harmonic Electromagnetic wave equation that I want to solve reads:
f * v dx = (mu^1).(curl E) * (curl v) dx + k^2. (E * v) dx
f = i * omega * J
k^2 = (i * sigma * omega)  (eps * omega^2)
I wrote the variation form which reads :
varf vEMwaveEqn([Ex,Ey,Ez],[Vx,Vy,Vz])
= int3d(Team7)(Nu * Curl(Ex,Ey,Ez)’ * Curl(Vx,Vy,Vz) )
 int3d(Team7)([Ex, Ey, Ez]’ * [Vx, Vy, Vz])
 int3d(Team7)([fx,fy,fz]’ * [Vx,Vy,Vz])
 on(4,Ex =0 ,Ey = 0,Ez =0);
//…
matrix EmWaveEqn = vEMwaveEqn(FeEdge,FeEdge,solver = GMRES) ;
complex[int] RHS = vEMwaveEqn(0,FeEdge) ;
I solved the values of :
FeEdge [Jx,Jy,Jz]; // Edge elements has 3 components !
[Jx,Jy,Jz] = Grad(u) * sigma; //*(region == 1) ;
FE1 k2 ; // FE1 is P2 element
FeEdge[Ex,Ey,Ez] ; // FeEdge Edge13d element .
After using iterative solvers, or direct solvers my job is getting “Killed” after few minutes in the terminal .
So I am trying to get some possible answers to the problem.

load "MUMPS" I think it is a plugin and I am trying to install along with PETSc and link it to Freefem++
3.
macro def(i)[i, i#y, i#z]// EOM . Can you please clarify the meaning of" # "in Freefem++.
I wanted to write you after I ran your code but I am having some issues with the PETSc installation due to 2 times installation of MPI (I will fix it soon !)
I would be grateful if you could get back to me.
Regards
Saptarshi Mandal
Bergische University of Wuppertal
Germany