On equivalence of weak and viscosity solutions to nonlocal double phase problems with nonhomogeneous data

This work focuses on the nonhomogeneous nonlocal double phase problem \begin{align*} L_au(x)=f(x,u,D_s^p u, D_{a,t}^q u) $\text{ in }$ Ω, \end{align*} where $\mathrm{\Omega }\subset {\mathbb{R}}^N$ is a bounded domain with Lipschitz boundary, $0<s,t<1<p\le q<\mathrm{\infty }$ with $tq\le sp$ and the operator $L_a$ is defined as \begin{align*} L_a u(x)&=2\operatorname{P.V.}\int_{$\mathbb{R}$^N}|u(x)-u(y)|^{p-2}(u(x)-u(y))K_{s,p}(x,y) &\ \ \ +2\operatorname{P.V.}\int_{$\mathbb{R}$^N}a(x,y)|u(x)-u(y)|^{q-2}(u(x)-u(y))K_{t,q}(x,y)dy. \end{align*} We establish the equivalence between weak and viscosity solutions under boundedness and continuity assumptions. In addition, the local boundedness of weak solutions in some special cases on $f$ is also obtained using the notion of De Giorgi classes.