De Huang

Bojin Chen, De Huang, Xiangyuan Li

We present novel self-similar finite-time blowup scenarios for the 1D Hou--Luo model. We numerically demonstrate that solutions that initially satisfy certain derivative degeneracy condition can develop asymptotically self-similar finite-time blowups with singular self-similar profiles that are unbounded at some point. Moreover, this blowup phenomenon exhibits a two-stage feature: the solution first undergoes a local $L^{\infty}$ blowup at some time $\tilde{T}$, then continues in the weak sense beyond $\tilde{T}$ and develops a local $L^p$ blowup at a later time $T>\tilde{T}$ for some $p>0$. A further numerical investigation indicates that both stages are asymptotically self-similar. Finally, we extend our numerical study to the 2D Boussinesq equations and discover similar self-similar finite-time blowups with singular profiles that also exhibit a two-stage feature.

De Huang, Jonathan Niles-Weed, Rachel Ward

We analyze Oja's algorithm for streaming $k$-PCA and prove that it achieves performance nearly matching that of an optimal offline algorithm. Given access to a sequence of i.i.d. $d \times d$ symmetric matrices, we show that Oja's algorithm can obtain an accurate approximation to the subspace of the top $k$ eigenvectors of their expectation using a number of samples that scales polylogarithmically with $d$. Previously, such a result was only known in the case where the updates have rank one. Our analysis is based on recently developed matrix concentration tools, which allow us to prove strong bounds on the tails of the random matrices which arise in the course of the algorithm's execution.

Thomas Y. Hou, De Huang, Ka Chun Lam et al.

In this paper we propose a new iterative method to hierarchically compute a relatively large number of leftmost eigenpairs of a sparse symmetric positive matrix under the multiresolution operator compression framework. We exploit the well-conditioned property of every decomposition components by integrating the multiresolution framework into the Implicitly restarted Lanczos method. We achieve this combination by proposing an extension-refinement iterative scheme, in which the intrinsic idea is to decompose the target spectrum into several segments such that the corresponding eigenproblem in each segment is well-conditioned. Theoretical analysis and numerical illustration are also reported to illustrate the efficiency and effectiveness of this algorithm.