Towards Coding for Human and Machine Vision A Scalable Image Coding Approach. The past decades have witnessed the rapid development of image and video coding techniques in the era of big data. However, the signal fidelity driven coding pipeline design limits the capability of the existing image video coding frameworks to fulfill the needs of both machine and human vision. In this paper, we come up with a novel image coding framework by leveraging both the compressive and the generative models, to support machine vision and human perception tasks jointly. Given an input image, the feature analysis is first applied, and then the generative model is employed to perform image reconstruction with features and additional reference pixels, in which compact edge maps are extracted in this work to connect both kinds of vision in a scalable way. The compact edge map serves as the basic layer for machine vision tasks, and the reference pixels act as a sort of enhanced layer to guarantee signal fidelity for human vision. By introducing advanced generative models, we train a flexible network to reconstruct images from compact feature representations and the reference pixels. Experimental results demonstrate the superiority of our framework in both human visual quality and facial landmark detection, which provide useful evidence on the emerging standardization efforts on MPEG VCM (Video Coding for Machine).


Predictive Accuracy of Wall Modelled Large Eddy Simulation on Unstructured Grids. The predictive accuracy of wall modelled LES is influenced by a combination of the subgrid model, the wall model, the numerical dissipation induced primarily by the convective numerical scheme, and also by the density and topology of the computational grid. The latter factor is of particular importance for industrial flow problems, where unstructured grids are typically employed due to the necessity to handle complex geometries. Here, a systematic simulation based study is presented, investigating the effect of grid cell type on the predictive accuracy of wall modelled LES in the framework of a general purpose finite volume solver. Following standard practice for meshing near wall regions, it is proposed to use prismatic cells. Three candidate shapes for the base of the prisms are considered a triangle, a quadrilateral, and an arbitrary polygon. The cell centre distance is proposed as a metric to determine the spatial resolution of grids with different cell types. The simulation campaign covers two test cases with attached boundary layers fully developed turbulent channel flow, and a zero pressure gradient flat plate turbulent boundary layer. A grid construction strategy is employed, which adapts the grid metric to the outer length scale of the boundary layer. The results are compared with DNS data concerning mean wall shear stress and profiles of flow statistics. The principle outcome is that unstructured simulations may provide the same accuracy as simulations on structured orthogonal hexahedral grids. The choice of base shape of the near wall cells has a significant impact on the computational cost, but in terms of accuracy appears to be a factor of secondary importance.


MuonTrap Preventing Cross Domain Spectre Like Attacks by Capturing Speculative State. The disclosure of the Spectre speculative execution attacks in January 2018 has left a severe vulnerability that systems are still struggling with how to patch. The solutions that currently exist tend to have incomplete coverage, perform badly, or have highly undesirable edge cases that cause application domains to break. MuonTrap allows processors to continue to speculate, avoiding significant reductions in performance, without impacting security. We instead prevent the propagation of any state based on speculative execution, by placing the results of speculative cache accesses into a small, fast L0 filter cache, that is non inclusive, non exclusive with the rest of the cache hierarchy. This isolates all parts of the system that can't be quickly cleared on any change in threat domain. MuonTrap uses these speculative filter caches, which are cleared on context and protection domain switches, along with a series of extensions to the cache coherence protocol and prefetcher. This renders systems immune to cross domain information leakage via Spectre and a host of similar attacks based on speculative execution, with low performance impact and few changes to the CPU design.


Existence of radial bounded solutions for some quasilinear elliptic equations in R N. We study the quasilinear equation (P) qquad rm div (A(x,u) nabla u p 2 nabla u) frac1p A t(x,u) nabla u p u p 2 u g(x,u) qquad hbox in mathbb R N , with N ge 3 , p 1 , where A(x,t) , A t(x,t) frac partial A partial t (x,t) and g(x,t) are Carath 'eodory functions on mathbb R N times mathbb R . Suitable assumptions on A(x,t) and g(x,t) set off the variational structure of (P) and its related functional cal J is C 1 on the Banach space X W 1,p ( mathbb R N) cap L infty( mathbb R N) . In order to overcome the lack of compactness, we assume that the problem has radial symmetry, then we look for critical points of cal J restricted to X r , subspace of the radial functions in X . Following an approach which exploits the interaction between cdot X and the norm on W 1,p ( mathbb R N) , we prove the existence of at least one weak bounded radial solution of (P) by applying a generalized version of the Ambrosetti Rabinowitz Mountain Pass Theorem.


KAM tori for the generalized Bejamin Bona Mahony equation. A generalized Benjamin Bona Mahony (gBBM) equation subject to the periodic boundary condition is studied in this paper. Based on a new infinite dimensional Kolomogorov Arnold Moser (KAM) theorem with normal frequencies of finite limit points, it is shown that the gBBM equation admits plenty of time quasi periodic solutions with two frequencies of high modes.


Generating and Detecting High Frequency Liquid Based Sound Resonances with Nanoplasmonics. We use metal nanostructures (nanoplasmonics) excited with dual frequency lasers to generate and detect high frequency ( 10 GHz) sound wave resonances in water. The difference frequency between the two lasers causes beating in the intensity, which results in a drop in the transmission through the nanostructure when an acoustic resonance is excited. By observing the resonance frequency shifts with changing nanostructure size, the transition from slow to fast sound in water is inferred, which has been measured by inelastic scattering methods in the past. The observed behavior shows remarkable similarities to a simple Debye model (without fitting parameters). The ability to directly excite high frequency sound waves in water may unlock the secret of how the nanofluidic environment, that is typically considered to be extremely viscous, can efficiently support the energetic dynamics of life via protein vibrations at the nanometer scale.


Extremals for a series of sub Finsler problems with 2 dimensional control via convex trigonometry. We consider a series of optimal control problems with 2 dimensional control lying in an arbitrary convex compact set Omega . The considered problems are well studied for the case when Omega is a unit disc, but barely studied for arbitrary Omega . We derive extremals to these problems in general case by using machinery of convex trigonometry, which allows us to do this identically and independently on the shape of Omega . The paper describes geodesics in (i) the Finsler problem on the Lobachevsky hyperbolic plane (ii) left invariant sub Finsler problems on all unimodular 3D Lie groups (SU(2), SL(2), SE(2), SH(2)) (iii) the problem of rolling ball on a plane with distance function given by Omega (iv) a series of "yacht problems" generalizing Euler's elastic problem, Markov Dubins problem, Reeds Shepp problem and a new sub Riemannian problem on SE(2) and (v) the plane dynamic motion problem.


Well posedness for good Boussinesq equations subject to quasi periodic initial data. This paper concerns the local well posedness for the "good" Boussinesq equation subject to quasi periodic initial conditions. By constructing a delicately and subtly iterative process together with an explicit combinatorial analysis, we show that there exists a unique solution for such a model in a small region of time. The size of this region depends on both the given data and the frequency vector involved. Moreover the local solution has an expansion with exponentially decaying Fourier coefficients.


Hyperplane Arrangements of Trained ConvNets Are Biased. We investigate the geometric properties of the functions learned by trained ConvNets in the preactivation space of their convolutional layers, by performing an empirical study of hyperplane arrangements induced by a convolutional layer. We introduce statistics over the weights of a trained network to study local arrangements and relate them to the training dynamics. We observe that trained ConvNets show a significant statistical bias towards regular hyperplane configurations. Furthermore, we find that layers showing biased configurations are critical to validation performance for the architectures considered, trained on CIFAR10, CIFAR100 and ImageNet.


Representing Point Clouds with Generative Conditional Invertible Flow Networks. In this paper, we propose a simple yet effective method to represent point clouds as sets of samples drawn from a cloud specific probability distribution. This interpretation matches intrinsic characteristics of point clouds the number of points and their ordering within a cloud is not important as all points are drawn from the proximity of the object boundary. We postulate to represent each cloud as a parameterized probability distribution defined by a generative neural network. Once trained, such a model provides a natural framework for point cloud manipulation operations, such as aligning a new cloud into a default spatial orientation. To exploit similarities between same class objects and to improve model performance, we turn to weight sharing networks that model densities of points belonging to objects in the same family share all parameters with the exception of a small, object specific embedding vector. We show that these embedding vectors capture semantic relationships between objects. Our method leverages generative invertible flow networks to learn embeddings as well as to generate point clouds. Thanks to this formulation and contrary to similar approaches, we are able to train our model in an end to end fashion. As a result, our model offers competitive or superior quantitative results on benchmark datasets, while enabling unprecedented capabilities to perform cloud manipulation tasks, such as point cloud registration and regeneration, by a generative network.


Searching and Sorting with O(n 2) processors in O(1) time. The proliferation of number of processing elements (PEs) in parallel computer systems, along with the use of more extensive parallelization of algorithms causes the interprocessor communications dominate VLSI chip space. This paper proposes a new architecture to overcome this issue by using simple crosspoint switches to pair PEs instead of a complex interconnection network. Based on the cyclic permutation wiring idea described in cite oruc2016self , this pairing leads to a linear crosspoint array of n(n 1) 2 processing elements and as many crosspoints. We demonstrate the versatility of this new parallel architecture by designing fast searching and sorting algorithms for it. In particular, we show that finding a minimum, maximum, and searching a list of n elements can all be performed in O(1) time with elementary logic gates with O(n) fan in, and in O( lg n) time with O(1) fan in. We further show that sorting a list of n elements can also be carried out in O(1) time using elementary logic gates with O(n) fan in and threshold logic gates. The sorting time increases to O( lg n lg lg n) if only elementary logic gates with O(1) fan in are used. The algorithm can find the maximum among n elements in O(1) time, and sort n elements in O( lg n ( lg lg n)) time. In addition, we show how other fundamental queries can be handled within the same order of time complexities.


