这是我博士课题中的一个小算法，用于解析形如 \( f(X,Y,Z) \) 以空间坐标为变量的文本函数并进行表达式计算。当时的目的是为了能够统一地设置流场的初始状态，由于流场的初始状态的界面往往可以通过解析函数来表示，因此设计了以空间坐标 \( (X,Y,Z) \) 为变量的函数表达式解析。通过文本函数解析避免了初始状态文件的导出，并且在后续做多相流的时候还可以拿来计算相应的初始体积分数。

针对数据量特别大的对象实例，运算过程中都要特别关注对其的拷贝操作，例如大型矩阵的乘法运算\(A_{r \times s}\cdot B_{s \times t}\) 等等。出于代码的复用性考虑，我们会将这些运算封装成函数调用。在C++新标准前，通常的做法是将对象实例通过形参列表将引用或者指针传入以避免大数据的拷贝操作，在C++17新标准之后，我们可以采用C++编译器的具名返回值优化(NRVO)功能，并搭配浅拷贝的移动构造函数来避免大数据的拷贝操作。

Congratulations to ourself for marriage certificate!