При использовании формулировки Галеркина для уравнений Максвелла первого порядка, уравнение

$rotH=J$

(для проводящей среды) переходит в

$(rotH, \phi)=(J, \phi)$, где $\phi \in H(div, \Omega)$

Т.е. в дискретном смысле $\phi - базисная функция из $H^h(div, \Omega)$, которое натянуто на face-конформный базис.

Так вот, при стремлении h к 0, $rotH$, как выяснилось, не сходится к $J$. Не сходится. Ни в норме, ни поточечно, никак.

Пиздец, пойду повешусь.

Есть, конечно, еще способы, но все они требуют серьезных ресурсов, в первую очередь времени, которого катастрофически нет.