David Cohen

David Cohen, Ludwig Gauckler, Ernst Hairer et al.

For trigonometric and modified trigonometric integrators applied to oscillatory Hamiltonian differential equations with one or several constant high frequencies, near-conservation of the total and oscillatory energies are shown over time scales that cover arbitrary negative powers of the step size. This requires non-resonance conditions between the step size and the frequencies, but in contrast to previous results the results do not require any non-resonance conditions among the frequencies. The proof uses modulated Fourier expansions with appropriately modified frequencies.

Rikard Anton, David Cohen, Stig Larsson et al.

A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size restriction as in the often used Störmer-Verlet-leap-frog scheme. Furthermore, it satisfies an almost trace formula (i.e., a linear drift of the expected value of the energy of the problem). Numerical experiments are presented and confirm the theoretical results.

Rikard Anton, David Cohen, Lluis Quer-Sardanyons

A fully discrete approximation of the one-dimensional stochastic heat equation driven by multiplicative space-time white noise is presented. The standard finite difference approximation is used in space and a stochastic exponential method is used for the temporal approximation. Observe that the proposed exponential scheme does not suffer from any kind of CFL-type step size restriction. When the drift term and the diffusion coefficient are assumed to be globally Lipschitz, this explicit time integrator allows for error bounds in $L^q(Ω)$, for all $q\geq2$, improving some existing results in the literature. On top of this, we also prove almost sure convergence of the numerical scheme. In the case of non-globally Lipschitz coefficients, we provide sufficient conditions under which the numerical solution converges in probability to the exact solution. Numerical experiments are presented to illustrate the theoretical results.

David Cohen, Tal Shnitzer, Yuval Kluger et al.

In this paper, we present a new method for few-sample supervised feature selection (FS). Our method first learns the manifold of the feature space of each class using kernels capturing multi-feature associations. Then, based on Riemannian geometry, a composite kernel is computed, extracting the differences between the learned feature associations. Finally, a FS score based on spectral analysis is proposed. Considering multi-feature associations makes our method multivariate by design. This in turn allows for the extraction of the hidden manifold underlying the features and avoids overfitting, facilitating few-sample FS. We showcase the efficacy of our method on illustrative examples and several benchmarks, where our method demonstrates higher accuracy in selecting the informative features compared to competing methods. In addition, we show that our FS leads to improved classification and better generalization when applied to test data.

Rikard Anton, David Cohen

We study an explicit exponential scheme for the time discretisation of stochastic Schrödinger equations driven by additive or multiplicative Ito noise. The numerical scheme is shown to converge with strong order $1$ if the noise is additive and with strong order $1/2$ for multiplicative noise. In addition, if the noise is additive, we show that the exact solutions of our problems satisfy trace formulas for the expected mass, energy, and momentum (i.e., linear drifts in these quantities). Furthermore, we inspect the behaviour of the numerical solutions with respect to these trace formulas. Several numerical simulations are presented and confirm our theoretical results.

David Cohen, Xavier Raynaud

We propose a fully discretised numerical scheme for the hyperelastic rod wave equation on the line. The convergence of the method is established. Moreover, the scheme can handle the blow-up of the derivative which naturally occurs for this equation. By using a time splitting integrator which preserves the invariants of the problem, we can also show that the scheme preserves the positivity of the energy density.

Yuto Miyatake, David Cohen, Daisuke Furihata et al.

We present novel geometric numerical integrators for Hunter--Saxton-like equations by means of new multi-symplectic formulations and known Hamiltonian structures of the problems. We consider the Hunter--Saxton equation, the modified Hunter--Saxton equation, and the two-component Hunter--Saxton equation. Multi-symplectic discretisations based on these new formulations of the problems are exemplified by means of the explicit Euler box scheme, and Hamiltonian-preserving discretisations are exemplified by means of the discrete variational derivative method. We explain and justify the correct treatment of boundary conditions in a unified manner. This is necessary for a proper numerical implementation of these equations and was never explicitly clarified in the literature before, to the best of our knowledge. Finally, numerical experiments demonstrate the favourable behaviour of the proposed numerical integrators.

Chuchu Chen, David Cohen, Jialin Hong

This paper proposes a novel conservative method for numerical computation of general stochastic differential equations in the Stratonovich sense with a conserved quantity. We show that the mean-square order of the method is $1$ if noises are commutative and that the weak order is also $1$. Since the proposed method may need the computation of a deterministic integral, we analyse the effect of the use of quadrature formulas on the convergence orders. Furthermore, based on the splitting technique of stochastic vector fields, we construct conservative composition methods with similar orders as the above method. Finally, numerical experiments are presented to support our theoretical results.