ACCURATE COMPUTATION FOR IRREDUCIBLE NONSINGULAR M-MATRIX

Computing the eigenvalue of smallest modulus and its corresponding eigneveclor of an irreducible nonsingular M-matrix A is considered, It is shown that if the entries of A are known with high relative accuracy, its eigenvalue of smallest modulus and each component of the corresponding eigenvector will be determined to much higher accuracy than the standard perturbation theory suggests. An algorithm is presented to compute them with a small componentwise backward error, which is consistent with the perturbation results.

Alternative Methods of Regular and Singular Perturbation Problems

Making exact approximations to solve equations distinguishes applied mathematicians from pure mathematicians, physicists, and engineers. Perturbation problems, both regular and singular, are pervasive in diverse fields of applied mathematics and engineering. This research paper provides a comprehensive overview of algebraic methods for solving perturbation problems, featuring a comparative analysis of their strengths and limitations. Serving as a valuable resource for researchers and practitioners, it offers insights and guidance for tackling perturbation problems in various disciplines, facilitating the advancement of applied mathematics and engineering.

An Initial Perturbation Method for the Multiscale Singular Vector in Global Ensemble Prediction

Ensemble prediction is widely used to represent the uncertainty of single deterministic Numerical Weather Prediction(NWP) caused by errors in initial conditions(ICs). The traditional Singular Vector(SV) initial perturbation method tends only to capture synoptic scale initial uncertainty rather than mesoscale uncertainty in global ensemble prediction. To address this issue, a multiscale SV initial perturbation method based on the China Meteorological Administration Global Ensemble Prediction System(CMA-GEPS) is proposed to quantify multiscale initial uncertainty. The multiscale SV initial perturbation approach entails calculating multiscale SVs at different resolutions with multiple linearized physical processes to capture fast-growing perturbations from mesoscale to synoptic scale in target areas and combining these SVs by using a Gaussian sampling method with amplitude coefficients to generate initial perturbations. Following that, the energy norm,energy spectrum, and structure of multiscale SVs and their impact on GEPS are analyzed based on a batch experiment in different seasons. The results show that the multiscale SV initial perturbations can possess more energy and capture more mesoscale uncertainties than the traditional single-SV method. Meanwhile, multiscale SV initial perturbations can reflect the strongest dynamical instability in target areas. Their performances in global ensemble prediction when compared to single-scale SVs are shown to(i) improve the relationship between the ensemble spread and the root-mean-square error and(ii) provide a better probability forecast skill for atmospheric circulation during the late forecast period and for short-to medium-range precipitation. This study provides scientific evidence and application foundations for the design and development of a multiscale SV initial perturbation method for the GEPS.

THE SUPERCLOSENESS OF THE FINITE ELEMENT METHOD FOR A SINGULARLY PERTURBED CONVECTION-DIFFUSION PROBLEM ON A BAKHVALOV-TYPE MESH IN 2D

For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.

Periodic signal extraction of GNSS height time series based on adaptive singular spectrum analysis

Singular spectrum analysis is widely used in geodetic time series analysis.However,when extracting time-varying periodic signals from a large number of Global Navigation Satellite System(GNSS)time series,the selection of appropriate embedding window size and principal components makes this method cumbersome and inefficient.To improve the efficiency and accuracy of singular spectrum analysis,this paper proposes an adaptive singular spectrum analysis method by combining spectrum analysis with a new trace matrix.The running time and correlation analysis indicate that the proposed method can adaptively set the embedding window size to extract the time-varying periodic signals from GNSS time series,and the extraction efficiency of a single time series is six times that of singular spectrum analysis.The method is also accurate and more suitable for time-varying periodic signal analysis of global GNSS sites.

THE ABSENCE OF SINGULAR CONTINUOUS SPECTRUM FOR PERTURBED JACOBI OPERATORS

