0

我在 fenics 中有一个函数,它是在产品函数空间(用于实部和虚部)上定义的。现在我想计算这个函数在域上的积分。我没有胶水应该如何工作。而且我无法在文档中找到某些内容。我认为本周文档发生了变化,因为谷歌链接不再有效。函数空间定义如下:

FuncSpace = FiniteElement(“CG”,mesh.ufl_cell(),1)

FuncSpace2 = FunctionSpace(mesh, MixedElement([FuncSpace, FuncSpace])

函数定义为:

Psi = 函数(FuncSpace2)

我正在解决一些部分积分问题,然后想计算这个函数的积分。任何人都可以帮我解决这个问题,我很困惑吗?

或者我可以将函数转换为一些 python 数组,因为这样我的集成就很简单了

4

1 回答 1

0

所以我现在找到了答案。首先,我将结果函数 Psi 拆分为实数和复数部分

实数,图像 = 分裂(Psi)

接下来我通过以下方式定义函数的 abs 值:

absPsi = sqrt( 实 * 实 + 图像 * 图像 )

现在我必须在使用的网格上定义一个度量:

ds = 测量(“dx”,域 = 网格)

其中 mesh 包含由 fenics 定义的网格。最后积分通过以下方式求解:

规范=组装(d * ds)

于 2018-01-20T10:16:57.263 回答