Learning Local Forward Models on Unforgiving Games. This paper examines learning approaches for forward models based on local cell transition functions. We provide a formal definition of local forward models for which we propose two basic learning approaches. Our analysis is based on the game Sokoban, where a wrong action can lead to an unsolvable game state. Therefore, an accurate prediction of an action's resulting state is necessary to avoid this scenario. In contrast to learning the complete state transition function, local forward models allow extracting multiple training examples from a single state transition. In this way, the Hash Set model, as well as the Decision Tree model, quickly learn to predict upcoming state transitions of both the training and the test set. Applying the model using a statistical forward planner showed that the best models can be used to satisfying degree even in cases in which the test levels have not yet been seen. Our evaluation includes an analysis of various local neighbourhood patterns and sizes to test the learners' capabilities in case too few or too many attributes are extracted, of which the latter has shown do degrade the performance of the model learner.


Change detection in SAR time series based on the coefficient of variation. This paper discusses change detection in SAR time series. Firstly, several statistical properties of the coefficient of variation highlight its pertinence for change detection. Then several criteria are proposed. The coefficient of variation is suggested to detect any kind of change. Then other criteria based on ratios of coefficients of variations are proposed to detect long events such as construction test sites, or point event such as vehicles. These detection methods are evaluated first on theoretical statistical simulations to determine the scenarios where they can deliver the best results. Then detection performance is assessed on real data for different types of scenes and sensors (Sentinel 1, UAVSAR). In particular, a quantitative evaluation is performed with a comparison of our solutions with state of the art methods.


Models of dynamic damage and phase field fracture, and their various time discretisations. Several variants of models of damage in viscoelastic continua under small strains in the Kelvin Voigt rheology are presented and analyzed by using the Galerkin method. The particular case, known as a phase field fracture approximation of cracks, is discussed in detail. All these models are dynamic (i.e. involve inertia to model vibrations or waves possibly emitted during fast damage fracture or induced by fast varying forcing) and consider viscosity which is also damageable. Then various options for time discretisation are devised. Eventually, extensions to more complex rheologies or a modification for large strains are briefly exposed, too.


Comparing the Roller and B(X) boundaries of CAT(0) cube complexes. The Roller boundary is a well known compactification of a CAT(0) cube complex X. When X is locally finite, essential, irreducible, non Euclidean and admits a cocompact action by a group G, Nevo Sageev show that a subset, B(X), of the Roller boundary is the realization of the Poisson boundary and that the action of G on B(X) is minimal and strongly proximal. Additionally, these authors show B(X) satisfies many other desirable dynamical and topological properties. In this article we give several equivalent characterizations for when B(X) is equal to the entire Roller boundary. As an application we show, under mild hypotheses, that if X is also 2 dimensional then X is G equivariantly quasi isometric to a CAT(0) cube complex X' whose Roller boundary is equal to B(X'). Additionally, we use our characterization to show that the usual CAT(0) cube complex for which an infinite right angled Coxeter Artin group acts on geometrically has Roller boundary equal to B(X), as long as the corresponding group does not decompose as a direct product.


Fast diffusion on noncompact manifolds well posedness theory and connections with semilinear elliptic equations. We investigate the well posedness of the fast diffusion equation (FDE) in a wide class of noncompact Riemannian manifolds. Existence and uniqueness of solutions for globally integrable initial data was established in 5 . However, in the Euclidean space, it is known from Herrero and Pierre 20 that the Cauchy problem associated with the FDE is well posed for initial data that are merely in L 1 mathrm loc . We establish here that such data still give rise to global solutions on general Riemannian manifolds. If, in addition, the radial Ricci curvature satisfies a suitable pointwise bound from below (possibly diverging to infty at spatial infinity), we prove that also uniqueness holds, for the same type of data, in the class of strong solutions. Besides, under the further assumption that the initial datum is in L 2 mathrm loc and nonnegative, a minimal solution is shown to exist, and we are able to establish uniqueness of purely (nonnegative) distributional solutions, which to our knowledge was not known before even in the Euclidean space. The required curvature bound is in fact sharp, since on model manifolds it turns out to be equivalent to stochastic completeness, and it was shown in 13 that uniqueness for the FDE fails even in the class of bounded solutions on manifolds that are not stochastically complete. Qualitatively this amounts to asking that the curvature diverges at most quadratically at infinity. A crucial ingredient of the uniqueness result is the proof of nonexistence of distributional subsolutions to certain semilinear elliptic equations with power nonlinearities, of independent interest.


Local limit theorems for occupancy models. We present a rather general method for proving local limit theorems, with a good rate of convergence, for sums of dependent random variables. The method is applicable when a Stein coupling can be exhibited. Our approach involves both Stein's method for distributional approximation and Stein's method for concentration. As applications, we prove local central limit theorems with rate of convergence for the number of germs with d neighbours in a germ grain model, and the number of degree d vertices in an Erd H o s R 'enyi random graph. In both cases, the error rate is optimal, up to logarithmic factors.


Central extensions of filiform Zinbiel algebras. In this paper we describe central extensions (up to isomorphism) of all complex null filiform and filiform Zinbiel algebras. It is proven that every non split central extension of an n dimensional null filiform Zinbiel algebra is isomorphic to an (n 1) dimensional null filiform Zinbiel algebra. Moreover, we obtain all pairwise non isomorphic quasi filiform Zinbiel algebras.


Small Footprint Open Vocabulary Keyword Spotting with Quantized LSTM Networks. We explore a keyword based spoken language understanding system, in which the intent of the user can directly be derived from the detection of a sequence of keywords in the query. In this paper, we focus on an open vocabulary keyword spotting method, allowing the user to define their own keywords without having to retrain the whole model. We describe the different design choices leading to a fast and small footprint system, able to run on tiny devices, for any arbitrary set of user defined keywords, without training data specific to those keywords. The model, based on a quantized long short term memory (LSTM) neural network, trained with connectionist temporal classification (CTC), weighs less than 500KB. Our approach takes advantage of some properties of the predictions of CTC trained networks to calibrate the confidence scores and implement a fast detection algorithm. The proposed system outperforms a standard keyword filler model approach.


CDL Curriculum Dual Learning for Emotion Controllable Response Generation. Emotion controllable response generation is an attractive and valuable task that aims to make open domain conversations more empathetic and engaging. Existing methods mainly enhance the emotion expression by adding regularization terms to standard cross entropy loss and thus influence the training process. However, due to the lack of further consideration of content consistency, the common problem of response generation tasks, safe response, is intensified. Besides, query emotions that can help model the relationship between query and response are simply ignored in previous models, which would further hurt the coherence. To alleviate these problems, we propose a novel framework named Curriculum Dual Learning (CDL) which extends the emotion controllable response generation to a dual task to generate emotional responses and emotional queries alternatively. CDL utilizes two rewards focusing on emotion and content to improve the duality. Additionally, it applies curriculum learning to gradually generate high quality responses based on the difficulties of expressing various emotions. Experimental results show that CDL significantly outperforms the baselines in terms of coherence, diversity, and relation to emotion factors.


Inquisitive Question Generation for High Level Text Comprehension. Inquisitive probing questions come naturally to humans in a variety of settings, but is a challenging task for automatic systems. One natural type of question to ask tries to fill a gap in knowledge during text comprehension, like reading a news article we might ask about background information, deeper reasons behind things occurring, or more. Despite recent progress with data driven approaches, generating such questions is beyond the range of models trained on existing datasets. We introduce INQUISITIVE, a dataset of 19K questions that are elicited while a person is reading through a document. Compared to existing datasets, INQUISITIVE questions target more towards high level (semantic and discourse) comprehension of text. We show that readers engage in a series of pragmatic strategies to seek information. Finally, we evaluate question generation models based on GPT 2 and show that our model is able to generate reasonable questions although the task is challenging, and highlight the importance of context to generate INQUISITIVE questions.


Zero noise limit for multidimensional SDEs driven by a pointy gradient. The purpose of the article is to address the limiting behavior of the solutions of stochastic differential equations driven by a pointy d dimensional gradient as the intensity of the underlying Brownian motion tends to 0 . By pointy gradient, we here mean that the drift derives from a potential that is mathcal C 1,1 on any compact subset that does not contain the origin. As a matter of fact, the corresponding deterministic version of the differential equation may have an infinite number of solutions when initialized from 0 mathbb R d , in which case the limit theorem proved in the paper reads as a selection theorem of the solutions to the zero noise system. Generally speaking, our result says that, under suitable conditions, the probability that the particle leaves the origin by going through regions of higher potential tends to 1 as the intensity of the noise tends to 0 . In particular, our result extends the earlier one due to Bafico and Baldi for the zero noise limit of one dimensional stochastic differential equations.


DIII clan combinatorics for the orthogonal Grassmannian. Borel subgroup orbits of the classical symmetric space SO 2n GL n are parametrized by DIII (n,n) clans. We group the clans into "sects" corresponding to Schubert cells of the orthogonal Grassmannian, thus providing a cell decomposition for SO 2n GL n . We also compute a recurrence for the rank polynomial of the weak order poset on DIII clans, and then describe explicit bijections between such clans, diagonally symmetric rook placements, certain pairs of minimally intersecting set partitions, and a class of weighted Delannoy paths. Clans of the largest sect are in bijection with fixed point free partial involutions.


