Bence, Merriman, and Osher proposed an algorithm for computing the motion a hypersurface by mean curvature in terms of solutions of the usual heat equation, continually reinitialized after short time steps. In this pa...
详细信息
Bence, Merriman, and Osher proposed an algorithm for computing the motion a hypersurface by mean curvature in terms of solutions of the usual heat equation, continually reinitialized after short time steps. In this paper, applying some techniques of asymptotic analysis for the Allen-Cahn equation, we give a rate of convergence of their algorithm for the motion of a smooth and compact hypersurface by mean curvature. We also consider the special case of a circle evolving by curvature and show that our rate is optimal.
Electrical impedance tomography (EIT) is a medical imaging technique used to reconstruct images inside the domain of interest. EIT collects data on the boundary of the domain to infer the conductivity distribution ins...
详细信息
Electrical impedance tomography (EIT) is a medical imaging technique used to reconstruct images inside the domain of interest. EIT collects data on the boundary of the domain to infer the conductivity distribution inside the domain. The conductivity distribution will then be used to produce a tomographic image of the inside of the domain. This paper aims to recover geometric properties of a spherical perturbation in the conductivity inside a domain using sensitivity values of the electric potential on the boundary of the domain. The continuum model for EIT is first considered, as it holds more boundary information compared to other models of EIT. A change on the conductivity inside the domain is applied, and the impact on the electric potential is studied. The inverse EIT problem is then solved by formulating relations between the sensitivity values on the boundary and the geometric properties of the spherical perturbation: the radius and the projection onto the boundary and depth of its center. A reconstruction method using these relations is proposed and the method is examined by performing numerical simulations on different domains to model the head and the thorax. Lastly, the proposed method is applied to the complete electrode model of the EIT problem to analyze the performance of the method when the boundary data is limited on the electrodes.
The general problem considered in the paper is partitioning of a matrix operation between processors of a parallel system in an optimum load-balanced way without potential memory contention. The considered parallel sy...
详细信息
The general problem considered in the paper is partitioning of a matrix operation between processors of a parallel system in an optimum load-balanced way without potential memory contention. The considered parallel system is defined by several features the main of which is availability of a virtual shared memory divided into segments. if partitioning of a matrix operation causes parallel access to the same memory segment with writing data to the segment by al least one processor, then contention between processors arises which implies performance degradation. To eliminate such situation, a restriction is imposed on a class of possible partitionings, so that no two processors would write data to the same segment. On the resulting class of contention-free partitionings, a load-balanced optimum partitioning is defined as satisfying independent minimax criteria. The main result of the paper is an algorithm for finding the optimum partitioning by means of analytical solution of respective minimax problems. The paper also discusses implementation and performance issues related to the algorithm, on the basis of experience at Kendall Square Research Corporation, where the partitioning algorithm was used for creating high-performance parallel matrix libraries.
The kinematic shakedown theorem is formulated for some deformation processes as a kinematic extremum problem based on a polyhedral load domain and a deformation mode domain. An eigen-mode method is used to construct t...
详细信息
The kinematic shakedown theorem is formulated for some deformation processes as a kinematic extremum problem based on a polyhedral load domain and a deformation mode domain. An eigen-mode method is used to construct the deformation mode domain. Every kinematically admissible strain field within some time interval can be derived from the deformation domain and used in the proposed shakedown formulation to determine the safety load factor. Several simple shakedown problems are examined by using this proposed methodology. These numerical results are discussed and compared with available analytical and other numerical results. (C) 2002 Elsevier Science Ltd. All rights reserved.
This paper considers the problem of factorizing a matrix with missing components into a product of two smaller matrices, also known as principal component analysis with missing data (PCAMD). The Wiberg algorithm is a ...
详细信息
This paper considers the problem of factorizing a matrix with missing components into a product of two smaller matrices, also known as principal component analysis with missing data (PCAMD). The Wiberg algorithm is a numerical algorithm developed for the problem in the community of applied mathematics. We argue that the algorithm has not been correctly understood in the computer vision community. Although there are many studies in our community, almost every one of which refers to the Wiberg study, as far as we know, there is no literature in which the performance of the Wiberg algorithm is investigated or the detail of the algorithm is presented. In this paper, we present derivation of the algorithm along with a problem in its implementation that needs to be carefully considered, and then examine its performance. The experimental results demonstrate that the Wiberg algorithm shows a considerably good performance, which should contradict the conventional view in our community, namely that minimization-based algorithms tend to fail to converge to a global minimum relatively frequently. The performance of the Wiberg algorithm is such that even starting with random initial values, it converges in most cases to a correct solution, even when the matrix has many missing components and the data are contaminated with very strong noise. Our conclusion is that the Wiberg algorithm can also be used as a standard algorithm for the problems of computer vision.
This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of...
详细信息
This paper focuses on the dynamic thermo-mechanical coupled response of random particulate composite materials. Both the inertia term and coupling term are considered in the dynamic coupled problem. The formulation of the problem by a statistical second-order two-scale (SSOTS) analysis method and the algorithm procedure based on the finite-element difference method are presented. numerical results of coupled cases are compared with those of uncoupled cases. It shows that the coupling effects on temperature, thermal flux, displacement, and stresses are very distinct, and the micro- characteristics of particles affect the coupling effect of the random composites. Furthermore, the coupling effect causes a lag in the variations of temperature, thermal flux, displacement, and stresses.
Under appropriate technical assumptions, the simple-loop theory allows to derive various types of asymptotic expansions for the eigenvalues of Toeplitz matrices T-n (f) generated by a function f. Unfortunately, such a...
详细信息
Under appropriate technical assumptions, the simple-loop theory allows to derive various types of asymptotic expansions for the eigenvalues of Toeplitz matrices T-n (f) generated by a function f. Unfortunately, such a theory is not available in the preconditioning setting, that is for matrices of the formT(n)(-1)(g)T-n(l) with g,l real-valued, g nonnnegative and not identically zero almost everywhere. Independently and under the milder hypothesis that f = l/g is even and monotonic over [0, pi], matrix-less algorithms have been developed for the fast eigenvalue computation of large preconditioned matrices of the type above, within a linear complexity in the matrix order: behind the high efficiency of such algorithms there are the expansions as in the case g = 1, combined with the extrapolation idea, and hence we conjecture that the simple-loop theory has to be extended in such a new setting, as the numerics strongly suggest. Here we focus our attention on a change of variable, followed by the asymptotic expansion of the new variable, and we consider new matrix-less algorithms ad hoc for the current case. numerical experiments show a much higher accuracy till machine precision and the same linear computational cost, when compared with the matrix-less procedures already proposed in the literature.
Here we investigate the optimal harvesting problem for some periodic age-dependent population dynamics;namely, we consider the linear Lotka-McKendrick model with periodic vital rates and a periodic forcing term that s...
详细信息
Here we investigate the optimal harvesting problem for some periodic age-dependent population dynamics;namely, we consider the linear Lotka-McKendrick model with periodic vital rates and a periodic forcing term that sustains oscillations. Existence and uniqueness of a positive periodic solution are demonstrated and the existence and uniqueness of the optimal control are established. We also state necessary optimality conditions. A numerical algorithm is developed to approximate the optimal control and the optimal harvest. Some numerical results are presented.
We consider a problem of optimizing multi-cell downlink throughput in multiple-input single-output (MISO) beamforming with single user per sub-channel in the wireless communication system. Previous work based on the g...
详细信息
We consider a problem of optimizing multi-cell downlink throughput in multiple-input single-output (MISO) beamforming with single user per sub-channel in the wireless communication system. Previous work based on the generalization of uplink-downlink duality has already reformulated the maximum achievable downlink throughput into dual uplink throughput maximization problem. Since the dual uplink problem is nonconvex, it is difficult to find its optimal solution. The main contribution of this paper is a novel practical algorithm based on heuristic to find the solution of beamforer design satisfying the necessary optimality conditions of the dual uplink problem. Meanwhile the converged beamforming vectors can in turn improve the system sum rate significantly. As the dual problem is a mixed optimization, we also provide algorithms for its two sub-optimal problems. Simulation results validate the convergence and the efficiency of proposed algorithms.
Given n points in Euclidean space E (d) , we propose an algebraic algorithm to compute the best fitting (d-1)-cylinder. This algorithm computes the unknown direction of the axis of the cylinder. The location of the ax...
详细信息
Given n points in Euclidean space E (d) , we propose an algebraic algorithm to compute the best fitting (d-1)-cylinder. This algorithm computes the unknown direction of the axis of the cylinder. The location of the axis and the radius of the cylinder are deduced analytically from this direction. Special attention is paid to the case d = 3 when n = 4 and n = 5. For the former, the minimal radius enclosing cylinder is computed algebrically from constrained minimization of a quartic form of the unknown direction of the axis. For the latter, an analytical condition of existence of the circumscribed cylinder is given, and the algorithm reduces to find the zeroes of an one unknown polynomial of degree at most 6. In both cases, the other parameters of the cylinder are deduced analytically. The minimal radius enclosing cylinder is computed analytically for the regular tetrahedron and for a trigonal bipyramids family with a symmetry axis of order 3.
暂无评论