An analysis on local convergence of inexact newton-gauss method for solving singular systems of equations.
Ontology highlight
ABSTRACT:
SUBMITTER: Zhou F
PROVIDER: S-EPMC3984777 | biostudies-other | 2014
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We present a fixed-point iterative method for solving systems of nonlinear equations. The convergence theorem of the proposed method is proved under suitable conditions. In addition, some numerical results are also reported in the paper, which confirm the good theoretical properties of our approach.
| S-EPMC3982483 | biostudies-other
Project description:We give an operator-algebraic interpretation of the notion of an ideal generated by the unbounded operators associated with the elements of the Lie algebra of a Lie group that implements the symmetries of a quantum system. We use this interpretation to establish a link between Rieffel induction and the implementation of a local Gauss law in lattice gauge theories similar to the method discussed by Kijowski and Rudolph (J Math Phys 43:1796-1808, 2002; J Math Phys 46:032303, 2004).
| S-EPMC6182777 | biostudies-other
Project description:Elastic materials are ubiquitous in nature and indispensable components in man-made devices and equipments. When a device or equipment involves composite or multiple elastic materials, elasticity interface problems come into play. The solution of three dimensional (3D) elasticity interface problems is significantly more difficult than that of elliptic counterparts due to the coupled vector components and cross derivatives in the governing elasticity equation. This work introduces the matched interface and boundary (MIB) method for solving 3D elasticity interface problems. The proposed MIB elasticity interface scheme utilizes fictitious values on irregular grid points near the material interface to replace function values in the discretization so that the elasticity equation can be discretized using the standard finite difference schemes as if there were no material interface. The interface jump conditions are rigorously enforced on the intersecting points between the interface and the mesh lines. Such an enforcement determines the fictitious values. A number of new techniques has been developed to construct efficient MIB elasticity interface schemes for dealing with cross derivative in coupled governing equations. The proposed method is extensively validated over both weak and strong discontinuity of the solution, both piecewise constant and position-dependent material parameters, both smooth and nonsmooth interface geometries, and both small and large contrasts in the Poisson's ratio and shear modulus across the interface. Numerical experiments indicate that the present MIB method is of second order convergence in both L∞ and L2 error norms for handling arbitrarily complex interfaces, including biomolecular surfaces. To our best knowledge, this is the first elasticity interface method that is able to deliver the second convergence for the molecular surfaces of proteins..
| S-EPMC4404754 | biostudies-literature
Project description:This paper proposes a new method for solving fuzzy system of linear equations with crisp coefficients matrix and fuzzy or interval right hand side. Some conditions for the existence of a fuzzy or interval solution of m × n linear system are derived and also a practical algorithm is introduced in detail. The method is based on linear programming problem. Finally the applicability of the proposed method is illustrated by some numerical examples.
| S-EPMC4152988 | biostudies-other
Project description:In this paper, we present a numerical scheme used to solve the nonlinear time fractional Navier-Stokes equations in two dimensions. We first employ the meshless local Petrov-Galerkin (MLPG) method based on a local weak formulation to form the system of discretized equations and then we will approximate the time fractional derivative interpreted in the sense of Caputo by a simple quadrature formula. The moving Kriging interpolation which possesses the Kronecker delta property is applied to construct shape functions. This research aims to extend and develop further the applicability of the truly MLPG method to the generalized incompressible Navier-Stokes equations. Two numerical examples are provided to illustrate the accuracy and efficiency of the proposed algorithm. Very good agreement between the numerically and analytically computed solutions can be observed in the verification. The present MLPG method has proved its efficiency and reliability for solving the two-dimensional time fractional Navier-Stokes equations arising in fluid dynamics as well as several other problems in science and engineering.
| S-EPMC4823229 | biostudies-other
Project description:In this paper, we derive fully implementable first-order time-stepping schemes for McKean-Vlasov stochastic differential equations, allowing for a drift term with super-linear growth in the state component. We propose Milstein schemes for a time-discretized interacting particle system associated with the McKean-Vlasov equation and prove strong convergence of order 1 and moment stability, taming the drift if only a one-sided Lipschitz condition holds. To derive our main results on strong convergence rates, we make use of calculus on the space of probability measures with finite second-order moments. In addition, numerical examples are presented which support our theoretical findings.
| S-EPMC7897642 | biostudies-literature
Project description:The question of finite-time blowup of the 3D incompressible Euler equations is numerically investigated in a periodic cylinder with solid boundaries. Using rotational symmetry, the equations are discretized in the (2D) meridian plane on an adaptive (moving) mesh and is integrated in time with adaptively chosen time steps. The vorticity is observed to develop a ring-singularity on the solid boundary with a growth proportional to ∼(ts - t)(-2.46), where ts ∼ 0.0035056 is the estimated singularity time. A local analysis also suggests the existence of a self-similar blowup. The simulations stop at τ(2) = 0.003505 at which time the vorticity amplifies by more than (3 × 10(8))-fold and the maximum mesh resolution exceeds (3 × 10(12))(2). The vorticity vector is observed to maintain four significant digits throughout the computations.
| S-EPMC4246962 | biostudies-literature
Project description:The local behavior of the lowest order boundary element method on quasi-uniform meshes for Symm's integral equation and the stabilized hyper-singular integral equation on polygonal/polyhedral Lipschitz domains is analyzed. We prove local a priori estimates in L2 for Symm's integral equation and in H1 for the hyper-singular equation. The local rate of convergence is limited by the local regularity of the sought solution and the sum of the rates given by the global regularity and additional regularity provided by the shift theorem for a dual problem.
| S-EPMC6154049 | biostudies-other
Project description:To combine a feedforward neural network (FNN) and Lie group (symmetry) theory of differential equations (DEs), an alternative artificial NN approach is proposed to solve the initial value problems (IVPs) of ordinary DEs (ODEs). Introducing the Lie group expressions of the solution, the trial solution of ODEs is split into two parts. The first part is a solution of other ODEs with initial values of original IVP. This is easily solved using the Lie group and known symbolic or numerical methods without any network parameters (weights and biases). The second part consists of an FNN with adjustable parameters. This is trained using the error back propagation method by minimizing an error (loss) function and updating the parameters. The method significantly reduces the number of the trainable parameters and can more quickly and accurately learn the real solution, compared to the existing similar methods. The numerical method is applied to several cases, including physical oscillation problems. The results have been graphically represented, and some conclusions have been made.
| S-EPMC8985966 | biostudies-literature
Project description:In this research, we have investigated doubly singular ordinary differential equations and a real application problem of studying the temperature profile in a porous fin model. We have suggested a novel soft computing strategy for the training of unknown weights involved in the feed-forward artificial neural networks (ANNs). Our neuroevolutionary approach is used to suggest approximate solutions to a highly nonlinear doubly singular type of differential equations. We have considered a real application from thermodynamics, which analyses the temperature profile in porous fins. For this purpose, we have used the optimizer, namely, the fractional-order particle swarm optimization technique (FO-DPSO), to minimize errors in solutions through fitness functions. ANNs are used to design the approximate series of solutions to problems considered in this paper. We find the values of unknown weights such that the approximate solutions to these problems have a minimum residual error. For global search in the domain, we have initialized FO-DPSO with random solutions, and it collects best so far solutions in each generation/ iteration. In the second phase, we have fine-tuned our algorithm by initializing FO-DPSO with the collection of best so far solutions. It is graphically illustrated that this strategy is very efficient in terms of convergence and minimum mean squared error in its best solutions. We can use this strategy for the higher-order system of differential equations modeling different important real applications.
| S-EPMC7347205 | biostudies-literature