Huayi Wei

p 和 u 全部是分片线性连续空间吗? 如果都是, 你可以 p[:] = 2*u**2.

这就是一个插值的操作, 当然你也可以做投影. 无论是插值还是投影, 目前这种做法都会存在误差的. 我建议下次你问关于编程的问题,最好提供一段最小的测试代码, 方便大家直接理解你的意思并提供帮助.

这是一个线性插值, 这个逼近阶是 h^2 的, 就看你非线迭代的要求是什么?   我的问题是  p 这个中间变量是必须的吗? 它如果只是你双线性形的系数, 你根本不需要插值. 只接当成系数用就可以. 比如  (p*u, v) 这样的项,  你就把 p 当成一个系数函数即可, 根本不需要做插值.  ``` from fealpy.decorator import barycentric @barycentric # 声明定义的函数传入的参数是重心坐标  def p(bc):  ...

前面邮件回复的, 格式有点乱了 ``` from fealpy.decorator import barycentric @barycentric # 声明定义的函数传入的参数是重心坐标 def p(bc): return 2*u(bc)**2 ```

质量应该是只有单元内部有啊.     ------------------ Original ------------------ From: ***@***.***>; Date:  Wed, Mar 24, 2021 09:47 AM To: ***@***.***>; Cc: ***@***.***>; ***@***.***>; Subject:  Re: [weihuayi/fealpy] 插值函数 (#33)   谢谢魏老师，我还有个关于WG的问题，就是我看程序包里面都是关于其刚度矩阵的组装，没有WG的质量矩阵组装，质量矩阵组装只有单元内函数积分块有值，其他例如单元函数和边界函数，边界函数之间的积分值为0吗？ — You are receiving this...

你看的是哪一段代码? 把相应的文件告诉我, 我不确定你看的是哪一个文件.      ------------------ Original ------------------ From: ***@***.***>; Date:  Wed, Mar 24, 2021 10:49 AM To: ***@***.***>; Cc: ***@***.***>; ***@***.***>; Subject:  Re: [weihuayi/fealpy] 插值函数 (#33)   魏老师，WG针对于自由度，我看您写的Ndof=NT+NE，但是其他是这么写的Ndof=3*NT+NE，其中3代表一个单元内部有三个基函数，是因为您这个解是在节点和边上的系数，其它的是内部和边上的系数 — You are...

我以为你说的是 FEALPy 中的代码, 在 FEALPy 中你可以看这个空间: https://github.com/weihuayi/fealpy/blob/master/fealpy/functionspace/WeakGalerkinSpace2d.py

你说的这个系数函数一个函数，也是一个数组，在FEALPy 当中， 假设 C 就是一个 np.ndarray, 就是你系数函数的值向量, 下面的代码就把 C 封装成一个有限元函数，你可以像函数一样调用它。 space = LagrangeFiniteElementSpace(mesh, p=1) F = space.function(array=C)

这些例子需要你输入参数啊。 sys 是用来管理输入参数的 另外， 要用 petsc 的话，也要安装才行, 只是现在在有些平台上安装有问题所以就没有说明 ``` pip install petsc petsc4py ```

你可以把 built 的错误信息贴出来。也许你需要安装 visual studio.