---
title: 'Quasi-Random PINNs: Enhanced Sampling for PDEs'
url: https://www.emergentmind.com/topics/quasi-random-physics-informed-neural-networks-qrpinns
type: topic
---

# Quasi-Random PINNs: Enhanced Sampling for PDEs

Searching arXiv for the primary QRPINN paper and closely related PINN collocation literature.
Physics-informed neural networks (PINNs) approximate solutions of differential equations by minimizing losses that penalize violations of governing equations, initial conditions, and boundary conditions at sampled points in the domain. Quasi-Random Physics-informed Neural Networks (QRPINNs) modify this standard PINN workflow at the sampling stage: instead of drawing collocation points directly from the domain by ordinary random sampling, they use low-discrepancy sequences, specifically quasi-Monte Carlo constructions such as Halton and Sobol’, or random minibatches drawn from pools generated by those sequences [2507.08121]. The core motivation is that the PINN loss is an empirical quadrature of a continuous residual-based objective, so the distribution of collocation points directly affects how accurately the discrete loss approximates its integral counterpart [2507.08121]. In the broader PINN literature, the baseline formulation is explicitly built around enforcing physics on “random points in the space-time domain,” which makes QRPINNs a sampling-level refinement rather than a change to the PINN architecture or residual definition [2205.08332].

## 1. Definition and conceptual position

QRPINNs are a variant of PINNs in which “the usual random sampling of collocation/training points is replaced by sampling from low-discrepancy, quasi-random point sets” [2507.08121]. In ordinary PINNs, the approximation $\boldsymbol{u}^\theta$ is trained by minimizing residual, initial-condition, and boundary-condition losses evaluated at sampled points from the corresponding domain subsets. QRPINNs “keep the PINN architecture and objective essentially unchanged; the key modification is the generation of training/collocation points” [2507.08121].

This places QRPINNs within the standard mesh-free PINN paradigm rather than in alternative physics-informed learning frameworks. The review literature describes classical PINNs as methods that “can solve a partial differential equation (PDE) by directly incorporating the PDE into the loss function of the neural network (NN) and employing automatic differentiation to represent all the differential operators,” with the physics enforced by penalizing residuals at random collocation points [2205.08332]. QRPINNs preserve this formulation and intervene only in the construction of the sampled point sets. This suggests that QRPINNs are best understood as a collocation-design modification of PINNs, analogous in role to changing the quadrature rule used to discretize a continuous residual norm.

A frequent source of confusion is the abbreviation “QR.” In QRPINNs, it refers to “quasi-random,” whereas in the adaptive collocation method based on QR Discrete Empirical Interpolation Method (QR-DEIM), “QR” refers to QR factorization with column pivoting rather than low-discrepancy sampling [2501.07700]. The two lines of work address the same broad question—how to place collocation points in PINNs—but they are not the same method.

## 2. Mathematical formulation

The standard time-dependent PDE form used in the QRPINN paper is
$$
\begin{aligned}
& \partial_{t}\boldsymbol{u}+\mathcal{N}[\boldsymbol{u}]=0, \quad t \in[0, T],\ \boldsymbol{x} \in \Omega, \\
& \boldsymbol{u}(0, \boldsymbol{x})=\boldsymbol{g}(\boldsymbol{x}), \quad \boldsymbol{x} \in \Omega, \\
& \mathcal{B}[\boldsymbol{u}]=0, \quad t \in[0, T],\ \boldsymbol{x} \in \partial \Omega,
\end{aligned}
$$
with neural approximation $\boldsymbol{u}^\theta$ and composite loss
$$
\mathcal{L}(\theta)=\mathcal{L}_{i c}(\theta)+\mathcal{L}_{b c}(\theta)+\mathcal{L}_r(\theta),
$$
where
$$
\begin{aligned}
& \mathcal{L}_r(\theta)=\frac{1}{N_r} \sum_{i=1}^{N_r}\varepsilon_r(\theta,t_r^i, \boldsymbol{x}_r^i), \quad \varepsilon_r(\theta,t,\boldsymbol{x})=\left|\partial_{t}\boldsymbol{u}^\theta\left(t,\boldsymbol{x}\right)+\mathcal{N}\left[\boldsymbol{u}^\theta\right]\left(t,\boldsymbol{x}\right)\right|^2,\\
& \mathcal{L}_{i c}(\theta)=\frac{1}{N_{i c}} \sum_{i=1}^{N_{i c}}\varepsilon_{ic}\left(\boldsymbol{x}_{i c}^i\right), \quad \varepsilon_{ic}\left(\boldsymbol{x}\right)=\left|\boldsymbol{u}^\theta\left(0, \boldsymbol{x}\right)-\boldsymbol{g}\left(\boldsymbol{x}\right)\right|^2,\\
& \mathcal{L}_{b c}(\theta)=\frac{1}{N_{b c}} \sum_{i=1}^{N_{b c}}\varepsilon_{bc}\left(t_{b c}^i, \boldsymbol{x}_{b c}^i\right), \quad \varepsilon_{bc}\left(t, \boldsymbol{x}\right)=\left|\mathcal{B}\left[\boldsymbol{u}^\theta\right]\left(t, \boldsymbol{x}\right)\right|^2.
\end{aligned}
$$
In standard PINNs, the point sets $\boldsymbol{x}_{ic}, \boldsymbol{x}_{bc}, \boldsymbol{x}_r$ are typically sampled uniformly at random from the corresponding domain pieces [2507.08121].

