5,366 results
Search Results
2. Comment on the paper 'On conservation laws by Lie symmetry analysis for (2+1)-dimensional Bogoyavlensky–Konopelchenko equation in wave propagation' by S. Saha Ray
- Author
-
Muhammad Alim Abdulwahhab
- Subjects
Conservation law ,Wave propagation ,010102 general mathematics ,One-dimensional space ,010103 numerical & computational mathematics ,System of linear equations ,01 natural sciences ,Symmetry (physics) ,Computational Mathematics ,symbols.namesake ,Computational Theory and Mathematics ,Modeling and Simulation ,symbols ,0101 mathematics ,Noether's theorem ,Mathematics ,Mathematical physics - Abstract
In a recent paper referred to in the title, the author used the concept of quasi self-adjointness to obtain conservation laws for a system of the ( 2 + 1 ) -dimensional Bogoyavlensky–Konopelchenko equation. Apart from the adjoint system of equations, all the results on the quasi self-adjointness and the subsequent conservation laws obtained are inaccurate. In this comment, we clarify these inaccuracies and also generate conservation laws for a potential form of the underlying equation through Noether theorem.
- Published
- 2018
3. Comment on the paper 'Convection from an inverted cone in a porous medium with cross-diffusion effects, F.G. Awad, P. Sibanda, S.S. Motsa, O.D. Makinde, Comput. Math. Appl. 61 (2011) 1431–1441'
- Author
-
Asterios Pantokratoras
- Subjects
Convection ,Cross diffusion ,Mathematical analysis ,Thermodynamics ,010103 numerical & computational mathematics ,01 natural sciences ,010305 fluids & plasmas ,Computational Mathematics ,Computational Theory and Mathematics ,Cone (topology) ,Modeling and Simulation ,0103 physical sciences ,0101 mathematics ,Porous medium ,Mathematics - Published
- 2017
4. A note on the paper 'Combination of interval-valued fuzzy set and soft set' [Comput. Math. Appl. 58 (2009) 521–527]
- Author
-
Young Bae Jun and Xibei Yang
- Subjects
Discrete mathematics ,Fuzzy classification ,Mathematics::General Mathematics ,Astrophysics::High Energy Astrophysical Phenomena ,Fuzzy set ,Fuzzy subalgebra ,Type-2 fuzzy sets and systems ,Computational Mathematics ,Computational Theory and Mathematics ,Modeling and Simulation ,Modelling and Simulation ,Fuzzy number ,Fuzzy set operations ,Membership function ,Soft set ,Mathematics - Abstract
Yang (the second author of this paper) et al. (2009) [3] established the distributive law of interval-valued fuzzy soft sets. In this article, it is pointed that it is not true in general. But, in order to establish the distributive law of interval-valued fuzzy soft sets, more general notions than interval-valued fuzzy soft sets and interval-valued fuzzy soft equal are considered. Using such notions, the generalized distributive law of interval-valued fuzzy soft sets is established.
- Published
- 2011
- Full Text
- View/download PDF
5. A comment on the papers 'A study on controllability of semilinear integrodifferential systems in banach spaces' and 'controllability of neutral functional integrodifferential systems in banach spaces'
- Author
-
E. Hernández, P. Táboas, and Michelle Pierri
- Subjects
Controllability ,Neutral differential equations ,Unbounded operator ,Pure mathematics ,Differential equations in abstract spaces ,Approximation property ,Mathematical analysis ,Banach space ,Banach manifold ,Finite-rank operator ,Computational Mathematics ,Computational Theory and Mathematics ,Distributed parameter system ,Modelling and Simulation ,Modeling and Simulation ,Semigroups of linear operators ,C0-semigroup ,Mathematics - Published
- 2005
6. A note correcting the proof of a lemma in a recent paper
- Author
-
Mingshu Peng
- Subjects
Matrix difference equation ,Discrete mathematics ,Pure mathematics ,Lemma (mathematics) ,Differential equation ,Difference equation ,Positive solution ,Computational Mathematics ,Computational Theory and Mathematics ,Modeling and Simulation ,Modelling and Simulation ,Linear difference equation ,Nonoscillating ,Mathematics - Abstract
A nonoscillation criterion for a second-order linear difference equation is established correcting a result in [1].
- Published
- 2003
- Full Text
- View/download PDF
7. Brief comment on C.A. Los' papers
- Author
-
Arnold Zellner
- Subjects
Computational Mathematics ,Computational Theory and Mathematics ,Modeling and Simulation ,Modelling and Simulation ,Library science ,Mathematics - Published
- 1992
- Full Text
- View/download PDF
8. Note on a paper by H. S. Qin
- Author
-
André Vanderbauwhede
- Subjects
Scalar (mathematics) ,Mathematical analysis ,Banach space ,Codimension ,Combinatorics ,Simple eigenvalue ,Computational Mathematics ,Bifurcation theory ,Computational Theory and Mathematics ,Modeling and Simulation ,Modelling and Simulation ,Bifurcation ,Eigenvalues and eigenvectors ,Subspace topology ,Mathematics - Abstract
In a recent paper entitled "Some sufficient conditions for occurring bifurcation" [1], H. S. Qin claims to have obtained sufficient conditions for bifurcation from multiple eigenvalues. The setting is that of an equation F(x,/~) = 0, (1) where F: R" x R--, R" is smooth, with F(0,/~) = 0, V# e R. It is assumed that for some Poe R we have rank DxF(O, ~) = n - rn, with 1 l, is misleading and gives a false impression that the results have to do with "higher multiplicities". What happens is that one looks for solutions in a subspace with codimension m - 1, while assuming that: (a) m - 1 of the n scalar equations in equation (1) are automatically satisfied for x in this subspace; (b) zero is a simple eigenvalue of DxF(O, ~o) when the nonrelevant dimensions are taken out of the picture. Condition (a) is especially severe when m > 1, except when it is a consequence of symmetry considerations (e.g. Ref. [3]). To be more precise let us formulate the following "generalized" version of the Crandall- Rabinowitz theorem. Theorem I Let X and Y be Banach spaces, X 0 a closed subspace of X and Y0 a closed subspace of Y. Let F: Y x R --* Y be a smooth mapping such that and
- Full Text
- View/download PDF
9. A note correcting the proof of a lemma in a recent paper
- Author
-
Peng, Mingshu
- Subjects
- *
OSCILLATION theory of differential equations , *LINEAR differential equations , *LINEAR systems , *EQUATIONS , *MATHEMATICS - Abstract
A nonoscillation criterion for a second-order linear difference equation is established correcting a result in [1]. [Copyright &y& Elsevier]
- Published
- 2003
- Full Text
- View/download PDF
10. Preface A collection of papers to commemorate the Cornelius Lanczos centennial
- Author
-
E.L. Ortiz and T.J. Rivlin
- Subjects
Lanczos resampling ,Computational Mathematics ,Centennial ,Computational Theory and Mathematics ,Modeling and Simulation ,Modelling and Simulation ,Engineering physics ,Classics ,Mathematics - Full Text
- View/download PDF
11. On augmentation block triangular preconditioners for regularized saddle point problems
- Author
-
Shuqian Shen, Wen-Di Bao, and Ling Jian
- Subjects
Combinatorics ,Computational Mathematics ,Computational Theory and Mathematics ,Preconditioner ,Modeling and Simulation ,Saddle point ,Block (telecommunications) ,Short paper ,Applied mathematics ,Spectral analysis ,Computer Science::Numerical Analysis ,Mathematics::Numerical Analysis ,Mathematics - Abstract
Two augmentation block triangular preconditioners were introduced by Shen et?al. (2012) for the regularized saddle point problem. However, the spectral analysis of the preconditioner based on the augmentation of the ( 2 , 2 ) block was not throughly derived there. In this short paper, we give a detailed spectral analysis on this preconditioner. A numerical example is employed to validate the efficiency of the presented theoretical results.
- Published
- 2015
12. Scalable DPG multigrid solver for Helmholtz problems: A study on convergence.
- Author
-
Badger, Jacob, Henneking, Stefan, Petrides, Socratis, and Demkowicz, Leszek
- Subjects
- *
MULTIGRID methods (Numerical analysis) , *BENCHMARK problems (Computer science) , *DEGREES of freedom , *THEORY of wave motion , *HELMHOLTZ equation , *MATHEMATICS - Abstract
This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from O (10 7) degrees of freedom (DOFs) to O (10 9) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. The paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast. [ABSTRACT FROM AUTHOR]
- Published
- 2023
- Full Text
- View/download PDF
13. New Lagrangian function for nonconvex primal-dual decomposition
- Author
-
H. Mukai and Akio Tanikawa
- Subjects
0209 industrial biotechnology ,Mathematical optimization ,Optimization problem ,Short paper ,Structure (category theory) ,MathematicsofComputing_NUMERICALANALYSIS ,0211 other engineering and technologies ,Mathematics::Optimization and Control ,010103 numerical & computational mathematics ,02 engineering and technology ,01 natural sciences ,Separable space ,symbols.namesake ,020901 industrial engineering & automation ,TheoryofComputation_ANALYSISOFALGORITHMSANDPROBLEMCOMPLEXITY ,Modelling and Simulation ,Decomposition (computer science) ,0101 mathematics ,Mathematics ,021103 operations research ,Primal dual ,Computational Mathematics ,Computational Theory and Mathematics ,Lagrangian relaxation ,Modeling and Simulation ,symbols ,Lagrangian - Abstract
In this paper, a new Lagrangian function is reported which is particularly suited for large-scale nonconvex optimization problems with separable structure. Our modification convexifies the standard Lagrangian function without destroying its separable structure so that the primal-dual decomposition technique can be applied even to nonconvex optimization problems. Furthermore, the proposed Lagrangian results in two levels of iterative optimization as compared with the three levels needed for techniques recently proposed for nonconvex primal-dual decomposition.
- Full Text
- View/download PDF
14. On the existence of a global minimum in inverse parameters identification by Self-Optimizing inverse analysis method.
- Author
-
Yun, Gun Jin and Shang, Shen
- Subjects
- *
INVERSE problems , *PARAMETERS (Statistics) , *PARTIAL differential equations , *MATHEMATICS , *MATHEMATICAL optimization - Abstract
Abstract In this paper, a mathematical proof of the existence of a global minimum of Self-Optim (Self-Optimizing Inverse Analysis Method) cost functional is presented based upon weak-solution theory of partial differential equations. The Self-Optim provides single global minimum rather than having multiple global minima corresponding to unrealistic solutions of the inverse problem. Furthermore, discrete approximation of the inverse problem and computational methods for the cost functional are proposed and the proof is numerically verified. This paper provides a rigorous mathematical foundation for applications of the Self-Optim method to various inverse problems in mechanics. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
15. Virtual element approximation of two-dimensional parabolic variational inequalities
- Author
-
Sundararajan Natarajan, Dibyendu Adak, and Gianmarco Manzini
- Subjects
Polynomial ,Degrees of freedom (statistics) ,010103 numerical & computational mathematics ,01 natural sciences ,Upper and lower bounds ,Projection (linear algebra) ,010101 applied mathematics ,Computational Mathematics ,Quadratic equation ,Computational Theory and Mathematics ,Rate of convergence ,Modeling and Simulation ,Variational inequality ,Applied mathematics ,0101 mathematics ,Voronoi diagram ,Mathematics - Abstract
We design a virtual element method for the numerical treatment of the two-dimensional parabolic variational inequality problem on unstructured polygonal meshes. Due to the expected low regularity of the exact solution, the virtual element method is based on the lowest-order virtual element space that contains the subspace of the linear polynomials defined on each element. The connection between the nonnegativity of the virtual element functions and the nonnegativity of the degrees of freedom, i.e., the values at the mesh vertices, is established by applying the Maximum and Minimum Principle Theorem. The mass matrix is computed through an approximate L 2 polynomial projection, whose properties are carefully investigated in the paper. We prove the well-posedness of the resulting scheme in two different ways that reveal the contractive nature of the VEM and its connection with the minimization of quadratic functionals. The convergence analysis requires the existence of a nonnegative quasi-interpolation operator, whose construction is also discussed in the paper. The variational crime introduced by the virtual element setting produces five error terms that we control by estimating a suitable upper bound. Numerical experiments confirm the theoretical convergence rate for the refinement in space and time on three different mesh families including distorted squares, nonconvex elements, and Voronoi tesselations.
- Published
- 2022
16. Asymptotical profiles of a viral infection model with multi-target cells and spatial diffusion.
- Author
-
Wang, Xiaoyan, Yang, Junyuan, and Luo, Xiaofeng
- Subjects
- *
VIRUS diseases , *COMMUNICABLE diseases , *DYNAMICS , *MATHEMATICS , *ANALYTICAL mechanics - Abstract
Abstract In this paper, we propose a viral infection model with multi-target cells on a heterogeneous environment. Then we use the semigroup theory and a variational characterization to compute the local reproduction number R (x) and the basic reproduction number R 0. Furthermore, the model exhibits a threshold dynamics: the virus-free steady state E 0 is globally asymptotically stable if R 0 < 1 ; otherwise, the endemic steady state E ∗ is globally asymptotically stable. Finally, we perform some numerical examples to illustrate the theoretical results. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
17. Equivalent a posteriori error estimates for elliptic optimal control problems with [formula omitted]-control cost.
- Author
-
Leng, Haitao, Chen, Yanping, and Huang, Yunqing
- Subjects
- *
DISCRETIZATION methods , *MATHEMATICS , *COMPUTATIONAL fluid dynamics , *NUMERICAL solutions to differential equations , *FINITE element method - Abstract
Abstract An elliptic optimal control problem involving the L 1 norm of the control in the cost functional is considered in this paper. We use the full discretization and the variational discretization to approximate the control problem and the efficient and reliable a posteriori error estimates are obtained for the two cases. For the variational discretization, we also analyze the convergence of adaptive finite element methods. In the end, some examples are provided to validate our analysis. [ABSTRACT FROM AUTHOR]
- Published
- 2019
- Full Text
- View/download PDF
18. Some new solitonary solutions of the modified Benjamin–Bona–Mahony equation
- Author
-
Noor, Muhammad Aslam, Noor, Khalida Inayat, Waheed, Asif, and Al-Said, Eisa A.
- Subjects
- *
NONLINEAR evolution equations , *SOLITONS , *MATHEMATICAL physics , *MATHEMATICAL analysis , *NUMERICAL analysis , *MATHEMATICS - Abstract
Abstract: In this paper, we use the exp-function method to construct some new soliton solutions of the Benjamin–Bona–Mahony and modified Benjamin–Bona–Mahony equations. These equations have important and fundamental applications in mathematical physics and engineering sciences. The exp-function method is used to find the soliton solution of a wide class of nonlinear evolution equations with symbolic computation. This method provides the concise and straightforward solution in a very easier way. The results obtained in this paper can be viewed as a refinement and improvement of the previously known results. [Copyright &y& Elsevier]
- Published
- 2011
- Full Text
- View/download PDF
19. On generalized -summability factors
- Author
-
Savaş, Ekrem
- Subjects
- *
SUMMABILITY theory , *INFINITE series (Mathematics) , *MATHEMATICAL sequences , *MATHEMATICS - Abstract
Abstract: The paper deals with absolute summability factors for infinite series. The main result obtained in this paper generalizes a recent paper of Rhoades and Savas [B. Rhoades, E. Savas, A summability factor theorem for generalized absolute summability, Real Anal. Exchange 31 (2) (2005–2006) 355–363]. [Copyright &y& Elsevier]
- Published
- 2010
- Full Text
- View/download PDF
20. New bounds for the solutions of second order linear differential equations
- Author
-
Almenar, Pedro and Jódar, Lucas
- Subjects
- *
NUMERICAL solutions to linear differential equations , *MATHEMATICAL inequalities , *OSCILLATION theory of differential equations , *DIFFERENTIAL equations , *ASYMPTOTIC theory of algebraic ideals , *CALCULUS , *MATHEMATICS - Abstract
Abstract: This paper presents a new method to construct explicit bounds for the solutions of second order linear differential equations of the type ,, , which complements the method described in recent papers of the authors. The method lies in the use of Beesack’s version of Opial’s inequality. Some results on asymptotics of the solution for the case are also included. [Copyright &y& Elsevier]
- Published
- 2010
- Full Text
- View/download PDF
21. Pseudo-isochronicity in a class of septic differential systems
- Author
-
Wu, Yusen and Li, Peiluan
- Subjects
- *
SET theory , *DIFFERENTIAL equations , *MATHEMATICAL transformations , *MATHEMATICAL analysis , *MATHEMATICS , *CALCULUS - Abstract
Abstract: This paper deals with the pseudo-isochronicity for a class of septic differential systems. In this paper, we transform infinity into the origin so that the properties of infinity can be investigated with the methods developed for finite critical points. By calculating the singular point quantities and period constants of the origin, the problem of infinity being a pseudo-isochronous center has been solved in this case. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
22. Some remarks on sequences of binomial type
- Author
-
Sun, Xie-Hua
- Subjects
- *
MATHEMATICAL sequences , *BINOMIAL theorem , *CALCULUS , *FUNCTION algebras , *MATHEMATICS - Abstract
Abstract: A paper by D.L. Reiner researched function sequences of binomial type and established many interesting theorems. But because of an oversight in a lemma, some theorems in his paper need to be revised. This paper will revise these theorems and establish some new results on sequences of binomial type. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
23. A new method for roots of monic quaternionic quadratic polynomial
- Author
-
Jia, Zhigang, Cheng, Xuehan, and Zhao, Meixiang
- Subjects
- *
POLYNOMIALS , *QUADRATIC forms , *QUATERNIONS , *EXISTENCE theorems , *NUMERICAL roots , *MATHEMATICS - Abstract
Abstract: The purpose of this paper is to show how the problem of finding roots (or zeros) of the monic quaternionic quadratic polynomials (MQQP) can be solved by its equivalent real quadratic form. The real quadratic form matrices, firstly defined in this paper, are used to form a simple equivalent real quadratic form of MQQP. Some necessary and sufficient conditions for the existence of roots of MQQP are also presented. The main idea of the practical method proposed in this work can be summarized in two steps: translating MQQP into its equivalent real quadratic form, and giving directly the quaternionic roots of MQQP by solving its equivalent real quadratic form. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
24. Comments on “Notes on triple I method of fuzzy reasoning”
- Author
-
Yuan, Xue-hai and Stanley Lee, E.
- Subjects
- *
FUZZY logic , *REASONING , *MULTIPLE integrals , *GENERALIZATION , *SYSTEMS theory , *MATHEMATICS - Abstract
Abstract: In this article, we will show that (1) the results about the triple I method for fuzzy reasoning obtained in the paper “Computers and Mathematics with Applications 44 (2002) 1567–1579” are correct; (2) Example 2.1 in the paper “Notes on triple I method of fuzzy reasoning” presented by Hua-wen Liu is incorrect; (3) -triple I method presented by Liu is a generalization of -triple I method. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
25. Erratum to: “Notes on Triple I method of fuzzy reasoning” [Comput. Math. Appl. 44 (2002) 1567–1579]
- Author
-
Liu, Hua-Wen
- Subjects
- *
REASONING , *FUZZY logic , *MULTIPLE integrals , *SYSTEMS theory , *MATHEMATICS - Abstract
Abstract: The aim of this paper is to correct the results about triple I method for fuzzy reasoning obtained in paper [S. Song, C. Feng, E.S. Lee, Triple I method of fuzzy reasoning, Computers and Mathematics with Applications 44 (2002) 1567–1579]. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
26. On some new operations in soft set theory
- Author
-
Ali, M. Irfan, Feng, Feng, Liu, Xiaoyan, Min, Won Keun, and Shabir, M.
- Subjects
- *
SET theory , *UNCERTAINTY , *MATHEMATICS , *BOOLEAN algebra - Abstract
Abstract: Molodtsov introduced the theory of soft sets, which can be seen as a new mathematical approach to vagueness. In this paper, we first point out that several assertions (Proposition 2.3 (iv)–(vi), Proposition 2.4 and Proposition 2.6 (iii), (iv)) in a previous paper by Maji et al. [P.K. Maji, R. Biswas, A.R. Roy, Soft set theory, Comput. Math. Appl. 45 (2003) 555–562] are not true in general, by counterexamples. Furthermore, based on the analysis of several operations on soft sets introduced in the same paper, we give some new notions such as the restricted intersection, the restricted union, the restricted difference and the extended intersection of two soft sets. Moreover, we improve the notion of complement of a soft set, and prove that certain De Morgan’s laws hold in soft set theory with respect to these new definitions. [Copyright &y& Elsevier]
- Published
- 2009
- Full Text
- View/download PDF
27. Eigenvalue of fourth-order -point boundary value problem with derivatives
- Author
-
Zhang, Xinguang and Liu, Lishan
- Subjects
- *
DIFFERENTIAL equations , *BOUNDARY value problems , *CALCULUS , *MATHEMATICS , *MATHEMATICAL programming - Abstract
Abstract: This paper deals with the existence of positive solutions for the fourth-order nonlinear ordinary differential equation subject to the boundary conditions: where are constants such that , and . By means of a fixed-point theorem due to Krasnaselskii, some new existence results of positive solutions for the above multi-point boundary value problem are obtained, which improve the main results of Graef et al. [J.R. Graef, C. Qian, B. Yang, A three-point boundary value problem for nonlinear fourth-order differential equations, J. Math. Anal. Appl. 287 (2003) 217–233]. An example is given to demonstrate the main results of this paper. [Copyright &y& Elsevier]
- Published
- 2008
- Full Text
- View/download PDF
28. Newton-Type method for a class of mathematical programs with complementarity constraints
- Author
-
Tao, Yan
- Subjects
- *
METHODOLOGY , *MATHEMATICS , *EQUATIONS , *SIMULATION methods & models , *MATHEMATICAL models , *PROBLEM solving , *TRUTHFULNESS & falsehood , *RESEARCH - Abstract
In this paper, we present a Newton-type method for a class of mathematical programs with complementarity constraints. Under the MPEC-LICQ, we use the definition of B-stationary point to construct a constrained equations model, and apply the Newton method to solve the problem. At the end of this paper, numerical results are reported to show our method's validity. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
29. Level Crossings of an Oscillating Marked Random Walk
- Author
-
Dshalalow, J.H. and Liew, A.
- Subjects
- *
RANDOM walks , *FINANCE , *MATHEMATICS , *OSCILLATIONS , *MATHEMATICAL models - Abstract
Abstract: This paper deals with a class of real-valued random-walk processes, observed over random epochs of time, that forms a delayed renewal process. The present model does not restrict this class to a merely monotone random walk, which is easier to analyze and find explicit form functional. The objective is to find the first passage of the process exiting a rectangular set and registering the value of the process at this time, thus generalizing past models where either the observed process was monotone or the first passage time reduced to the moment of the first drop. The joint transformation of the named random characteristics of the process are derived in a closed form. The paper concludes with examples, including numerical examples, demonstrating the use of the results as well as practical applications to finance. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
30. Asymptotic behavior of solutions for a cooperation-diffusion model with a saturating interaction
- Author
-
Wang, Yuan-Ming
- Subjects
- *
BOUNDARY value problems , *NONNEGATIVE matrices , *DIFFUSION , *ASYMPTOTES , *MATHEMATICS - Abstract
Abstract: This paper is concerned with a Lotka-Volterra cooperation-diffusion model with a saturating interaction term for one species. The goal of the paper is to investigate the asymptotic behavior of the time- dependent solution in relation to the corresponding steady-state solutions under homogeneous Neumann boundary condition. Some simple and easily verifiable conditions are given to the rate constants so that for every nontrivial nonnegative initial function the corresponding time-dependent solution converges to one of the nonnegative constant steady-state solutions as time tends to infinity. This convergence result leads to the existence and uniqueness of a positive (or nonnegative) steady-state solution and the global asymptotic stability of a given nonnegative constant steady-state solution. In terms of ecological dynamics, it also gives some coexistence, permanence and extinction results for the model. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
31. Reflections on optimality and dynamic programming
- Author
-
Galperin, E.A.
- Subjects
- *
MATHEMATICAL optimization , *MATHEMATICAL programming , *DYNAMIC programming , *MATHEMATICAL analysis , *MATHEMATICS - Abstract
Abstract: Discrete models and continuous control systems are considered in regard to optimality of their trajectories. Some aspects of the principle of optimality [1, p. 83] are analyzed, and it is shown to imply total optimality, that is, the optimality of every part of an optimal trajectory. Certain autonomous systems with free admissible variations possess this property. Nonautonomous optimal systems are not, in general, totally optimal, in which case the principle of optimality is not valid. A modification is proposed for the derivation of the main functional equation to demonstrate that dynamic programming and its functional equations are valid also in the case of nonoptimal remaining trajectories under a certain contiguity condition that is defined and analyzed in the paper. Control systems with incomplete information or structural limitations on controls do not, in general, satisfy the contiguity condition. Control problems for such systems may have optimal solutions which, however, cannot be obtained by dynamic programming. This fact is shown in an example of a widely used engineering system for which an optimal trajectory has all its remaining parts nonoptimal and noncontiguous to the optimal trajectory. The paper presents theoretical justification of dynamic programming for contiguous systems that do not conform to the principle of optimality. Examples are presented to illustrate the results which open new avenues in modeling and optimization of general (not totally optimal) control systems. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
32. Fuzzy set based multiobjective allocation of resources and its applications
- Author
-
Ekel, P.Ya., Martins, C.A.P.S., Pereira, J.G., Palhares, R.M., and Canha, L.N.
- Subjects
- *
FUZZY mathematics , *MATHEMATICAL optimization , *MULTIPLE criteria decision making , *ALGORITHMS , *MATHEMATICS , *DECISION making - Abstract
Abstract: This paper presents results of research into the use of the Bellman-Zadeh approach to decision making in a fuzzy environment for solving multiobjective optimization problems. Its application conforms to the principle of guaranteed result and provides constructive lines in obtaining harmonious solutions on the basis of analyzing associated maxmin problems. The use of the Bellman-Zadeh approach has served as a basis for solving a problem of multiobjective allocation of resources (or their shortages) and developing a corresponding adaptive interactive decision-making system (AIDMS1). Its calculating kernel permits one to solve maxmin problems using an algorithm based on a nonlocal search (modification of the Gelfand''s and Tsetlin''s “long valley” method). The AIDMS1 includes procedures for considering linguistic variables to reflect conditions that are difficult to formalize as well as procedures for constructing and correcting vectors of importance factors for goals. The use of these procedures permits one to realize an adaptive approach to processing information of a decision maker to provide successive improving of the solution quality. C++ windows of the AIDMS1 are presented for input, output, and special possibilities related to considering linguistic variables and constructing and correcting vectors of importance factors. The results of the paper are universally applicable and are already being used to solve power engineering problems. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
33. Sensitivity analysis in multiobjective differential programming
- Author
-
Jiménez Guerra, P., Melguizo, M.A., and Munoz-Bouzo, M.J.
- Subjects
- *
MATHEMATICAL programming , *MATHEMATICAL optimization , *MULTIPLE criteria decision making , *MATHEMATICS , *DECISION making - Abstract
Abstract: In this paper we deal with sensitivity analysis in multiobjective differential programs with equality constraints. We analyze the quantitative behavior of the optimal solutions according to changes of right-hand side values included in the original optimization problem. One of the difficulties lies in the fact that the efficient solution in multiobjective optimization in general becomes a set. If the preference of the decision maker is represented by a scalar utility which transforms optimal solutions in optima, we may apply existing methods of sensitivity analysis. However, when dealing with a subset of optimal points, the existence of a Fréchet differentiable selection of such a set-valued map is usually assumed. The aim of the paper is to investigate the derivative of certain set-valued maps of efficient points. We show that the sensitivity depends on a set-valued map associated to the T-Lagrange multipliers and a projection of its sensitivity. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
34. The M/G/1 processor-sharing queue with disasters
- Author
-
Li, Quan-Lin and Lin, Chuang
- Subjects
- *
LAPLACE transformation , *DIFFERENTIAL equations , *MATHEMATICAL transformations , *MATHEMATICS , *DISASTERS - Abstract
Abstract: n this paper, the M/G/1 processor-sharing queue with disasters is given a detailed analysis by means of extending the supplementary variable method. The transient and steady-state distributions of the queue length are expressed as a simple and computable form, the Laplace-Stieltjes transform of the sojourn time is derived, and the Laplace transform of the busy period and its mean are obtained. Also, the approach developed in this paper is shown to be able to study more complicated M/G/1 processor-sharing models. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
35. Exponential periodicity and stability of neural networks with reaction-diffusion terms and both variable and unbounded delays
- Author
-
Zhao, Z.J., Song, Q.K., and Zhang, J.Y.
- Subjects
- *
ARTIFICIAL neural networks , *LIPSCHITZ spaces , *CONTINUOUS functions , *MATHEMATICAL variables , *MATHEMATICS - Abstract
Abstract: In this paper, the exponential periodicity and stability of neural networks with Lipschitz continuous activation functions are investigated, without assuming the boundedness of the activation functions and the differentiability of time-varying delays, as needed in most other papers. The neural networks contain reaction-diffusion terms and both variable and unbounded delays. Some sufficient conditions ensuring the existence and uniqueness of periodic solution and stability of neural networks with reaction-diffusion terms and both variable and unbounded delays are obtained by analytic methods and inequality technique. Furthermore, the exponential converging index is also estimated. The methods, which does not make use of Lyapunov functional, is simple and valid for the periodicity and stability analysis of neural networks with variable and/or unbounded delays. The results extend some previous results. Two examples are given to show the effectiveness of the obtained results. [Copyright &y& Elsevier]
- Published
- 2006
- Full Text
- View/download PDF
36. A new existence theory for positive periodic solutions to functional differential equations
- Author
-
Wan, Aying, Jiang, Daqing, and Xu, Xiaojie
- Subjects
- *
DIFFERENTIAL equations , *BESSEL functions , *CALCULUS , *EQUATIONS , *MATHEMATICS - Abstract
This paper deals with a new existence theory for positive periodic solutions to a kind of nonautonomous functional differential equation by employing the fixed-point theorem in cones. Applying the general theorems established to several biomathematical models, the paper improves some previous results and obtains some new results. [Copyright &y& Elsevier]
- Published
- 2004
- Full Text
- View/download PDF
37. Infinitely many positive solutions of the diophantine equation <F>x2 − kxy + y2 + x = 0</F>
- Author
-
Marlewski, A. and Zarzycki, P.
- Subjects
- *
ALGEBRA , *ASYMPTOTIC expansions , *ASYMPTOTES , *EQUATIONS , *MATHEMATICAL analysis , *MATHEMATICS - Abstract
We prove that the equation
x2 − kxy + y2 + x = 0 withk ∊ N+ has an infinite number of positive integer solutionsx andy if and only ifk = 3 . Fork = 3 the quotientSHAPE="SOL" ALIGN="C" STYLE="S"> is asymptotically equal to (3 + √5)/2 or (3 − √5)/2. Results of the paper are based on data obtained via Computer Algebra System (derive 5). Some derive procedures presented in the paper made it possible to discover interesting regularities concerning simple continued fractions of certain numbers. [Copyright &y& Elsevier]x y - Published
- 2004
- Full Text
- View/download PDF
38. Wavelets and regularization of the sideways heat equation
- Author
-
Qiu, Chun-Yu, Fu, Chu-Li, and Zhu, You-Bin
- Subjects
- *
PERTURBATION theory , *WAVELETS (Mathematics) , *MATHEMATICS , *MATHEMATICAL analysis , *DYNAMICS , *MATHEMATICAL physics , *FUNCTIONAL analysis , *APPROXIMATION theory - Abstract
In this paper, the following inverse heat conduction problem: ut = uxx, x≥0, t≥0, u(x,0) = 0, x≥0 u(1,t) = g(t), t≥0, u&z.sfnc;x→∞ is bounded, is considered again. This problem is severely ill-posed: its solution (if it exists) does not depend continuously on the data; a small perturbation in the data may cause a dramatically large error in the solution for
0 < x < 1 . In this paper, a new wavelet regularization method for this problem is given. Moreover, we can easily find the regularization parameter J such that some sharp stable estimates between the exact solution and the approximate one inHr(R)-norm meaning is given. [Copyright &y& Elsevier]- Published
- 2003
- Full Text
- View/download PDF
39. Conservative upwind difference schemes for the Euler equations
- Author
-
Glaister, P.
- Subjects
- *
EULER polynomials , *MATHEMATICAL sequences , *POLYNOMIALS , *EQUATIONS , *ALGEBRA , *MATHEMATICS - Abstract
In a recent paper [1] a number of numerical schemes for the shallow water equations based on a conservative linearization are analyzed. In particular, it is established that the schemes are related through the use of a source term. In this paper this technique is applied to the Euler equations, and further analysis suggests a new formulation of an existing scheme having the same key properties. [Copyright &y& Elsevier]
- Published
- 2003
- Full Text
- View/download PDF
40. Stability of the modified projected dynamical systems
- Author
-
Noor, M.A.
- Subjects
- *
RANDOM dynamical systems , *MATHEMATICS - Abstract
In this paper, we propose and analyse a modified projection-type dynamical system associated with variational inequalities by using the technique of updating the solution. We prove that the globally asymptotic stability of this dynamical system requires only the pseudomonotonicity of the underlying operator, which is a weaker condition than monotonicity. The results obtained in this paper represent a significant improvement of the previously known results. [Copyright &y& Elsevier]
- Published
- 2002
- Full Text
- View/download PDF
41. High-performance computation of pricing two-asset American options under the Merton jump-diffusion model on a GPU
- Author
-
Chittaranjan Mishra and Abhijit Ghosh
- Subjects
Computational Mathematics ,Alternating direction implicit method ,Computational Theory and Mathematics ,Complementarity theory ,Modeling and Simulation ,Fast Fourier transform ,Jump diffusion ,Graphics processing unit ,Algorithm ,Toeplitz matrix ,Cyclic reduction ,Mathematics ,Block (data storage) - Abstract
This paper is concerned with fast, parallel and numerically accurate pricing of two-asset American options under the Merton jump-diffusion model, which gives rise to a two-dimensional partial integro-differential complementarity problem (PIDCP) with a nonlocal two-dimensional integral term. Following method-of-lines approach, the solution to the PIDCP can be computed quite accurately by a robust numerical technique that combines Ikonen–Toivanen splitting with an alternating direction implicit scheme. However, we observed that computing the numerical solution with this technique becomes extremely time consuming, mainly due to the handling of the integral term. In this paper we parallelize this technique by applying a parallel fast Fourier transformation algorithm to all matrix-vector multiplications involving the huge and dense integral approximation matrix by exploiting its block Toeplitz with Toeplitz block structure. We also parallelize other computationally intensive steps of this technique by applying a recently developed parallel cyclic reduction algorithm for pentadiagonal systems. Our solutions computed on a graphics processing unit (GPU) using CUDA® platform are compared for accuracy with those available in the literature. It is observed that by solving the PIDCP parallelly we could bring down the computational times from several hours to a few seconds in certain cases in our experiments.
- Published
- 2022
42. Sixth order compact finite difference schemes for Poisson interface problems with singular sources
- Author
-
Peter D. Minev, Qiwei Feng, and Bin Han
- Subjects
Constant coefficients ,Weak solution ,Mathematical analysis ,Compact finite difference ,Dirac delta function ,010103 numerical & computational mathematics ,01 natural sciences ,010101 applied mathematics ,Computational Mathematics ,Matrix (mathematics) ,symbols.namesake ,Maximum principle ,Computational Theory and Mathematics ,Modeling and Simulation ,Dirichlet boundary condition ,symbols ,0101 mathematics ,Coefficient matrix ,Mathematics - Abstract
Let Γ be a smooth curve inside a two-dimensional rectangular region Ω. In this paper, we consider the Poisson interface problem − ∇ 2 u = f in Ω ∖ Γ with Dirichlet boundary condition such that f is smooth in Ω ∖ Γ and the jump functions [ u ] and [ ∇ u ⋅ n → ] across Γ are smooth along Γ. This Poisson interface problem includes the weak solution of − ∇ 2 u = f + g δ Γ in Ω as a special case. Because the source term f is possibly discontinuous across the interface curve Γ and contains a delta function singularity along the curve Γ, both the solution u of the Poisson interface problem and its flux ∇ u ⋅ n → are often discontinuous across the interface. To solve the Poisson interface problem with singular sources, in this paper we propose a sixth order compact finite difference scheme on uniform Cartesian grids. Our proposed compact finite difference scheme with explicitly given stencils extends the immersed interface method (IIM) to the highest possible accuracy order six for compact finite difference schemes on uniform Cartesian grids, but without the need to change coordinates into the local coordinates as in most papers on IIM in the literature. Also in contrast with most published papers on IIM, we explicitly provide the formulas for all involved stencils and therefore, our proposed scheme can be easily implemented and is of interest to practitioners dealing with Poisson interface problems. Note that the curve Γ splits Ω into two disjoint subregions Ω + and Ω − . The coefficient matrix A in the resulting linear system A x = b , following from the proposed scheme, is independent of any source term f, jump condition g δ Γ , interface curve Γ and Dirichlet boundary conditions, while only b depends on these factors and is explicitly given, according to the configuration of the nine stencil points in Ω + or Ω − . The constant coefficient matrix A facilitates the parallel implementation of the algorithm in case of a large size matrix and only requires the update of the right hand side vector b for different Poisson interface problems. Due to the flexibility and explicitness of the proposed scheme, it can be generalized to obtain the highest order compact finite difference scheme for non-uniform grids as well. We prove the order 6 convergence for the proposed scheme using the discrete maximum principle. Our numerical experiments confirm the sixth accuracy order of the proposed compact finite difference scheme on uniform meshes for the Poisson interface problems with various singular sources.
- Published
- 2021
43. Recovering the aqueous concentration in a multi-layer porous media
- Author
-
Quan Pham Hoang
- Subjects
Computational Mathematics ,Exact solutions in general relativity ,Computational Theory and Mathematics ,Rate of convergence ,Modeling and Simulation ,Applied mathematics ,Point (geometry) ,Filter (signal processing) ,Inverse problem ,Layer (object-oriented design) ,Porous medium ,Regularization (mathematics) ,Mathematics - Abstract
In this paper, we consider an inverse problem for time-fractional advection-dispersion equation in a multi-layer composite medium. The main goal of our paper is to approximate the initial information, which is inaccessible for measurement, from the observation data at a certain point in second layer by constructing a regularized solution using a filter regularization method. Under appropriate regularity assumptions of the exact solution, we give convergence rate of the error between the reconstructed solution and the exact one. A numerical example is provided to illustrate the main results.
- Published
- 2021
44. On convergence of a structure preserving difference scheme for two-dimensional space-fractional nonlinear Schrödinger equation and its fast implementation
- Author
-
Yushun Wang, Yuezheng Gong, and Dongdong Hu
- Subjects
Discretization ,Preconditioner ,Fast Fourier transform ,Finite difference ,Krylov subspace ,Solver ,Computational Mathematics ,symbols.namesake ,Computational Theory and Mathematics ,Modeling and Simulation ,Norm (mathematics) ,symbols ,Applied mathematics ,Nonlinear Schrödinger equation ,Mathematics - Abstract
In this paper we intend to construct a structure preserving difference scheme for two-dimensional space-fractional nonlinear Schrodinger (2D SFNS) equation with the integral fractional Laplacian. The temporal direction is discretized by the modified Crank-Nicolson method, and the spatial variable is approximated by a novel fractional central difference method. The mass and energy conservations and the convergence are rigorously proved for the proposed scheme. For 1D SFNS equation, the convergence relies heavily on the L ∞ -norm boundness of the numerical solution of the proposed scheme. However, we cannot obtain the L ∞ -norm boundness of the numerical solution by using the similar process for the 2D SFNS equation. One of the major significance of this paper is that we first obtain the L ∞ -norm boundness of the numerical solution and L 2 -norm error estimate via the popular “cut-off” function for the 2D SFNS equation. Further, we reveal that the spatial discretization generates a block-Toeplitz coefficient matrix, and it will be ill-conditioned as the spatial grid mesh number M and the fractional order α increase. Thus, we exploit an linearized iteration algorithm for the nonlinear system, so that it can be efficiently solved by the Krylov subspace solver with a suitable preconditioner, where the 2D fast Fourier transform (2D FFT) is applied in the solver to accelerate the matrix-vector product, and the standard orthogonal projection approach is used to eliminate the drift of mass and energy. Extensive numerical results are reported to confirm the theoretical analysis and high efficiency of the proposed algorithm.
- Published
- 2021
45. Decay for solutions of a nonlinear damped wave equation with variable-exponent nonlinearities.
- Author
-
Messaoudi, Salim A., Al-Smail, Jamal H., and Talahmeh, Ala A.
- Subjects
- *
NONLINEAR equations , *EXPONENTS , *EQUATIONS , *ALGEBRA , *MATHEMATICS - Abstract
Abstract In this paper, we consider the following nonlinear waveequation with variable exponents: u t t − div (| ∇ u | r (⋅) − 2 ∇ u) + | u t | m (⋅) − 2 u t = 0. By using a lemma by Komornik, we prove the decay estimates for the solution under suitable assumptions on the variable exponents m , r and the initial data. We also give two numerical applications to illustrate our theoretical results. [ABSTRACT FROM AUTHOR]
- Published
- 2018
- Full Text
- View/download PDF
46. Unconditionally optimal error estimates of two linearized Galerkin FEMs for the two-dimensional nonlinear fractional Rayleigh–Stokes problem
- Author
-
Zhen Guan, Jungang Wang, and Yufeng Nie
- Subjects
Discrete system ,Computational Mathematics ,Nonlinear system ,Elliptic curve ,Exact solutions in general relativity ,Operator (computer programming) ,Computational Theory and Mathematics ,Modeling and Simulation ,Norm (mathematics) ,Applied mathematics ,Galerkin method ,Finite element method ,Mathematics - Abstract
In this paper, two linearized Galerkin finite element methods, which are based on the L1 approximation and the WSGD operator, respectively, are proposed to solve the nonlinear fractional Rayleigh-Stokes problem. In order to obtain the unconditionally optimal error estimate, we firstly introduce a time-discrete elliptic equation, and derive the unconditional error estimate between the exact solution and the solution of the time-discrete system in H 2 -norm. Secondly, we obtain the boundedness of the fully discrete finite element solution in L ∞ -norm through the more detailed study of the error equation. Then, the optimal L 2 -norm error estimate is derived for the fully discrete system without any restriction conditions on the time step size. Finally, some numerical experiments are presented to confirm the theoretical results, showing that the two linearized schemes given in this paper are efficient and reliable.
- Published
- 2021
47. Polar differentiation matrices for the Laplace equation in the disk under nonhomogeneous Dirichlet, Neumann and Robin boundary conditions and the biharmonic equation under nonhomogeneous Dirichlet conditions
- Author
-
Marcela Molina Meyer and Frank Richard Prieto Medina
- Subjects
Laplace's equation ,Dirichlet conditions ,010103 numerical & computational mathematics ,01 natural sciences ,Dirichlet distribution ,Robin boundary condition ,Quadrature (mathematics) ,010101 applied mathematics ,Computational Mathematics ,symbols.namesake ,Computational Theory and Mathematics ,Modeling and Simulation ,Ordinary differential equation ,symbols ,Biharmonic equation ,Applied mathematics ,Pseudo-spectral method ,0101 mathematics ,Mathematics - Abstract
In this paper we present a pseudospectral method in the disk. Unlike the methods already known, the disk is not duplicated. Moreover, we solve the Laplace equation under nonhomogeneous Dirichlet, Neumann and Robin boundary conditions, as well as the biharmonic equation subject to nonhomogeneous Dirichlet conditions, by only using the elements of the corresponding differentiation matrices. It is worth mentioning that we do not use any quadrature, nor need to solve any decoupled system of ordinary differential equations, nor use any pole condition, nor require any lifting. We also solve several numerical examples to show the spectral convergence. The pseudospectral method developed in this paper is applied to estimate Sherwood numbers integrating the mass flux to the disk, and it can be implemented to solve Lotka–Volterra systems and nonlinear diffusion problems involving chemical reactions.
- Published
- 2021
48. Generalization of the Multipoint meshless FDM application to the nonlinear analysis
- Author
-
Irena Jaworska
- Subjects
Geometrically nonlinear ,Generalization ,MathematicsofComputing_NUMERICALANALYSIS ,Finite difference method ,010103 numerical & computational mathematics ,Computer Science::Numerical Analysis ,01 natural sciences ,Collatz conjecture ,Computer Science::Performance ,010101 applied mathematics ,Computational Mathematics ,Nonlinear system ,Computational Theory and Mathematics ,Modeling and Simulation ,Computer Science::Networking and Internet Architecture ,Order (group theory) ,Applied mathematics ,Computer Science::Symbolic Computation ,Boundary value problem ,0101 mathematics ,Mathematics - Abstract
The paper focuses on the new Multipoint meshless finite difference method, following the original Collatz higher order multipoint concept and the essential ideas of the Meshless FDM. The method was formulated, developed, and tested for various boundary value problems. Generalization of the multipoint method application to nonlinear analysis is the purpose of this research. The first attempt of the multipoint technique application to the geometrically nonlinear problems was successfully done recently. The case of physically nonlinear problem is considered in this paper. Several benefits of the proposed approach are highlighted, numerical algorithm and selected results are presented, and application of the multipoint method to nonlinear analysis is summarized.
- Published
- 2021
49. On the sources placement in the method of fundamental solutions for time-dependent heat conduction problems
- Author
-
Jakub Krzysztof Grabski
- Subjects
Boundary (topology) ,010103 numerical & computational mathematics ,Thermal conduction ,Space (mathematics) ,01 natural sciences ,Domain (mathematical analysis) ,010101 applied mathematics ,Computational Mathematics ,Distribution (mathematics) ,Computational Theory and Mathematics ,Modeling and Simulation ,Method of fundamental solutions ,Applied mathematics ,Transient (oscillation) ,0101 mathematics ,Value (mathematics) ,Mathematics - Abstract
The method of fundamental solutions is a more and more popular meshless method for solving boundary or initial–boundary value problems. The most important issue in this method is the determination of the positions of the source points. The accuracy of the method depends strongly on the distribution of the source points. In this paper placement of these points for transient heat conduction problems is studied. The problems are initial–boundary value problems and they are considered in a time-space domain. Because of that, the placement of the source points differs from the classical distribution of the source points for boundary values problems. In the paper, four different possible sources distributions are considered for 1D, 2D and 3D transient heat conduction problems. The results show very good accuracy in case of the source points placed in a space much bigger than the considered region, additionally with the negative time coordinate.
- Published
- 2021
50. Contraction operator transformation for the complex heterogeneous Helmholtz equation
- Author
-
N. Yavich, Michael S. Zhdanov, Nikolay I. Khokhlov, and M. Malovichko
- Subjects
Helmholtz equation ,Discretization ,Preconditioner ,Fast Fourier transform ,010103 numerical & computational mathematics ,Solver ,Computer Science::Numerical Analysis ,01 natural sciences ,010101 applied mathematics ,Computational Mathematics ,Matrix (mathematics) ,Transformation (function) ,Computational Theory and Mathematics ,Modeling and Simulation ,Convergence (routing) ,Applied mathematics ,0101 mathematics ,Mathematics - Abstract
An efficient solution of the three-dimensional Helmholtz equation is known to be crucial in many applications, especially geophysics. In this paper, we present and test two preconditioning approaches for the discrete problem resulting from the second order finite-difference discretization of this equation. The first approach combines shifted-Laplacian preconditioner with inversion of a separable matrix, corresponding to the horizontally-layered velocity model, using fast Fourier based transforms. The second approach is novel and involves a special transformation resulting in a preconditioner with a contraction operator (CO preconditioner). The two approaches have near the same arithmetical complexity; however, the second approach, developed in this paper, provides a faster convergence of an iterative solver as illustrated by numerical experiments and analysis of the spectral properties of the preconditioned matrices. Our numerical experiments involve parallel modeling of highly heterogeneous lossy and lossless media at different frequencies. We show that the CO-based solver can tackle problems with hundreds of millions of unknowns on a conventional cluster node. The CO preconditioned solver demonstrates a very moderate increase of iteration count with the frequency. We have conducted a comparison of the performance of the developed method versus open-source parallel sweeping preconditioner. The results indicate that, the CO solver is several times faster with respect to the wall-clock time and consumes substantially less memory than the code based on the sweeping preconditioner at least in the example we tested.
- Published
- 2021
Discovery Service for Jio Institute Digital Library
For full access to our library's resources, please sign in.