Stokeslet和stresslet核通常用于边界元模拟和慢粘性流动的奇异性方法。通过直接求和计算Stokeslet和Stresslet集合所诱导的速度需要$\mathcal{O}(N^2)$操作,其中$N$是系统大小。目前的工作开发了一种用于3D Stokeslets和Stresslet的树代码算法,该算法将成本降低到$\mathcal{O}(N\log N)$。粒子被划分为簇的层次,并通过远场笛卡尔-泰勒近似计算分离良好的粒子-簇相互作用。为了提高计算效率,对近似值中的项进行了收缩。串行和并行结果显示了树码在几个测试用例中的性能。特别是,该方法具有相对简单的结构和较低的内存使用率,这提高了大型系统的并行效率。