This paper is mainly about the spectral properties of a class of Jacobi operators(H_(c,b)u)(n)=c_(n)u(n+1)+c_(n-1)u(n-1)+b_(n)u(n),.where∣c_(n)−1∣=O(n^(−α))and b_(n)=O(n^(−1)).We will show that,forα≥1,the singular continuous spectrum of the operator is empty.

Structural Foundation and Geometry of the Material Singularity (and Its Quantum Entanglement)

In this paper we develop and study, as the second part of one more general development, the energy transmutation equation for the material singularity, previously obtained through the symmetrisation of a wave packet, that is, we develop the correlation between the terms of this equation, which accounts for the formation of matter from a previous vibrational state, and the different possible energy species. These energetic species are ascribed, in a simplified form, to the equation E¯ω=E¯k+E¯f, which allows us, through its associated phase factor, to gain an insight into the wave character of the kinetic energy and thus to attain the basis of the matter-wave, and all sorts of related phenomenologies, including that concerning quantum entanglement. The formation of the matter was previously identified as an energetic process, analogous to the kinetic one, in which finally the inertial mass is consolidated as a mass in a different phase, now, in addition, the mass of the material singularity is identified as a volumetric density of waves of toroidal geometry created in the process of singularisation or energy transfer between species, which makes it possible to establish the real relation or correspondence between the corpuscular and photonic energy equation (E=mc2=hν), i.e. to explain through m the intimate sense of the first equivalence, which explains what νis in the second one.

Chien-physics-informed neural networks for solving singularly perturbed boundary-layer problems

A physics-informed neural network(PINN)is a powerful tool for solving differential equations in solid and fluid mechanics.However,it suffers from singularly perturbed boundary-layer problems in which there exist sharp changes caused by a small perturbation parameter multiplying the highest-order derivatives.In this paper,we introduce Chien's composite expansion method into PINNs,and propose a novel architecture for the PINNs,namely,the Chien-PINN(C-PINN)method.This novel PINN method is validated by singularly perturbed differential equations,and successfully solves the wellknown thin plate bending problems.In particular,no cumbersome matching conditions are needed for the C-PINN method,compared with the previous studies based on matched asymptotic expansions.

A reweighted damped singular spectrum analysis method for robust seismic noise suppression

