我正在尝试解决$Ax=b$线性方程组。矩阵$A=(a_{ij})$构造为

$$a_{ij}=\int_0^1\int_0^1 K(x,y)h_i(x)h_j(y) dxdy.$$ 积分的求值时间比较长，整体计算速度很慢。我可以使用 GPU 内核并行评估这些积分吗？请建议我如何同时计算这些积分？

谢谢