This paper is devoted to the study of a discrepancy-type characteristic -- the fixed volume discrepancy -- of the Korobov point sets in the unit cube. It was observed recently that this new characteristic allows us to obtain optimal rate of dispersion from numerical integration results. This observation motivates us to thoroughly study this new version of discrepancy, which seems to be interesting by itself. This paper extends recent results by V. Temlyakov and M. Ullrich on the fixed volume discrepancy of the Fibonacci point sets.
The smooth fixed volume discrepancy in the periodic case is studied here. It is proved that the Frolov point sets adjusted to the periodic case have optimal in a certain sense order of decay of the smooth periodic discrepancy. The upper bounds for the $r$-smooth fixed volume periodic discrepancy for these sets are established.
Solving linear systems is a ubiquitous task in science and engineering. Because directly inverting a large-scale linear system can be computationally expensive, iterative algorithms are often used to numerically find the inverse. To accommodate the dynamic range and precision requirements, these iterative algorithms are often carried out on floating-point processing units. Low-precision, fixed-point processors require only a fraction of the energy per operation consumed by their floating-point counterparts, yet their current usages exclude iterative solvers due to the computational errors arising from fixed-point arithmetic. In this work, we show that for a simple iterative algorithm, such as Richardson iteration, using a fixed-point processor can provide the same rate of convergence and achieve high-precision solutions beyond its native precision limit when combined with residual iteration. These results indicate that power-efficient computing platform consisting of analog computing devices can be used to solve a broad range of problems without compromising the speed or precision.
It is proved that the Fibonacci and the Frolov point sets, which are known to be very good for numerical integration, have optimal rate of decay of dispersion with respect to the cardinality of sets. This implies that the Fibonacci and the Frolov point sets provide universal discretization of the uniform norm for natural collections of subspaces of the multivariate trigonometric polynomials. It is shown how the optimal upper bounds for dispersion can be derived from the upper bounds for a new characteristic -- the smooth fixed volume discrepancy. It is proved that the Fibonacci point sets provide the universal discretization of all integral norms.
In this paper, we combine the nonlinear HWENO reconstruction in cite{newhwenozq} and the fixed-point iteration with Gauss-Seidel fast sweeping strategy, to solve the static Hamilton-Jacobi equations in a novel HWENO framework recently developed in cite{mehweno1}. The proposed HWENO frameworks enjoys several advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function $phi$. They are now computed from the current value of $phi$ and the previous spatial derivatives of $phi$. This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO scheme. In addition, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed-point fast sweeping method is used to update the numerical approximation. It is an explicit method and does not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton-Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some known issues, including that the iterative number is very sensitive to the parameter $varepsilon$ used in the nonlinear weights, as observed in previous studies. Finally, in order to further reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two-dimensional problems, which demonstrate the good performance of the proposed fixed-point fast sweeping HWENO methods.
Recently developed concept of dissipative measure-valued solution for compressible flows is a suitable tool to describe oscillations and singularities possibly developed in solutions of multidimensional Euler equations. In this paper we study the convergence of the first-order finite volume method based on the exact Riemann solver for the complete compressible Euler equations. Specifically, we derive entropy inequality and prove the consistency of numerical method. Passing to the limit, we show the weak and strong convergence of numerical solutions and identify their limit. The numerical results presented for the spiral, Kelvin-Helmholtz and the Richtmyer-Meshkov problem are consistent with our theoretical analysis.