Coherent structures and secondary flow in turbulent square duct. The aim of the present work is to investigate the role of coherent structures in the generation of the secondary flow in a turbulent square duct. The coherent structures are defined as connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long time turbulence statistics, in the spirit of the three dimensional quadrant analysis proposed by Lozano Dur 'an et al. (J. Fluid Mech., vol. 694, 2012, pp. 100 130). We consider both the direct contribution of the structures to the mean in plane velocity components and their geometrical properties. The instantaneous phenomena taking place in the turbulent duct are compared with turbulent channel flow at Reynolds numbers of Re tau 180 and 360 , based on friction velocity at the center plane and channel half height. In the core region of the duct, the fractional contribution of intense events to the wall normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three dimensional objects do not differ significantly in both flows. On the other hand, in the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity, which is less relevant than that of the complementary portion of the flow not included in such objects. Our results show however that strong Reynolds shear stress events, despite the differences observed between channel and duct, do not contribute directly to the secondary motion, and thus other phenomena need to be considered instead.


Augmenting reality On the shared history of perceptual illusion and video projection mapping. Perceptual illusions based on the spatial correspondence between objects and displayed images have been pursued by artists and scientists since the 15th century, mastering optics to create crucial techniques as the linear perspective and devices as the Magic Lantern. Contemporary video projection mapping inherits and further extends this drive to produce perceptual illusions in space by incorporating the required real time capabilities for dynamically superposing the imaginary onto physical objects under fluid real world conditions. A critical milestone has been reached in the creation of the technical possibilities for all encompassing, untethered synthetic reality experiences available to the plain senses, where every surface may act as a screen and the relation to everyday objects is open to alterations.


Simple and robust speckle detection method for fire and heat detection in harsh environments. Standard laser based fire detection systems are often based on measuring variation of optical signal amplitude. However, mechanical noise interference and loss from dust and steam can obscure the detection signal, resulting in faulty results or inability to detect a potential fire. The presented fire detection technology will allow the detection of fire in harsh and dusty areas, which are prone to fires, where current systems show limited performance or are unable to operate. It is not the amount of light nor its wavelength that is used for detecting fire, but how the refractive index randomly fluctuates due to the heat convection from the fire. In practical terms this means that light obstruction from ambient dust particles will not be a problem as long as a small fraction of the light is detected and that fires without visible flames can still be detected. The standalone laser system consists of a Linux based Red Pitaya system, a cheap 650 nm laser diode, and a PIN photo detector. Laser light propagates through the monitored area and reflects off a retroreflector generating a speckle pattern. Every 3 seconds time traces and frequency noise spectra are measured and 8 descriptors are deduced to identify a potential fire. Both laboratory and factory acceptance tests have been performed with success.


Validation and noise robustness assessment of microscopic anisotropy estimation with clinically feasible double diffusion encoding MRI. Purpose Double diffusion encoding (DDE) MRI enables the estimation of microscopic diffusion anisotropy, yielding valuable information on tissue microstructure. A recent study proposed that the acquisition of rotationally invariant DDE metrics, typically obtained using a spherical "5 design", could be greatly simplified by assuming Gaussian diffusion, facilitating reduced acquisition times that are more compatible with clinical settings. Here, we aim to validate the new minimal acquisition scheme against the standard DDE 5 design, and to quantify the proposed method's noise robustness to facilitate future clinical use. Methods DDE MRI experiments were performed on both ex vivo and in vivo rat brains at 9.4 T using the 5 design and the proposed minimal design and taking into account the difference in the number of acquisitions. The ensuing microscopic fractional anisotropy ( mu FA) maps were compared over a range of b values up to 5000 s mm2. Noise robustness was studied using analytical calculations and numerical simulations. Results The minimal protocol quantified mu FA at an accuracy comparable to the estimates obtained via the more theoretically robust DDE 5 design. mu FA's sensitivity to noise was found to strongly depend on compartment anisotropy and tensor magnitude in a non linear fashion. When mu FA 0.75 or when mean diffusivity is particularly low, very high signal to noise ratio (SNR) is required for precise quantification of mu FA. Conclusion Our work supports using DDE for quantifying microscopic diffusion anisotropy in clinical settings but raises hitherto overlooked precision issues when measuring mu FA with DDE and typical clinical SNR.


Role of rotational coherence in femtosecond pulse driven nitrogen ion lasing. We experimentally investigated the rotationally resolved polarization characteristics of N 2 lasing at 391 and 428 nm using a pump seed scheme. By varying the relative angle between the linear polarizations of the pump and seed, it is found that the polarizations of the P and R branches of 391 nm lasing are counter rotated. By contrast, both branches of 428 nm lasing remain polarized along the pump. The origin of the puzzled abnormal polarization characteristics is found based on a complete physical model that simultaneously includes the transient photoionization and the subsequent coupling among the electronic, vibrational and rotational quantum states of ions.It suggests that the cascaded resonant Raman processes following ionization create negative coherence between the rotational states of J and J 2 in the ionic ground state X 2 Sigma g ( nu 0) , which leads to mirror symmetrical polarization for the P and R branches of 391 nm lasing. Both the experiment and theory indicate that the demonstrated rotational coherence plays an extremely pivotal role in clarifying the gain mechanism of N 2 lasing and opens up the route toward quantum optics under ultrafast strong fields.


New Ideas and Trends in Deep Multimodal Content Understanding A Review. The focus of this survey is on the analysis of two modalities of multimodal deep learning image and text. Unlike classic reviews of deep learning where monomodal image classifiers such as VGG, ResNet and Inception module are central topics, this paper will examine recent multimodal deep models and structures, including auto encoders, generative adversarial nets and their variants. These models go beyond the simple image classifiers in which they can do uni directional (e.g. image captioning, image generation) and bi directional (e.g. cross modal retrieval, visual question answering) multimodal tasks. Besides, we analyze two aspects of the challenge in terms of better content understanding in deep multimodal applications. We then introduce current ideas and trends in deep multimodal feature learning, such as feature embedding approaches and objective function design, which are crucial in overcoming the aforementioned challenges. Finally, we include several promising directions for future research.


Curvature homogeneous manifolds in dimension 4. We classify complete curvature homogeneous metrics on simply connected four dimensional manifolds which are invariant under a cohomogeneity one action. We show that they are either isometric to a symmetric space with one of its cohomogeneity one actions, or to a complete example by Tsukada on the normal bundle of the Veronese surface in CP 2. Along the way we show (in any dimension) that via an equivariant diffeomorphism the functions describing the metric can be partially diagonalized, a fact that may be useful for other problems as well


Exponential Lower Bounds for Batch Reinforcement Learning Batch RL can be Exponentially Harder than Online RL. Several practical applications of reinforcement learning involve an agent learning from past data without the possibility of further exploration. Often these applications require us to 1) identify a near optimal policy or to 2) estimate the value of a target policy. For both tasks we derive emph exponential information theoretic lower bounds in discounted infinite horizon MDPs with a linear function representation for the action value function even if 1) emph realizability holds, 2) the batch algorithm observes the exact reward and transition emph functions , and 3) the batch algorithm is given the emph best a priori data distribution for the problem class. Our work introduces a new oracle batch algorithm' framework to prove lower bounds that hold for every distribution. The work shows an exponential separation between batch and online reinforcement learning.


Hardware Implementation of Adaptive Watermarking Based on Local Spatial Disorder Analysis. With the increasing use of the internet and the ease of exchange of multimedia content, the protection of ownership rights has become a significant concern. Watermarking is an efficient means for this purpose. In many applications, real time watermarking is required, which demands hardware implementation of low complexity and robust algorithm. In this paper, an adaptive watermarking is presented, which uses embedding in different bit planes to achieve transparency and robustness. Local disorder of pixels is analyzed to control the strength of the watermark. A new low complexity method for disorder analysis is proposed, and its hardware implantation is presented. An embedding method is proposed, which causes lower degradation in the watermarked image. Also, the performance of proposed watermarking architecture is improved by a pipe line structure and is tested on an FPGA device. Results show that the algorithm produces transparent and robust watermarked images. The synthesis report from FPGA implementation illustrates a low complexity hardware structure.


Validation Frameworks for Self Driving Vehicles A Survey. As a part of the digital transformation, we interact with more and more intelligent gadgets. Today, these gadgets are often mobile devices, but in the advent of smart cities, more and more infrastructure such as traffic and buildings in our surroundings becomes intelligent. The intelligence, however, does not emerge by itself. Instead, we need both design techniques to create intelligent systems, as well as approaches to validate their correct behavior. An example of intelligent systems that could benefit smart cities are self driving vehicles. Self driving vehicles are continuously becoming both commercially available and common on roads. Accidents involving self driving vehicles, however, have raised concerns about their reliability. Due to these concerns, the safety of self driving vehicles should be thoroughly tested before they can be released into traffic. To ensure that self driving vehicles encounter all possible scenarios, several millions of hours of testing must be carried out therefore, testing self driving vehicles in the real world is impractical. There is also the issue that testing self driving vehicles directly in the traffic poses a potential safety hazard to human drivers. To tackle this challenge, validation frameworks for testing self driving vehicles in simulated scenarios are being developed by academia and industry. In this chapter, we briefly introduce self driving vehicles and give an overview of validation frameworks for testing them in a simulated environment. We conclude by discussing what an ideal validation framework at the state of the art should be and what could benefit validation frameworks for self driving vehicles in the future.


