亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

In this paper, we introduce a randomized algorithm for solving the non-symmetric eigenvalue problem, referred to as randomized Implicitly Restarted Arnoldi (rIRA). This method relies on using a sketch-orthogonal basis during the Arnoldi process while maintaining the Arnoldi relation and exploiting a restarting scheme to focus on a specific part of the spectrum. We analyze this method and show that it retains useful properties of the Implicitly Restarted Arnoldi (IRA) method, such as restarting without adding errors to the Ritz pairs and implicitly applying polynomial filtering. Experiments are presented to validate the numerical efficiency of the proposed randomized eigenvalue solver.

In this paper we propose and analyse a new formulation and pointwise divergence-free mixed finite element methods for the numerical approximation of Darcy--Brinkman equations in vorticity--velocity--pressure form, coupled with a transport equation for thermal energy with viscous dissipative effect and mixed Navier-type boundary conditions. The solvability analysis of the continuous and discrete problems is significantly more involved than usual as it hinges on Banach spaces needed to properly control the advective and dissipative terms in the non-isothermal energy balance equation. We proceed by decoupling the set of equations and use the Banach fixed-point theorem in combination with the abstract theory for perturbed saddle-point problems. Some of the necessary estimates are straightforward modifications of well-known results, while other technical tools require a more elaborated analysis. The velocity is approximated by Raviart--Thomas elements, the vorticity uses N\'ed\'elec spaces of the first kind, the pressure is approximated by piecewise polynomials, and the temperature by continuous and piecewise polynomials of one degree higher than pressure. Special care is needed to establish discrete inf-sup conditions since the curl of the discrete vorticity is not necessarily contained in the discrete velocity space, therefore suggesting to use two different Raviart--Thomas interpolants. A discrete fixed-point argument is used to show well-posedness of the Galerkin scheme. Error estimates in appropriate norms are derived, and a few representative numerical examples in 2D and 3D and with mixed boundary conditions are provided.

In this paper, we investigate a nonlocal traffic flow model based on a scalar conservation law, where a stochastic velocity function is assumed. In addition to the modeling, theoretical properties of the stochastic nonlocal model are provided, also addressing the question of well-posedness. A detailed numerical analysis offers insights how the stochasticity affects the evolution of densities. Finally, numerical examples illustrate the mean behavior of solutions and the influence of parameters for a large number of realizations.

In this paper, we propose a novel shape optimization approach for the source identification of elliptic equations. This identification problem arises from two application backgrounds: actuator placement in PDE-constrained optimal controls and the regularized least-squares formulation of source identifications. The optimization problem seeks both the source strength and its support. By eliminating the variable associated with the source strength, we reduce the problem to a shape optimization problem for a coupled elliptic system, known as the first-order optimality system. As a model problem, we derive the shape derivative for the regularized least-squares formulation of the inverse source problem and propose a gradient descent shape optimization algorithm, implemented using the level-set method. Several numerical experiments are presented to demonstrate the efficiency of our proposed algorithms.

We show that the limiting variance of a sequence of estimators for a structured covariance matrix has a general form that appears as the variance of a scaled projection of a random matrix that is of radial type and a similar result is obtained for the corresponding sequence of estimators for the vector of variance components. These results are illustrated by the limiting behavior of estimators for a linear covariance structure in a variety of multivariate statistical models. We also derive a characterization for the influence function of corresponding functionals. Furthermore, we derive the limiting distribution and influence function of scale invariant mappings of such estimators and their corresponding functionals. As a consequence, the asymptotic relative efficiency of different estimators for the shape component of a structured covariance matrix can be compared by means of a single scalar and the gross error sensitivity of the corresponding influence functions can be compared by means of a single index. Similar results are obtained for estimators of the normalized vector of variance components. We apply our results to investigate how the efficiency, gross error sensitivity, and breakdown point of S-estimators for the normalized variance components are affected simultaneously by varying their cutoff value.

In this work, we study and extend a class of semi-Lagrangian exponential methods, which combine exponential time integration techniques, suitable for integrating stiff linear terms, with a semi-Lagrangian treatment of nonlinear advection terms. Partial differential equations involving both processes arise for instance in atmospheric circulation models. Through a truncation error analysis, we show that previously formulated semi-Lagrangian exponential schemes are limited to first-order accuracy due to the discretization of the linear term; we then formulate a new discretization leading to second-order accuracy. Also, a detailed stability study is conducted to compare several Eulerian and semi-Lagrangian exponential schemes, as well as a well-established semi-Lagrangian semi-implicit method, which is used in operational atmospheric models. Numerical simulations of the shallow-water equations on the rotating sphere are performed to assess the orders of convergence, stability properties, and computational cost of each method. The proposed second-order semi-Lagrangian exponential method was shown to be more stable and accurate than the previously formulated schemes of the same class at the expense of larger wall-clock times; however, the method is more stable and has a similar cost compared to the well-established semi-Lagrangian semi-implicit method; therefore, it is a competitive candidate for potential operational applications in atmospheric circulation modeling.

The paper analyzes how the enlarging of the sample affects to the mitigation of collinearity concluding that it may mitigate the consequences of collinearity related to statistical analysis but not necessarily the numerical instability. The problem that is addressed is of importance in the teaching of social sciences since it discusses one of the solutions proposed almost unanimously to solve the problem of multicollinearity. For a better understanding and illustration of the contribution of this paper, two empirical examples are presented and not highly technical developments are used.

In this paper, we introduce some new notions of aging based on geometric, harmonic means of failure rate and aging intensity function. We define a generalized version of aging functions called specific interval-average geometric hazard rate, specific interval-average harmonic hazard rate. We focus on some characterization results and their inter-relationships among the resulting non-parametric classes of distributions. Monotonic nature of so defined aging classes are exhibited by some well known probability distributions. Probabilistic orders based on these functions are taken up for further study. The work is illustrated through case studies and a simulated data having applications in reliability/survival analysis.

In this article, we investigate both forward and backward problems for coupled systems of time-fractional diffusion equations, encompassing scenarios of strong coupling. For the forward problem, we establish the well-posedness of the system, leveraging the eigensystem of the corresponding elliptic system as the foundation. When considering the backward problem, specifically the determination of initial values through final time observations, we demonstrate a Lipschitz stability estimate, which is consist with the stability observed in the case of a single equation. To numerically address this backward problem, we refer to the explicit formulation of Tikhonov regularization to devise a multi-channel neural network architecture. This innovative architecture offers a versatile approach, exhibiting its efficacy in multidimensional settings through numerical examples and its robustness in handling initial values that have not been trained.

In this paper, we present a rigorous proof of the convergence of first order and second order exponential time differencing (ETD) schemes for solving the nonlocal Cahn-Hilliard (NCH) equation. The spatial discretization employs the Fourier spectral collocation method, while the time discretization is implemented using ETD-based multistep schemes. The absence of a higher-order diffusion term in the NCH equation poses a significant challenge to its convergence analysis. To tackle this, we introduce new error decomposition formulas and employ the higher-order consistency analysis. These techniques enable us to establish the $\ell^\infty$ bound of numerical solutions under some natural constraints. By treating the numerical solution as a perturbation of the exact solution, we derive optimal convergence rates in $\ell^\infty(0,T;H_h^{-1})\cap \ell^2(0,T; \ell^2)$. We conduct several numerical experiments to validate the accuracy and efficiency of the proposed schemes, including convergence tests and the observation of long-term coarsening dynamics.