0

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

    0 <= x <= 3.2895, 0 <= y <= 2

    gh = movemesh(gh,[x, ((x == 0.)*(0.) + (0. < x <= 1.2622)*(((- 8.04576683e-07) 
    + y*(2. - (- 8.04576683e-07 + 1.82453782e-04*x + 5.34970465e-03*(x)^2 
    + 4.81687917e-02*(x)^3 + 1.03840134e-03*(x)^4 - 2.05216562e-03*(x)^5))/2. 
    + 1.82453782e-04*x + 5.34970465e-03*(x)^2 + 4.81687917e-02*(x)^3 
    + 1.03840134e-03*(x)^4 - 2.05216562e-03*(x)^5)) + (1.2622 < x <= 2.5372)*
    (((- 1.21682548) + y(2. - (- 1.21682548 + 3.46330834*x - 3.84020385*(x)^2 
    + 2.11549888*(x)^3 - 0.53002281*(x)^4 + 0.04873513*(x)^5))/2. + 3.46330834*x 
    - 3.84020385*(x)^2 + 2.11549888*(x)^3 - 0.53002281*(x)^4 + 0.04873513*(x)^5)) 
    + (2.5372 < x <= 3.1748)*(((30.84120707) + y(2. - (30.84120707 - 55.44215109*x 
    + 39.41337914*(x)^2 - 13.738279*(x)^3 + 2.36768666*(x)^4 - 0.16226731*(x)^5))/2. 
    - 55.44215109*x + 39.41337914*(x)^2 - 13.738279*(x)^3 + 2.36768666*(x)^4 
    - 0.16226731*(x)^5)) + (3.1748 < x <= 3.2895)*(((- 45.55263021) 
    + y(2. - (- 45.55263021 + 71.59032958*x - 45.04785632*(x)^2 + 14.32888663*(x)^3 
    - 2.29397151*(x)^4 + 0.14731438*(x)^5))/2. + 71.59032958*x - 45.04785632*(x)^2 
    + 14.32888663*(x)^3 - 2.29397151*(x)^4 + 0.14731438*(x)^5)))]);

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

4

0 回答 0