C++ Codes of Implicit Lu Algorithms for Absdll01

Xing Li,Ying Liu,Antonino Del Popolo

DOI: https://doi.org/10.48550/arXiv.math/0105167

2001-05-21

Abstract:This report is devoted to some C++ codes implementing the implicit LU class algorithms for solving linear determined, and undetermined systems with $n$ variables and $m$ equations. A main program used in part of the numerical test is given in the last section.

Numerical Analysis,Astrophysics

What problem does this paper attempt to address?

Nonlinear Structural Analysis Algorithm Based on the Null Space Basis

陈务军,董石麟,付功义,关富玲

DOI: https://doi.org/10.3969/j.issn.1007-4708.2001.04.006

2001-01-01

Abstract:Based on the basis of null space and M-P inverse theorem, a unified incremental solution formulae, expressed as the null space basis of the augmented matrix, were firstly formulated herein for the incremental equilibrium equation at the limit point and non-critical points. Three numerical analysis techniques, based on the M-P inverse, singular value decomposition (SVD) and pseudo up-triangular decomposition (PUTD), respectively, were developed to compute the null space basis for the augmented equilibrium matrix, as well as how to obtain the numerical incremental solution. Numerical implementations are presented. A specialized nonlinear analysis program were developed with Legend PⅡ300 in Matlab integrated development environment (IDE) using Matlab language. Finally, two classical numerical examples are used to verify the accuracy of the method aforementioned. Efficiency and accuracy, i.e. computation time and buckling load differences, have been compared in detail for various numerical approaches. The algorithm and program were proved to be efficient. This method can be used for medium scale engineering project.
Abs Algorithms for Linear Equations and Abspack

E. Spedicato,Z. Xia,E. Bodon,A. Del Popolo

DOI: https://doi.org/10.48550/arXiv.math/0105074

2001-05-09

Numerical Analysis

Abstract:We present the main results obtained during a research on ABS methods in the framework of the project Analisi Numerica e Matematica Computazionale.
Operational Flexibility Analysis of High-Dimensional Systems via Cylindrical Algebraic Decomposition

Chenglin Zheng,Fei Zhao,Lingyu Zhu,Xi Chen

DOI: https://doi.org/10.1021/acs.iecr.9b06061

2020-01-01

Industrial & Engineering Chemistry Research

Abstract:The cylindrical algebraic decomposition (CAD) method has been proposed for flexibility analysis to derive analytical expressions of a feasible region. Due to the heavy computational burden caused by symbolic computation, this method can only handle small-scale problems currently. To overcome this limitation, a novel method is proposed for high-dimensional systems with a number of equalities and limited inequalities. A surrogate model is first built to correlate the inequality constraints based on an initial sample set. Then, the flexibility region is obtained with explicit expressions via the CAD method. Next, for any violation, a refinement will be activated by taking an iterative process of boundary check, surrogate modeling, region deriving, and underestimation check, until the termination condition is satisfied. The case studies show the proposed method can effectively describe the flexibility region for both the convex and nonconvex systems.
A New Algorithm for Linear Systems of the Pascal Type

Xiao-Guang Lv,Ting-Zhu Huang,Zhi-Gang Ren

DOI: https://doi.org/10.1016/j.cam.2008.07.045

IF: 2.872

2009-01-01

Journal of Computational and Applied Mathematics

Abstract:In this paper, we give an algorithm for solving linear systems of the Pascal matrices. The method is based on the explicit factorization of the Pascal matrices. The algorithm costs no multiplications and O(n2) additions. The linear systems of the generalized Pascal matrices are also considered. Some examples are given.
A Complete Algorithm for Automated Discovering of a Class of Inequality-Type Theorems

Yang Lu,Hou Xiaorong,Xia Bican

DOI: https://doi.org/10.1007/bf02713938

2001-01-01

Abstract:>Making use of the discriminant sequence for polynomials, WR algorithm, Wu’ s elimination and a partial cylindrical algebraic decomposition, we present here a practical algorithm for automated inequality discovering which can discover new inequalities automatically without requiring to put forward any conjectures beforehand. That is complete for an extensive class of inequality-type theorems. Also this algorithm is applied to the classification of the real physical solutions of geometric constraint problems. Many inequalities with various backgrounds have been discovered or rediscovered by our program, DISCOVERER, which implements the algorithm in Maple.
A New Class of Methods for Solving Nonlinear Equations

