A note on computing the generalized inverse of a matrix .
Necessity of computing large sparse Hessian matrices gave birth to many methods for their effective approximation by differences of gradients. We adopt the so-called direct methods for this problem that we faced when developing programs for nonlinear optimization. A new approach used in the frame of symmetric sequential coloring is described. Numerical results illustrate the differences between this method and the popular Powell-Toint method.
The complex (bio)chemical reaction systems, frequently possess fast/slow phenomena, represent both difficulties and challenges for numerical simulation. We develop and test an enhancement of the classical QSSA (quasi-steady-state approximation) model reduction method applied to a system of chemical reactions. The novel model reduction method, the so-called delayed quasi-steady-state approximation method, proposed by Vejchodský (2014) and further developed by Papáček (2021) and Matonoha (2022), is...
Based on QR-like decomposition with column pivoting, a new and efficient numerical method for solving symmetric matrix inverse eigenvalue problems is proposed, which is suitable for both the distinct and multiple eigenvalue cases. A locally quadratic convergence analysis is given. Some numerical experiments are presented to illustrate our results.