Associating Uncertainty to Extended Poses for on Lie Group IMU Preintegration with Rotating Earth. The recently introduced matrix group SE2(3) provides a 5x5 matrix representation for the orientation, velocity and position of an object in the 3 D space, a triplet we call "extended pose". In this paper we build on this group to develop a theory to associate uncertainty with extended poses represented by 5x5 matrices. Our approach is particularly suited to describe how uncertainty propagates when the extended pose represents the state of an Inertial Measurement Unit (IMU). In particular it allows revisiting the theory of IMU preintegration on manifold and reaching a further theoretic level in this field. Exact preintegration formulas that account for rotating Earth, that is, centrifugal force and Coriolis force, are derived as a byproduct, and the factors are shown to be more accurate. The approach is validated through extensive simulations and applied to sensor fusion where a loosely coupled fixed lag smoother fuses IMU and LiDAR on one hour long experiments using our experimental car. It shows how handling rotating Earth may be beneficial for long term navigation within incremental smoothing algorithms.


Theory of reflectionless scattering modes. We develop the theory of a special type of scattering state in which a set of asymptotic channels are chosen as inputs and the complementary set as outputs, and there is zero reflection back into the input channels. In general an infinite number of such solutions exist at discrete complex frequencies. Our results apply to linear electromagnetic and acoustic wave scattering and also to quantum scattering, in all dimensions, for arbitrary geometries including scatterers in free space, and for any choice of the input output sets. We refer to such a state as reflection zero (R zero) when it occurs off the real frequency axis and as an Reflectionless Scattering Mode (RSM) when it is tuned to a real frequency as a steady state solution. Such reflectionless behavior requires a specific monochromatic input wavefront, given by the eigenvector of a filtered scattering matrix with eigenvalue zero. Steady state RSMs may be realized by index tuning which do not break flux conservation or by gain loss tuning. RSMs of flux conserving cavities are bidirectional while those of non flux conserving cavities are generically unidirectional. Cavities with cal PT symmetry have unidirectional R zeros in complex conjugate pairs, implying that reflectionless states naturally arise at real frequencies for small gain loss parameter but move into the complex frequency plane after a spontaneous cal PT breaking transition. Numerical examples of RSMs are given for one dimensional cavities with and without gain loss, a cal PT cavity, a two dimensional multiwaveguide junction, and a two dimensional deformed dielectric cavity in free space. We outline and implement a general technique for solving such problems, which shows promise for designing photonic structures which are perfectly impedance matched for specific inputs, or can perfectly convert inputs from one set of modes to a complementary set.


Periodicity in fields of elongating dunes. Dune fields are commonly associated with periodic patterns that are among the most recognizable landscapes on Earth and other planetary bodies. However, in zones of limited sediment supply, where periodic dunes elongate and align in the direction of the resultant sand flux, there has been no attempt to explain the emergence of such a regular pattern. Here, we show, by means of numerical simulations, that the elongation growth mechanism does not produce a pattern with a specific wavelength. Periodic elongating dunes appear to be a juxtaposition of individual structures, the arrangement of which is due to regular landforms at the border of the field acting as boundary conditions. This includes, among others, dune patterns resulting from bed instability, or the crestline reorganization induced by dune migration. The wavelength selection in fields of elongating dunes therefore reflects the interdependence of dune patterns over the course of their evolution.


The Hessian Discriminant. We express the Hessian discriminant of a cubic surface in terms of fundamental invariants. This answers Question 15 from the emph 27 questions on the cubic surface . We also explain how to compute the fundamental invariants for smooth cubics of rank 6.


Alternating links with totally geodesic checkerboard surfaces. We prove that alternating links with two totally geodesic checkerboard surfaces are three links with projection the 1 skeleton of the octahedron, the cuboctahedron and the icosidodecahedron. We also characterize these links as right angled completely realisable links and show that all hyperbolic weaving knots with two exceptions have both checkerboard surfaces not totally geodesic.


Kinetic diffusion asymptotic preserving Monte Carlo algorithm for Boltzmann BGK in the diffusive scaling. We develop a novel Monte Carlo strategy for the simulation of the Boltzmann BGK model with both low collisional and high collisional regimes present. The presented solution to maintain accuracy in low collisional regimes and remove exploding simulation costs in high collisional regimes uses hybridized particles that exhibit both kinetic behaviour and diffusive behaviour depending on the local collisionality. In this work, we develop such a method that maintains the correct mean, variance, and correlation of the positional increments over multiple time steps of fixed step size for all values of the collisionality, under the condition of spatial homogeneity during the time step. In the low collisional regime, the method reverts to the standard velocity jump process. In the high collisional regime, the method collapses to a standard random walk process. We analyze the error of the presented scheme in the low collisional regime for which we obtain the order of convergence in the time step size. We furthermore provide an analysis in the high collisional regime that demonstrates the asymptotic preserving property.


Refining the Uluzzian through a new lithic assemblage from Roccia San Sebastiano (Mondragone, southern Italy). Roccia San Sebastiano is a tectonic karstic cave located at the foot of the southern slope of Mt. Massico, in the territory of Mondragone (Caserta) in Campania (southern Italy). Systematic excavation has been carried out since 2001, leading to the partial exploration of an important Pleistocene deposit, extraordinarily rich in lithic and faunal remains. The aim of this paper is to (1) present the stratigraphic sequence of Roccia San Sebastiano, and (2) technologically describe the lithic materials of squares F14 t18, t19, t20 E16 t16, t17, t18 recently recognised as Uluzzian. The stratigraphic sequence is more than 3 metres thick and dates from the Middle to the Upper Palaeolithic. It contains different techno complexes Gravettian, Aurignacian, Uluzzian and Mousterian. In the Uluzzian lithic assemblage mostly local pebbles of chert were used in order to produce small sized objects. The concept of debitage mainly deals with unidirectional debitage with absent or fairly accurate management of the convexities and angles the striking platforms are usually natural or made by one stroke. It is attested the use of both direct freehand percussion and bipolar technique on anvil in the same reduction sequence. Amongst the retouched tools the presence of two lunates is of note. This study of the Roccia San Sebastiano Uluzzian lithic complexes is significant for understanding the dynamics of the transition from Middle to Upper Palaeolithic in the Tyrrhenian margin of southern Italy.


Cumulative Message Authentication Codes for Resource Constrained IoT Networks. In resource constrained IoT networks, the use of conventional message authentication codes (MACs) to provide message authentication and integrity is not possible due to the large size of the MAC output. A straightforward yet naive solution to this problem is to employ a truncated MAC which undesirably sacrifices cryptographic strength in exchange for reduced communication overhead. In this paper, we address this problem by proposing a novel approach for message authentication called textit Cumulative Message Authentication Code (CuMAC), which consists of two distinctive procedures textit aggregation and textit accumulation . In aggregation, a sender generates compact authentication tags from segments of multiple MACs by using a systematic encoding procedure. In accumulation, a receiver accumulates the cryptographic strength of the underlying MAC by collecting and verifying the authentication tags. Embodied with these two procedures, CuMAC enables the receiver to achieve an advantageous trade off between the cryptographic strength and the latency in processing of the authentication tags. Furthermore, for some latency sensitive messages where this trade off may be unacceptable, we propose a variant of CuMAC that we refer to as textit CuMAC with Speculation (CuMAC S). In addition to the aggregation and accumulation procedures, CuMAC S enables the sender and receiver to employ a speculation procedure for predicting future message values and pre computing the corresponding MAC segments. For the messages which can be reliably speculated, CuMAC S significantly reduces the MAC verification latency without compromising the cryptographic strength. We have carried out comprehensive evaluation of CuMAC and CuMAC S through simulation and a prototype implementation on a real car.


Notes on constructions of knots with the same trace. The m trace of a knot is the 4 manifold obtained from mathbf B 4 by attaching a 2 handle along the knot with m framing. In 2015, Abe, Jong, Luecke and Osoinach introduced a technique to construct infinitely many knots with the same m trace, which is called the operation ( ast m) . In this paper, we prove that their technique can be explained in terms of Gompf and Miyazaki's dualizable pattern. In addition, we show that the family of knots admitting the same 4 surgery given by Teragaito can be explained by the operation ( ast m) .


Envelope of intermediate lines of a plane curve. For a pair of points in a smooth closed convex planar curve gamma , its mid line is the line containing its mid point and the intersection point of the corresponding pair of tangent lines. It is well known that the envelope of the mid lines ( EML ) is formed by the union of three affine invariants sets Affine Envelope Symmetry Sets ( AESS ) Mid Parallel Tangent Locus ( MPTL ) and Affine Evolute of gamma . In this paper, we generalized these concepts by considering the envelope of the intermediate lines. For a pair of points of gamma , its intermediate line is the line containing an intermediate point and the intersection point of the corresponding pair of tangent lines. Here, we present the envelope of intermediate lines ( EIL ) of the curve gamma and prove that this set is formed by three disconnected sets when the intermediate point is different from the mid point Affine Envelope of Intermediate Lines ( AEIL ) the curve gamma itself and the Intermediate Parallel Tangent Locus ( IPTL ). When the intermediate point coincides with the mid point, the EIL coincides with the EML , and thus these sets are connected. Moreover, we introduce some standard techniques of singularity theory and use them to explain the local behavior of this set.


A deformation of Robert Wagner foam evaluation and link homology. We consider a deformation of the Robert Wagner foam evaluation formula, with an eye toward a relation to formal groups. Integrality of the deformed evaluation is established, giving rise to state spaces for planar GL(N) MOY graphs (Murakami Ohtsuki Yamada graphs). Skein relations for the deformation are worked out in details in the GL(2) case. These skein relations deform GL(2) foam relations of Beliakova, Hogancamp, Putyra and Wehrli. We establish the Reidemeister move invariance of the resulting chain complexes assigned to link diagrams, giving us a link homology theory.


