Нормализовать функцию в фенике

Question

У меня есть функция в фенике, которая определена в пространстве функций произведения (для вещественных и мнимых частей). Теперь я хочу вычислить интеграл этой функции по области. У меня нет клея, как это должно работать. И я не могу найти что-то в документации. Я думаю, что документация изменилась на этой неделе, так как ни одна из ссылок Google больше не работает. Функциональное пространство определяется следующим образом:

FuncSpace = FiniteElement ("CG", mesh.ufl_cell (), 1)

FuncSpace2 = FunctionSpace (сетка, MixedElement( [ FuncSpace, FuncSpace ])

И функция определяется как:

Psi = Функция ( FuncSpace2)

Я решаю некоторую проблему частичной интеграции, а затем хочу вычислить интеграл этой функции. Может кто-нибудь помочь мне с этой проблемой, я застрял?

Или я могу преобразовать функцию в некоторый массив Python, потому что тогда интеграция будет простой

Answer 1

Так что я нашел ответ сейчас. Сначала я разбил результирующую функцию Psi на вещественную и сложную части.

Реальный, Imag = раскол ( Psi)

Затем я определяю значение абс функции:

absPsi = sqrt (Реальный * Реальный + Imag * Imag)

Теперь я должен определить меру на используемой сетке:

ds = Measure ("dx", домен = сетка)

где mesh содержит сетку, определенную fenics. И последний интеграл решается путем:

норма = сборка (д * дс)