B HAN,DL ZHANG,JQ LIU

DOI: https://doi.org/10.1016/0377-0427(94)00009-3

IF: 2.872

1994-01-01

Journal of Computational and Applied Mathematics

Abstract:A new class of implicit methods for solving nonlinear equations is proposed in this paper. Comparisons of computational results are made with other well-known methods on a number of difficult problems. The implicit methods are accurate, robust and economical.
SDPNAL+: A Matlab software for semidefinite programming with bound constraints (version 1.0)

Defeng Sun,Kim-Chuan Toh,Yancheng Yuan,Xin-Yuan Zhao

DOI: https://doi.org/10.1080/10556788.2019.1576176

2019-02-22

Optimization Methods and Software

Abstract:Sdpnal+ is a MATLAB software package that implements an augmented Lagrangian based method to solve large scale semidefinite programming problems with bound constraints. The implementation was initially based on a majorized semismooth Newton-CG augmented Lagrangian method, here we designed it within an inexact symmetric Gauss-Seidel based semi-proximal ADMM/ALM (alternating direction method of multipliers/augmented Lagrangian method) framework for the purpose of deriving simpler stopping conditions and closing the gap between the practical implementation of the algorithm and the theoretical algorithm. The basic code is written in MATLAB, but some subroutines in C language are incorporated via Mex files. We also design a convenient interface for users to input their SDP models into the solver. Numerous problems arising from combinatorial optimization and binary integer quadratic programming problems have been tested to evaluate the performance of the solver. Extensive numerical experiments conducted in [L.Q. Yang, D.F. Sun, and K.C. Toh, SDPNAL+: A majorized semismooth Newton-CG augmented Lagrangian method for semidefinite programming with nonnegative constraints, Math. Program. Comput. 7 (2015), pp. 331–366] show that the proposed method is quite efficient and robust, in that it is able to solve 98.9% of the 745 test instances of SDP problems arising from various applications to the accuracy of in the relative KKT residual.

operations research & management science,mathematics, applied,computer science, software engineering
Parallelization of LAPACK Linear Equation Functions Based on Loongson 3A

Bin-bin LIU,Nai-jie GU,Kai-xin REN,Jie ZHANG

DOI: https://doi.org/10.3969/j.issn.1000-1220.2014.05.028

2014-01-01

Abstract:According to the characteristics of Loongson 3A architecture,this paper achieves the parallelized version of all linear equation functions based on the level of the original library structure through parallelization of BLAS functions and LAPACK functions,by methods such as substitution or elimination of variables,simplifing dependence,increasing copies of hot data and parallel flow method. Users can use it by simply setting the number of cores and recompiling the multi-core LAPACK library. By running the LAPACK Timing Programs,experimental results are obtained and it shows that most functions have a speedup of 2 or more,the arithmetic mean of speedups of all linear equation functions reached 2. 04.
ABS Algorithms for Linear Systems and Optimization

Emilio Spedicato,Elena Bodon,Antonino Del Popolo,Zunquan Xia

DOI: https://doi.org/10.48550/arXiv.math/0105056

2001-05-10

Abstract:We present a review and bibliography of the main results obtained during a research on ABS (Abaffy, Broyden, Spedicato) methods.

Numerical Analysis,Astrophysics
Detection of soft errors in LU decomposition with partial pivoting using algorithm-based fault tolerance

Erlin Yao,Jiutian Zhang,Mingyu Chen,Guangming Tan,Ninghui Sun

DOI: https://doi.org/10.1177/1094342015578487

2015-01-01

