Leszek Marcinkowski

Konrad Sakowski, Leszek Marcinkowski, Pawel Strak et al.

In this paper, a variant of discretization of the van Roosbroeck equations in the equilibrium state with the Composite Discontinuous Galerkin Method for the rectangular domain is discussed. It is based on Symmetric Interior Penalty Galerkin (SIPG) method. The proposed method accounts for lower regularity of the solution on the interfaces of devices' layers. It is shown that the discrete problem is well-defined and that discrete solution is unique. Error estimates are derived. Finally, numerical simulations are presented.

Erik Eikeland, Leszek Marcinkowski, Talal Rahman

We propose two variants of the overlapping additive Schwarz method for the finite element discretiza- tion of the elliptic problem in 3D with highly heterogeneous coefficients. The methods are efficient and simple to construct using the abstract framework of the additive Schwarz method, and an idea of adaptive coarse spaces. In one variant, the coarse space consists of finite element functions associated with the wire basket nodes and functions based on solving some generalized eigenvalue problem on the faces, and in the other variant, it contains functions associated with the vertex nodes with functions based on solving some generalized eigenvalue problems on subdomain faces and on subdomain edges. The functions that are used to build the coarse spaces are chosen adaptively, they correspond to the eigenvalues that are smaller than a given threshold. The convergence rate of the preconditioned conjugate gradients method in both cases, is shown to be independent of the variations in the coefficients for sufficient number of eigenfunctions in the coarse space. Numerical results are given to support the theory.

Leszek Marcinkowski, Talal Rahman

We present an analysis of the additive average Schwarz preconditioner with two newly proposed adaptively enriched coarse spaces which was presented at the 23rd International conference on domain decomposition methods in Korea, for solving second order elliptic problems with highly varying and discontinuous coefficients. It is shown that the condition number of the preconditioned system is bounded independently of the variations and the jumps in the coefficient, and depends linearly on the mesh parameter ratio H/h, that is the ratio between the subdomain size and the mesh size, thereby retaining the same optimality and scalablity of the original additive average Schwarz preconditioner.

Atle Loneland, Leszek Marcinkowski, Talal Rahman

In this paper we introduce an additive Schwarz method for a Crouzeix-Raviart Finite Volume Element (CRFVE) discretization of a second order elliptic problem with discontinuous coefficients, where the discontinuities are both inside the subdomains and across and along the subdomain boundaries. We show that, depending on the distribution of the coefficient in the model problem, the parameters describing the GMRES convergence rate of the proposed method depend linearly or quadratically on the mesh parameters $H/h$. Also, under certain restrictions on the distribution of the coefficient, the convergence of the GMRES method is independent of jumps in the coefficient.

Bin Wu, Leszek Marcinkowski, Xue-Cheng Tai et al.

We propose a set of iterative regularization algorithms for the TV-Stokes model to restore images from noisy images with Gaussian noise. These are some extensions of the iterative regularization algorithm proposed for the classical Rudin-Osher-Fatemi (ROF) model for image reconstruction, a single step model involving a scalar field smoothing, to the TV-Stokes model for image reconstruction, a two steps model involving a vector field smoothing in the first and a scalar field smoothing in the second. The iterative regularization algorithms proposed here are Richardson's iteration like. We have experimental results that show improvement over the original method in the quality of the restored image. Convergence analysis and numerical experiments are presented.

Erik Eikeland, Leszek Marcinkowski, Talal Rahman

In this paper, we propose a two-level overlapping additive Schwarz domain decomposition preconditioner for the symmetric interior penalty discontinuous Galerkin method for the second order elliptic boundary value problem with highly heterogeneous coefficients. A specific feature of this preconditioner is that it is based on adaptively enriching its coarse space with functions created by solving generalized eigenvalue problems on thin patches covering the subdomain interfaces. It is shown that the condition number of the underlined preconditioned system is independent of the contrast if an adequate number of functions are used to enrich the coarse space. Numerical results are provided to confirm this claim.

Leszek Marcinkowski, Talal Rahman, Atle Loneland et al.

A symmetric and a nonsymmetric variant of the additive Schwarz preconditioner are proposed for the solution of a nonsymmetric system of algebraic equations arising from a general finite volume element discretization of symmetric elliptic problems with large jumps in the entries of the coefficient matrices across subdomains. It is shown in the analysis, that the convergence of the preconditioned GMRES iteration with the proposed preconditioners, depends polylogarithmically on the mesh parameters, in other words, the convergence is only weakly dependent on the mesh parameters, and it is robust with respect to the jumps in the coefficients.