[Test]

Goal: $$ Q_{k}=\sum_{l=-n}^{n} \hat{f}{l} \hat{f}{k-l} \hat{\beta}(l, k-l),\quad k = -n,\dots,n $$

​

3d example

[G. Dimarco and L. Pareschi 2014] Section 5.1 $$ \hat{Q}{k}=\sum{l=-N }^{N} \hat{f}{l} \hat{f}{k-l} \hat{\beta}(l, k-l), \quad k=-n, \ldots, n $$ where the Boltzmann kernel modes $\hat{\beta}(l, m)=\hat{B}(l, m)-\hat{B}(m, m)$ are now given by $$ \hat{B}(l, m)=\int_{\mathcal{B}{0}(2 \lambda \pi)} \int{\mathbb{S}^{2}}|q| \sigma(|q|, \cos \theta) e^{-i\left(l \cdot q^{+}+m \cdot q^{-}\right)} d \omega d q $$ In the VHS case, $|q| \sigma(|q|, \cos \theta)=C_{\alpha}|q|^{\alpha}$ where we can chose $C_{\alpha}=\left((4 \pi)^{2}(2 \lambda \pi)^{3+\alpha}\right)^{-1}$. Furthermore, the integral can be reduces to a one-dimensional integral (Pareschi and Russo 2000b) $$ \hat{B}(l, m)=\int_{0}^{1} r^{2+\alpha} \operatorname{Sinc}(\xi r) \operatorname{Sinc}(\eta r) d r=F_{\alpha}(\xi, \eta) $$ where $\xi=|l+m| \lambda \pi, \eta=|l-m| \lambda \pi$, $\lambda = 2/(3+\sqrt{2})$. We report the expressions for the case of Maxwell molecules $\alpha = 0$ and hard spheres $\alpha = 1$ $$ \begin{aligned} F_{0}(\xi, \eta) &=\frac{p \sin (q)-q \sin (p)}{2 \xi \eta p q} \ F_{1}(\xi, \eta) &=\frac{q \sin (q)+\cos (q)}{2 \xi \eta q^{2}}-\frac{p \sin (p)+\cos (p)}{2 \xi \eta p^{2}}-\frac{2}{p^{2} q^{2}} \end{aligned} $$ where $p=(\xi+\eta), q=(\xi-\eta)$.

Special cases $$ \begin{aligned} F_0(\xi,0) &= \frac{-\xi \cos (\xi) + \sin (\xi)}{\xi^3}\ F_0(0,0) &= \frac{1}{3}\ F_0(\xi,\xi) &= \frac{\xi-\cos(\xi)\sin(\xi)}{2\xi^3}\ F_1(\xi,0) &= \frac{-2-(-2+\xi^2)\cos(\xi) + 2\xi \sin (\xi)}{\xi^4}\ F_1(0,0) &= \frac{1}{4}\ F_1(\xi,\xi) &= -\frac{-1-2\xi^2+\cos(2\xi) + 2\xi \sin (2\xi)}{8\xi^4} \end{aligned} $$

1. initial guess
2. fftw3 interface
3. 1rk-4
4. nofiltering: for smooth case, it enjoys higher accuracy
5. change to array for faster evaluation

Remarks

1. $J_{max} = N/2 + c, c < N/2$

 ./testBKW3D 17 0.000855 0.5 490 

• entropy tu
• 比较误差 dt 取非常小，
• 不加entropy fix

Setup