Abstract:Soft errors in scientific computing applications are becoming inevitable with the ever-increasing system scale and execution time, and new technologies that feature increased transistor density and lower voltage. Soft errors can be mainly classified into two categories: bit-flipping error e.g. 1 becomes −1 in random access memory; and computation error e.g. 1+1 = 3 in floating point units. Traditionally, bit-flipping error is handled by the Error Correcting Code ECC technique, and computation error is dealt with the Triple Modular Redundancy TMR method. Note that, ECC cannot handle computation error, while TMR cannot deal with bit-flipping error and is not efficient on handling computation error. To uniformly and efficiently handle both computation and bit-flipping errors in matrix operations, the Algorithm-Based Fault Tolerance ABFT method is developed. This paper focuses on the detection of soft errors in the LU Decomposition with Partial Pivoting LUPP algorithm, which is widely used in scientific computing applications. First, this paper notes that existing ABFT methods are not adequate to detect soft errors in LUPP in terms of time or space. Then we propose a new ABFT algorithm which can detect soft errors in LUPP both flexible in time and comprehensive in space. Flexible in time means that soft errors can be detected flexibly during the execution instead of only at the end of LUPP, while comprehensive in space indicates that all of the elements in data matrices L and U will be covered for detecting soft errors. To show the feasibility and efficiency of the proposed algorithm, this paper has incorporated it into the implementation of LUPP in the widely used benchmark High Performance Linpack HPL. Experiment results verify the feasibility of this algorithm: for soft errors injected at various timings and to different elements in LUPP, this algorithm has detected most of the injected errors, which have covered all of the errors that cannot pass the residual check of HPL. Both theoretical overhead analysis and experiments demonstrate that this ABFT algorithm is also very efficient at detecting soft errors in LUPP.
New iterative methods for solving linear systems

joshua du,baodong zheng,liancheng wang

DOI: https://doi.org/10.11948/2011024

2011-01-01

Abstract:In this paper, we introduce some new iterative methods to solve linear systems Ax = b. We show that these methods, comparing to the classical Jacobi or Gauss-Seidel method, can be applied to more systems and have faster convergence.
On the Performance of a Novel Class of Linear System Solvers and Comparison with State-of-The-Art Algorithms

Chun Lau,Bahman Kalantari

2023-04-24

Abstract:We present a comprehensive computational study of a class of linear system solvers, called {\it Triangle Algorithm} (TA) and {\it Centering Triangle Algorithm} (CTA), developed by Kalantari \cite{kalantari23}. The algorithms compute an approximate solution or minimum-norm solution to $Ax=b$ or $A^TAx=A^Tb$, where $A$ is an $m \times n$ real matrix of arbitrary rank. The algorithms specialize when $A$ is symmetric positive semi-definite. Based on the description and theoretical properties of TA and CTA from \cite{kalantari23}, we give an implementation of the algorithms that is easy-to-use for practitioners, versatile for a wide range of problems, and robust in that our implementation does not necessitate any constraints on $A$. Next, we make computational comparisons of our implementation with the Matlab implementations of two state-of-the-art algorithms, GMRES and ``lsqminnorm". We consider square and rectangular matrices, for $m$ up to $10000$ and $n$ up to $1000000$, encompassing a variety of applications. These results indicate that our implementation outperforms GMRES and ``lsqminnorm" both in runtime and quality of residuals. Moreover, the relative residuals of CTA decrease considerably faster and more consistently than GMRES, and our implementation provides high precision approximation, faster than GMRES reports lack of convergence. With respect to ``lsqminnorm", our implementation runs faster, producing better solutions. Additionally, we present a theoretical study in the dynamics of iterations of residuals in CTA and complement it with revealing visualizations. Lastly, we extend TA for LP feasibility problems, handling non-negativity constraints. Computational results show that our implementation for this extension is on par with those of TA and CTA, suggesting applicability in linear programming and related problems.

Numerical Analysis
Alternately linearized implicit iteration methods for the minimal nonnegative solutions of the nonsymmetric algebraic Riccati equations

Zhong-Zhi Bai,Xiao-Xia Guo,Shu-Fang Xu

DOI: https://doi.org/10.1002/nla.500

2006-01-01

Numerical Linear Algebra with Applications

