Source author record

David Cohen

David Cohen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

ManiFeSt: Manifold-based Feature Selection for Small Data Sets

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.

preprint2022arXiv

Splitting schemes for FitzHugh--Nagumo stochastic partial differential equations

We design and study splitting integrators for the temporal discretization of the stochastic FitzHugh--Nagumo system. This system is a model for signal propagation in nerve cells where the voltage variable is solution of a one-dimensional parabolic PDE with a cubic nonlinearity driven by additive space-time white noise. We first show that the numerical solutions have finite moments. We then prove that the splitting schemes have, at least, the strong rate of convergence $1/4$. Finally, numerical experiments illustrating the performance of the splitting schemes are provided.

preprint2020arXiv

Exponential integrators for the stochastic Manakov equation

This article presents and analyses an exponential integrator for the stochastic Manakov equation, a system arising in the study of pulse propagation in randomly birefringent optical fibers. We first prove that the strong order of the numerical approximation is $1/2$ if the nonlinear term in the system is globally Lipschitz-continuous. Then, we use this fact to prove that the exponential integrator has convergence order $1/2$ in probability and almost sure order $1/2$, in the case of the cubic nonlinear coupling which is relevant in optical fibers. Finally, we present several numerical experiments in order to support our theoretical findings and to illustrate the efficiency of the exponential integrator as well as a modified version of it.

preprint2020arXiv

High order numerical integrators for single integrand Stratonovich SDEs

We show that applying any deterministic B-series method of order $p_d$ with a random step size to single integrand SDEs gives a numerical method converging in the mean-square and weak sense with order $\lfloor p_d/2\rfloor$.As an application, we derive high order energy-preserving methods for stochastic Poisson systems as well as further geometric numerical schemes for this wide class of Stratonovich SDEs.

preprint2019arXiv

Exponential Integrators for Stochastic Maxwell's Equations Driven by Itô Noise

This article presents explicit exponential integrators for stochastic Maxwell's equations driven by both multiplicative and additive noises. By utilizing the regularity estimate of the mild solution, we first prove that the strong order of the numerical approximation is $\frac 12$ for general multiplicative noise. Combing a proper decomposition with the stochastic Fubini's theorem, the strong order of the proposed scheme is shown to be $1$ for additive noise. Moreover, for linear stochastic Maxwell's equation with additive noise, the proposed time integrator is shown to preserve exactly the symplectic structure, the evolution of the energy as well as the evolution of the divergence in the sense of expectation. Several numerical experiments are presented in order to verify our theoretical findings.

preprint2016arXiv

Cognitive Sub-Nyquist Hardware Prototype of a Collocated MIMO Radar

We present the design and hardware implementation of a radar prototype that demonstrates the principle of a sub-Nyquist collocated multiple-input multiple-output (MIMO) radar. The setup allows sampling in both spatial and spectral domains at rates much lower than dictated by the Nyquist sampling theorem. Our prototype realizes an X-band MIMO radar that can be configured to have a maximum of 8 transmit and 10 receive antenna elements. We use frequency division multiplexing (FDM) to achieve the orthogonality of MIMO waveforms and apply the Xampling framework for signal recovery. The prototype also implements a cognitive transmission scheme where each transmit waveform is restricted to those pre-determined subbands of the full signal bandwidth that the receiver samples and processes. Real-time experiments show reasonable recovery performance while operating as a 4x5 thinned random array wherein the combined spatial and spectral sampling factor reduction is 87.5% of that of a filled 8x10 array.

preprint2016arXiv

Exponential integrators for stochastic Schrödinger equations driven by Ito noise

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.

preprint2016arXiv

On the Existence of Synchrostates in Multichannel EEG Signals during Face-perception Tasks

Phase synchronisation in multichannel EEG is known as the manifestation of functional brain connectivity. Traditional phase synchronisation studies are mostly based on time average synchrony measures hence do not preserve the temporal evolution of the phase difference. Here we propose a new method to show the existence of a small set of unique phase synchronised patterns or "states" in multi-channel EEG recordings, each "state" being stable of the order of ms, from typical and pathological subjects during face perception tasks. The proposed methodology bridges the concepts of EEG microstates and phase synchronisation in time and frequency domain respectively. The analysis is reported for four groups of children including typical, Autism Spectrum Disorder (ASD), low and high anxiety subjects - a total of 44 subjects. In all cases, we observe consistent existence of these states - termed as synchrostates - within specific cognition related frequency bands (beta and gamma bands), though the topographies of these synchrostates differ for different subject groups with different pathological conditions. The inter-synchrostate switching follows a well-defined sequence capturing the underlying inter-electrode phase relation dynamics in stimulus- and person-centric manner. Our study is motivated from the well-known EEG microstate exhibiting stable potential maps over the scalp. However, here we report a similar observation of quasi-stable phase synchronised states in multichannel EEG. The existence of the synchrostates coupled with their unique switching sequence characteristics could be considered as a potentially new field over contemporary EEG phase synchronisation studies.

