看一下 FEniCS 的 stokes-iterative 演示,在使用求解线性系统之后
U = Function(W)
solver.solve(U.vector(), bb)
解决方案U分为:
u, p = U.split()
是否有一个函数可以让我获得U包含u解决方案和p解决方案的索引?例如,如果U =[u1 p1 u2 p2,...]那么我想要数组
indices_u = [0,2,4,...]
indices_p = [1,3,5,...]
我猜它在 FEniCS 某处,因为它应该用于拆分解决方案?我想知道如何在 C++ 和 python 中做到这一点。谢谢!