Using Improved Directions of Negative Curvature for the Solution of Bound-Constrained Nonconvex Problems

In this work, we describe the efficient use of improved directions of negative curvature for the solution of bound-constrained nonconvex problems. We follow an interior-point framework, in which the key point is the inclusion of computational low-cost procedures to improve directions of negative curvature obtained from a factorisation of the KKT matrix. From a theoretical point of view, it is well known that these directions ensure convergence to second-order KKT points. As a novelty, we consider the convergence rate of the algorithm with exploitation of negative curvature information. Finally, we test the performance of our proposal on both CUTEr/st and simulated problems, showing empirically that the enhanced directions affect positively the practical performance of the procedure.     

Cleft extensions of a family of braided Hopf algebras

We introduce a family of braided Hopf algebras that (in characteristic zero) generalizes the rank 1 Hopf algebras introduced by Krop and Radford and we study its cleft extensions.     

Bergman Subspaces and Subkernels: Degenerate $$L^p$$ Mapping and Zeroes

Regularity and irregularity of the Bergman projection on \(L^p\) spaces is established on a natural family of bounded, pseudoconvex domains. The family is parameterized by a real variable \(\gamma \). A surprising consequence of the analysis is that, whenever \(\gamma \) is irrational, the Bergman projection is bounded only for \(p=2\).     

A fractional step method for 2D parabolic convection-diffusion singularly perturbed problems: uniform convergence and order reduction

In this work, we are concerned with the efficient resolution of two dimensional parabolic singularly perturbed problems of convection-diffusion type. The numerical method combines the fractional implicit Euler method to discretize in time on a uniform mesh and the classical upwind finite difference scheme, defined on a Shishkin mesh, to discretize in space. We consider general time-dependent Dirichlet boundary conditions, and we show that classical evaluations of the boundary conditions cause an order reduction in the consistency of the time integrator. An appropriate correction for the evaluations of the boundary data permits to remove such order reduction. Using this correction, we prove that the fully discrete scheme is uniformly convergent of first order in time and of almost first order in space. Some numerical experiments, which corroborate in practice the robustness and the efficiency of the proposed numerical algorithm, are shown; from them, we bring to light the influence in practice of the two options for the boundary data considered here, which is in agreement with the theoretical results.     

Generalizing the ENO-DB2<Emphasis Type="Italic">p</Emphasis> transform using the inverse wavelet transform

The essentially non-oscillatory (ENO)-wavelet transform developed by Chan and Zhou (SIAM J. Numer. Anal. 40(4), 1369–1404, 2002) is based on a combination of the Daubechies-2p wavelet transform and the ENO technique. It uses extrapolation methods to compute the scaling coefficients without differencing function values across jumps and obtains a multiresolution framework (essentially) free of edge artifacts. In this work, we present a different way to compute the ENO-DB2p wavelet transform of Chan and Zhou which allows us to simplify the process and easily generalize it to other families of orthonormal wavelets.     

ThrEEBoost: Thresholded Boosting for Variable Selection and Prediction via Estimating Equations

Most variable selection techniques for high-dimensional models are designed to be used in settings, where observations are independent and completely observed. At the same time, there is a rich literature on approaches to estimation of low-dimensional parameters in the presence of correlation, missingness, measurement error, selection bias, and other characteristics of real data. In this article, we present ThrEEBoost (Thresholded EEBoost), a general-purpose variable selection technique which can accommodate such problem characteristics by replacing the gradient of the loss by an estimating function. ThrEEBoost generalizes the previously proposed EEBoost algorithm (Wolfson 2011 Wolfson, J. (2011), “EEBoost: A General Method for Prediction and Variable Selection Based on Estimating Equations,” Journal of the American Statistical Association, 106, 296–305.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) by allowing the number of regression coefficients updated at each step to be controlled by a thresholding parameter. Different thresholding parameter values yield different variable selection paths, greatly diversifying the set of models that can be explored; the optimal degree of thresholding can be chosen by cross-validation. ThrEEBoost was evaluated using simulation studies to assess the effects of different threshold values on prediction error, sensitivity, specificity, and the number of iterations to identify minimum prediction error under both sparse and nonsparse true models with correlated continuous outcomes. We show that when the true model is sparse, ThrEEBoost achieves similar prediction error to EEBoost while requiring fewer iterations to locate the set of coefficients yielding the minimum error. When the true model is less sparse, ThrEEBoost has lower prediction error than EEBoost and also finds the point yielding the minimum error more quickly. The technique is illustrated by applying it to the problem of identifying predictors of weight change in a longitudinal nutrition study. Supplementary materials are available online.     

Precision Matrix Estimation With ROPE

It is known that the accuracy of the maximum likelihood-based covariance and precision matrix estimates can be improved by penalized log-likelihood estimation. In this article, we propose a ridge-type operator for the precision matrix estimation, ROPE for short, to maximize a penalized likelihood function where the Frobenius norm is used as the penalty function. We show that there is an explicit closed form representation of a shrinkage estimator for the precision matrix when using a penalized log-likelihood, which is analogous to ridge regression in a regression context. The performance of the proposed method is illustrated by a simulation study and real data applications. Computer code used in the example analyses as well as other supplementary materials for this article are available online.     

On distributional solutions of local and nonlocal problems of porous medium type

We present a theory of well-posedness and a priori estimates for bounded distributional (or very weak) solutions of

(0.1)${?}_{t}u?{\mathfrak{L}}^{\sigma ,\mu }[\phi (u)]=g(x,t)\text{in}{\mathbb{R}}^N\times (0,T),$

where φ is merely continuous and nondecreasing, and ${\mathfrak{L}}^{\sigma ,\mu }$ is the generator of a general symmetric Lévy process. This means that ${\mathfrak{L}}^{\sigma ,\mu }$ can have both local and nonlocal parts like, e.g., ${\mathfrak{L}}^{\sigma ,\mu }=\mathrm{\Delta }?{(?\mathrm{\Delta })}^{\frac{1}{2}}$. New uniqueness results for bounded distributional solutions to this problem and the corresponding elliptic equation are presented and proven. A key role is played by a new Liouville type result for ${\mathfrak{L}}^{\sigma ,\mu }$. Existence and a priori estimates are deduced from a numerical approximation, and energy-type estimates are also obtained.     

Solid-Phase Synthesis of an Insect Pyrokinin Analog Incorporating an Imidazoline Ring as Isosteric Replacement of a trans Peptide Bond

A facile solid-phase synthetic method for incorporating the imidazoline ring motif, a surrogate for a trans peptide bond, into bioactive peptides is reported. The example described is the synthesis of an imidazoline peptidomimetic analog of an insect pyrokinin neuropeptide via a cyclization reaction of an iminium salt generated from the preceding amino acid and 2,4-diaminopropanoic acid (Dap).