On Julia limiting directions in higher dimensions. In this paper we study, for the first time, Julia limiting directions of quasiregular mappings in mathbb R n of transcendental type. First, we give conditions under which every direction is a Julia limiting direction. Along the way, our methods show that if a quasi Fatou component contains a sectorial domain, then there is a polynomial bound on the growth in the sector. Second, we give a contribution to the inverse problem in mathbb R 3 of determining which compact subsets of S 2 can give rise to Julia limiting directions. The methods here will require showing that certain sectorial domains in mathbb R 3 are ambient quasiballs, which is a contribution to the notoriously hard problem of determining which domains are the image of the unit ball mathbb B 3 under an ambient quasiconformal map of mathbb R 3 to itself.


Constructing Depth Optimum Circuits for Adders and AND OR Paths. We examine the fundamental problem of constructing depth optimum circuits for binary addition. More precisely, as in literature, we consider the following problem Given auxiliary inputs t 0, dotsc, t m 1 , so called generate and propagate signals, construct a depth optimum circuit over the basis AND2, OR2 computing all n carry bits of an n bit adder, where m 2n 1 . In fact, carry bits are AND OR paths, i.e., Boolean functions of the form t 0 lor ( t 1 land (t 2 lor ( dots t m 1 ) dots )) . Classical approaches construct so called prefix circuits which do not achieve a competitive depth. For instance, the popular construction by Kogge and Stone is only a 2 approximation. A lower bound on the depth of any prefix circuit is 1.44 log 2 m const, while recent non prefix circuits have a depth of log 2 m log 2 log 2 m const. However, it is unknown whether any of these polynomial time approaches achieves the optimum depth for all m . We present a new exponential time algorithm solving the problem optimally. The previously best exact algorithm with a running time of mathcal O(2.45 m) is viable only for m leq 29 . Our algorithm is significantly faster We achieve a running time of mathcal O(2.02 m) and apply sophisticated pruning strategies to improve practical running times dramatically. This allows us to compute optimum circuits for all m leq 64 . Combining these computational results with new theoretical insights, we derive the optimum depths of 2 k bit adder circuits for all k leq 13 , previously known only for k leq 4 . In fact, we solve a more general problem occurring in VLSI design delay optimization of a generalization of AND OR paths where AND and OR do not necessarily alternate. Our algorithm arises from our new structure theorem which characterizes delay optimum generalized AND OR path circuits.


The Switch Point Algorithm. The Switch Point Algorithm is a new approach for solving optimal control problems whose solutions are either singular or bang bang or both singular and bang bang, and which possess a finite number of jump discontinuities in an optimal control at the points in time where the solution structure changes. Problems in this class can often be reduced to an optimization over the switching points. Formulas are derived for the derivative of the objective with respect to the switch points, the initial costate, and the terminal time. All these derivatives can be computed simultaneously in just one integration of the state and costate dynamics. Hence, gradient based unconstrained optimization techniques, including the conjugate gradient method or quasi Newton methods, can be used to compute an optimal control. The performance of the algorithm is illustrated using test problems with known solutions and comparisons with other algorithms from the literature.


Hopf algebras with enough quotients. A family of algebra maps H to A i whose common domain is a Hopf algebra is said to be jointly inner faithful if it does not factor simultaneously through a proper Hopf quotient of H . We show that tensor and free products of jointly inner faithful maps of Hopf algebras are again jointly inner faithful, generalizing a number of results in the literature on torus generation of compact quantum groups.


Measuring prominence of scientific work in online news as a proxy for impact. The impact made by a scientific paper on the work of other academics has many established metrics, including metrics based on citation counts and social media commenting. However, determination of the impact of a scientific paper on the wider society is less well established. For example, is it important for scientific work to be newsworthy? Here we present a new corpus of newspaper articles linked to the scientific papers that they describe. We find that Impact Case studies submitted to the UK Research Excellence Framework (REF) 2014 that refer to scientific papers mentioned in newspaper articles were awarded a higher score in the REF assessment. The papers associated with these case studies also feature prominently in the newspaper articles. We hypothesise that such prominence can be a useful proxy for societal impact. We therefore provide a novel baseline approach for measuring the prominence of scientific papers mentioned within news articles. Our measurement of prominence is based on semantic similarity through a graph based ranking algorithm. We find that scientific papers with an associated REF case study are more likely to have a stronger prominence score. This supports our hypothesis that linguistic prominence in news can be used to suggest the wider non academic impact of scientific work.


Persuading Voters It's Easy to Whisper, It's Hard to Speak Loud. We focus on the following natural question is it possible to influence the outcome of a voting process through the strategic provision of information to voters who update their beliefs rationally? We investigate whether it is computationally tractable to design a signaling scheme maximizing the probability with which the sender's preferred candidate is elected. We focus on the model recently introduced by Arieli and Babichenko (2019) (i.e., without inter agent externalities), and consider, as explanatory examples, k voting rule and plurality voting. There is a sharp contrast between the case in which private signals are allowed and the more restrictive setting in which only public signals are allowed. In the former, we show that an optimal signaling scheme can be computed efficiently both under a k voting rule and plurality voting. In establishing these results, we provide two general (i.e., applicable to settings beyond voting) contributions. Specifically, we extend a well known result by Dughmi and Xu (2017) to more general settings, and prove that, when the sender's utility function is anonymous, computing an optimal signaling scheme is fixed parameter tractable w.r.t. the number of receivers' actions. In the public signaling case, we show that the sender's optimal expected return cannot be approximated to within any factor under a k voting rule. This negative result easily extends to plurality voting and problems where utility functions are anonymous.


On the Convergence of the JKO scheme and Blow up of solutions for a Multi species Chemotaxis System with no Mass Preservation. This work considers a chemotaxis system for multi species that includes birth or death rate terms, which implies no mass preservation of the populations. We aim to show the convergence to a L infty weak solutions, that is local in time, of the JKO scheme arising from the Optimal Transport Theory, in the spirit of 35,14 . Currently, L infty solutions have shown to be important in order to get uniqueness. Since death rate case does not ensure global solutions, for arbitrary initial data, in this framework, it could be interest to analyze the Blowing up phenomenon of this system. Therefore, in the last section, we get sufficient conditions that implies blowing up phenomenon in finite time and we draw several stages where this occurs. This last part can be seen as a partial generalization of the blowing up results in 16 .


Progressive DNN Compression A Key to Achieve Ultra High Weight Pruning and Quantization Rates using ADMM. Weight pruning and weight quantization are two important categories of DNN model compression. Prior work on these techniques are mainly based on heuristics. A recent work developed a systematic frame work of DNN weight pruning using the advanced optimization technique ADMM (Alternating Direction Methods of Multipliers), achieving one of state of art in weight pruning results. In this work, we first extend such one shot ADMM based framework to guarantee solution feasibility and provide fast convergence rate, and generalize to weight quantization as well. We have further developed a multi step, progressive DNN weight pruning and quantization framework, with dual benefits of (i) achieving further weight pruning quantization thanks to the special property of ADMM regularization, and (ii) reducing the search space within each step. Extensive experimental results demonstrate the superior performance compared with prior work. Some highlights (i) we achieve 246x,36x, and 8x weight pruning on LeNet 5, AlexNet, and ResNet 50 models, respectively, with (almost) zero accuracy loss (ii) even a significant 61x weight pruning in AlexNet (ImageNet) results in only minor degradation in actual accuracy compared with prior work (iii) we are among the first to derive notable weight pruning results for ResNet and MobileNet models (iv) we derive the first lossless, fully binarized (for all layers) LeNet 5 for MNIST and VGG 16 for CIFAR 10 and (v) we derive the first fully binarized (for all layers) ResNet for ImageNet with reasonable accuracy loss.


Towards Large Intelligent Surface (LIS) based Communications. The concept of large intelligent surface (LIS) based communication has recently raised research attention, in which a LIS is regarded as an antenna array whose entire surface area can be used for radio signal transmission and reception. To provide a fundamental understanding of LIS based communication, this paper studies the uplink (UL) performance of LIS based communication with matched filtering. We first investigate the new properties introduced by LIS. In particular, the array gain, spatial resolution, and the capability of interference suppression are theoretically presented and characterized. Then, we study two possible LIS system layouts in terms of UL, i.e., centralized LIS (C LIS) and distributed LIS (D LIS). Our analysis showcases that a centralized system has strong capability of interference suppression in fact, interference can nearly be eliminated if the surface area is sufficient large or the frequency band is sufficient high. For D LIS, we propose a series of resource allocation algorithms, including user association scheme, orientation control, and power control, to extend the coverage area of a distributed system. Simulation results show that the proposed algorithms significantly improve the system performance, and even more importantly, we observe that D LIS outperforms CLIS in microwave bands, while C LIS is superior to D LIS in mmWave bands. These observations serve as useful guidelines for practical LIS deployments.


Automaton semigroup free products revisited. An improvement on earlier results on free products of automaton semigroups showing that a free product of two automaton semigroups is again an automaton semigroup providing there exists a homomorphism from one of the base semigroups to the other. The result is extended by induction to give a condition for a free product of finitely many automaton semigroups to be an automaton semigroup.


