skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Accelerating the DC algorithm for smooth functions

; ;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Published Article
Journal Name:
Mathematical Programming
Additional Journal Information:
Related Information: CHORUS Timestamp: 2017-07-17 14:18:59; Journal ID: ISSN 0025-5610
Springer Science + Business Media
Country of Publication:

Citation Formats

Aragón Artacho, Francisco J., Fleming, Ronan M. T., and Vuong, Phan T.. Accelerating the DC algorithm for smooth functions. Netherlands: N. p., 2017. Web. doi:10.1007/s10107-017-1180-1.
Aragón Artacho, Francisco J., Fleming, Ronan M. T., & Vuong, Phan T.. Accelerating the DC algorithm for smooth functions. Netherlands. doi:10.1007/s10107-017-1180-1.
Aragón Artacho, Francisco J., Fleming, Ronan M. T., and Vuong, Phan T.. 2017. "Accelerating the DC algorithm for smooth functions". Netherlands. doi:10.1007/s10107-017-1180-1.
title = {Accelerating the DC algorithm for smooth functions},
author = {Aragón Artacho, Francisco J. and Fleming, Ronan M. T. and Vuong, Phan T.},
abstractNote = {},
doi = {10.1007/s10107-017-1180-1},
journal = {Mathematical Programming},
number = ,
volume = ,
place = {Netherlands},
year = 2017,
month = 7

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1007/s10107-017-1180-1

Save / Share:
  • We discuss the possibility to represent smooth nonnegative matrix-valued functions as finite linear combinations of fixed matrices with positive real-valued coefficients whose square roots are Lipschitz continuous. This issue is reduced to a similar problem for smooth functions with values in a polyhedron.
  • A criterion describing all functions with finitely many critical points on two-dimensional surfaces that can be height functions corresponding to some immersions of the surface in three-dimensional Euclidean space is established. It is proved that each smooth deformation of a Morse function on the surface can be realized as the deformation of the height function induced by a suitable deformation of the immersion of the surface in R{sup 3}. A new proof of the well-known result on the path connectedness of the space of all smooth immersions of a two-dimensional sphere in R{sup 3} obtained. A new description of anmore » eversion of a two-dimensional sphere in R{sup 3} is given. Generalizations of S. Matveev's result on the connectedness of the space of Morse functions with fixed numbers of minima and maxima on a closed surface are established.« less
  • In this paper an algorithm for determining pressure-derivative functions from measured well-test data is presented. The pressure data are fit with a spline function that can account for curve shape constraints. The first and second derivatives of the resulting spline function are calculated analytically and can be used for type-curve matching. Two example illustrate the performance of the algorithm.
  • When applying the standard {delta}f particle-in-cell simulation method to simulate linear and nonlinear collective instabilities with coherent structures, wave-particle interaction may result in large weight growth for resonant or nearly resonant simulation particles. In this paper, we demonstrate that the large noise associated with the large weight of nearly resonant simulation particles can produce significant error fields at the nonlinear stage of the instability. To overcome this deleterious effect, we have developed a modified {delta}f method that contains a smooth switching algorithm between the {delta}f and total-f methods. Before the switch, the simulation effectively makes use of the desirable low-noisemore » feature of the {delta}f method for small weight to accurately follow unstable mode structures. When the weight function becomes large during the nonlinear phase, the low-noise advantage of the {delta}f method ceases to be significant and the simulation is switched to the total-f method to avoid the large noise induced by nearly resonant simulation particles. This algorithm has been successfully applied to simulation studies of the electrostatic Harris instability driven by large temperature anisotropy in high-intensity charged particle beams typical of applications in high current accelerators, including high-energy density physics and heavy ion fusion.« less
  • The calculation of form factors for the analysis of radiation heat transfer is computationally expensive when surfaces are not in full view of one another. Because of its lower time complexity, the hemi-cube method is efficient for problems with large numbers of surfaces. In this note, methods to improve the hemi-cube algorithm by taking advantage of graphics hardware and two types of geometric coherence are presented. The implementation of the methods on a computer workstation is described, and timings demonstrating speedups by factors of up to 6.7 are given. The acceleration of the calculation of form factors is not anmore » isolated computational problem. The ability to calculate form factors very quickly can alter the overall computational methodology used to analyze radiant heat transfer, and can greatly extend the complexity of geometries that can be studied.« less