Abstract:For the non-symmetric algebraic Riccati equations, we establish a class of alternately linearized implicit (ALI) iteration methods for computing its minimal non-negative solutions by technical combination of alternate splitting and successive approximating of the algebraic Riccati operators. These methods include one iteration parameter, and suitable choices of this parameter may result in fast convergent iteration methods. Under suitable conditions, we prove the monotone convergence and estimate the asymptotic convergence factor of the ALI iteration matrix sequences. Numerical experiments show that the ALI iteration methods are feasible and effective, and can outperform the Newton iteration method and the fixed-point iteration methods. Besides, we further generalize the known fixed-point iterations, obtaining an extensive class of relaxed splitting iteration methods for solving the non-symmetric algebraic Riccati equations. Copyright (C) 2006 John Wiley & Sons, Ltd.
Modified Picard-like Method for Solving Absolute Value Equations

Yuan Liang,Chaoqian Li

DOI: https://doi.org/10.3390/math11040848

IF: 2.4

2023-02-07

Mathematics

Abstract:We present a modified Picard-like method to solve absolute value equations by equivalently expressing the implicit fixed-point equation form of the absolute value equations as a two-by-two block nonlinear equation. This unifies some existing matrix splitting algorithms and improves the efficiency of the algorithm by introducing the parameter ω. For the choice of ω in the new method, we give a way to determine the quasi-optimal values. Numerical examples are given to show the feasibility of the proposed method. It is also shown that the new method is better than those proposed by Ke and Ma in 2017 and Dehghan and Shirilord in 2020.

mathematics
Subproper and Regular Splittings for Restricted Rectangular Linear System

YM Wei,XZ Li,HB Wu

DOI: https://doi.org/10.1016/s0096-3003(02)00078-4

IF: 4.397

2003-01-01

Applied Mathematics and Computation

Abstract:Iterative schemes for approximating a solution to restricted rectangular but consistent linear system of equations Ax=b, x∈T, are considered. The methods are based upon so-called subproper splitting A=M−N, which is a generalization of the concept of subproper splitting introduced by Neumann and studied further by Berman and Neumann and others. We give a necessary and sufficient condition on the splitting such that the iterative sequence converges to a solution of Ax=b in the case of b∈AT for every x0, where A∈Cm×n and T is a subspace of Cn. Monotonicity and the concept of regular subproper splitting are used to determine a necessary and a sufficient condition for the convergence of the iterative scheme. Finally, we present two numerical examples to verify our conclusions.
New High Performance Sparse Static Solver in Finite Element Analysis with Loop-Unrolling

P Chen,SL Sun

DOI: https://doi.org/10.1007/s10338-005-0531-2

2005-01-01

Acta Mechanica Solida Sinica

Abstract:In the previous papers, a high performance sparse static solver with two-level unrolling based on a cell-sparse storage scheme was reported. Although the solver reaches quite a high efficiency for a big percentage of finite element analysis benchmark tests, the MFLOPS (million floating operations per second) of LDLT factorization of benchmark tests vary on a Dell Pentium IV 850 MHz machine from 100 to 456 depending on the average size of the super-equations, i.e., on the average depth of unrolling. In this paper, a new sparse static solver with two-level unrolling that employs the concept of master-equations and searches for an appropriate depths of unrolling is proposed. The new solver provides higher MFLOPS for LDLT factorization of benchmark tests, and therefore speeds up the solution process.
A general alternating direction implicit iteration method for solving complex symmetric linear systems

Juan Zhang,Wenlu Xun

2024-04-18

Abstract:We have introduced the generalized alternating direction implicit iteration (GADI) method for solving large sparse complex symmetric linear systems and proved its convergence properties. Additionally, some numerical results have demonstrated the effectiveness of this algorithm. Furthermore, as an application of the GADI method in solving complex symmetric linear systems, we utilized the flattening operator and Kronecker product properties to solve Lyapunov and Riccati equations with complex coefficients using the GADI method. In solving the Riccati equation, we combined inner and outer iterations, first simplifying the Riccati equation into a Lyapunov equation using the Newton method, and then applying the GADI method for solution. Finally, we provided convergence analysis of the method and corresponding numerical results.

Numerical Analysis
A Class of Smoothing Modulus-Based Iterative Method for Solving Implicit Complementarity Problems

Cong Guo,Chenliang Li,Tao Luo

