No abstract uploaded!

The matrix completion problem is to recover a low-rank matrix from a subset of its entries. The main solution strategy for this problem has been based on nuclear-norm minimization which requires computing singular value decompositionsa task that is increasingly costly as matrix sizes and ranks increase. To improve the capacity of solving large-scale problems, we propose a low-rank factorization model and construct a nonlinear successive over-relaxation (SOR) algorithm that only requires solving a linear least squares problem per iteration. Extensive numerical experiments show that the algorithm can reliably solve a wide range of problems at a speed at least several times faster than many nuclear-norm minimization algorithms. In addition, convergence of this nonlinear SOR algorithm to a stationary point is analyzed.

No abstract uploaded!

In this paper, we propose and analyze a novel stabilized finite element method (FEM) for the system of generalized Stokes equations arising from the time-discretization of transient Stokes problem. The system involves a small viscosity, which is proportional to the inverse of large Reynolds number, and a large reaction coefficient, which is the inverse of small time step. The proposed stabilized FEM employs the $C^0$ piecewise linear elements for both velocity field and pressure on the same mesh and uses the residuals of the momentum equation and the divergence-free equation to define the stabilization terms. The stabilization parameters are fixed and element-independent, without a comparison of the viscosity, the reaction coefficient and the mesh size. Using the finite element solution of an auxiliary boundary value problem as the interpolating function for velocity and the $H^1$-seminorm projection for pressure, instead of the usual nodal interpolants, we derive error estimates for the stabilized finite element approximations to velocity and pressure in the $L^2$ and $H^1$ norms and most importantly, we explicitly establish the dependence of error bounds on the viscosity, the reaction coefficient and the mesh size. Our analysis reveals that this stabilized FEM is particularly suitable for the generalized Stokes system with a small viscosity and a large reaction coefficient, which has never been achieved before in the error analysis of other stabilization methods in the literature. We then numerically confirm the effectiveness of the proposed stabilized FEM. Comparisons made with other existing stabilization methods show that the newly proposed method can attain better accuracy and stability.

In this paper, we apply discontinuous Galerkin (DG) methods to solve two model equations: Krause's consensus models and pressureless Euler equations. These two models are used to describe the collisions of particles, and the distributions can be identified as density functions. If the particles are placed at a single point, then the density function turns out to be a $\delta$-function and is difficult to be well approximated numerically. In this paper, we use DG method to approximate such a singularity and demonstrate the good performance of the scheme. Since the density functions are always positive, we apply a positivity-preserving limiter to them. Moreover, for pressureless Euler equations, the velocity satisfies the maximum principle. We also construct special limiters to fulfill this requirement.