Source author record

Zeyao Mo

Zeyao Mo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Computation Computational Geometry cond-mat.stat-mech Numerical Analysis

Catalog footprint

What is connected

3works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A New Subspace Iteration Algorithm for Solving Generalized Eigenvalue Problems

It is needed to solve generalized eigenvalue problems (GEP) in many applications, such as the numerical simulation of vibration analysis, quantum mechanics, electronic structure, etc. The subspace iteration is a kind of widely used algorithm to solve eigenvalue problems. To solve the generalized eigenvalue problem, one kind of subspace iteration method, Chebyshev-Davidson algorithm, is proposed recently. In Chebyshev-Davidson algorithm, the Chebyshev polynomial filter technique is incorporated in the subspace iteration. In this paper, based on Chebyshev-Davidson algorithm, a new subspace iteration algorithm is constructed. In the new algorithm, the Chebyshev filter and inexact Rayleigh quotient iteration techniques are combined together to enlarge the subspace in the iteration. Numerical results of a vibration analysis problem show that the number of iteration and computing time of the proposed algorithm is much less than that of the Chebyshev-Davidson algorithm and some typical GEP solution algorithms. Furthermore, the new algorithm is more stable and reliable than the Chebyshev-Davidson algorithm in the numerical results.

preprint2016arXiv

Information Splitting for Big Data Analytics

Many statistical models require an estimation of unknown (co)-variance parameter(s) in a model. The estimation usually obtained by maximizing a log-likelihood which involves log determinant terms. In principle, one requires the \emph{observed information}--the negative Hessian matrix or the second derivative of the log-likelihood---to obtain an accurate maximum likelihood estimator according to the Newton method. When one uses the \emph{Fisher information}, the expect value of the observed information, a simpler algorithm than the Newton method is obtained as the Fisher scoring algorithm. With the advance in high-throughput technologies in the biological sciences, recommendation systems and social networks, the sizes of data sets---and the corresponding statistical models---have suddenly increased by several orders of magnitude. Neither the observed information nor the Fisher information is easy to obtained for these big data sets. This paper introduces an information splitting technique to simplify the computation. After splitting the mean of the observed information and the Fisher information, an simpler approximate Hessian matrix for the log-likelihood can be obtained. This approximated Hessian matrix can significantly reduce computations, and makes the linear mixed model applicable for big data sets. Such a spitting and simpler formulas heavily depends on matrix algebra transforms, and applicable to large scale breeding model, genetics wide association analysis.

preprint2010arXiv

Cluster Identification and Characterization of Physical Fields

The description of complex configuration is a difficult issue. We present a powerful technique for cluster identification and characterization. The scheme is designed to treat with and analyze the experimental and/or simulation data from various methods. Main steps are as follows. We first divide the space using face or volume elements from discrete points. Then, combine the elements with the same and/or similar properties to construct clusters with special physical characterizations. In the algorithm, we adopt administrative structure of hierarchy-tree for spatial bodies such as points, lines, faces, blocks, and clusters. Two fast search algorithms with the complexity are realized. The establishing of the hierarchy-tree and the fast searching of spatial bodies are general, which are independent of spatial dimensions. Therefore, it is easy to extend the skill to other fields. As a verification and validation, we treated with and analyzed some two-dimensional and three-dimensional random data.