Step 81 - tutorial - Maxwell Equations

[Pratyusa Kar]

Hi, 

I have a question in code mentioned in step - 81 tutorial. 

I see the discretized form is written as (attached image).

Aij=∫Ω(μ−1r∇×φj)⋅(∇×φ¯i)dx−∫Ωεrφj⋅φ¯idx−i∫Σ(σΣrφjT)⋅φ¯iTdox−i∫∂Ω(μ−1rε−−−−√φjT)⋅(∇×φ¯iT)dox,

Fi=i∫ΩJa⋅φi¯dx−∫Ωμ−1rMa⋅(∇×φi¯)dx.

But in the code, 

the magnetization term in the cell_rhs is absent. 
This part of the rhs accounts for interior domain.

My question is: The reason of omitting that term is PML boundary condition? or assumption Ma = 0? 

I am not able to understand the flow here. 

I would request some help with this. 

const auto rhs_value =

                  (imag * scalar_product(J_a, phi_i)) *

fe_values.JxW(q_point);

                cell_rhs(i) += rhs_value.real();

Thanks & Regards, 

Pratyusa Kar