问题标签 [freefem++]

我正在使用 FreeFem++ 解决复杂几何上的泊松方程。我想在异构区域上解决它。ΔΦ=f中的f值在多个子区域之间变化。有没有办法做到这一点？到目前为止，我只发现了同质区域的例子。我对有限元分析比较陌生。

我正在试验FreeFem++，它基本上是一个c++带有一些附加库的编译器。

正如您在此处看到的：https ://doc.freefem.org/references/types.html一个 5 元素数组可以定义并打印为：

（请注意，real 是一种自定义类型，实际上是双精度类型）。

作为一个（新手）c++开发人员，我会real使用模板定义类，并且使用正确的构造函数，我可以通过一些调用来做同样的事情，比如real<int>.

问题

是否可以编写一个c++库来使上述FreeFem++代码工作（具体来说，我想要方[]括号，而不是<>）？

我的尝试

作为FreeFem++开源，我试着看一下代码。

https://github.com/FreeFem/FreeFem-sources/blob/master/src/fflib/array_real.cpp

我无法真正破译它，但是我的猜测是这是不可能的，实际上是FreeFem++编译器[int]将<int>.

奖金

有时我还会看到部分代码，例如

它在哪里打印一个向量，但是要访问我必须做的向量的元素v[][5]。那是什么法术？重载operator[](void)永远行不通，对吧？

我正在与您联系以尝试获得有关与非线性问题相关的问题的一些答案。

事实上，我试图将方程线性化，所以为此我取了 c1 = c10 + dc1，其中 c1 是我的变量，c10 是初始常数，dc1 是增量。

对应的代码是：

为了确定增量 dc1，我使用了一个循环：

当我编译代码时，c1 的值完全不正确，并且增加得非常快。

有人可以帮助我吗？或者给我一些关于这段代码的提示？

我希望我已经提供了足够的信息来让您理解代码，尤其是我必须编写的错误。

先感谢您。

我构建了用于将矩形网格转换为凹凸形状的离散解析表达式。我在 python 上使用样条曲线导出了表达式，我试图在 FreeFem++ 中使用相同的表达式，但它没有为我提供所需的网格转换几何（根据需要）。下面是 FreeFem++ 中使用的一段代码。解析表达式与 Python 完美配合，为网格提供所需的几何形状。gh 是指从

任何形式的帮助将不胜感激。

我在使用 Freefem++ 语法时遇到问题。我明白那个

是创建数组的方法。但是，我无法构造一个返回这个数组的函数。

任何帮助表示赞赏，

谢谢

在此处输入图像描述

其中 u : Ω → R n 是位移矢量场，C 是弹性介质的刚度张量，(u) 是位移的线性化应变张量，ρ 是质量密度，f 是外力。