Weak adversarial tensor neural network for high-dimensional Fokker-Planck equations

Solving high-dimensional stationary Fokker-Planck (SFP) equations remains challenging. Recently, deep learning methods have shown great potential for solving high-dimensional SFP equations due to their universal approximation and mesh-free training. However, the high-dimensional integrals involved in loss computation and normalization can be expensive to approximate accurately using Monte Carlo integration. To address this issue, we propose a weak adversarial tensor neural network (WATNN) method for solving high-dimensional SFP equations, which combines a weak adversarial (WA) formulation with the tensor neural network (TNN). The WA formulation avoids the cross terms introduced by squaring the strong residual, reducing the number of integrals to be computed. When both the primal network and the adversarial test network are set as TNNs, the TNN quadrature scheme provides accurate and efficient normalization and loss computation with relatively few fixed quadrature points, avoiding the statistical sampling error that arises when using Monte Carlo integration. We demonstrate the effectiveness numerically on four high-dimensional SFP equations. For the 6D unimodal, 6D multi-modal, 10D multi-modal and 20D multi-modal problems, WATNN achieves the mean relative errors of 1.50%, 1.13%, 1.20% and 7.93% respectively, in only a few minutes.