DOI: https://doi.org/10.4236/ajcm.2022.122011

2023-06-08

Abstract:In this paper, a class of smoothing modulus-based iterative method was presented for solving implicit complementarity problems. The main idea was to transform the implicit complementarity problem into an equivalent implicit fixed-point equation, then introduces a smoothing function to obtain its approximation solutions. The convergence analysis of the algorithm was given, and the efficiency of the algorithms was verified by numerical experiment

Numerical Analysis
Comparison of Subspace Iteration, Iterative Ritz Vector Method and Iterative Lanczos Method

GONG Yu-cai,ZHOU Hong-Wei,CHEN Pu,YUAN Ming-wu

DOI: https://doi.org/10.3969/j.issn.1004-4523.2005.02.019

2005-01-01

Abstract:Based on the cell sparse fast solver and loop-unrolling, this paper implements three efficient eigenvalue algorithms-subspace iteration, iterative Ritz method and iterative Lanczos method. Slight modifications are made for iterative Ritz method and iterative Lanczos method. These eigenvalue algorithms are examined under the mode error, i.e, the ratio of out-of-balance nodal point forces that is the difference of maximum elastic nodal point forces and maximum inertia nodal point forces, and the maximum elastic nodal point. Averagely, iterative Ritz method is the most effcient one among them. Engineering projects are used as examples to verify the methods. Under the mode error convergence criterion the eigenvalue extracting processes are more stable than eigenvalue convergence criteria. Compared with ANSYS's subspace iteration and block Lanczos approaches, the subspace iteration of this paper is much more efficient, and Lanczos approach has almost equal efficiency. The methods proposed are of industrial strength and efficient. Large scale tests show that the improvement in terms of CPU time and storage request is tremendous.
ITERATIVE ALGORITHMS FOR THE LEAST-SQUARES SYMMETRIC SOLUTION OF AXB = C WITH A SUBMATRIX CONSTRAINT

Minghui Wang,Yan Feng

2009-01-01

Abstract:Iterative algorithms are proposed for the least-squares symmetric solution of AXB = E with a submatrix constraint. We characterize the linear mappings from their independent element space to the constrained solution sets, study their properties and use these properties to propose two matrix iterative algorithms that can find the minimum and quasi-minimum norm solution based on the classical LSQR algorithm for solving the unconstrained LS problem. Numerical results are provided that show the efficiency of the proposed methods.

C++ Codes of Implicit Lu Algorithms for Absdll01

Nonlinear Structural Analysis Algorithm Based on the Null Space Basis

Abs Algorithms for Linear Equations and Abspack

Operational Flexibility Analysis of High-Dimensional Systems via Cylindrical Algebraic Decomposition

A New Algorithm for Linear Systems of the Pascal Type

A Complete Algorithm for Automated Discovering of a Class of Inequality-Type Theorems

A New Class of Methods for Solving Nonlinear Equations

SDPNAL+: A Matlab software for semidefinite programming with bound constraints (version 1.0)

Parallelization of LAPACK Linear Equation Functions Based on Loongson 3A

ABS Algorithms for Linear Systems and Optimization

Detection of soft errors in LU decomposition with partial pivoting using algorithm-based fault tolerance

New iterative methods for solving linear systems

On the Performance of a Novel Class of Linear System Solvers and Comparison with State-of-The-Art Algorithms

Alternately linearized implicit iteration methods for the minimal nonnegative solutions of the nonsymmetric algebraic Riccati equations

Modified Picard-like Method for Solving Absolute Value Equations

Subproper and Regular Splittings for Restricted Rectangular Linear System

New High Performance Sparse Static Solver in Finite Element Analysis with Loop-Unrolling

A general alternating direction implicit iteration method for solving complex symmetric linear systems

A Class of Smoothing Modulus-Based Iterative Method for Solving Implicit Complementarity Problems

Comparison of Subspace Iteration, Iterative Ritz Vector Method and Iterative Lanczos Method

ITERATIVE ALGORITHMS FOR THE LEAST-SQUARES SYMMETRIC SOLUTION OF AXB = C WITH A SUBMATRIX CONSTRAINT