lamb_util.h

一些使用广义闭合解求解 Lamb 问题过程中可能用到的辅助函数

Author

Zhu Dengda (zhudengda@mail.iggcas.ac.cn)

Date

2025-11

Functions

void grt_rayleigh1_roots(real_t nu, cplx_t y3[3])

求解如下一元三次形式的 Rayleigh 方程的根, 其中 \( \nu \) 为泊松比

\[ x^3 - \dfrac{2\nu^2 + 1}{2(1 - \nu)} x^2 + \dfrac{4\nu^3 - 4\nu^2 + 4\nu - 1}{4(1 - \nu)^2} x - \dfrac{\nu^4}{8(1-\nu)^3} = 0 \]

参数:
  • nu[in] 泊松比, (0, 0.5)

  • y3[out] 三个根,其中 y3[2] 为正根

void grt_rayleigh2_roots(real_t m, cplx_t y3[3])

求解如下一元三次形式的 Rayleigh 方程的根, 其中 \( m=\dfrac{1}{2}\dfrac{1-2\nu}{1-\nu}, \nu \) 为泊松比

\[ x^3 + \dfrac{2m - 3}{2(1 - m)} x^2 + \dfrac{1}{2(1-m)} x - \dfrac{1}{16(1-m)} = 0 \]

参数:
  • m[in] 系数 m

  • y3[out] 三个根,其中 y3[2] 为正根

cplx_t grt_evalpoly2(const cplx_t *C, const int n, const cplx_t y, const int offset)

做如下多项式求值, \( \sum_{m=0}^n C_{2m+o} y^m \)

参数:
  • C[in] 数组 C

  • n[in] 最高幂次 n

  • y[in] 自变量 y

  • o[in] 偏移量

返回:

多项式结果