An integral equation method for the simulation of doubly periodic suspensions of rigid bodies in a shearing viscous flow. With rheology applications in mind, we present a fast solver for the time dependent effective viscosity of an infinite lattice containing one or more neutrally buoyant smooth rigid particles per unit cell, in a two dimensional Stokes fluid with given shear rate. At each time, the mobility problem is reformulated as a 2nd kind boundary integral equation, then discretized to spectral accuracy by the Nystrom method and solved iteratively, giving typically 10 digits of accuracy. Its solution controls the evolution of particle locations and angles in a first order system of ordinary differential equations. The formulation is placed on a rigorous footing by defining a generalized periodic Green's function for the skew lattice. Numerically, the periodized integral operator is split into a near image sum applied in linear time via the fast multipole method plus a correction field solved cheaply via proxy Stokeslets. We use barycentric quadratures to evaluate particle interactions and velocity fields accurately, even at distances much closer than the node spacing. Using first order time stepping we simulate, eg, 25 ellipses per unit cell to 3 digit accuracy on a desktop in 1 hour per shear time. Our examples show equilibration at long times, force chains, and two types of blow ups (jamming) whose power laws match lubrication theory asymptotics.


Action Sensitive Phonological Dependencies. This paper defines a subregular class of functions called the tier based synchronized strictly local (TSSL) functions. These functions are similar to the the tier based input output strictly local (TIOSL) functions, except that the locality condition is enforced not on the input and output streams, but on the computation history of the minimal subsequential finite state transducer. We show that TSSL functions naturally describe rhythmic syncope while TIOSL functions cannot, and we argue that TSSL functions provide a more restricted characterization of rhythmic syncope than existing treatments within Optimality Theory.


Adversarial Examples and Metrics. Adversarial examples are a type of attack on machine learning (ML) systems which cause misclassification of inputs. Achieving robustness against adversarial examples is crucial to apply ML in the real world. While most prior work on adversarial examples is empirical, a recent line of work establishes fundamental limitations of robust classification based on cryptographic hardness. Most positive and negative results in this field however assume that there is a fixed target metric which constrains the adversary, and we argue that this is often an unrealistic assumption. In this work we study the limitations of robust classification if the target metric is uncertain. Concretely, we construct a classification problem, which admits robust classification by a small classifier if the target metric is known at the time the model is trained, but for which robust classification is impossible for small classifiers if the target metric is chosen after the fact. In the process, we explore a novel connection between hardness of robust classification and bounded storage model cryptography.


SeismiQB a novel framework for deep learning with seismic data. In recent years, Deep Neural Networks were successfully adopted in numerous domains to solve various image related tasks, ranging from simple classification to fine borders annotation. Naturally, many researches proposed to use it to solve geological problems. Unfortunately, many of the seismic processing tools were developed years before the era of machine learning, including the most popular SEG Y data format for storing seismic cubes. Its slow loading speed heavily hampers experimentation speed, which is essential for getting acceptable results. Worse yet, there is no widely used format for storing surfaces inside the volume (for example, seismic horizons). To address these problems, we've developed an open sourced Python framework with emphasis on working with neural networks, that provides convenient tools for (i) fast loading seismic cubes in multiple data formats and converting between them, (ii) generating crops of desired shape and augmenting them with various transformations, and (iii) pairing cube data with labeled horizons or other types of geobodies.


A descriptive Main Gap Theorem. Answering one of the main questions of FHK14, Chapter 7 , we show that there is a tight connection between the depth of a classifiable shallow theory T and the Borel rank of the isomorphism relation cong kappa T on its models of size kappa , for kappa any cardinal satisfying kappa kappa kappa 2 aleph 0 . This is achieved by establishing a link between said rank and the mathcal L infty kappa Scott height of the kappa sized models of T , and yields to the following descriptive set theoretical analogue of Shelah's Main Gap Theorem Given a countable complete first order theory T , either cong kappa T is Borel with a countable Borel rank (i.e. very simple, given that the length of the relevant Borel hierarchy is kappa aleph 1 ), or it is not Borel at all. The dividing line between the two situations is the same as in Shelah's theorem, namely that of classifiable shallow theories. We also provide a Borel reducibility version of the above theorem, discuss some limitations to the possible (Borel) complexities of cong kappa T , and provide a characterization of categoricity of T in terms of the descriptive set theoretical complexity of cong kappa T .


Neural Knowledge Extraction From Cloud Service Incidents. In the last decade, two paradigm shifts have reshaped the software industry the move from boxed products to services and the widespread adoption of cloud computing. This has had a huge impact on the software development life cycle and the DevOps processes. Particularly, incident management has become critical for developing and operating large scale services. Incidents are created to ensure timely communication of service issues and, also, their resolution. Prior work on incident management has been heavily focused on the challenges with incident triaging and de duplication. In this work, we address the fundamental problem of structured knowledge extraction from service incidents. We have built SoftNER, a framework for unsupervised knowledge extraction from service incidents. We frame the knowledge extraction problem as a Named entity Recognition task for extracting factual information. SoftNER leverages structural patterns like key,value pairs and tables for bootstrapping the training data. Further, we build a novel multi task learning based BiLSTM CRF model which leverages not just the semantic context but also the data types for named entity extraction. We have deployed SoftNER at Microsoft, a major cloud service provider and have evaluated it on more than 2 months of cloud incidents. We show that the unsupervised machine learning based approach has a high precision of 0.96. Our multi task learning based deep learning model also outperforms the state of the art NER models. Lastly, using the knowledge extracted by SoftNER we are able to build significantly more accurate models for important downstream tasks like incident triaging.


Hardness of Network Satisfaction for Relation Algebras with Normal Representations. We study the computational complexity of the general network satisfaction problem for a finite relation algebra A with a normal representation B . If B contains a non trivial equivalence relation with a finite number of equivalence classes, then the network satisfaction problem for A is NP hard. As a second result, we prove hardness if B has domain size at least three and contains no non trivial equivalence relations but a symmetric atom a with a forbidden triple (a,a,a) , that is, a not leq a circ a . We illustrate how to apply our conditions on two small relation algebras.


Discontinuous Constituency Parsing with a Stack Free Transition System and a Dynamic Oracle. We introduce a novel transition system for discontinuous constituency parsing. Instead of storing subtrees in a stack i.e. a data structure with linear time sequential access the proposed system uses a set of parsing items, with constant time random access. This change makes it possible to construct any discontinuous constituency tree in exactly 4n 2 transitions for a sentence of length n . At each parsing step, the parser considers every item in the set to be combined with a focus item and to construct a new constituent in a bottom up fashion. The parsing strategy is based on the assumption that most syntactic structures can be parsed incrementally and that the set the memory of the parser remains reasonably small on average. Moreover, we introduce a provably correct dynamic oracle for the new transition system, and present the first experiments in discontinuous constituency parsing using a dynamic oracle. Our parser obtains state of the art results on three English and German discontinuous treebanks.


Error estimates for the Scaled Boundary Finite Element Method. The Scaled Boundary Finite Element Method (SBFEM) is a technique in which approximation spaces are constructed using a semi analytical approach. They are based on partitions of the computational domain by polygonal polyhedral subregions, where the shape functions approximate local Dirichlet problems with piecewise polynomial trace data. Using this operator adaptation approach, and by imposing a starlike scaling requirement on the subregions, the representation of local SBFEM shape functions in radial and surface directions are obtained from eigenvalues and eigenfunctions of an ODE system, whose coefficients are determined by the element geometry and the trace polynomial spaces. The aim of this paper is to derive a priori error estimates for SBFEM's solutions of harmonic test problems. For that, the SBFEM spaces are characterized in the context of Duffy's approximations for which a gradient orthogonality constraint is imposed. As a consequence, the scaled boundary functions are gradient orthogonal to any function in Duffy's spaces vanishing at the mesh skeleton, a mimetic version of a well known property valid for harmonic functions. This orthogonality property is applied to provide a priori SBFEM error estimates in terms of known finite element interpolant errors of the exact solution. Similarities with virtual harmonic approximations are also explored for the understanding of SBFEM convergence properties. Numerical experiments with 2D and 3D polytopal meshes confirm optimal SBFEM convergence rates for two test problems with smooth solutions. Attention is also paid to the approximation of a point singular solution by using SBFEM close to the singularity and finite element approximations elsewhere, revealing optimal accuracy rates of standard regular contexts.


Counterfactual Off Policy Training for Neural Response Generation. Open domain dialogue generation suffers from the data insufficiency problem due to the vast size of potential responses. In this paper, we propose to explore potential responses by counterfactual reasoning. Given an observed response, the counterfactual reasoning model automatically infers the outcome of an alternative policy that could have been taken. The resulting counterfactual response synthesized in hindsight is of higher quality than the response synthesized from scratch. Training on the counterfactual responses under the adversarial learning framework helps to explore the high reward area of the potential response space. An empirical study on the DailyDialog dataset shows that our approach significantly outperforms the HRED model as well as the conventional adversarial learning approaches.


Seeing the Un Scene Learning Amodal Semantic Maps for Room Navigation. We introduce a learning based approach for room navigation using semantic maps. Our proposed architecture learns to predict top down belief maps of regions that lie beyond the agent's field of view while modeling architectural and stylistic regularities in houses. First, we train a model to generate amodal semantic top down maps indicating beliefs of location, size, and shape of rooms by learning the underlying architectural patterns in houses. Next, we use these maps to predict a point that lies in the target room and train a policy to navigate to the point. We empirically demonstrate that by predicting semantic maps, the model learns common correlations found in houses and generalizes to novel environments. We also demonstrate that reducing the task of room navigation to point navigation improves the performance further.


Automatic Generation of Chatbots for Conversational Web Browsing. In this paper, we describe the foundations for generating a chatbot out of a website equipped with simple, bot specific HTML annotations. The approach is part of what we call conversational web browsing, i.e., a dialog based, natural language interaction with websites. The goal is to enable users to use content and functionality accessible through rendered UIs by "talking to websites" instead of by operating the graphical UI using keyboard and mouse. The chatbot mediates between the user and the website, operates its graphical UI on behalf of the user, and informs the user about the state of interaction. We describe the conceptual vocabulary and annotation format, the supporting conversational middleware and techniques, and the implementation of a demo able to deliver conversational web browsing experiences through Amazon Alexa.