The QRPINN paper interprets these finite sums as empirical approximations to an “infinite-sample” objective
$$
\begin{aligned}
\mathcal{L}_{int}(\theta):= & \int_t\int_{\Omega} \left|\partial_{t}u^\theta\left(t, \boldsymbol{x}\right)+\mathcal{N}\left[u^\theta\right]\left(t, \boldsymbol{x}\right)\right|^2\mathrm{d}\boldsymbol{x}\mathrm{d}t \\
& +\int_\Omega \left|u^\theta\left(0, \boldsymbol{x}\right)-g\left(\boldsymbol{x}\right)\right|^2\mathrm{d}\boldsymbol{x} \\
& + \int_t\int_{\partial\Omega}\left|\mathcal{B}\left[u^\theta\right]\left(t, \boldsymbol{x}\right)\right|^2\mathrm{d}\boldsymbol{x}\mathrm{d}t,
\end{aligned}
$$
with
$$
\theta^*_\infty := \arg \min_\theta \mathcal{L}_{int}(\theta).
$$
For the residual term, the paper writes the empirical average as converging to the corresponding space-time integral:
$$
\frac{1}{N_r} \sum_{i=1}^{N_r}\left|\partial_{t}u\left(t_r^i, \boldsymbol{x}_r^i\right)+\mathcal{N}\left[u\right]\left(t_r^i, \boldsymbol{x}_r^i\right)\right|^2 \rightarrow \int_t\int_{\Omega} \left|\partial_{t}u\left(t, \boldsymbol{x}\right)+\mathcal{N}\left[u\right]\left(t, \boldsymbol{x}\right)\right|^2\mathrm{d}\boldsymbol{x}\mathrm{d}t.
$$
This loss-as-quadrature perspective is the formal entry point for QRPINNs [2507.08121].

The broader PINN review presents the same structure in a generic PDE form on $\Omega$, with residual points $\mathcal{T}_f$ and boundary points $\mathcal{T}_b$ and total loss
$$
\mathcal{L}(\boldsymbol{\theta} ; \mathcal{T})=w_{f} \mathcal{L}_{f}\left(\boldsymbol{\theta} ; \mathcal{T}_{f}\right)+w_{b} \mathcal{L}_{b}\left(\boldsymbol{\theta} ; \mathcal{T}_{b}\right),
$$
where the paper states explicitly: “To evaluate the residual of the PDE and the boundary/initial conditions, we introduce a set of random points inside the domain ($\mathcal{T}_f$) and another random set of points on the boundary ($\mathcal{T}_b$)” [2205.08332]. In a QRPINN, the loss decomposition is not changed; only the generation of $\mathcal{T}_f$ and $\mathcal{T}_b$ is changed.

## 3. Quasi-random sampling and theoretical rationale

The specific low-discrepancy sequences implemented in the QRPINN paper are Halton and Sobol’ [2507.08121]. For Halton, with prime bases $p_1,\dots,p_d$, the points are
$$
\boldsymbol{x}_i=\left(\phi_{p_1}(i), \phi_{p_2}(i), \ldots, \phi_{p_d}(i)\right), \quad i=0,1, \ldots,
$$
with
$$
\phi_b(i):=\sum_{a=1}^{\infty}\frac{i_a}{b^a}.
$$
For Sobol’, the paper gives
$$
t_{i, j}=i_0 v_{1, j} \oplus i_1 v_{2, j} \oplus \cdots \oplus i_{r-1} v_{r, j},
$$
where $i=i_0+2 i_1+\cdots+2^{r-1} i_{r-1}$ is the dyadic expansion and
$$
v_{k,j}:=\frac{m_{k,j}}{2^k},
$$
with $m_{k,j}$ defined recursively from primitive polynomials over $\mathbb{Z}_2[x]$ [2507.08121].