preprint2016arXiv

SUMMeR: Sub-Nyquist MIMO Radar

Multiple input multiple output (MIMO) radar exhibits several advantages with respect to traditional radar array systems in terms of flexibility and performance. However, MIMO radar poses new challenges for both hardware design and digital processing. In particular, achieving high azimuth resolution requires a large number of transmit and receive antennas. In addition, the digital processing is performed on samples of the received signal, from each transmitter to each receiver, at its Nyquist rate, which can be prohibitively large when high resolution is needed. Overcoming the rate bottleneck, sub-Nyquist sampling methods have been proposed that break the link between radar signal bandwidth and sampling rate. In this work, we extend these methods to MIMO configurations and propose a sub-Nyquist MIMO radar (SUMMeR) system that performs both time and spatial compression. We present a range-azimuth-Doppler recovery algorithm from sub-Nyquist samples obtained from a reduced number of transmitters and receivers, that exploits the sparsity of the recovered targets' parameters. This allows us to achieve reduction in the number of deployed antennas and the number of samples per receiver, without degrading the time and spatial resolutions. Simulations illustrate the detection performance of SUMMeR for different compression levels and shows that both time and spatial resolution are preserved, with respect to classic Nyquist MIMO configurations. We also examine the impact of design parameters, such as antennas' locations and carrier frequencies, on the detection performance, and provide guidelines for their choice.

preprint2016arXiv

X-ray emission from magnetic massive stars

Magnetically confined winds of early-type stars are expected to be sources of bright and hard X-rays. To clarify the systematics of the observed X-ray properties, we have analyzed a large series of Chandra and XMM observations, corresponding to all available exposures of known massive magnetic stars (over 100 exposures covering ~60% of stars compiled in the catalog of Petit et al. 2013). We show that the X-ray luminosity is strongly correlated with the stellar wind mass-loss-rate, with a power-law form that is slightly steeper than linear for the majority of the less luminous, lower-Mdot B stars and flattens for the more luminous, higher-Mdot O stars. As the winds are radiatively driven, these scalings can be equivalently written as relations with the bolometric luminosity. The observed X-ray luminosities, and their trend with mass-loss rates, are well reproduced by new MHD models, although a few overluminous stars (mostly rapidly rotating objects) exist. No relation is found between other X-ray properties (plasma temperature, absorption) and stellar or magnetic parameters, contrary to expectations (e.g. higher temperature for stronger mass-loss rate). This suggests that the main driver for the plasma properties is different from the main determinant of the X-ray luminosity. Finally, variations of the X-ray hardnesses and luminosities, in phase with the stellar rotation period, are detected for some objects and they suggest some temperature stratification to exist in massive stars' magnetospheres.

preprint2015arXiv

Confirming HD 23478 as a new magnetic B star hosting an H$α$-bright centrifugal magnetosphere

In this paper we report 23 magnetic field measurements of the B3IV star HD 23478: 12 obtained from high resolution Stokes $V$ spectra using the ESPaDOnS (CFHT) and Narval (TBL) spectropolarimeters, and 11 from medium resolution Stokes $V$ spectra obtained with the DimaPol spectropolarimeter (DAO). HD 23478 was one of two rapidly rotating stars identified as potential "centrifugal magnetosphere" hosts based on IR observations from the Apache Point Observatory Galactic Evolution Experiment survey. We derive basic physical properties of this star including its mass ($M=6.1^{+0.8}_{-0.7}\,M_\odot$), effective temperature ($T_{\rm eff}=20\pm2\,$kK), radius ($R=2.7^{+1.6}_{-0.9}\,R_\odot$), and age ($τ_{\rm age}=3^{+37}_{-1}\,$Myr). We repeatedly detect weakly-variable Zeeman signatures in metal, He and H lines in all our observations corresponding to a longitudinal magnetic field of $\langle B_z\rangle\approx-2.0\,$kG. The rotational period is inferred from Hipparcos photometry ($P_{\rm rot}=1.0498(4)\,$d). Under the assumption of the Oblique Rotator Model, our obsevations yield a surface dipole magnetic field of strength $B_d\geq9.5\,$kG that is approximately aligned with the stellar rotation axis. We confirm the presence of strong and broad H$α$ emission and gauge the volume of this star's centrifugal magnetosphere to be consistent with those of other H$α$ emitting centrifugal magnetosphere stars based on the large inferred Alfvén to Kepler radius ratio.

preprint2015arXiv

Dilations of matricies

