• Journal of applied mathematics & informatics
• /
• v.9 no.1
• /
• pp.101-124
• /
• 2002

The convergence rate of a numerical procedure barred on Schwarz Alternating Method (SAM) for solving elliptic boundary value problems (BVP's) depends on the selection of the interface conditions applied on the interior boundaries of the overlapping subdomains. It hee been observed that the Robin condition(mixed interface condition), controlled by a parameter, can optimize SAM's convergence rate. Since the convergence rate is very sensitive to the parameter, Tang[17] suggested another interface condition called over-determined interface condition. Based on the over-determined interface condition, we formulate the two-layer multi-parameterized SAM. For the SAM and the one-dimensional elliptic model BVP's, we determine analytically the optimal values of the parameters. For the two-dimensional elliptic BVP's , we also formulate the two-layer multi-parameterized SAM and suggest a choice of multi-parameter to produce good convergence rate .

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.32 no.10
• /
• pp.335-340
• /
• 1983

The analysis of static power converter circuit using state-space method is presented. Semiconductors are modeled in two-state resistors depending on their ON or OFF states. Then the modes of circuit are determined according to the conducting states of semiconductors and different describing matrices are given automatically for each mode. Newton-Raphson algorithm is used as an iterative method for obtaining steady-state solution and an adjoint network is introduced for the efficient and accurate evaluation of the Jacobi matrix in the algorithm. Using the porogram exploited from the above algorithm, it is shown through examples that the results are in good agreement with the analytic solutions and computation time is considerably reduced for obtaining the steady-state solutions.

• Journal of applied mathematics & informatics
• /
• v.25 no.1_2
• /
• pp.435-443
• /
• 2007

Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

• Journal of the Korea Society of Computer and Information
• /
• v.17 no.9
• /
• pp.9-16
• /
• 2012

In this paper, we propose the two dimensional variable-length vector storage format which can be used for efficient storage of sparse matrix in the FEM (finite element method). The proposed storage format is the method storing only actual needed non-zero values of each row on upper triangular matrix with the total rows N, by using two dimensional variable-length vector instead of $N{\times}N$ large sparse matrix of entire equation of finite elements. This method only needs storage spaces of the number of minimum 1 to maximum 5 in 2D grid structure and the number of minimum 1 to maximum 14 in 3D grid structure of analysis target. The number doesn't excess two times although involving index number. From the experimental result, we can find out that the proposed storage format can reduce the memory space more effectively, as the total number of nodes increases, than the existing skyline storage format storing maximum column height.