The paper’s theoretical motivation begins from the Monte Carlo approximation
$$
I\left(f\right):=\int_\Omega f(\boldsymbol{x})\mathrm{d}\boldsymbol{x}\approx I_{\text{MC}}(f):=\frac{1}{N}\sum_{i=1}^N f(\boldsymbol{x}_i),
$$
for which the PINN paper emphasizes the usual Monte Carlo convergence rate $\mathcal{O}(N^{-1/2})$ [2507.08121]. It also gives
$$
\operatorname{Var}\left[I_{\text{MC}}(f)\right]=\frac{I(f^2)-I\left(f\right)^2}{N}.
$$
For quasi-Monte Carlo, the paper states a better rate, generally
$$
\mathcal{O}(N^{-(1-\epsilon)}), \qquad \epsilon\in(0,1),
$$
and in weighted settings
$$
\mathcal{O}(N^{-1/p}), \qquad p\in[1,2].
$$
For Halton and Sobol’, it mentions the classical form
$$
\mathcal{O}\!\left(\frac{(\log N)^d}{N}\right).
$$

The paper reproduces a discrepancy-based inequality. In one dimension,
$$
E_{\text{QMC}} =\left|\int_0^1f(x)\mathrm{d}x-\frac{1}{N}\sum_{i=1}^N f(x_i)\right| \leq \|\Delta p\|_{L_q} \|f^{\prime}\|_{L_q},
$$
and for $p=\infty, q=1$,
$$
E_{\text{QMC}}\leq V(f)D_N^*\left(x_1,\cdots x_N\right),
$$
where $D_N^*$ is the star discrepancy. The star discrepancy for a point set $P=\{\boldsymbol{x}_i\}$ is defined by
$$
D^*_N(P)=\sup_{u_1,\cdots,u_d\in[0,1]}\left|\frac{A(J;P)}{N}-\lambda_d(J)\right|,
$$
where $J=\prod_{i=1}^d[0,u_i)$, $\lambda_d(J)=\prod_{i=1}^d u_i$, and
$$
A(J,P)=\sum_{n=1}^N \mathbbm{1}_{J}(\boldsymbol{x}_n).
$$
The intended conclusion is that low-discrepancy point sets reduce the mismatch between empirical point coverage and uniform volume, thereby improving the approximation of the continuous PINN objective by the discrete loss [2507.08121].

This argument is tied to a theorem in the QRPINN paper:
$$
\textbf{Theorem 1.}\quad \text{If the PDE system is well defined and the solution } u \in \mathcal{H} \text{ is a smooth function, then there exists a constant } r,
$$
such that for $u^\theta \in B_r(u^{\theta^*})$,
$$
\mathbb{E}\left[\|u^{\theta^*_\infty}-u^{\theta^*}\|_{L^2(\boldsymbol{\Omega})}\right] = \mathcal{O}\left(\mathbb{E}\left[\left|\mathcal{L}_{int}(\theta)-\mathcal{L}(\theta)\right|\right]\right).
$$
The paper uses this to link PINN approximation quality to quadrature error: once network approximation error is controlled, the gap between the empirical minimizer and the ideal integral-loss minimizer is of the same order as the loss-approximation error [2507.08121]. This suggests that improving the sampling rule may improve the overall solution by improving the quadrature stage.

## 4. Randomized quasi-Monte Carlo training

The QRPINN paper does not advocate training on one fixed deterministic low-discrepancy set at every epoch. Instead, it introduces Random Quasi Monte Carlo (RQMC): a deterministic low-discrepancy pool is generated first, and minibatches are then sampled randomly from that pool during training [2507.08121]. The authors summarize the method by stating that QRPINNs “randomly sample points from low-discrepancy sequences rather than the input domain” [2507.08121].

This design reflects a specifically machine-learning-oriented compromise. Fully deterministic QMC preserves discrepancy properties, but minibatch stochasticity is often useful in neural optimization. The paper therefore frames QRPINNs as standard PINNs trained with an RQMC sampling scheme rather than as fully deterministic quadrature solvers.