Social Biases in NLP Models as Barriers for Persons with Disabilities. Building equitable and inclusive NLP technologies demands consideration of whether and how social attitudes are represented in ML models. In particular, representations encoded in models often inadvertently perpetuate undesirable social biases from the data on which they are trained. In this paper, we present evidence of such undesirable biases towards mentions of disability in two different English language models toxicity prediction and sentiment analysis. Next, we demonstrate that the neural embeddings that are the critical first step in most NLP pipelines similarly contain undesirable biases towards mentions of disability. We end by highlighting topical biases in the discourse about disability which may contribute to the observed model biases for instance, gun violence, homelessness, and drug addiction are over represented in texts discussing mental illness.


Seamless design of information system architecture based on adaptive clustering method. The paper considers the concept of building the architecture of an information system that provides a seamless connection between architectural representations of various levels of abstraction. The concept is based on the application of the adaptive clustering method of information systems developed by the author. Seamless connection is understood as the presence of connections between elements of architectural models related to architectural representations of various levels of abstraction.


Navigating Planar Topologies in Near Optimal Space and Time. We show that any embedding of a planar graph can be encoded succinctly while efficiently answering a number of topological queries near optimally. More precisely, we build on a succinct representation that encodes an embedding of m edges within 4m bits, which is close to the information theoretic lower bound of about 3.58m . With 4m o(m) bits of space, we show how to answer a number of topological queries relating nodes, edges, and faces, most of them in any time in omega(1) . Further, we show that with O(m) bits of space we can solve all those operations in O(1) time.


On k antichains in the unit n cube. A emph chain in the unit n cube is a set C subset 0,1 n such that for every mathbf x (x 1, ldots,x n) and mathbf y (y 1, ldots,y n) in C we either have x i le y i for all i in n , or x i ge y i for all i in n . We consider subsets, A , of the unit n cube 0,1 n that satisfy text card (A cap C) le k, , text for all chains , C subset 0,1 n , , where k is a fixed positive integer. We refer to such a set A as a k antichain. We show that the (n 1) dimensional Hausdorff measure of a k antichain in 0,1 n is at most kn and that the bound is asymptotically sharp. Moreover, we conjecture that there exist k antichains in 0,1 n whose (n 1) dimensional Hausdorff measure equals kn and we verify the validity of this conjecture when n 2 .


Can robust optimization offer improved portfolio performance? An empirical study of Indian market. The emergence of robust optimization has been driven primarily by the necessity to address the demerits of the Markowitz model. There has been a noteworthy debate regarding consideration of robust approaches as superior or at par with the Markowitz model, in terms of portfolio performance. In order to address this skepticism, we perform empirical analysis of three robust optimization models, namely the ones based on box, ellipsoidal and separable uncertainty sets. We conclude that robust approaches can be considered as a viable alternative to the Markowitz model, not only in simulated data but also in a real market setup, involving the Indian indices of S P BSE 30 and S P BSE 100. Finally, we offer qualitative and quantitative justification regarding the practical usefulness of robust optimization approaches from the point of view of number of stocks, sample size and types of data.


A stochastic spatial model for the sterile insect control strategy. In the system we study, 1's and 0's represent occupied and vacant sites in the contact process with births at rate lambda and deaths at rate 1. 1 's are sterile individuals that do not reproduce but appear spontaneously on vacant sites at rate alpha and die at rate theta alpha . We show that the system (which is attractive but has no dual) dies out at the critical value and has a nontrivial stationary distribution when it is supercritical. Our most interesting results concern the asymptotics when alpha to 0 . In this regime the process resembles the contact process in a random environment.


Multi pulse Photoacoustic Microscopy Based on k space Pseudospectral Method Simulation Study. Photoacoustic imaging has emerged in the past decades. Compared with the traditional medical imaging mode, it has better imaging performance and has great development potential in the field of biological imaging. In traditional photoacoustic microscopy, a single laser pulse is generally used to irradiate the sample to produce photoacoustic signal. And signal to noise ratio (SNR) is a very important indicator for photoacoustic imaging. In order to obtain the image with high SNR, multiple acquisition or increasing laser pulse energy is usually adopted. The former will lead to slower imaging speed, and the latter will lead to photobleaching or phototoxicity. Here, we proposes multi pulse photoacoustic microscopy, the photoacoustic signals were stimulated sequentially using multiple laser pulses in each A line data acquisition. In order to verify the feasibility of this method, a multi pulse photoacoustic imaging simulation platform is established using k Wave toolbox. The performance of multi pulse photoacoustic imaging is verified through the three scanning modes of photoacoustic microscopy A scan, B scan, and C scan. The results indicate that the SNR is proportion to the number of laser pulses used, high SNR can be achieved by low energy laser pulse. This work will help to expand the application of photoacoustic imaging.


Vanishing viscosity limit of a conservation law regularised by a Riesz Feller operator. We study a nonlocal regularisation of a scalar conservation law given by a fractional derivative of order between one and two. The nonlocal operator is of Riesz Feller type with skewness two minus its order. This equation describes the internal structure of hydraulic jumps in a shallow water model. The main purpose of the paper is the study of the vanishing viscosity limit of the Cauchy problem for this equation. First, we study the properties of the solution of the regularised problem and then we show that solutions converge to the entropy solution of the scalar conservation law in this limit in C( 0,T L 1 loc ( mathbb R )) for initial data in L infty( mathbb R ) , and in C( 0,T L 1( mathbb R )) for initial data in L infty( mathbb R ) cap BV( mathbb R ) . In order to prove these results we use weak entropy inequalities and the double scale technique of Kruzhkov. Such techniques also allow to show the L 1( mathbb R ) contraction of the regularised problem. For completeness, we study the behaviour in the tail of travelling wave solutions for genuinely nonlinear fluxes. These waves converge to shock waves in the vanishing viscosity limit, but decay algebraically as x ct to infty , rather than exponentially, the latter being a behaviour that they exhibit as x ct to infty , however. Finally, we generalise the results concerning the vanishing viscosity limit to Riesz Feller operators.


Creating a Physicist The Impact of Informal Programs on University Student Development. Physics outreach programs provide a critical context for informal experiences that promote the transition from new student to contributing physicist. Prior studies have suggested a positive link between participation in informal physics outreach programs and the development of a student's physics identity. In this study, we adopt a student focused investigation to explore the effects of informal programs on dimensions of physics identity, sense of community, 21st century skill development, and motivation. We employed a mixed methods study combining a survey instrument (117 responses) and interviews (35) with current and former undergraduate and graduate students who participated in five programs through a physics and astronomy department at a large land grant university. To examine interviews, we employed a framework based on situated learning theory, transformative learning theory, and the Dynamic Systems Model of Role Identity. Our findings, based on self reported data, show that students who facilitated informal physics programs positively developed their physics identity, experienced increased sense of belonging to the physics community, and developed 21st century career skills. Specifically, students reported positive benefits to their communication, teamwork and networking, and design skills. The benefits of these programs can be achieved by departments of any size without significant commitment of funds or changes to curriculum.


Congruences for critical values of higher derivatives of twisted Hasse Weil L functions, III. Let A be an abelian variety defined over a number field k , let p be an odd prime number and let F k be a cyclic extension of p power degree. Under not too stringent hypotheses we give an interpretation of the p component of the relevant case of the equivariant Tamagawa number conjecture in terms of integral congruence relations involving the evaluation on appropriate points of A of the rm Gal (F k) valued height pairing of Mazur and Tate. We then discuss the numerical computation of this pairing, and in particular obtain the first numerical verifications of this conjecture in situations in which the p completion of the Mordell Weil group of A over F is not a projective Galois module.


Flooding and Diameter in General Weighted Random Graphs. We study in this paper, the first passage percolation on a random graph model, the configuration model. We first introduce, the notions of weighted diameter, which is the maximum of the weighted lengths of all optimal paths between any two vertices in the graph, and the flooding time, which represents the time (weighted length) needed to reach all the vertices in the graph starting from a uniformly chosen vertex. Our result consists of describing the asymptotic behavior of the diameter and the flooding time, as the number of vertices n tends to infinity, in the case where the weight distribution G has an exponential tail behavior, and proving that this category of distributions is the largest possible for which the asymptotic behavior holds.


On Minor Left Prime Factorization Problem for Multivariate Polynomial Matrices. A new necessary and sufficient condition for the existence of minor left prime factorizations of multivariate polynomial matrices without full row rank is presented. The key idea is to establish a relationship between a matrix and its full row rank submatrix. Based on the new result, we propose an algorithm for factorizing matrices and have implemented it on the computer algebra system Maple. Two examples are given to illustrate the effectiveness of the algorithm, and experimental data shows that the algorithm is efficient.