We explore aspects of dilation theory in the finite dimensional case and show that for a commuting $n$-tuple of operators $T=(T_1,...,T_n) $ acting on some finite dimensional Hilbert space $H$ and a compact set $X\subset \mathbb{C}^n$ the following are equivalent: 1. $T$ has a normal $ X$-dilation. 2. For any $m\in \mathbb{N}$ there exists some finite dimensional Hilbert space $K$ containing $H$ and a tuple of commuting normal operators $N=(N_1,...,N_n)$ acting on $K$ such that $$ q(T)=P_Hq(N)|_H$$ for all polynomials $q$ of degree at most $m$ and such that the joint spectrum of $N$ is contained in $X$ (where $P_H$ is the projection from $K$ to $H$).

preprint2015arXiv

Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise

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.

preprint2015arXiv

Multi-symplectic discretisation of wave map equations

We present a new multi-symplectic formulation of constrained Hamiltonian partial differential equations, and we study the associated local conservation laws. A multi-symplectic discretisation based on this new formulation is exemplified by means of the Euler box scheme. When applied to the wave map equation, this numerical scheme is explicit, preserves the constraint and can be seen as a generalisation of the Shake algorithm for constrained mechanical systems. Furthermore, numerical experiments show excellent conservation properties of the numerical solutions.

preprint2014arXiv

Conservative methods for stochastic differential equations with a conserved quantity

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.

preprint2014arXiv

Iterative Plan Construction for the Workflow Satisfiability Problem

The \emph{Workflow Satisfiability Problem (WSP)} is a problem of practical interest that arises whenever tasks need to be performed by authorized users, subject to constraints defined by business rules. We are required to decide whether there exists a \emph{plan} -- an assignment of tasks to authorized users -- such that all constraints are satisfied. Several bespoke algorithms have been constructed for solving the WSP, optimised to deal with constraints (business rules) of particular types. It is natural to see the WSP as a subclass of the {\em Constraint Satisfaction Problem (CSP)} in which the variables are tasks and the domain is the set of users. What makes the WSP distinctive as a CSP is that we can assume that the number of tasks is very small compared to the number of users. This is in sharp contrast with traditional CSP models where the domain is small and the number of variables is very large. As such, it is appropriate to ask for which constraint languages the WSP is fixed-parameter tractable (FPT), parameterized by the number of tasks. We have identified a new FPT constraint language, user-independent constraint, that includes many of the constraints of interest in business processing systems. We are also able to prove that the union of FPT languages remains FPT if they satisfy a simple compatibility condition. In this paper we present our generic algorithm, in which plans are grouped into equivalence classes, each class being associated with a \emph{pattern}. We demonstrate that our generic algorithm has running time $O^*(2^{k\log k})$, where $k$ is the number of tasks, for the language of user-independent constraints. We also show that there is no algorithm of running time $O^*(2^{o(k\log k)})$ for user-independent constraints unless the Exponential Time Hypothesis fails.

preprint2014arXiv

X-rays from Magnetically Confined Wind Shocks: Effect of Cooling-Regulated Shock Retreat

We use 2D MHD simulations to examine the effects of radiative cooling and inverse Compton (IC) cooling on X-ray emission from magnetically confined wind shocks (MCWS) in magnetic massive stars with radiatively driven stellar winds. For the standard dependence of mass loss rate on luminosity $\Mdot \sim L^{1.7} $, the scaling of IC cooling with $L$ and radiative cooling with $\Mdot$ means that IC cooling become formally more important for lower luminosity stars. However, because the sense of the trends is similar, we find the overall effect of including IC cooling is quite modest. More significantly, for stars with high enough mass loss to keep the shocks radiative, the MHD simulations indicate a linear scaling of X-ray luminosity with mass loss rate; but for lower luminosity stars with weak winds, X-ray emission is reduced and softened by a {\em shock retreat} resulting from the larger post-shock cooling length, which within the fixed length of a closed magnetic loop forces the shock back to lower pre-shock wind speeds. A semi-analytic scaling analysis that accounts both for the wind magnetic confinement and this shock retreat yields X-ray luminosities that have a similar scaling trend, but a factor few higher values, compared to time-averages computed from the MHD simulations. The simulation and scaling results here thus provide a good basis for interpreting available X-ray observations from the growing list of massive stars with confirmed large-scale magnetic fields.

preprint2011arXiv

Convergent Numerical Schemes for the Compressible Hyperelastic Rod Wave Equation

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.

preprint2011arXiv

The Physical Basis of the Lx-Lbol Empirical Law for O-star X-rays

X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^-7 Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^1.7 if the shocks are radiative, or with Lx ~ Lx ~ Mdot^2 ~ Lbol^3.4 if they are adiabatic. We present here a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we then propose that the associated softening and weakening of the X-ray emission can be parametrized by the cooling length ratio raised to a power m, the "mixing exponent." For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^0.6 ~ Lbol that matches the empirical scaling. We conclude by noting that such thin-shell mixing may also be important for X-rays from colliding wind binaries, and that future numerical simulation studies will be needed to test this thin-shell mixing ansatz for X-ray emission.