python-2.7 - 单位球面上的曲面积分

Question

我正在尝试求解球体表面上的方程，这与球体表面上的 Cahn-Hilliard 问题非常相似。我已经按照 Cahn-Hilliard 问题中的描述创建了网格。但是如何计算表面上的总概率，因为数据不是按顺序排列的。Fipy中是否有任何直接的方法可以计算总概率？谁能帮我解决这个问题？

标签： python-2.7fipy

Answer 1

可以通过以下方式实现域上的集成

var.cellVolumeAverage * mesh.celVolumes.sum()

或者

(var * mesh.celVolumes).sum()