Network Slicing for Vehicular Communication. Ultra reliable vehicle to everything (V2X) communication is essential for enabling the next generation of intelligent vehicles. V2X communication is a growing area of communication, that connects vehicles to neighboring vehicles (V2V), infrastructure (V2I) and pedestrians (V2P). Network slicing is one of the promising technology for connectivity of the next generation devices, creating several logical networks on a common and programmable physical infrastructure. Network slicing offers an efficient way to satisfy the diverse use case requirements by exploiting the benefits of shared physical infrastructure. In this regard, we propose a network slicing based communication solution for vehicular networks. In this work, we model a highway scenario with vehicles having heterogeneous traffic demands. The autonomous driving slice (safety messages) and the infotainment slice (video stream) are the two logical slices created on a common infrastructure. We formulated a network clustering and slicing algorithm to partition the vehicles into different clusters and allocate slice leaders to each cluster. Slice leaders serve its clustered vehicles with high quality V2V links and forwards safety information with low latency. On the other hand, road side unit provides infotainment service using high quality V2I links. An extensive Long Term Evolution Advanced (LTE A) system level simulator with enhancement of cellular V2X (C V2X) standard is used to evaluate the performance of the proposed method, in which it is shown that the proposed network slicing technique achieves low latency and high reliability communication.


Delone sets generated by square roots. Delone sets are locally finite point sets, such that (a) any two points are separated by a given minimum distance, and (b) there is a given radius so that every ball of that radius contains at least one point. Important examples include the vertex set of Penrose tilings and other regular model sets, which serve as a mathematical model for quasicrystals. In this note we show that the point set given by the values sqrt n e 2 pi i alpha sqrt n with n 1,2,3, ldots is a Delone set in the complex plane, for any alpha 0 . This complements Akiyama's recent observation (see arXiv 1904.10815) that sqrt n e 2 pi i alpha n with n 1,2,3, ldots forms a Delone set, if and only if alpha is badly approximated by rationals. A key difference is that our setting does not require Diophantine conditions on alpha .


Integrated pulse scope for tunable generation and intrinsic characterization of structured femtosecond laser. Numerous techniques have been demonstrated for effective generation of orbital angular momentum carrying radiation, but intracavity generation of continuously tunable pulses in the femtosecond regime remains challenging. Even if such a creation was realized, the generated pulses, like all pulses in reality, are complex and transitory objects that can only be comprehensively characterized via multidimensional spaces. An integrated lasing system that generates pulses while simultaneously quantifies them can achieve adaptive pulse tailoring. Here, we report a femtosecond pulse scope that unifies vector vortex mode locked lasing and vectorial quantification. With intracavity controlled Pancharatnam Berry phase modulation, continuous and ergodic generation of spirally polarized states along a broadband higher order Poincare sphere was realized. By intrinsically coupling a two dimensional polarization sensitive time scanning interferometer to the laser, multidimensional spatiotemporal features of the pulse were further visualized. The proposed methodology paves the way for design optimization of ultrafast optics by integrating complex femtosecond pulse generation and structural customization, facilitating its applications in optical physics research and laser based manufacturing.


Progress and Prospects in Weather and Climate Modelling. This popular article provides a short summary of the progress and prospects in Weather and Climate Modelling for the benefit of high school and undergraduate college students and early career researchers. Although this is not a comprehensive scientific article, the basic information provided here is intended to introduce students and researchers to the topic of Weather and Climate Modelling which comes under the broad discipline of Atmospheric Oceanic Climate Earth Sciences. This article briefly summarizes the historical developments, progress, scientific challenges in weather and climate modelling and career opportunities.


Hardy spaces associated with the Dunkl setting on the real line. We do research on the real method of Hardy spaces associated with the Dunkl setting on the real line for the range of 0 p 1.


A Delay Equation Model for the Atlantic Multidecadal Oscillation. A new technique to derive delay models from systems of partial differential equations, based on the Mori Zwanzig formalism, is used to derive a delay difference equation model for the Atlantic Multidecadal Oscillation. The Mori Zwanzig formalism gives a rewriting of the original system of equations which contains a memory term. This memory term can be related to a delay term in a resulting delay equation. Here the technique is applied to an idealized, but spatially extended, model of the Atlantic Multidecadal Oscillation. The resulting delay difference model is of a different type than the delay differential model which has been used to describe the El Ni no Southern Oscillation. In addition to this model, which can also be obtained by integration along characteristics, error terms for a smoothing approximation of the model have been derived from the Mori Zwanzig formalism. Our new method of deriving delay models from spatially extended models has a large potential to use delay models to study a range of climate variability phenomena.


Towards a Theory of Multi Parameter Geometrical Variational Problems Fibre Bundles, Differential Forms and Riemannian Quasiconvexity. We are concerned with the relaxation and existence theories of a general class of geometrical minimisation problems, with action integrals defined via differential forms over fibre bundles. We find natural algebraic and analytic conditions which give rise to a relaxation theory. Moreover, we propose the notion of Riemannian quasiconvexity'' for cost functions whose variables are differential forms on Riemannian manifolds, which extends the classical quasiconvexity condition in the Euclidean settings. The existence of minimisers under the Riemannian quasiconvexity condition has been established. This work may serve as a tentative generalisation of the framework developed in the recent paper B. Dacorogna and W. Gangbo, Quasiconvexity and relaxation in optimal transportation of closed differential forms, textit Arch. Ration. Mech. Anal. (2019), to appear. DOI texttt https doi.org 10.1007 s00205 019 01390 9 .


COMETA A Corpus for Medical Entity Linking in the Social Media. Whilst there has been growing progress in Entity Linking (EL) for general language, existing datasets fail to address the complex nature of health terminology in layman's language. Meanwhile, there is a growing need for applications that can understand the public's voice in the health domain. To address this we introduce a new corpus called COMETA, consisting of 20k English biomedical entity mentions from Reddit expert annotated with links to SNOMED CT, a widely used medical knowledge graph. Our corpus satisfies a combination of desirable properties, from scale and coverage to diversity and quality, that to the best of our knowledge has not been met by any of the existing resources in the field. Through benchmark experiments on 20 EL baselines from string to neural based models we shed light on the ability of these systems to perform complex inference on entities and concepts under 2 challenging evaluation scenarios. Our experimental results on COMETA illustrate that no golden bullet exists and even the best mainstream techniques still have a significant performance gap to fill, while the best solution relies on combining different views of data.


Discrete to continuum limits of planar disclinations. In materials science, wedge disclinations are defects caused by angular mismatches in the crystallographic lattice. To describe such disclinations, we introduce an atomistic model in planar domains. This model is given by a nearest neighbor type energy for the atomic bonds with an additional term to penalize change in volume. We enforce the appearance of disclinations by means of a special boundary condition. Our main result is the discrete to continuum limit of this energy as the lattice size tends to zero. Our proof method is relaxation of the energy. The main mathematical novelty of our proof is a density theorem for the special boundary condition. In addition to our limit theorem, we construct examples of planar disclinations as solutions to numerical minimization of the model and show that classical results for wedge disclinations are recovered by our analysis.


Rational Kernel on Pricing Models of Inflation Derivatives. The aim of this thesis is to analyze and renovate few main stream models on inflation derivatives. In the first chapter of the thesis, concepts of financial instruments and fundamental terms are introduced, such as coupon bond, inflation indexed bond, swap. In the second chapter of the thesis, classic models along the history of developing quantified interest rate models are introduced and analyzed. Moreover, the classification of interest rate models is introduced to help audiences understand the intrinsic ideology behind each type of models. In the third chapter of the thesis, the related mathematical knowledge is introduced. This part has the contribution on understanding the terms and relation among terms in each model introduced previously. In the fourth part of the thesis, the renovation of HJM frame work is introduced and analysis has been initiated.


FakeLocator Robust Localization of GAN Based Face Manipulations. Full face synthesis and partial face manipulation by virtue of the generative adversarial networks (GANs) and its variants have raised wide public concerns. In the multi media forensics area, detecting and ultimately locating the image forgery has become an imperative task. In this work, we investigate the architecture of existing GAN based face manipulation methods and observe that the imperfection of upsampling methods therewithin could be served as an important asset for GAN synthesized fake image detection and forgery localization. Based on this basic observation, we have proposed a novel approach, termed FakeLocator, to obtain high localization accuracy, at full resolution, on manipulated facial images. To the best of our knowledge, this is the very first attempt to solve the GAN based fake localization problem with a gray scale fakeness map that preserves more information of fake regions. To improve the universality of FakeLocator across multifarious facial attributes, we introduce an attention mechanism to guide the training of the model. To improve the universality of FakeLocator across different DeepFake methods, we propose partial data augmentation and single sample clustering on the training images. Experimental results on popular FaceForensics , DFFD datasets and seven different state of the art GAN based face generation methods have shown the effectiveness of our method. Compared with the baselines, our method performs better on various metrics. Moreover, the proposed method is robust against various real world facial image degradations such as JPEG compression, low resolution, noise, and blur.


On the p Laplacian Lichnerowicz equation on compact Riemannian manifolds. In this paper, we deal with a singular quasilinear critical elliptic equation of Lichnerowicz type involving the p Laplacian operator. With the help of the subcritical approach from variational method, we obtain the non existence, existence, and multiplicity results under some given assumptions.


Polynomial invariants on matrices and partition, Brauer algebra. We identify the dimension of the centralizer of the symmetric group mathfrak S d in the partition algebra mathcal A d( delta) and in the Brauer algebra mathcal B d( delta) with the number of multidigraphs with d arrows and the number of disjoint union of directed cycles with d arrows, respectively. Using Schur Weyl duality as a fundamental theory, we conclude that each centralizer is related with the G invariant space P d(M n( mathbf k )) G of degree d homogeneous polynomials on n times n matrices, where G is the orthogonal group and the group of permutation matrices, respectively. Our approach gives a uniform way to show that the dimensions of P d(M n( mathbf k )) G are stable for sufficiently large n .


math.NA cs.NA physics.comp-ph physics.data-an physics.flu-dyn physics.geo-ph

Continuation of Probability Density Functions using a Generalized Lyapunov Approach. Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation.