{"title":"Edge importance in complex networks","authors":"Silvia Noschese, Lothar Reichel","doi":"10.1007/s11075-024-01881-1","DOIUrl":"https://doi.org/10.1007/s11075-024-01881-1","url":null,"abstract":"<p>Complex networks are made up of vertices and edges. The latter connect the vertices. There are several ways to measure the importance of the vertices, e.g., by counting the number of edges that start or end at each vertex, or by using the subgraph centrality of the vertices. It is more difficult to assess the importance of the edges. One approach is to consider the line graph associated with the given network and determine the importance of the vertices of the line graph, but this is fairly complicated except for small networks. This paper compares two approaches to estimate the importance of edges of medium-sized to large networks. One approach computes partial derivatives of the total communicability of the weights of the edges, where a partial derivative of large magnitude indicates that the corresponding edge may be important. Our second approach computes the Perron sensitivity of the edges. A high sensitivity signals that the edge may be important. The performance of these methods and some computational aspects are discussed. Applications of interest include to determine whether a network can be replaced by a network with fewer edges with about the same communicability.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"31 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935625","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Erik Burman, Rebecca Durst, Miguel A. Fernández, Johnny Guzmán, Sijing Liu
{"title":"Estimates of discrete time derivatives for the parabolic-parabolic Robin-Robin coupling method","authors":"Erik Burman, Rebecca Durst, Miguel A. Fernández, Johnny Guzmán, Sijing Liu","doi":"10.1007/s11075-024-01902-z","DOIUrl":"https://doi.org/10.1007/s11075-024-01902-z","url":null,"abstract":"<p>We consider a loosely coupled, non-iterative Robin-Robin coupling method proposed and analyzed in Burman et al. (J. Numer. Math. <b>31</b>(1):59–77, 2023) for a parabolic-parabolic interface problem and prove estimates for the discrete time derivatives of the scalar field in different norms. When the interface is flat and perpendicular to two of the edges of the domain we prove error estimates in the <span>(H^2)</span>-norm. Such estimates are key ingredients to analyze a defect correction method for the parabolic-parabolic interface problem. Numerical results are shown to support our findings.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"22 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935624","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Discrete non-commutative hungry Toda lattice and its application in matrix computation","authors":"Zheng Wang, Shi-Hao Li, Kang-Ya Lu, Jian-Qing Sun","doi":"10.1007/s11075-024-01915-8","DOIUrl":"https://doi.org/10.1007/s11075-024-01915-8","url":null,"abstract":"<p>In this paper, we plan to show an eigenvalue algorithm for block Hessenberg matrices by using the idea of non-commutative integrable systems and matrix-valued orthogonal polynomials. We introduce adjacent families of matrix-valued <span>(theta )</span>-deformed bi-orthogonal polynomials, and derive corresponding discrete non-commutative hungry Toda lattice from discrete spectral transformations for polynomials. It is shown that this discrete system can be used as a pre-precessing algorithm for block Hessenberg matrices. Besides, some convergence analysis and numerical examples of this algorithm are presented.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"41 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935553","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Approximate moment functions for logistic stochastic differentialequations","authors":"Coşkun Çetin, Jasmina Đorđević","doi":"10.1007/s11075-024-01911-y","DOIUrl":"https://doi.org/10.1007/s11075-024-01911-y","url":null,"abstract":"<p>In this paper, we introduce a method of successive approximations for moment functions of logistic stochastic differential equations. We first reduce the system of the corresponding moment functions to an infinite system of linear ordinary differential equations. Then, we determine certain upper and lower bounds on the moment functions, and utilize these bounds to solve the resulting systems approximately via suitable truncations, iterations and a local improvement step. After obtaining some general theoretical results on the error norms and describing a general algorithm for logistic SDE, we focus on stochastic Verhulst systems in numerical implementations. We compare their moment approximations with numerical solutions via simulation-based methods that include discretizations of the pathwise solutions as well as other convergent numerical procedures like semi-implicit split-step Euler methods.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"59 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A novel higher-order efficient computational method for pricing European and Asian options","authors":"Saurabh Bansal, Srinivasan Natesan","doi":"10.1007/s11075-024-01909-6","DOIUrl":"https://doi.org/10.1007/s11075-024-01909-6","url":null,"abstract":"<p>In this article, we present a fourth-order accurate numerical method for solving generalized Black-Scholes PDE describing European and Asian options. Initially, we discretize the time derivative by the Crank-Nicolson scheme, and then the resultant semi-discrete problem by the central difference scheme on uniform meshes. In order to enhance the order of convergence of the proposed scheme, we employ the Richardson extrapolation method, by using two different meshes to solve the fully discrete problem. The stability and convergence are studied. To validate the proposed technique, several numerical experiments are carried out.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"2 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935555","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Convergence analysis of the augmented Lagrangian method for $$ell _{p}$$ -norm cone optimization problems with $$p ge 2$$","authors":"Benqi Liu, Kai Gong, Liwei Zhang","doi":"10.1007/s11075-024-01912-x","DOIUrl":"https://doi.org/10.1007/s11075-024-01912-x","url":null,"abstract":"<p>This paper focuses on the convergence analysis of the augmented Lagrangian method (ALM) for <span>(varvec{ell }_{varvec{p}})</span>-norm cone optimization problems. We investigate some properties of the augmented Lagrangian function and <span>(varvec{ell }_{varvec{p}})</span>-norm cone. Moreover, under the Jacobian uniqueness conditions, we prove that the local convergence rate of ALM for solving <span>(varvec{ell }_{varvec{p}})</span>-norm cone optimization problems with <span>(varvec{p} varvec{ge } varvec{2})</span> is proportional to <span>(varvec{1}varvec{/}varvec{r})</span>, where the penalty parameter <span>(varvec{r})</span> is not less than a threshold <span>(varvec{hat{r}})</span>. In numerical simulations, we successfully validate the effectiveness and convergence properties of ALM.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"6 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935556","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
James Baglama, Jonathan A. Chávez-Casillas, Vasilije Perović
{"title":"A hybrid algorithm for computing a partial singular value decomposition satisfying a given threshold","authors":"James Baglama, Jonathan A. Chávez-Casillas, Vasilije Perović","doi":"10.1007/s11075-024-01906-9","DOIUrl":"https://doi.org/10.1007/s11075-024-01906-9","url":null,"abstract":"<p>In this paper, we describe a new hybrid algorithm for computing all singular triplets above a given threshold and provide its implementation in MATLAB/Octave and <span>R</span>. The high performance of our codes and ease at which they can be used, either independently or within a larger numerical scheme, are illustrated through several numerical examples with applications to matrix completion and image compression. Well-documented MATLAB and <span>R</span> codes are provided for public use.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"16 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935622","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Stochastic theta methods for random periodic solution of stochastic differential equations under non-globally Lipschitz conditions","authors":"Ziheng Chen, Liangmin Cao, Lin Chen","doi":"10.1007/s11075-024-01892-y","DOIUrl":"https://doi.org/10.1007/s11075-024-01892-y","url":null,"abstract":"<p>This work focuses on the numerical approximations of random periodic solutions of stochastic differential equations (SDEs). Under non-globally Lipschitz conditions, we prove the existence and uniqueness of random periodic solutions for the considered equations and its numerical approximations generated by the stochastic theta (ST) methods with <span>(theta in (1/2,1])</span>. It is shown that the random periodic solution of each ST method converges strongly in the mean square sense to that of SDEs. More precisely, the mean square convergence order is 1/2 for SDEs with multiplicative noise and 1 for SDEs with additive noise, respectively. Numerical results are finally reported to confirm these theoretical findings.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"2 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141935557","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Optimization algorithms for stabilization of multi-input vibration system with time delay using eigenvalues assignment technique","authors":"Peizhao Yu, Fuheng Zhao, Haoming Xin","doi":"10.1007/s11075-024-01899-5","DOIUrl":"https://doi.org/10.1007/s11075-024-01899-5","url":null,"abstract":"<p>The study considers the robust and minimum norm problems for stabilization using partial eigenvalue assignment technique in nonsingular vibration system with time delay via the acceleration-velocity-displacement active controller. The new gains expressions of active controller are derived by orthogonality relations, which keeps the no spill-over property of the vibration system. To discuss the stabilization problem using eigenvalues assignment technique, the linear equation is solved by constructing a special matrix which is proved to be nonsingular. Solving algorithm is proposed to obtain the parametric expressions of active controller. A new gradient-based optimization method is proposed to discuss the robust and minimum norm controller design by establishing the gradient formulas of cost functions. The optimization algorithm is proposed to discuss the robust and minimum norm stabilization of closed-loop eigenvalues in vibration system with time delay. The presented algorithms are feasible to the case of time delay between measurements of state and actuation of control. Numerical examples show the effectiveness of the method.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"97 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141881205","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Structure-preserving joint Lanczos bidiagonalization with thick-restart for the partial quaternion GSVD","authors":"Zhe-Han Hu, Si-Tao Ling, Zhi-Gang Jia","doi":"10.1007/s11075-024-01900-1","DOIUrl":"https://doi.org/10.1007/s11075-024-01900-1","url":null,"abstract":"<p>A new Krylov subspace method is designed in the computation of partial quaternion generalized singular value decomposition (QGSVD) of a large-scale quaternion matrix pair <span>({textbf{A}, textbf{B}})</span>. Explicitly, we present the structure-preserving joint Lanczos bidiagonalization method to reduce <span>(textbf{A})</span> and <span>(textbf{B})</span> to lower and upper real bidiagonal matrices, respectively. We carry out the thick-restarted technique with the combination of a robust selective reorthogonalization strategy in the structure-preserving joint Lanczos bidiagonalization process. In the iteration process we avoid performing the explicit QR decomposition of the quaternion matrix pair. Numerical experiments illustrate the effectiveness of the proposed method.</p>","PeriodicalId":54709,"journal":{"name":"Numerical Algorithms","volume":"77 1","pages":""},"PeriodicalIF":2.1,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141870103","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}