(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression as a low-rank reconstruction problem.However,in some cases the seismic geophones receive some erratic disturbances and the amplitudes are dramatically larger than other receivers.The presence of this kind of noise,called erratic noise,makes singular spectrum analysis(SSA)reconstruction unstable and has undesirable effects on the final results.We robustify the low-rank reconstruction of seismic data by a reweighted damped SSA(RD-SSA)method.It incorporates the damped SSA,an improved version of SSA,into a reweighted framework.The damping operator is used to weaken the artificial disturbance introduced by the low-rank projection of both erratic and random noise.The central idea of the RD-SSA method is to iteratively approximate the observed data with the quadratic norm for the first iteration and the Tukeys bisquare norm for the rest iterations.The RD-SSA method can suppress seismic incoherent noise and keep the reconstruction process robust to the erratic disturbance.The feasibility of RD-SSA is validated via both synthetic and field data examples.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Adaptive Optimal Output Regulation of Interconnected Singularly Perturbed Systems With Application to Power Systems

This article studies the adaptive optimal output regulation problem for a class of interconnected singularly perturbed systems(SPSs) with unknown dynamics based on reinforcement learning(RL).Taking into account the slow and fast characteristics among system states,the interconnected SPS is decomposed into the slow time-scale dynamics and the fast timescale dynamics through singular perturbation theory.For the fast time-scale dynamics with interconnections,we devise a decentralized optimal control strategy by selecting appropriate weight matrices in the cost function.For the slow time-scale dynamics with unknown system parameters,an off-policy RL algorithm with convergence guarantee is given to learn the optimal control strategy in terms of measurement data.By combining the slow and fast controllers,we establish the composite decentralized adaptive optimal output regulator,and rigorously analyze the stability and optimality of the closed-loop system.The proposed decomposition design not only bypasses the numerical stiffness but also alleviates the high-dimensionality.The efficacy of the proposed methodology is validated by a load-frequency control application of a two-area power system.

A SINGULAR DIRICHLET PROBLEM FOR THE MONGE-AMPÈRE TYPE EQUATION

We consider the singular Dirichlet problem for the Monge-Ampère type equation{\rm det}\D^2 u=b(x)g(-u)(1+|\nabla u|^2)^{q/2},\u<0,\x\in\Omega,\u|_{\partial\Omega}=0,whereΩis a strictly convex and bounded smooth domain inℝn,q∈[0,n+1),g∈C∞(0,∞)is positive and strictly decreasing in(0,∞)with\lim\limits_{s\rightarrow 0^+}g(s)=\infty,and b∈C∞(Ω)is positive inΩ.We obtain the existence,nonexistence and global asymptotic behavior of the convex solution to such a problem for more general b and g.Our approach is based on the Karamata regular variation theory and the construction of suitable sub-and super-solutions.

New Raw Materials Ignite the Singularity of Cosmetic Innovation

As traced by the Big Bang theory,the starting point of the universe,is called the“Singularity”in Da Ci Hai,an unabridged,comprehensive dictionary.According to cosmological reasoning,the singularity has an infinite density of matter,an infinite curvature of space and time,and it is invisible and infinite.These characteristics are analogous to the human imagination at the level of innovation.For the innovation of cosmetic raw materials,there is also the possibility of infinite evolution.For example,in recent years,the scientific research in cosmetic industry the for promoting upgrade in raw materials is quite proactive.From the raw material enterprises,down to the brand company,the investment in raw material innovation is also strengthened at a visible rate.

On the Vacuum Hydrodynamics of Moving Bodies—The Theory of General Singularity

The Theory of General Singularity is presented, unifying quantum field theory, general relativity, and the standard model. This theory posits phonons as fundamental excitations in a quantum vacuum, modeled as a Bose-Einstein condensate. Through key equations, the role of phonons as intermediaries between matter, energy, and spacetime geometry is demonstrated. The theory expands Einsteins field equations to differentiate between visible and dark matter, and revises the standard model by incorporating phonons. It addresses dark matter, dark energy, gravity, and phase transitions, while making testable predictions. The theory proposes that singularities, the essence of particles and black holes, are quantum entities ubiquitous in nature, constituting the very essence of elementary particles, seen as micro black holes or quantum fractal structures of spacetime. As the theory is refined with increasing mathematical rigor, it builds upon the foundation of initial physical intuition, connecting the spacetime continuum of general relativity with the hydrodynamics of the quantum vacuum. Inspired by the insights of Tesla and Majorana, who believed that physical intuition justifies the infringement of mathematical rigor in the early stages of theory development, this work aims to advance the understanding of the fundamental laws of the universe and the perception of reality.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Global Existence and Decay of Solutions for a Class of a Pseudo-Parabolic Equation with Singular Potential and Logarithmic Nonlocal Source

This article investigates the well posedness and asymptotic behavior of Neumann initial boundary value problems for a class of pseudo-parabolic equations with singular potential and logarithmic nonlinearity. By utilizing cut-off techniques and combining with the Faedo Galerkin approximation method, local solvability was established. Based on the potential well method and Hardy Sobolev inequality, derive the global existence of the solution. In addition, we also obtained the results of decay.

基于Singularity的超算资源容器化

深度学习往往伴随着大量的计算,属于典型的高性能计算(HPC)应用。基于Python的深度学习中往往使用PyTorch或者TensorFlow库,环境的配置与管理极为复杂,可重复性低,导致代码的升级、管理、测试、迁移成为大难题。为此,诞生了容器化技术。容器化技术以Docker为代表,通过在迁移代码的同时迁移代码运行环境,从而避免了代码在移植的时候存在繁杂的环境配置过程。但是超算服务器与普通的服务器又有很大的差别,如开发人员无root权限,测试代码时需要进行资源管理(内存,CPU)。这些问题都让Docker技术在HPC环境的应用受限,针对以上问题,使用Singularity代替Docker,在神威蓝光II型超级计算机上尝试生成了Python3.7+Numpy+Pandas+Scipy数据科学计算环境的Singularity镜像文件,在阿里云服务器上无root权限创建容器并成功运行了距离判别Python计算程序。结果表明,使用Singularity可以解决Docker在服务器上面的root权限问题,有助于深度学习环境在不同HPC服务器上面的部署。

Discrete vortex bound states with a van Hove singularity in the vicinity of the Fermi level

A theoretical study on discrete vortex bound states is carried out near a vortex core in the presence of a van Hove singularity(VHS) near the Fermi level by solving Bogoliubov–de Gennes(Bd G) equations. When the VHS lies exactly at the Fermi level and also at the middle of the band, a zero-energy state and other higher-energy states whose energy ratios follow integer numbers emerge. These discrete vortex bound state peaks undergo a splitting behavior when the VHS or Fermi level moves away from the middle of the band. Such splitting behavior will eventually lead to a new arrangement of quantized vortex core states whose energy ratios follow half-odd-integer numbers.

High impedance fault detection in distribution network based on S-transform and average singular entropy

When a high impedance fault(HIF)occurs in a distribution network,the detection efficiency of traditional protection devices is strongly limited by the weak fault information.In this study,a method based on S-transform(ST)and average singular entropy(ASE)is proposed to identify HIFs.First,a wavelet packet transform(WPT)was applied to extract the feature frequency band.Thereafter,the ST was investigated in each half cycle.Afterwards,the obtained time-frequency matrix was denoised by singular value decomposition(SVD),followed by the calculation of the ASE index.Finally,an appropriate threshold was selected to detect the HIFs.The advantages of this method are the ability of fine band division,adaptive time-frequency transformation,and quantitative expression of signal complexity.The performance of the proposed method was verified by simulated and field data,and further analysis revealed that it could still achieve good results under different conditions.

Local singularity and S–A methods for analyzing ore-producing anomalies in the Jianbiannongchang area of Heilongjiang,China

The Heilongjiang Jianbiannongchang area is located at the confluence of the Great and Lesser Xing’an Ranges.This area has a complex magmatic and tectonic evolutionary history that has resulted in a complex and diverse geological background for mineralization.In this study,isometric logarithmic ratio(ILR)transformations of Au,Cu,Pb,Zn,and Sb contents were performed in the1:50,000 soil geochemical data of the Jianbiannongchang area.Robust principal component analysis(RPCA)was conducted based on ILR transformation.The local singularity and spectrum-area(S-A)methods were used to extract information on mineralogic anomalies.The results showed that:(1)the transformed data eliminated the influence of the original data closure effect,and the PC1and PC2 information obtained by applying RPCA reflected ore-producing element anomalies dominated by Au and Cu.(2)The local singularity method can enhance the information of the local strong and weak slow anomalies.After performing local singularity analysis on PC1 and PC2,the obtained local anomalies reflected the local singularity spatial anomaly patterns related to Cu and Au mineralization in this area,which is an effective method for trapping ore-producing anomalies.(3)Furthermore,the composite anomaly decomposition of PC1 and PC2 was performed using the S-A method,and the screened anomalous and background fields reflect the ore-producing anomalies related to Cu and Au mineralization.This information is in agreement with known Cu and Au mineralization.(4)The geochemical anomalies with mineralization potential were obtained outside the known mineralization sites by integrating the information of oreproducing anomalies extracted by the local singularity and S-A methods,providing the theoretical basis and exploration direction for future exploration in the study area.