The associated theorem is stated as:
$$
\textbf{Theorem 2.}\quad \text{For a given deterministic sequence } P_{N_{total}}=\{\boldsymbol{x}_i\}_{i=1}^{N_{total}},
$$
if $P_N=\{\boldsymbol{\tilde{x}}_j\}_{j=1}^N$ is a subset sampled from $P_{N_{total}}$, and if
$$
E_{\text{QMC}}=\left|\int f(\boldsymbol{x})\mathrm{d}\boldsymbol{x}-\frac{1}{N_{total}}\sum_{i=1}^{N_{total}} f(\boldsymbol{x}_i)\right|= \mathcal{O}(N_{total}^{-(1-\epsilon)})
$$
for some $\epsilon \in (0,1)$, and
$$
E_{\text{RQMC}}=\left|\int f(\boldsymbol{x})\mathrm{d}\boldsymbol{x}-\frac{1}{N}\sum_{j=1}^{N} f(\boldsymbol{\tilde{x}}_j)\right|,
$$
then
$$
E_{\text{RQMC}}\leq \frac{dD}{2k N_{total}}+dD \left(1-k\right)+C N_{total}^{-(1-\epsilon)},
$$
for constants $C,D>0$, $k=N/N_{total}$, and dimensionality $d$ [2507.08121]. The paper notes notation inconsistencies in the theorem statement, but the intended message is that random minibatching from a large low-discrepancy pool can preserve favorable asymptotic behavior.

The paper further adds a combinatorial coverage argument:
$$
P(\mathcal{X}= 0)=\frac{\sum_{i=0}^N(-1)^i\binom{N}{i}\left(\binom{N-i}{N_b}\right)^s}{\left(\binom{N}{N_b}\right)^s},
$$
and states that this is “almost 1” when $s=20$ and $N_b=0.1N$, which it uses to argue that $k=1$ is practical and that the error convergence of RQMC is almost the same as QMC [2507.08121].

A notable practical hyperparameter is the sampling-pool scale $N_{scale}$. In the main 100D experiments the QRPINN paper sets $N_{scale}=10$ and trains for 30 epochs [2507.08121]. It also reports that performance depends on $N_{scale}$, with Halton more sensitive than Sobol’, and emphasizes that using a larger quasi-random pool with stochastic minibatching can outperform fully deterministic training [2507.08121].

## 5. Empirical performance and dimensional regime

The QRPINN paper studies steady Poisson, steady Allen-Cahn, steady Sine-Gordon, and time-dependent high-dimensional PDEs coupled with the Stochastic Taylor Derivative Estimator (STDE), including semi-linear heat, time-dependent Allen-Cahn, and time-dependent Sine-Gordon [2507.08121]. The steady experiments consider $d=3,10,100$, while the STDE experiments extend to dimensions $10^2,10^3,10^4$ [2507.08121]. Baselines are Vanilla PINNs with uniform random sampling, RAD, ACLE, Halton-based QRPINN, and Sobol’-based QRPINN [2507.08121].

The empirical pattern is dimension-dependent. In lower dimensions, QRPINNs are competitive but not universally dominant. For example, in 3D Allen-Cahn the best result reported is Vanilla at $2.15\times 10^{-4}\pm 3.78\times 10^{-5}$, with Halton very close at $2.54\times 10^{-4}\pm 3.72\times 10^{-5}$; in 10D Allen-Cahn, Vanilla is again best at $2.07\times 10^{-3}\pm 3.40\times 10^{-5}$, while the other methods are described as very close [2507.08121]. In 3D Sine-Gordon, Sobol’ is best at $1.14\times 10^{-4}\pm 6.48\times 10^{-6}$, whereas in 10D Sine-Gordon, Vanilla is best at $2.09\times 10^{-3}\pm 4.08\times 10^{-5}$, with Halton and Sobol’ only slightly worse [2507.08121].

The paper identifies stronger evidence in 100D. Compared with the best result among Vanilla and RAD, QRPINNs improve by 22.6% for Allen-Cahn and 32.8% for Sine-Gordon [2507.08121]. The supporting values given in the appendix are:

| Problem | Setting | Comparison |
|---|---|---|
| Allen-Cahn | 100D, 1000 points | Vanilla $1.69\times 10^{-2}$, RAD $1.69\times 10^{-2}$, Halton $1.30\times 10^{-2}$ |
| Sine-Gordon | 100D, 1000 points | Vanilla $1.09\times 10^{-2}$, RAD $1.09\times 10^{-2}$, Halton $7.33\times 10^{-3}$ |

The STDE results push the method into much higher dimensions and are the strongest support for the claim that QRPINNs become more attractive as $d$ grows. The paper reports the following examples [2507.08121]:

| Problem | Dimension | Reported values |
|---|---:|---|
| Allen-Cahn | $10^2$ | Vanilla $8.12\times 10^{-3}$, Halton $2.54\times 10^{-3}$, Sobol $3.01\times 10^{-3}$ |
| Allen-Cahn | $10^4$ | Vanilla $1.55\times 10^{-2}$, Sobol $3.48\times 10^{-3}$ |
| Heat | $10^2$ | Vanilla $1.45\times 10^{-2}$, Sobol $1.05\times 10^{-2}$ |
| Sine-Gordon | $10^4$ | Vanilla $1.19\times 10^{-1}$, Sobol $3.28\times 10^{-2}$ |

The paper also records failure cases. For Sine-Gordon at $d=10^2$, Vanilla is best and QRPINNs do not improve it; for Heat at $d=10^3$ and $10^4$, neither Halton nor Sobol beats Vanilla [2507.08121]. This matters because it limits the scope of the empirical claim: QRPINNs are not presented as uniformly superior, but as increasingly favorable in high-dimensional regimes.

The evaluation metric throughout is relative error, reported as mean $\pm$ standard deviation over repeated runs, with experiments marked as failures when the relative error exceeds $10^{-1}$ [2507.08121].

## 6. Relationship to adaptive collocation methods

The literature surrounding QRPINNs repeatedly distinguishes static space-filling designs from adaptive refinement. The review article on scalable PINNs describes the classical baseline entirely in terms of random collocation points and discusses extensions such as self-adaptive weighting, gradient-enhanced PINNs, and domain decomposition, but “there is no explicit discussion of quasi-random sampling, low-discrepancy sequences, Latin hypercube sampling, stratified sampling, importance sampling, residual-based adaptive refinement, or quadrature-based collocation” [2205.08332]. This makes the review a baseline reference for standard random-point PINNs rather than a source on QRPINNs themselves.

The adaptive QR-DEIM paper is more directly relevant to sampling strategy. It distinguishes fixed methods—uniform random sampling, equispaced grids, Latin hypercube sampling, Sobol sequences, and Hammersley sequences—from adaptive ones, and argues that fixed methods are “agnostic to the evolving PDE residual” and can under-resolve regions with high gradients, phase boundaries, shocks, localized complex dynamics, and later times in time-dependent problems [2501.07700]. Its proposed method is not quasi-random; instead, it builds a residual snapshot matrix over a training window, compresses it with SVD, and selects informative points by QR with column pivoting [2501.07700].

The QR-DEIM paper’s core workflow is based on a training set $\mathcal{T}$ and a disjoint snapshot set $\mathcal{S}$, with residual snapshots collected over $P$ iterations,
$$
\mathbf{Y} = \big[ \mathbf{r}_1(\mathcal{S}),\dots,\mathbf{r}_P(\mathcal{S}) \big] \in \mathbb{R}^{N_{\text{snapshots}}\times P},
$$
followed by
$$
\mathbf{Y} = \mathbf{V}\Sigma\mathbf{W}^T
$$
and QR with column pivoting on the reduced basis transpose,
$$
\hat{\mathbf{V}}^T = \mathbf{Q} \begin{pmatrix} \mathbf{R} & \mathbf{0} \end{pmatrix} \Pi^T.
$$
The selected snapshot-point indices are then added to the training set [2501.07700]. The paper argues that this improves accuracy over uniform random sampling and representative residual-based adaptive methods on wave, Allen-Cahn, and Burgers’ equations, and reports that the method tends to place more points near final time, near phase boundaries in Allen-Cahn, and around the shock in Burgers while still maintaining some global coverage [2501.07700].

