Let a_(1),a_(2),a_(3)be nonzero integers with gcd(a_(1),a_(2),a_(3))=1,and let k be any positive integer,K=max[3,|a_(1)|,|a_(2)|,|a_(3)|,k].Suppose that l_(1),l_(2),l_(3)are integers each coprime to k.Suppose further ...Let a_(1),a_(2),a_(3)be nonzero integers with gcd(a_(1),a_(2),a_(3))=1,and let k be any positive integer,K=max[3,|a_(1)|,|a_(2)|,|a_(3)|,k].Suppose that l_(1),l_(2),l_(3)are integers each coprime to k.Suppose further that b is any integer satisfying some necessary congruent conditions.The solvability of linear equation a_(1)p_(1)+a_(2)p_(2)+a_(3)p_(3)=b(p_(j)=l_(j)(mod k),1≤j≤3)with prime variables pi,p_(2),ps is investigated.It is proved that if ai,a_(2),a_(3)are all positive,then the above equation is solvable whenever b≥K^(25);if a,a_(2),a_(3)are not all of the same sign,then the above equation has a solution p_(1),p_(2),p_(3)satisfying max(p_(1),p_(2),p_(3))≤3|b|+K^(25).展开更多
Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremend...Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.展开更多
The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditio...The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods. Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.展开更多
In this article, we study the complex oscillation problems of entire solutions to homogeneous and nonhomogeneous linear difference equations, and obtain some relations of the exponent of convergence of zeros and the o...In this article, we study the complex oscillation problems of entire solutions to homogeneous and nonhomogeneous linear difference equations, and obtain some relations of the exponent of convergence of zeros and the order of growth of entire solutions to complex linear difference equations.展开更多
In this paper,we mainly investigate the dynamical properties of entire solutions of complex differential equations.With some conditions on coefficients,we prove that the set of common limiting directions of Julia sets...In this paper,we mainly investigate the dynamical properties of entire solutions of complex differential equations.With some conditions on coefficients,we prove that the set of common limiting directions of Julia sets of solutions,their derivatives and their primitives must have a definite range of measure.展开更多
In this paper, we mainly investigate entire solutions of complex differential equations with coefficients involving exponential functions, and obtain the dynamical properties of the solutions, their derivatives and pr...In this paper, we mainly investigate entire solutions of complex differential equations with coefficients involving exponential functions, and obtain the dynamical properties of the solutions, their derivatives and primitives. With some conditions on coefficients, for the solutions, their derivatives and their primitives, we consider the common limiting directions of Julia set and the existence of Baker wandering domain.展开更多
In this article, the authors study the growth of certain second order linear differential equation f″＋A（z）f′＋B（z）f=0 and give precise estimates for the hyperorder of solutions of infinite order. Under similar ...In this article, the authors study the growth of certain second order linear differential equation f″＋A（z）f′＋B（z）f=0 and give precise estimates for the hyperorder of solutions of infinite order. Under similar conditions, higher order differential equations will be considered.展开更多
Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterativ...Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.展开更多
For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smoot...For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.展开更多
In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-m...In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.展开更多
Lyapunov equation is one of the most basic and important equations in control theory, which has various applications in, e.g., stability analysis and robust analysis of linear control systems. Inspired by the recent p...Lyapunov equation is one of the most basic and important equations in control theory, which has various applications in, e.g., stability analysis and robust analysis of linear control systems. Inspired by the recent progresses of quantum algorithms, we find that solving Lyapunov equation can be exponentially accelerated by quantum algorithms rather than traditional classical algorithms. Our algorithm is more efficient especially when the system matrix is sparse and has a low condition number. The results presented in this paper open up new dimensions of research in controlling classical system by quantum information processors, which has rarely been considered in the existing literature.展开更多
In this paper, we investigate the growth of solutions of higher order linear differential equations with meromorphic coefficients. Under certain conditions, we obtain precise estimation of growth order and hyper-order...In this paper, we investigate the growth of solutions of higher order linear differential equations with meromorphic coefficients. Under certain conditions, we obtain precise estimation of growth order and hyper-order of solutions of the equation.展开更多
A new direct method for solving unsymmetrical sparse linear systems（USLS） arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin （MLPG） method ...A new direct method for solving unsymmetrical sparse linear systems（USLS） arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin （MLPG） method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.展开更多
Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-ve...Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-vector product,matrix-matrix product,the sum of two matrices,and the calculation of determinant and inverse matrix.We encode the matrix entries into the probability amplitudes of the pure initial states of senders.After applying proper unitary transformation to the complete quantum system,the desired result can be found in certain blocks of the receiver’s density matrix.These quantum protocols can be used as subroutines in other quantum schemes.Furthermore,we present an alternative quantum algorithm for solving linear systems of equations.展开更多
Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing ...Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.展开更多
Based on linear interval equations, an accurate interval finite element method for solving structural static problems with uncertain parameters in terms of optimization is discussed. On the premise of ensuring the con...Based on linear interval equations, an accurate interval finite element method for solving structural static problems with uncertain parameters in terms of optimization is discussed. On the premise of ensuring the consistency of solution sets, the original interval equations are equivalently transformed into some deterministic inequations. On this basis, calculating the structural displacement response with interval parameters is predigested to a number of deterministic linear optimization problems. The results are proved to be accurate to the interval governing equations. Finally, a numerical example is given to demonstrate the feasibility and efficiency of the proposed method.展开更多
In this article, we will derive local elliptic type gradient estimates for positive solutions of linear parabolic equations （△-e/et）u（x,t）＋q（x,t）u^p（x,t）=0 and nonlinear parabolic equations （△-e/et）u（x,...In this article, we will derive local elliptic type gradient estimates for positive solutions of linear parabolic equations （△-e/et）u（x,t）＋q（x,t）u^p（x,t）=0 and nonlinear parabolic equations （△-e/et）u（x,t）＋h（x,t）u^p（x,t）=0（p 〉 1） on Riemannian manifolds.As applications, we obtain some theorems of Liouville type for positive ancient solutions of such equations. Our results generalize that of Souplet-Zhang （[1], Bull. London Math. Soc. 38（2006）, 1045-1053） and the author （[2], Nonlinear Anal. 74 （2011）, 5141-5146）.展开更多
Finite time blow up of the solutions to Boussinesq equation with linear restoring force and combined power nonlinearities is studied. Sufficient conditions on the initial data for nonexistence of global solutions are ...Finite time blow up of the solutions to Boussinesq equation with linear restoring force and combined power nonlinearities is studied. Sufficient conditions on the initial data for nonexistence of global solutions are derived. The results are valid for initial data with arbitrary high positive energy. The proofs are based on the concave method and new sign preserving functionals.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
文摘Let a_(1),a_(2),a_(3)be nonzero integers with gcd(a_(1),a_(2),a_(3))=1,and let k be any positive integer,K=max[3,|a_(1)|,|a_(2)|,|a_(3)|,k].Suppose that l_(1),l_(2),l_(3)are integers each coprime to k.Suppose further that b is any integer satisfying some necessary congruent conditions.The solvability of linear equation a_(1)p_(1)+a_(2)p_(2)+a_(3)p_(3)=b(p_(j)=l_(j)(mod k),1≤j≤3)with prime variables pi,p_(2),ps is investigated.It is proved that if ai,a_(2),a_(3)are all positive,then the above equation is solvable whenever b≥K^(25);if a,a_(2),a_(3)are not all of the same sign,then the above equation has a solution p_(1),p_(2),p_(3)satisfying max(p_(1),p_(2),p_(3))≤3|b|+K^(25).
文摘Numerical treatment of engineering application problems often eventually results in a solution of systems of linear or nonlinear equations.The solution process using digital computational devices usually takes tremendous time due to the extremely large size encountered in most real-world engineering applications.So,practical solvers for systems of linear and nonlinear equations based on multi graphic process units（GPUs）are proposed in order to accelerate the solving process.In the linear and nonlinear solvers,the preconditioned bi-conjugate gradient stable（PBi-CGstab）method and the Inexact Newton method are used to achieve the fast and stable convergence behavior.Multi-GPUs are utilized to obtain more data storage that large size problems need.
文摘The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods. Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.
基金supported by the National Natural Science Foundation of China (11171119 and 10871076)
文摘In this article, we study the complex oscillation problems of entire solutions to homogeneous and nonhomogeneous linear difference equations, and obtain some relations of the exponent of convergence of zeros and the order of growth of entire solutions to complex linear difference equations.
基金supported by Shanghai Center for Mathematical Sci-ences,China Scholarship Council(201206105015)National Science Foundation of China(11171119,11001057,11571049)Natural Science Foundation of Guangdong Province in China(2014A030313422)
文摘In this paper,we mainly investigate the dynamical properties of entire solutions of complex differential equations.With some conditions on coefficients,we prove that the set of common limiting directions of Julia sets of solutions,their derivatives and their primitives must have a definite range of measure.
基金supported by Shanghai Center for Mathematical Science China Scholarship Council(201206105015)the National Science Foundation of China(11171119,11001057,11571049)the Natural Science Foundation of Guangdong Province in China(2014A030313422)
文摘In this paper, we mainly investigate entire solutions of complex differential equations with coefficients involving exponential functions, and obtain the dynamical properties of the solutions, their derivatives and primitives. With some conditions on coefficients, for the solutions, their derivatives and their primitives, we consider the common limiting directions of Julia set and the existence of Baker wandering domain.
基金the National Natural Science Foundation of China(10161006,10571044)the Natural Science Foundation of Guangdong Prov(06025059)
文摘In this article, the authors study the growth of certain second order linear differential equation f″＋A（z）f′＋B（z）f=0 and give precise estimates for the hyperorder of solutions of infinite order. Under similar conditions, higher order differential equations will be considered.
基金Supported by the National Natural Science Foundation of China(61272300)
文摘Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.
基金Project supported by the National Science Foundation of China (10271097)
文摘For domains composed by balls in C^n, this paper studies the boundary behaviour of Cauchy type integrals with discrete holomorphic kernels and the corresponding linear singular integral equation on each piece of smooth lower dimensional edges on the boundary of the domain.
文摘In this paper, the asynchronous versions of classical iterative methods for solving linear systems of equations are considered. Sufficient conditions for convergence of asynchronous relaxed processes are given for H-matrix by which nor only the requirements of [3] on coefficient matrix are lowered, but also a larger region of convergence than that in [3] is obtained.
文摘Lyapunov equation is one of the most basic and important equations in control theory, which has various applications in, e.g., stability analysis and robust analysis of linear control systems. Inspired by the recent progresses of quantum algorithms, we find that solving Lyapunov equation can be exponentially accelerated by quantum algorithms rather than traditional classical algorithms. Our algorithm is more efficient especially when the system matrix is sparse and has a low condition number. The results presented in this paper open up new dimensions of research in controlling classical system by quantum information processors, which has rarely been considered in the existing literature.
文摘In this paper, we investigate the growth of solutions of higher order linear differential equations with meromorphic coefficients. Under certain conditions, we obtain precise estimation of growth order and hyper-order of solutions of the equation.
基金Project supported by the National Natural Science Foundation of China (Nos. 10232040, 10572002 and 10572003)
文摘A new direct method for solving unsymmetrical sparse linear systems（USLS） arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin （MLPG） method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.
基金supported by the National Natural Science Foundation of China(Grant No.12031004 and Grant No.12271474,61877054)the Fundamental Research Foundation for the Central Universities(Project No.K20210337)+1 种基金the Zhejiang University Global Partnership Fund,188170+194452119/003partially funded by a state task of Russian Fundamental Investigations(State Registration No.FFSG-2024-0002)。
文摘Fundamental matrix operations and solving linear systems of equations are ubiquitous in scientific investigations.Using the‘sender-receiver’model,we propose quantum algorithms for matrix operations such as matrix-vector product,matrix-matrix product,the sum of two matrices,and the calculation of determinant and inverse matrix.We encode the matrix entries into the probability amplitudes of the pure initial states of senders.After applying proper unitary transformation to the complete quantum system,the desired result can be found in certain blocks of the receiver’s density matrix.These quantum protocols can be used as subroutines in other quantum schemes.Furthermore,we present an alternative quantum algorithm for solving linear systems of equations.
基金Project supported by the National Natural Science Foundation of China(Nos.5130926141030747+3 种基金41102181and 51121005)the National Basic Research Program of China(973 Program)(No.2011CB013503)the Young Teachers’ Initial Funding Scheme of Sun Yat-sen University(No.39000-1188140)
文摘Fast solving large-scale linear equations in the finite element analysis is a classical subject in computational mechanics. It is a key technique in computer aided engineering （CAE） and computer aided manufacturing （CAM）. This paper presents a high-efficiency improved symmetric successive over-relaxation （ISSOR） preconditioned conjugate gradient （PCG） method, which maintains lelism consistent with the original form. Ideally, the by 50% as compared with the original algorithm. the convergence and inherent paralcomputation can It is suitable for be reduced nearly high-performance computing with its inherent basic high-efficiency operations. By comparing with the numerical results, it is shown that the proposed method has the best performance.
基金supported by the National Natural Science Foundation of China(Nos.90816024,10872017,and 10876100)the Defense Industrial Technology Development Program(Nos.A2120110001 and 2120110011)the 111 Project(No.B07009)
文摘Based on linear interval equations, an accurate interval finite element method for solving structural static problems with uncertain parameters in terms of optimization is discussed. On the premise of ensuring the consistency of solution sets, the original interval equations are equivalently transformed into some deterministic inequations. On this basis, calculating the structural displacement response with interval parameters is predigested to a number of deterministic linear optimization problems. The results are proved to be accurate to the interval governing equations. Finally, a numerical example is given to demonstrate the feasibility and efficiency of the proposed method.
基金supported by the National Science Foundation of China(41275063 and 11401575)
文摘In this article, we will derive local elliptic type gradient estimates for positive solutions of linear parabolic equations （△-e/et）u（x,t）＋q（x,t）u^p（x,t）=0 and nonlinear parabolic equations （△-e/et）u（x,t）＋h（x,t）u^p（x,t）=0（p 〉 1） on Riemannian manifolds.As applications, we obtain some theorems of Liouville type for positive ancient solutions of such equations. Our results generalize that of Souplet-Zhang （[1], Bull. London Math. Soc. 38（2006）, 1045-1053） and the author （[2], Nonlinear Anal. 74 （2011）, 5141-5146）.
基金partially supported by Grant No.DFNI I-02/9 of the Bulgarian Science Fund
文摘Finite time blow up of the solutions to Boussinesq equation with linear restoring force and combined power nonlinearities is studied. Sufficient conditions on the initial data for nonexistence of global solutions are derived. The results are valid for initial data with arbitrary high positive energy. The proofs are based on the concave method and new sign preserving functionals.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
