• Communications of the Korean Mathematical Society
• /
• v.34 no.2
• /
• pp.657-670
• /
• 2019

In this paper, we consider the preconditioned iterative methods for solving linear complementarity problem associated with an M-matrix. Based on the generalized Gunawardena's preconditioner, two preconditioned SSOR methods for solving the linear complementarity problem are proposed. The convergence of the proposed methods are analyzed, and the comparison results are derived. The comparison results showed that preconditioned SSOR methods accelerate the convergent rate of the original SSOR method. Numerical examples are used to illustrate the theoretical results.

• Journal of applied mathematics & informatics
• /
• v.31 no.5_6
• /
• pp.769-784
• /
• 2013

We consider the modulus-based successive overrelaxation method for the linear complementarity problems from the discretization of Black-Scholes American options model. The $H_+$-matrix property of the system matrix discretized from American option pricing which guarantees the convergence of the proposed method for the linear complementarity problem is analyzed. Numerical experiments confirm the theoretical analysis, and further show that the modulus-based successive overrelaxation method is superior to the classical projected successive overrelaxation method with optimal parameter.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2009.05a
• /
• pp.579-580
• /
• 2009

This paper presents a numerical simulation of quasi-brittle fracture in UHPFRC I-beam. A linear complementarity problem (LCP) is used to formulate the path-dependent hardening-softening behavior in non-holonomic rate form, and the PATH solver is employed to solve the LCP.

• Honam Mathematical Journal
• /
• v.19 no.1
• /
• pp.157-164
• /
• 1997

We are concerned with three classes of matrices that are relevant to the linear complementary problem. We prove that within the class of $P_0$-matrices, the Q-matrices are precisely the regular matrices and we show that the same characterizations hold for an L-matrix as well, and that the symmetric copositive-plus Q-matrices are precisely those which are strictly copositive.

• Communications of the Korean Mathematical Society
• /
• v.25 no.4
• /
• pp.655-669
• /
• 2010

In this paper we propose new primal-dual interior point methods (IPMs) for $P_*(\kappa)$ linear complementarity problems (LCPs) and analyze the iteration complexity of the algorithm. New search directions and proximity measures are defined based on a class of kernel functions, $\psi(t)=\frac{t^2-1}{2}-{\int}^t_1e{^{q(\frac{1}{\xi}-1)}d{\xi}$, $q\;{\geq}\;1$. If a strictly feasible starting point is available and the parameter $q\;=\;\log\;\(1+a{\sqrt{\frac{2{\tau}+2{\sqrt{2n{\tau}}+{\theta}n}}{1-{\theta}}\)$, where $a\;=\;1\;+\;\frac{1}{\sqrt{1+2{\kappa}}}$, then new large-update primal-dual interior point algorithms have $O((1\;+\;2{\kappa})\sqrt{n}log\;n\;log\;{\frac{n}{\varepsilon}})$ iteration complexity which is the best known result for this method. For small-update methods, we have $O((1\;+\;2{\kappa})q{\sqrt{qn}}log\;{\frac{n}{\varepsilon}})$ iteration complexity.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2026-2038
• /
• 1994

A numerical scheme is developed for the analysis of incipient sliding contact with orthotropic friction condition subjected to tangential load and twisting moment. The inherent nonlinearity in the orthotropic friction law has been treated by a polyhedral friction law. Then, a three-dimensional linear complementarity problem(LCP) formulation in an incremental form is obtained, and the existence of a solution is investigated. A Lemke's complementary pivoting algorithm is used for solving the LCP. The scheme is illustrated by spherical contact problems, and the effects of eccentricity of elliptical friction domain on the traction and stick region are discussed.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.3
• /
• pp.189-202
• /
• 2009

In this paper we propose new primal-dual interior point algorithms for $P_*(\kappa)$ linear complementarity problems based on a new class of kernel functions which contains the kernel function in [8] as a special case. We show that the iteration bounds are $O((1+2\kappa)n^{\frac{9}{14}}\;log\;\frac{n{\mu}^0}{\epsilon}$) for large-update and $O((1+2\kappa)\sqrt{n}log\frac{n{\mu}^0}{\epsilon}$) for small-update methods, respectively. This iteration complexity for large-update methods improves the iteration complexity with a factor $n^{\frac{5}{14}}$ when compared with the method based on the classical logarithmic kernel function. For small-update, the iteration complexity is the best known bound for such methods.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1511-1523
• /
• 2011

In this paper, we consider the smoothing Newton method for the nonlinear complementarity problems with $P_0$-function. The proposed algorithm is based on a new smoothing function and it needs only to solve one linear system of equations and perform one line search per iteration. Under the condition that the solution set is nonempty and bounded, the proposed algorithm is proved to be convergent globally. Furthermore, the local superlinearly(quadratic) convergence is established under suitable conditions. Preliminary numerical results show that the proposed algorithm is very promising.

• Journal of applied mathematics & informatics
• /
• v.32 no.1_2
• /
• pp.211-225
• /
• 2014

A new class of smoothing functions is introduced in this paper, which includes some important smoothing complementarity functions as its special cases. Based on this new smoothing function, we proposed a smoothing Newton method. Our algorithm needs only to solve one linear system of equations. Without requiring the nonemptyness and boundedness of the solution set, the proposed algorithm is proved to be globally convergent. Numerical results indicate that the smoothing Newton method based on the new proposed class of smoothing functions with ${\theta}{\in}(0,1)$ seems to have better numerical performance than those based on some other important smoothing functions, which also demonstrate that our algorithm is promising.

• Journal of applied mathematics & informatics
• /
• v.29 no.3_4
• /
• pp.653-668
• /
• 2011

This paper presents a new hybrid method for solving nonlinear complementarity problems with $P_0$-functions. It can be regarded as a combination of smoothing trust region method with ODE-based method and line search technique. A feature of the proposed method is that at each iteration, a linear system is only solved once to obtain a trial step, thus avoiding solving a trust region subproblem. Another is that when a trial step is not accepted, the method does not resolve the linear system but generates an iterative point whose step-length is defined by a line search. Under some conditions, the method is proven to be globally and superlinearly convergent. Preliminary numerical results indicate that the proposed method is promising.