For QRPINNs, the main significance of this literature is comparative rather than definitional. The QRPINN paper itself proposes combining low-discrepancy pools with RAD-style adaptive sampling and reports that, for 100D Poisson, replacing the RAD candidate pool by low-discrepancy points often helps. It gives, for example, $\alpha=0.1$: Vanilla+RAD $3.19\times 10^{-3}$, Halton+RAD $3.14\times 10^{-3}$, Sobol+RAD $3.11\times 10^{-3}$; and for $\alpha=1$: Vanilla+RAD $3.13\times 10^{-2}$, Halton+RAD $2.36\times 10^{-2}$, Sobol+RAD $3.15\times 10^{-2}$ [2507.08121]. It also states that low-discrepancy pools help RAD more consistently than ACLE, while for ACLE they can be invalid and decrease accuracy [2507.08121].

This suggests a division of labor. QRPINNs address the quality of global initial coverage through low discrepancy, whereas adaptive methods address evolving localized difficulty during training. The adaptive collocation paper makes this logic explicit by arguing that “static space-filling alone is not enough for difficult PDEs” and that collocation points should concentrate near shocks, interfaces, sharp fronts, localized layers, and late-time regions [2501.07700]. A plausible implication is that quasi-random collocation and residual-driven adaptivity are complementary rather than mutually exclusive.

## 7. Computational characteristics, limitations, and scope

The QRPINN paper argues that generating low-discrepancy pools is inexpensive relative to PINN training. It reports roughly linear runtime scaling in the number of points and linear memory scaling in both $N$ and $d$, with Sobol’ more sensitive to dimensionality in runtime than Halton at high dimension [2507.08121]. For example, at $N=1{,}024{,}000$, Halton generation time is about $147$ ms at $d=3$, $297$ ms at $d=10$, and $953$ ms at $d=50$; Sobol’ is about $12.7$ ms at $d=3$, $48.7$ ms at $d=10$, and $208$ ms at $d=50$ [2507.08121]. The overhead is therefore described as small but not zero.

The method’s limitations are also stated directly. First, only two classical low-discrepancy sequences, Halton and Sobol’, are implemented; modified Halton, Niederreiter, FNN, component-by-component constructions, random lattice rules, and weighted QMC are suggested as future possibilities but are not tested [2507.08121]. Second, generating very large quasi-random pools becomes increasingly costly in extremely high dimensions, so the paper does not claim efficiency indefinitely as $d\to\infty$ [2507.08121]. Third, QRPINNs are “not universally dominant”: in lower-dimensional problems, especially when locality is strong, adaptive samplers such as RAD or ACLE may outperform them [2507.08121]. Fourth, combining low-discrepancy pools with all adaptive samplers is not always beneficial; ACLE is given as a negative example [2507.08121]. Fifth, the paper notes that some theorem statements and bounds are “notationally rough,” so the theoretical claim is stronger at the level of intended mechanism than at the level of formal polish [2507.08121].

The broader PINN literature reinforces that sampling is only one component of a larger systems problem. The review article identifies “the large computational cost associated with the training of the neural networks, especially for forward multi-scale problems” as a major concern, and emphasizes multiscale difficulty, stiffness, loss-weight choice, and scalability through domain decomposition such as cPINN and XPINN [2205.08332]. This suggests that QRPINNs target one recognized bottleneck—collocation quality—without claiming to solve all PINN failure modes.

A concrete non-QRPINN example helps delimit the expected regime of benefit. In the Teukolsky-equation PINN for Kerr black-hole quasi-normal modes, the collocation scheme is a fixed uniform grid with $N_x=N_u=100$ over two one-dimensional domains, and the authors explicitly state that “a uniform grid is not a requirement” [2212.06103]. Because that problem involves one-dimensional radial and angular variables, the paper notes that uniform grids are already strong baselines; this suggests that the largest gains from quasi-random sampling may appear not in simple 1D settings, but in higher-dimensional or more stochastic PINN regimes [2212.06103]. The QRPINN paper’s own results, which are strongest at $d=100$ and beyond, are consistent with that pattern [2507.08121].

In the current literature, QRPINNs are therefore best characterized as a minimal modification of PINNs motivated by quasi-Monte Carlo integration: retain the standard residual-based objective and neural architecture, but replace direct random domain sampling with random minibatching from low-discrepancy pools such as Halton or Sobol’ [2507.08121]. The published evidence indicates that this modification is especially relevant for high-dimensional PDEs, where the mismatch between empirical collocation coverage and the continuous residual objective becomes increasingly consequential [2507.08121; 2205.08332].

Source: https://www.emergentmind.com/topics/quasi-random-physics-informed-neural-networks-qrpinns