A generalization for the expected value of the earth mover's distance. The earth mover's distance (EMD), also called the first Wasserstein distance, can be naturally extended to compare arbitrarily many probability distributions, rather than only two, on the set n 1, dots,n . We present the details for this generalization, along with a highly efficient algorithm inspired by combinatorics it turns out that in the special case of three distributions, the EMD is half the sum of the pairwise EMD's. Extending the methods of Bourn and Willenbring (arXiv 1903.03673), we compute the expected value of this generalized EMD on random d tuples of distributions, using a generating function which coincides with the Hilbert series of the Segre embedding. We then use the EMD to analyze a real world data set of grade distributions.


Lattices over Bass rings and graph agglomerations. We study direct sum decompositions of torsion free, finitely generated modules over a (commutative) Bass ring R through the factorization theory of the corresponding monoid T(R) . Results of Levy Wiegand and Levy Odenthal together with a study of the local case yield an explicit description of T(R) . The monoid is typically neither factorial nor cancellative. Nevertheless, we construct a transfer homomorphism to a monoid of graph agglomerations a natural class of monoids serving as combinatorial models for the factorization theory of T(R) . As a consequence, the monoid T(R) is transfer Krull of finite type and several finiteness results on arithmetical invariants apply. We also establish results on the elasticity of T(R) and characterize when T(R) is half factorial. (Factoriality, that is, torsion free Krull Remak Schmidt Azumaya, is characterized by a theorem of Levy Odenthal.) The monoids of graph agglomerations introduced here are also of independent interest.


Probabilistic potential theory and induction of dynamical systems. In this article, we outline a version of a balayage formula in probabilistic potential theory adapted to measure preserving dynamical systems. This balayage identity generalizes the property that induced maps preserve the restriction of the original invariant measure. As an application, we prove in some cases the invariance under induction of the Green Kubo formula, as well as the invariance of a new degree 3 invariant.


StackVault Protection from Untrusted Functions. Data exfiltration attacks have led to huge data breaches. Recently, the Equifax attack affected 147M users and a third party library Apache Struts was alleged to be responsible for it. These attacks often exploit the fact that sensitive data are stored unencrypted in process memory and can be accessed by any function executing within the same process, including untrusted third party library functions. This paper presents StackVault, a kernel based system to prevent sensitive stack based data from being accessed in an unauthorized manner by intra process functions. Stack based data includes data on stack as well as data pointed to by pointer variables on stack. StackVault consists of three components (1) a set of programming APIs to allow users to specify which data needs to be protected, (2) a kernel module which uses unforgeable function identities to reliably carry out the sensitive data protection, and (3) an LLVM compiler extension that enables transparent placement of stack protection operations. The StackVault system automatically enforces stack protection through spatial and temporal access monitoring and control over both sensitive stack data and untrusted functions. We implemented StackVault and evaluated it using a number of popular real world applications, including gRPC. The results show that StackVault is effective and efficient, incurring only up to 2.4 runtime overhead.


tACS Facilitates Flickering Driving by Boosting Steady State Visual Evoked Potentials. There has become of increasing interest in transcranial alternating current stimulation (tACS) since its inception nearly a decade ago. tACS in modulating brain state is an active area of research and has been demonstrated effective in various neuropsychological and clinical domains. In the visual domain, much effort has been dedicated to brain rhythms and rhythmic stimulation, i.e., tACS. However, little is known about the interplay between the rhythmic stimulation and visual stimulation. Here, we used steady state visual evoked potential (SSVEP), induced by flickering driving as a widely used technique for frequency tagging, to investigate the aftereffect of tACS in healthy human subjects. Seven blocks of 64 channel electroencephalogram were recorded before and after the administration of 20 min 10 Hz tACS, while subjects performed several blocks of SSVEP tasks. We characterized the physiological properties of tACS aftereffect by comparing and validating the temporal, spatial, spatiotemporal and signal to noise ratio (SNR) patterns between and within blocks in real tACS and sham tACS. Our result revealed that tACS boosted the 10 Hz SSVEP significantly. Besides, the aftereffect on SSVEP was mitigated with time and lasted up to 5 min. Our results demonstrate the feasibility of facilitating the flickering driving by external rhythmic stimulation and open a new possibility to alter the brain state in a direction by noninvasive transcranial brain stimulation.


LANCE Efficient Low Precision Quantized Winograd Convolution for Neural Networks Based on Graphics Processing Units. Accelerating deep convolutional neural networks has become an active topic and sparked an interest in academia and industry. In this paper, we propose an efficient low precision quantized Winograd convolution algorithm, called LANCE, which combines the advantages of fast convolution and quantization techniques. By embedding linear quantization operations into the Winograd domain, the fast convolution can be performed efficiently under low precision computation on graphics processing units. We test neural network models with LANCE on representative image classification datasets, including SVHN, CIFAR, and ImageNet. The experimental results show that our 8 bit quantized Winograd convolution improves the performance by up to 2.40x over the full precision convolution with trivial accuracy loss.


The Leading Edge Problem in Fluid Mechanics. The self similar momentum ordinary differential equation (MODE) and the self similar partial differential equation (MPDE) have been derived and the investigation of the integrability of the MODE and the MPDE has been done by performing Painlev 'e test. A detailed discussion of the leading order behavior of the MODE and the MPDE has been presented with the latter being analyzed for the cases in which terms of increasing orders of Reynolds number have been considered. We have provided a brief introduction to Lie point symmetries and have found the Lie infinitesimal operator which when acts on the MPDE to order mathcal O (R) satisfies the Lie symmetry condition. Explicit calculations and expressions for the Lie prolongation terms have been presented. We have also investigated the integrability of various self similar equations that arise from the generalized self similar equation for different values of constants alpha 1,2,3 . Foundational work on transitional boundary solutions has been presented and transition solutions have been found via application of a junction condition at the leading edge trailing edge boundary domain. A detailed discussion of semi analytical solutions via the homotopy perturbation method is presented. We find semi analytical solutions to the Falkner Skan equation and the MODE by considering a Taylor series expansion as the initial approximation. An algorithmic scheme that involves consideration of a multi dimensional Taylor expansion as the initial approximation to the MPDE has been presented.


True complexity of polynomial progressions in finite fields. The true complexity of a polynomial progression in finite fields corresponds to the smallest degree Gowers norm that controls the counting operator of the progression over finite fields of large characteristic. We give a conjecture that relates true complexity to algebraic relations between the terms of the progression, and we prove it for a number of progressions, including x, x y, x y 2, x y y 2 and x, x y, x 2y, x y 2 . As a corollary, we prove an asymptotic for the count of certain progressions of complexity 1 in subsets of finite fields. In the process, we obtain an equidistribution result for certain polynomial progressions, analogous to the counting lemma for systems of linear forms proved by Green and Tao.


Modification of Turbulent Boundary Layer Coherent Structures with Drag Reducing Polymer Solution. The modification of dominant coherent structures that extend through the log region of a drag reduced turbulent boundary layer is studied via examination of two point correlations from time resolved particle image velocimetry. Measurements were acquired in polymer oceans (uniform concentration) at drag reduction levels corresponding to the low drag reduction regime ( 40 ), the high drag reduction (HDR) regime ( 40 ), and at an intermediate level (46 ). The mean velocity profiles and two point correlations were compared with that of water (Newtonian, DR 0 ). These results show that with increasing drag reduction the inclination of these dominant coherent structures decrease, their streamwise extent increases, and the fluctuations in the correlations are suppressed (especially at HDR). These observations are examined in comparison with the coherent structure literature (Newtonian and polymeric).


Singularity formation in fractional Burgers' equations. The formation of singularities in finite time in non local Burgers' equations, with time fractional derivative, is studied in detail. The occurrence of finite time singularity is proved, revealing the underlying mechanism, and precise estimates on the blow up time are provided. The employment of the present equation to model a problem arising in job market is also analyzed.


Interactive Refinement of Cross Lingual Word Embeddings. Cross lingual word embeddings transfer knowledge between languages models trained on high resource languages can predict in low resource languages. We introduce CLIME, an interactive system to quickly refine cross lingual word embeddings for a given classification problem. First, CLIME ranks words by their salience to the downstream task. Then, users mark similarity between keywords and their nearest neighbors in the embedding space. Finally, CLIME updates the embeddings using the annotations. We evaluate CLIME on identifying health related text in four low resource languages Ilocano, Sinhalese, Tigrinya, and Uyghur. Embeddings refined by CLIME capture more nuanced word semantics and have higher test accuracy than the original embeddings. CLIME often improves accuracy faster than an active learning baseline and can be easily combined with active learning to improve results.


Run, Forest, Run? On Randomization and Reproducibility in Predictive Software Engineering. Machine learning (ML) has been widely used in the literature to automate software engineering tasks. However, ML outcomes may be sensitive to randomization in data sampling mechanisms and learning procedures. To understand whether and how researchers in SE address these threats, we surveyed 45 recent papers related to three predictive tasks defect prediction (DP), predictive mutation testing (PMT), and code smell detection (CSD). We found that less than 50 of the surveyed papers address the threats related to randomized data sampling (via multiple repetitions) only 8 of the papers address the random nature of ML and parameter values are rarely reported (only 18 of the papers). To assess the severity of these threats, we conducted an empirical study using 26 real world datasets commonly considered for the three predictive tasks of interest, considering eight common supervised ML classifiers. We show that different data resamplings for 10 fold cross validation lead to extreme variability in observed performance results. Furthermore, randomized ML methods also show non negligible variability for different choices of random seeds. More worryingly, performance and variability are inconsistent for different implementations of the conceptually same ML method in different libraries, as also shown through multi dataset pairwise comparison. To cope with these critical threats, we provide practical guidelines on how to validate, assess, and report the results of predictive methods.


A near linear time minimum Steiner cut algorithm for planar graphs. We consider the Minimum Steiner Cut problem on undirected planar graphs with non negative edge weights. This problem involves finding the minimum cut of the graph that separates a specified subset X of vertices (terminals) into two parts. This problem is of theoretical interest because it generalizes two classical optimization problems, Minimum s t Cut and Minimum Cut, and of practical importance because of its application to computing a lower bound for Steiner (Subset) TSP. Our algorithm has running time O(n log n log k ) where k is the number of terminals.


Bounds for the regularity of product of edge ideals. Let I and J be edge ideals in a polynomial ring R mathbb K x 1, ldots,x n with I subseteq J . In this paper, we obtain a general upper and lower bound for the Castelnuovo Mumford regularity of IJ in terms of certain invariants associated with I and J . Using these results, we explicitly compute the regularity of IJ for several classes of edge ideals. Let J 1, ldots,J d be edge ideals in a polynomial ring R with J 1 subseteq cdots subseteq J d . Finally, we compute the precise expression for the regularity of J 1 J 2 cdots J d when d in 3,4 and J d is the edge ideal of complete graph.


Large scale and high resolution analysis of food purchases and health outcomes. To complement traditional dietary surveys, which are costly and of limited scale, researchers have resorted to digital data to infer the impact of eating habits on people's health. However, online studies are limited in resolution they are carried out at regional level and do not capture precisely the composition of the food consumed. We study the association between food consumption (derived from the loyalty cards of the main grocery retailer in London) and health outcomes (derived from publicly available medical prescription records). The scale and granularity of our analysis is unprecedented we analyze 1.6B food item purchases and 1.1B medical prescriptions for the entire city of London over the course of one year. By studying food consumption down to the level of nutrients, we show that nutrient diversity and amount of calories are the strongest predictors of the prevalence of three diseases related to what is called the "metabolic syndrome" hypertension, high cholesterol, and diabetes. This syndrome is a cluster of symptoms generally associated with obesity, is common across the rich world, and affects one in four adults in the UK. Our linear regression models achieve an R2 of 0.6 when estimating the prevalence of diabetes in nearly 1000 census areas in London, and a classifier can identify (un)healthy areas with up to 91 accuracy. Interestingly, healthy areas are not necessarily well off (income matters less than what one would expect) and have distinctive features they tend to systematically eat less carbohydrates and sugar, diversify nutrients, and avoid large quantities. More generally, our study shows that analytics of digital records of grocery purchases can be used as a cheap and scalable tool for health surveillance and, upon these records, different stakeholders from governments to insurance companies to food companies could implement effective prevention strategies.


Sketch2code Generating a website from a paper mockup. An early stage of developing user facing applications is creating a wireframe to layout the interface. Once a wireframe has been created it is given to a developer to implement in code. Developing boiler plate user interface code is time consuming work but still requires an experienced developer. In this dissertation we present two approaches which automates this process, one using classical computer vision techniques, and another using a novel application of deep semantic segmentation networks. We release a dataset of websites which can be used to train and evaluate these approaches. Further, we have designed a novel evaluation framework which allows empirical evaluation by creating synthetic sketches. Our evaluation illustrates that our deep learning approach outperforms our classical computer vision approach and we conclude that deep learning is the most promising direction for future research.


MLANE Meta Learning Based Adaptive Network Embedding. Most existing random walk based network embedding methods often follow only one of two principles, homophily or structural equivalence. In real world networks, however, nodes exhibit a mixture of homophily and structural equivalence, which requires adaptive network embedding that can adaptively preserve both homophily and structural equivalence for different nodes in different down stream analysis tasks. In this paper, we propose a novel method called Meta Learning based Adaptive Network Embedding (MLANE), which can learn adaptive sampling strategy for different nodes in different tasks by incorporating sampling strategy learning with embedding learning into one optimization problem that can be solved via an end to end meta learning framework. In extensive experiments on real datasets, MLANE shows significant performance improvements over the baselines. The source code of MLANE and the datasets used in experiments and all the hyperparameter settings for baselines are available at https github.com 7733com MLANE.


Generalized Adversarial Examples Attacks and Defenses. Most of the works follow such definition of adversarial example that is imperceptible to humans but can fool the deep neural networks (DNNs). Some works find another interesting form of adversarial examples such as one which is unrecognizable to humans, but DNNs classify it as one class with high confidence and adversarial patch. Based on this phenomenon, in this paper, from the perspective of cognition of humans and machines, we propose a new definition of adversarial examples. We show that imperceptible adversarial examples, unrecognizable adversarial examples, and adversarial patches are derivates of generalized adversarial examples. Then, we propose three types of adversarial attacks based on the generalized definition. Finally, we propose a defence mechanism that achieves state of the art performance. We construct a lossy compression function to filter out the redundant features generated by the network. In this process, the perturbation produced by the attacker will be filtered out. Therefore, the defence mechanism can effectively improve the robustness of the model. The experiments show that our attack methods can effectively generate adversarial examples, and our defence method can significantly improve the adversarial robustness of DNNs compared with adversarial training. As far as we know, our defending method achieves the best performance even though we do not adopt adversarial training.


Asymptotic links between signal processing, acoustic metamaterials and biology. Biomimicry is a powerful science that takes advantage of nature's remarkable ability to devise innovative solutions to challenging problems. In this work, we use asymptotic methods to develop the mathematical foundations for the exchange of design inspiration and features between biological hearing systems, signal processing algorithms and acoustic metamaterials. Our starting point is a concise asymptotic analysis of high contrast acoustic metamaterials. We are able to fine tune this graded structure to mimic the biomechanical properties of the cochlea, at the same scale. We then turn our attention to developing a biomimetic signal processing algorithm. We use the response of the cochlea like metamaterial as an initial filtering layer and then add additional biomimetic processing stages, designed to mimic the human auditory system's ability to recognise the global properties of natural sounds. This demonstrates the three way exchange of ideas that, thanks to our analysis, is possible between signal processing, metamaterials and biology.


A new object oriented framework for solving multiphysics problems via combination of different numerical methods. Many interesting phenomena are characterized by the complex interaction of different physical processes, each often best modeled numerically via a specific approach. In this paper, we present the design and implementation of an object oriented framework for performing multiphysics simulations that allows for the monolithic coupling of different numerical schemes. In contrast, most of the currently available simulation tools are tailored towards a specific numerical model, so that one must resort to coupling different codes externally based on operator splitting. The current framework has been developed following the C 11 standard, and its main aim is to provide an environment that affords enough flexibility for developers to implement complex models while at the same time giving end users a maximum amount of control over finer details of the simulation without having to write additional code. The main challenges towards realizing these objectives are discussed in the paper, together with the manner in which they are addressed. Along with core objects representing the framework skeleton, we present the various polymorphic classes that may be utilized by developers to implement new formulations, material models and solution algorithms. The code architecture is designed to allow achievement of the aforementioned functionalities with a minimum level of inheritance in order to improve the learning curve for programmers who are not acquainted with the software. Key capabilities of the framework are demonstrated via the solution of numerical examples dealing on composite torsion, Biot poroelasticity (featuring a combined finite element finite volume formulation), and brittle crack propagation using a phase field approach.


An Overview of UPnP based IoT Security Threats, Vulnerabilities, and Prospective Solutions. Advances in the development and increased availability of smart devices ranging from small sensors to complex cloud infrastructures as well as various networking technologies and communication protocols have supported the rapid expansion of Internet of Things deployments. The Universal Plug and Play (UPnP) protocol has been widely accepted and used in the IoT domain to support interactions among heterogeneous IoT devices, in part due to zero configuration implementation which makes it feasible for use in large scale networks. The popularity and ubiquity of UPnP to support IoT systems necessitate an exploration of security risks associated with the use of the protocol for IoT deployments. In this work, we analyze security vulnerabilities of UPnP based IoT systems and identify attack opportunities by the adversaries leveraging the vulnerabilities. Finally, we propose prospective solutions to secure UPnP based IoT systems from adversarial operations.


Sturm type bounds for modular forms over functions fields. In this paper, we obtain two analogues of the Sturm bound for modular forms in the function field setting. In the case of mixed characteristic, we prove that any harmonic cochain is uniquely determined by an explicit finite number of its first Fourier coefficients where our bound is much smaller than the ones in the literature. A similar bound is derived for generators of the Hecke algebra on harmonic cochains. As an application, we present a computational criterion for checking whether two elliptic curves over the rational function field mathbb F q( theta) with same conductor are isogenous. In the case of equal characteristic, we also prove that any Drinfeld modular form is uniquely determined by an explicit finite number of its first coefficients in the t expansion.


The late start of the mean velocity overlap layer at y mathcal O (10 3) A generic feature of ducted turbulent flows. One of the key observations in the Princeton Superpipe was the late start of the logarithmic mean velocity overlap layer at a wall distance of the order of 10 3 inner units. Between y approx 150 , the start of the overlap layer in zero pressure gradient turbulent boundary layers, and y approx 10 3 , the Superpipe profile is modeled equally well by a power law or a log law with a larger slope than in the overlap layer. In this paper it is shown that the mean velocity profile in turbulent plane channel flow exhibits analogous characteristics, namely a sudden decrease of logarithmic slope (increase of kappa ) at a y mathrm break approx 600 , which marks the start of the actual overlap layer. This demonstration results from the first construction of the complete inner and outer asymptotic expansions up to order mathcal O ( mathrm Re tau) 1 from mean velocity profiles of direct numerical simulations (DNS) at moderate Reynolds numbers. A preliminary analysis of a Couette flow DNS, on the other hand, yields an increase of logarithmic slope (decrease of kappa ) at a y mathrm break approx 400 . The correlation between the sign of the slope change and the flow symmetry motivates the hypothesis that the breakpoint between the possibly universal short inner logarithmic region and the actual overlap log law corresponds to the penetration depth of large scale turbulent structures originating from the opposite wall.


Properties of mean dimension and metric mean dimension coming from the topological entropy. In the late 1990's, M. Gromov introduced the notion of mean dimension for a continuous map, which is, as well as the topological entropy, an invariant under topological conjugacy. The concept of metric mean dimension for a dynamical system was introduced by Lindenstrauss and Weiss in 2000. In this paper we will verify which properties coming from the topological entropy map are valid for both mean dimension and metric mean dimension. In particular, we will prove that the metric mean dimension map is not continuous anywhere on the set consisting of continuous maps on both the Cantor set, the interval or the circle. Finally we prove that the metric mean dimension on the set consisting of continuous map on the interval and on the circle is not lower semi continuous.


On the distribution of addition chains. In this paper we study the theory of addition chains producing any given number n geq 3 . With the goal of estimating the partial sums of an additive chain, we introduce the notion of the determiners and the regulators of an addition chain and prove the following identities begin align sum limits j 2 delta(n) 1 s j 2(n 1) ( delta(n) 1) a delta(n) r delta(n) 1 int limits 2 delta(n) 1 sum limits 2 leq j leq t r jdt. nonumber end align where begin align 2,s 3 a 3 r 3, ldots,s delta(n) a delta(n) r delta(n) ,s delta(n) 1 a delta(n) 1 r delta(n) 1 n nonumber end align are the associated generators of the chain 1,2, ldots,s delta(n) ,s delta(n) 1 n of length delta(n) . Also we obtain the identity begin align sum limits j 2 delta(n) 1 a j (n 1) ( delta(n) 1) a delta(n) r delta(n) 1 int limits 2 delta(n) 1 sum limits 2 leq j leq t r jdt. nonumber end align


Hardness of some variants of the graph coloring game. Very recently, a long standing open question proposed by Bodlaender in 1991 was answered the graph coloring game is PSPACE complete. In 2019, Andres and Lock proposed five variants of the graph coloring game and left open the question of PSPACE hardness related to them. In this paper, we prove that these variants are PSPACE complete for the graph coloring game and also for the greedy coloring game, even if the number of colors is the chromatic number. Finally, we also prove that a connected version of the graph coloring game, proposed by Charpentier et al. in 2019, is PSPACE complete.


Safe rendering of high impedance on a series elastic actuator with disturbance observer based torque control. An important performance metric for series elastic actuators is the range of impedance which they can safely render. Advanced torque control, using techniques such as the disturbance observer, improve torque tracking bandwidth and accuracy, but their impact on safe impedance range is not established. However, to define a safe impedance range requires a practical coupled stability condition. Here, passivity based conditions are proposed for two variants of DOB torque control, and validated experimentally in a high stiffness environment. While high gain PD torque control has been shown to reduce Z width, it is here shown that a DOB reduces the need for high gain PD feedback and allows a higher rendered impedance. A dynamic feedforward compensator is proposed which increases the maximum safe impedance of the DOB, validated in experimentally in collision with high stiffness environments and manual excitation.


Simulations of the spread of COVID 19 and control policies in Tunisia. We develop and analyze in this work an epidemiological model for COVID 19 using Tunisian data. Our aims are first to evaluate Tunisian control policies for COVID 19 and secondly to understand the effect of different screening, quarantine and containment strategies and the rule of the asymptomatic patients on the spread of the virus in the Tunisian population. With this work, we show that Tunisian control policies are efficient in screening infected and asymptomatic individuals and that if containment and curfew are maintained the epidemic will be quickly contained.


Unravelling the turbulent structures of temperature variations during a wind gust event a case study. The simultaneous observations from a Doppler weather radar and an instrumented micrometeorological tower, offer an opportunity to dissect the effects of a gust front on the surface layer turbulence in a tropical convective boundary layer. We present a case study where a sudden drop in temperature was noted at heights within the surface layer during the passage of a gust front in the afternoon time. Consequently, this temperature drop created an interface which separated two different turbulent regimes. In one regime the turbulent temperature fluctuations were large and energetic, whereas in the other regime they were weak and quiescent. Given its uniqueness, we investigated the size distribution and aggregation properties of the turbulent structures related to these two regimes. We found that, the size distributions of the turbulent structures for both of these regimes displayed a clear power law signature. Since power laws are synonymous with scale invariance, this indicated the passing of the gust front initiated a scale free response which governed the turbulent characteristics of the temperature fluctuations. We propose a hypothesis to link such behaviour with the self organized criticality as observed in the complex systems. However, the temporal organization of the turbulent structures, as indicated by their clustering tendencies, differed between these two regimes. For the regime, corresponding to large temperature fluctuations, the turbulent structures were significantly clustered, whose clustering properties changed with height. Contrarily, for the other regime where the temperature fluctuations were weak, the turbulent structures remained less clustered with no discernible change being observed with height.


Using Priming to Uncover the Organization of Syntactic Representations in Neural Language Models. Neural language models (LMs) perform well on tasks that require sensitivity to syntactic structure. Drawing on the syntactic priming paradigm from psycholinguistics, we propose a novel technique to analyze the representations that enable such success. By establishing a gradient similarity metric between structures, this technique allows us to reconstruct the organization of the LMs' syntactic representational space. We use this technique to demonstrate that LSTM LMs' representations of different types of sentences with relative clauses are organized hierarchically in a linguistically interpretable manner, suggesting that the LMs track abstract properties of the sentence.


Vortex axis tracking by iterative propagation (VATIP) a method for analysing three dimensional turbulent structures. Vortex is a central concept in the understanding of turbulent dynamics. Objective algorithms for the detection and extraction of vortex structures can facilitate the physical understanding of turbulence regeneration dynamics by enabling automated and quantitative analyses of these structures. Despite the wide availability of vortex identification criteria, they only label spatial regions belonging to vortices, without any information on the identity, topology and shape of individual vortices. This latter information is stored in the axis lines lining the contours of vortex tubes. In this study, a new tracking algorithm is proposed which propagates along the vortex axis lines and iteratively searches for new directions for growth. The method is validated in flow fields from transient simulations where vortices of different shapes are controllably generated. It is then applied to statistical turbulence for the analysis of vortex configurations and distributions. It is shown to reliably extract axis lines for complex three dimensional vortices generated from the walls. A new procedure is also proposed that classifies vortices into commonly observed shapes, including quasi streamwise vortices, hairpins, hooks and branches, based on their axis line topology. Clustering analysis is performed on the extracted axis lines to reveal vortex organization patterns and their potential connection to large scale motions in turbulence.


Infinite Horizon Impulse Control of Stochastic Functional Differential Equations. We consider impulse control of stochastic functional differential equations (SFDEs) driven by L 'evy processes under an additional L p Lipschitz condition on the coefficients. Our results, which are first derived for a general stochastic optimization problem over infinite horizon impulse controls and then applied to the case of a controlled SFDE, apply to the infinite horizon as well as the random horizon settings. The methodology employed to show existence of optimal controls is a probabilistic one based on the concept of Snell envelopes.


BIRADS Features Oriented Semi supervised Deep Learning for Breast Ultrasound Computer Aided Diagnosis. Breast ultrasound (US) is an effective imaging modality for breast cancer detection and diagnosis. US computer aided diagnosis (CAD) systems have been developed for decades and have employed either conventional hand crafted features or modern automatic deep learned features, the former relying on clinical experience and the latter demanding large datasets. In this paper, we have developed a novel BIRADS SDL network that integrates clinically approved breast lesion characteristics (BIRADS features) into semi supervised deep learning (SDL) to achieve accurate diagnoses with a small training dataset. Breast US images are converted to BIRADS oriented feature maps (BFMs) using a distance transformation coupled with a Gaussian filter. Then, the converted BFMs are used as the input of an SDL network, which performs unsupervised stacked convolutional auto encoder (SCAE) image reconstruction guided by lesion classification. We trained the BIRADS SDL network with an alternative learning strategy by balancing reconstruction error and classification label prediction error. We compared the performance of the BIRADS SDL network with conventional SCAE and SDL methods that use the original images as inputs, as well as with an SCAE that use BFMs as inputs. Experimental results on two breast US datasets show that BIRADS SDL ranked the best among the four networks, with classification accuracy around 92.00 and 83.90 on two datasets. These findings indicate that BIRADS SDL could be promising for effective breast US lesion CAD using small datasets.


A Tracer Based Algorithm for Automatic Generation of Seafloor Age Grids from Plate Tectonic Reconstructions. The age of the ocean floor and its time dependent age distribution control fundamental features of the Earth, such as bathymetry, sea level and mantle heat loss. Recently, the development of increasingly sophisticated reconstructions of past plate motions has provided models for plate kinematics and plate boundary evolution back in geological time. These models implicitly include the information necessary to determine the age of ocean floor that has since been lost to subduction. However, due to the lack of an automated and efficient method for generating global seafloor age grids, many tectonic models, most notably those extending back into the Paleozoic, are published without an accompanying set of age models for oceanic lithosphere. Here we present an automatic, tracer based algorithm that generates seafloor age grids from global plate tectonic reconstructions with defined plate boundaries. Our method enables us to produce the first seafloor age models for the Paleozoic's lost ocean basins. Estimated changes in sea level based on bathymetry inferred from our new age grids show good agreement with sea level record estimations from proxies, providing a possible explanation for the peak in sea level during the assembly phase of Pangea. This demonstrates how our seafloor age models can be directly compared with observables from the geologic record that extend further back in time than the constraints from preserved seafloor. Thus, our new algorithm may also aid the further development of plate tectonic reconstructions by strengthening the links between geological observations and tectonic reconstructions of deeper time.


Playing Text Based Games with Common Sense. Text based games are simulations in which an agent interacts with the world purely through natural language. They typically consist of a number of puzzles interspersed with interactions with common everyday objects and locations. Deep reinforcement learning agents can learn to solve these puzzles. However, the everyday interactions with the environment, while trivial for human players, present as additional puzzles to agents. We explore two techniques for incorporating commonsense knowledge into agents. Inferring possibly hidden aspects of the world state with either a commonsense inference model COMET, or a language model BERT. Biasing an agents exploration according to common patterns recognized by a language model. We test our technique in the 9to05 game, which is an extreme version of a text based game that requires numerous interactions with common, everyday objects in common, everyday scenarios. We conclude that agents that augment their beliefs about the world state with commonsense inferences are more robust to observational errors and omissions of common elements from text descriptions.


Fides Managing Data on Untrusted Infrastructure. Significant amounts of data are currently being stored and managed on third party servers. It is impractical for many small scale enterprises to own their private datacenters, hence renting third party servers is a viable solution for such businesses. But the increasing number of malicious attacks, both internal and external, as well as buggy software on third party servers is causing clients to lose their trust in these external infrastructures. While small enterprises cannot avoid using external infrastructures, they need the right set of protocols to manage their data on untrusted infrastructures. In this paper, we propose TFCommit, a novel atomic commitment protocol that executes transactions on data stored across multiple untrusted servers. To our knowledge, TFCommit is the first atomic commitment protocol to execute transactions in an untrusted environment without using expensive Byzantine replication. Using TFCommit, we propose an auditable data management system, Fides, residing completely on untrustworthy infrastructure. As an auditable system, Fides guarantees the detection of potentially malicious failures occurring on untrusted servers using tamper resistant logs with the support of cryptographic techniques. The experimental evaluation demonstrates the scalability and the relatively low overhead of our approach that allows executing transactions on untrusted infrastructure.


Analysis of Reference and Citation Copying in Evolving Bibliographic Networks. Extensive literature demonstrates how the copying of references (links) can lead to the emergence of various structural properties (e.g., power law degree distribution and bipartite cores) in bibliographic and other similar directed networks. However, it is also well known that the copying process is incapable of mimicking the number of directed triangles in such networks neither does it have the power to explain the obsolescence of older papers. In this paper, we propose RefOrCite, a new model that allows for copying of both the references from (i.e., out neighbors of) as well as the citations to (i.e., in neighbors of) an existing node. In contrast, the standard copying model (CP) only copies references. While retaining its spirit, RefOrCite differs from the Forest Fire (FF) model in ways that makes RefOrCite amenable to mean field analysis for degree distribution, triangle count, and densification. Empirically, RefOrCite gives the best overall agreement with observed degree distribution, triangle count, diameter, h index, and the growth of citations to newer papers.


A metric Suite for Systematic Quality Assessment of Linked Open Data. Abstract The vision of the Linked Open Data (LOD) initiative is to provide a distributed model for publishing and meaningfully interlinking open data. The realization of this goal depends strongly on the quality of the data that is published as a part of the LOD. This paper focuses on the systematic quality assessment of datasets prior to publication on the LOD cloud. To this end, we identify important quality deficiencies that need to be avoided and or resolved prior to the publication of a dataset. We then propose a set of metrics to measure these quality deficiencies in a dataset. This way, we enable the assessment and identification of undesirable quality characteristics of a dataset through our proposed metrics. This will help publishers to filter out low quality data based on the quality assessment results, which in turn enables data consumers to make better and more informed decisions when using the open datasets.


Using learning analytics to provide personalized recommendations for finding peers. This work aims to propose a method to support students in finding appropriate peers in collaborative and blended learning settings. The main goal of this research is to bridge the gap between pedagogical theory and data driven practice to provide personalized and adaptive guidance to students who engage in computer supported learning activities. The research hypothesis is that we can use Learning Analytics to model students' cognitive state and to assess whether the student is in the Zone of Proximal Development. Based on this assessment, we can plan how to provide scaffolding based on the principles of Contingent Tutoring and how to form study groups based on the principles of the Zone of Proximal Development.


Pulse delay time statistics in a superradiant laser with calcium atoms. Cold samples of calcium atoms are prepared in the metastable 3 P 1 state inside an optical cavity resonant with the narrow band (375 Hz) 1 S 0 rightarrow 3 P 1 intercombination line at 657 nm. We observe superradiant emission of hyperbolic secant shaped pulses into the cavity with an intensity proportional to the square of the atom number, a duration much shorter than the natural lifetime of the 3 P 1 state, and a delay time fluctuating from shot to shot in excellent agreement with theoretical predictions. Our incoherent pumping scheme to produce inversion on the 1 S 0 rightarrow 3 P 1 transition should be extendable to allow for continuous wave laser operation.


Local Plurisubharmonic Defining Functions on the Boundary. Necessary conditions for a domain Omega subset mathbb C n admitting a local plurisubharmonic defining function on the boundary are given. In tandem, we give an algorithm to construct a local plurisubharmonic defining function on the boundary when one exists. In some cases we show that the necessary conditions are also sufficient.


Quantum description of surface enhanced resonant Raman scattering within a hybrid optomechanical model. Surface Enhanced Raman Scattering (SERS) allows for detection and identification of molecular vibrational fingerprints in minute sample quantities. The SERS process can be also exploited for optical manipulation of molecular vibrations. We present a quantum description of Surface Enhanced Resonant Raman scattering (SERRS), in analogy to hybrid cavity optomechanics, and compare the resonant situation with the off resonant SERS. Our model predicts the existence of a regime of coherent interaction between electronic and vibrational degrees of freedom of a molecule, mediated by a plasmonic nanocavity. This coherent mechanism can be achieved by parametrically tuning the frequency and intensity of the incident pumping laser and is related to the optomechanical pumping of molecular vibrations. We find that vibrational pumping is able to selectively activate a particular vibrational mode, thus providing a mechanism to control its population and drive plasmon assisted chemistry.


Expected Linear Round Synchronization The Missing Link for Linear Byzantine SMR. State Machine Replication (SMR) solutions often divide time into rounds, with a designated leader driving decisions in each round. Progress is guaranteed once all correct processes synchronize to the same round, and the leader of that round is correct. Recently suggested Byzantine SMR solutions such as HotStuff, Tendermint, and LibraBFT achieve progress with a linear message complexity and a constant time complexity once such round synchronization occurs. But round synchronization itself incurs an additional cost. By Dolev and Reischuk's lower bound, any deterministic solution must have Omega(n 2) communication complexity. Yet the question of randomized round synchronization with an expected linear message complexity remained open. We present an algorithm that, for the first time, achieves round synchronization with expected linear message complexity and expected constant latency. Existing protocols can use our round synchronization algorithm to solve Byzantine SMR with the same asymptotic performance.


Metric Classification Network in Actual Face Recognition Scene. In order to make facial features more discriminative, some new models have recently been proposed. However, almost all of these models use the traditional face verification method, where the cosine operation is performed using the features of the bottleneck layer output. However, each of these models needs to change a threshold each time it is operated on a different test set. This is very inappropriate for application in real world scenarios. In this paper, we train a validation classifier to normalize the decision threshold, which means that the result can be obtained directly without replacing the threshold. We refer to our model as validation classifier, which achieves best result on the structure consisting of one convolution layer and six fully connected layers. To test our approach, we conduct extensive experiments on Labeled Face in the Wild (LFW) and Youtube Faces (YTF), and the relative error reduction is 25.37 and 26.60 than traditional method respectively. These experiments confirm the effectiveness of validation classifier on face recognition task.


Cycles of given lengths in unicyclic components in sparse random graphs. Let L be subset of 3,4, dots and let X n,M (L) be the number of cycles belonging to unicyclic components whose length is in L in the random graph G(n,M) . We find the limiting distribution of X n,M (L) in the subcritical regime M cn with c 1 2 and the critical regime M frac n 2 left(1 mu n 1 3 right) with mu O(1) . Depending on the regime and a condition involving the series sum l in L frac z l 2l , we obtain in the limit either a Poisson or a normal distribution as n to infty .


Parametric Flatten T Swish An Adaptive Non linear Activation Function For Deep Learning. Activation function is a key component in deep learning that performs non linear mappings between the inputs and outputs. Rectified Linear Unit (ReLU) has been the most popular activation function across the deep learning community. However, ReLU contains several shortcomings that can result in inefficient training of the deep neural networks, these are 1) the negative cancellation property of ReLU tends to treat negative inputs as unimportant information for the learning, resulting in a performance degradation 2) the inherent predefined nature of ReLU is unlikely to promote additional flexibility, expressivity, and robustness to the networks 3) the mean activation of ReLU is highly positive and leads to bias shift effect in network layers and 4) the multilinear structure of ReLU restricts the non linear approximation power of the networks. To tackle these shortcomings, this paper introduced Parametric Flatten T Swish (PFTS) as an alternative to ReLU. By taking ReLU as a baseline method, the experiments showed that PFTS improved classification accuracy on SVHN dataset by 0.31 , 0.98 , 2.16 , 17.72 , 1.35 , 0.97 , 39.99 , and 71.83 on DNN 3A, DNN 3B, DNN 4, DNN 5A, DNN 5B, DNN 5C, DNN 6, and DNN 7, respectively. Besides, PFTS also achieved the highest mean rank among the comparison methods. The proposed PFTS manifested higher non linear approximation power during training and thereby improved the predictive performance of the networks.


On geometric Goppa codes from Elementary Abelian p Extensions of mathbb F p s (x) . Let p be a prime number and s 0 an integer. In this short note, we investigate one point geometric Goppa codes associated with an elementary abelian p extension of mathbb F p s (x) . We determine their dimension and the exact minimum distance in a few cases. These codes are a special case of weak Castle codes. We also list the exact values of the second generalized Hamming weight of these codes in a few cases. Simple criteria for the self duality and the quasi self duality of these codes are also provided. Furthermore, we construct examples of quantum codes, convolutional codes, and locally recoverable codes on the function field.


ER index a referential index for encrypted genomic databases. Huge DBMSs storing genomic information are being created and engineerized for doing large scale, comprehensive and in depth analysis of human beings and their diseases. However, recent regulations like the GDPR require that sensitive data are stored and elaborated thanks to privacy by design methods and software. We designed and implemented ER index, a new full text index in minute space which was optimized for compressing and encrypting collections of genomic sequences, and for performing on them fast pattern search queries. Our new index complements the E2FM index, which was introduced to compress and encrypt collections of nucleotide sequences without relying on a reference sequence. When used on collections of highly similar sequences, the ER index allows to obtain compression ratios which are an order of magnitude smaller than those achieved with the E2FM index, but maintaining its very good search performance. Moreover, thanks to the ER index multi user and multiple keys encryption model, a single index can store the sequences related to a population of individuals so that users may perform search operations only on the sequences to which they were granted access. The ER index C source code plus scripts and data to assess the tool performance are available at https github.com EncryptedIndexes erindex.


Longterm existence of solutions of a reaction diffusion system with non local terms modeling an immune response an interpretation orientated proof. This paper shows the global existence and boundedness of solutions of a reaction diffusion system modeling liver infections. The existence proof is presented step by step and the focus lies on the interpretation of intermediate results in the context of liver infections which is modeled. Non local effects in the dynamics between the virus and the immune system cells coming from the immune response in the lymphs lead to an integro partial differential equation. While existence theorems for parabolic partial differential equations are textbook examples in the field, the additional integral term requires new approaches to proving the global existence of a solution. This allows to set up an existence proof with a focus on interpretation leading to more insight in the system and in the modeling perspective at the same time. We show the boundedness of the solution in the L 1( Omega) and the L 2( Omega) norms, and use these results to prove the global existence and boundedness of the solution. A core element of the proof is the handling of oppositely acting mechanisms in the reaction term, which occur in all population dynamics models and which results in reaction terms with opposite monotonicity behavior. In the context of modeling liver infections, the boundedness in the L infty( Omega) norm has practical relevance Large immune responses lead to strong inflammations of the liver tissue. Strong inflammations negatively impact the health of an infected person and lead to grave secondary diseases. The gained rough estimates are compared with numerical tests.


From DB nets to Coloured Petri Nets with Priorities (Extended Version). The recently introduced formalism of DB nets has brought in a new conceptual way of modelling complex dynamic systems that equally account for the process and data dimensions, considering local data as well as persistent, transactional data. DB nets combine a coloured variant of Petri nets with name creation and management (which we call nu CPN), with a relational database. The integration of these two components is realized by equipping the net with special view'' places that query the database and expose the resulting answers to the net, with actions that allow transitions to update the content of the database, and with special arcs capturing compensation in case of transaction failure. In this work, we study whether this sophisticated model can be encoded back into nu CPNs. In particular, we show that the meaningful fragment of DB nets where database queries are expressed using unions of conjunctive queries with inequalities can be faithfully encoded into nu CPNs with transition priorities. This allows us to directly exploit state of the art technologies such as CPN Tools to simulate and analyse this relevant class of DB nets.


Inventory Loops (i.e. Counting Sequences) have Pre period 2 max S 1 60 . An Inventory Sequence (S 0, S 1, S 2, ...) is the iteration of the map f defined roughly by taking an integer to its numericized description (e.g. f(1381) 211318 since " 1381 " has two 1 's, one 3 , and one 8 ). Our work analyzes the iteration under the infinite base. Any starting value of positive digits is known to be ultimately periodic 1 (e.g. S 0 1381 reaches the 1 cycle f(3122331418) 3122331418 ). Parametrizations of all possible cycles are also known 2,3 . We answer Bronstein and Fraenkel's open question of 26 years showing the pre period of any such starting value is no more than 2M 60 where M max S 1 . And oddly the period of the cycle can be determined after only O( log log M) iterations.


Self similar dynamics for the modified Korteweg de Vries equation. We prove a local well posedness result for the modified Korteweg de Vries equation in a critical space designed so that is contains self similar solutions. As a consequence, we can study the flow of this equation around self similar solutions in particular, we give an as ymptotic description of small solutions as t rightarrow infty and construct solutions with a prescribed blow up behavior as t rightarrow 0.


New approaches in modeling belt flesh pelvis interaction using obese GHBMC models. Obesity is associated with higher fatality risk and altered distribution of occupant injuries in automotive collisions partially because of the increased depth of abdominal soft tissue, which results in limited or delayed engagement of the lap belt with the pelvis and increases the risk of pelvis submarining under the lap belt exposing occupant abdomen to belt loading. Previous modeling studies have shown that pelvis submarining could not be replicated using existing human body models. The goal of this study is to perform model modifications and investigate whether they could lead to model submarining. By detaching the connections between the pelvis and surrounding flesh, submarining like belt kinematics were observed. By remeshing the flesh parts of the model, similar belt kinematics was observed but the pelvic wings were fractured. Finally, large shear deformation on the flesh together with submarining like kinematics were observed in the model with its flesh modeled using the meshless Smooth Particle Galerkin Method SPG method. The results of this study showed that SPG method has potential to simulate large deformations in soft tissue which may be necessary to improve the biofidelity of belt pelvis interaction.


EPIC An Energy Efficient, High Performance GPGPU Computing Research Infrastructure. The pursuit of many research questions requires massive computational resources. State of the art research in physical processes using simulations, the training of neural networks for deep learning, or the analysis of big data are all dependent on the availability of sufficient and performant computational resources. For such research, access to a high performance computing infrastructure is indispensable. Many scientific workloads from such research domains are inherently parallel and can benefit from the data parallel architecture of general purpose graphics processing units (GPGPUs). However, GPGPU resources are scarce at Norway's national infrastructure. EPIC is a GPGPU enabled computing research infrastructure at NTNU. It enables NTNU's researchers to perform experiments that otherwise would be impossible, as time to solution would simply take too long.


Second class particles and limit shapes of evacuation and sliding paths for random tableaux. We investigate two closely related setups. In the first one we consider a TASEP style system of particles with specified initial and final configurations. The probability of each history of the system is assumed to be equal. We show that the rescaled trajectory of the emph second class particle converges (as the size of the system tends to infinity) to a random arc of an ellipse. In the second setup we consider a uniformly random Young tableau of square shape and look for typical (in the sense of probability) sliding paths and evacuation paths in the asymptotic setting as the size of the square tends to infinity. We show that the probability distribution of such paths converges to a random meridian connecting the opposite corners of the square. We also discuss analogous results for non square Young tableaux.


Adversarial Graph Representation Adaptation for Cross Domain Facial Expression Recognition. Data inconsistency and bias are inevitable among different facial expression recognition (FER) datasets due to subjective annotating process and different collecting conditions. Recent works resort to adversarial mechanisms that learn domain invariant features to mitigate domain shift. However, most of these works focus on holistic feature adaptation, and they ignore local features that are more transferable across different datasets. Moreover, local features carry more detailed and discriminative content for expression recognition, and thus integrating local features may enable fine grained adaptation. In this work, we propose a novel Adversarial Graph Representation Adaptation (AGRA) framework that unifies graph representation propagation with adversarial learning for cross domain holistic local feature co adaptation. To achieve this, we first build a graph to correlate holistic and local regions within each domain and another graph to correlate these regions across different domains. Then, we learn the per class statistical distribution of each domain and extract holistic local features from the input image to initialize the corresponding graph nodes. Finally, we introduce two stacked graph convolution networks to propagate holistic local feature within each domain to explore their interaction and across different domains for holistic local feature co adaptation. In this way, the AGRA framework can adaptively learn fine grained domain invariant features and thus facilitate cross domain expression recognition. We conduct extensive and fair experiments on several popular benchmarks and show that the proposed AGRA framework achieves superior performance over previous state of the art methods.


Harnessing behavioral diversity to understand circuits for cognition. With the increasing acquisition of large scale neural recordings comes the challenge of inferring the computations they perform and understanding how these give rise to behavior. Here, we review emerging conceptual and technological advances that begin to address this challenge, garnering insights from both biological and artificial neural networks. We argue that neural data should be recorded during rich behavioral tasks, to model cognitive processes and estimate latent behavioral variables. Careful quantification of animal movements can also provide a more complete picture of how movements shape neural dynamics and reflect changes in brain state, such as arousal or stress. Artificial neural networks (ANNs) could serve as an important tool to connect neural dynamics and rich behavioral data. ANNs have already begun to reveal how particular behaviors can be optimally solved, generating hypotheses about how observed neural activity might drive behavior and explaining diversity in behavioral strategies.


The geometrical information encoded by the Euler obstruction of a map. In this work we investigate the topological information captured by the Euler obstruction of a map, f (X,0) to ( mathbb C 2 ,0) , with (X,0) a germ of a complex d equidimensional singular space, with d 2 , and its relation with the local Euler obstruction of the coordinate functions and, consequently, with the Brasselet number. Nevertheless, under some technical conditions on the departure variety we relate the Chern number of a special collection related to the map germ f at the origin with the number of cusps of a generic perturbation of f on a stabilization of (X,f) .


Electron FLASH Delivery at Treatment Room Isocenter for Efficient Reversible Conversion of a Clinical LINAC. Purpose In this study, procedures were developed to achieve efficient reversible conversion of a clinical linear accelerator (LINAC) and deliver electron FLASH (eFLASH) or conventional beams to the treatment room isocenter. Material Methods The LINAC was converted to deliver eFLASH beam within 20 minutes by retracting the x ray target from the beam's path, positioning the carousel on an empty port, and selecting 10 MV photon beam energy in the treatment console. Dose per pulse and average dose rate were measured in a solid water phantom at different depths with Gafchromic film and OSLD. A pulse controller counted the pulses via scattered radiation signal and gated the delivery for preset pulse count. A fast photomultiplier tube based Cherenkov detector measured per pulse beam output at 2 ns sampling rate. After conversion back to clinical mode, conventional beam output, flatness, symmetry, field size and energy were measured for all clinically commissioned energies. Results Dose per pulse of 0.86 0.01 Gy (310 7 Gy s average dose rate) were achieved at isocenter. The dose from simultaneous irradiation of film and OSLD were within 1 . The PMT showed the LINAC required about 5 pulses before the output stabilized and its long term stability was within 3 for measurements performed at 3 minutes intervals. The dose, flatness, symmetry, and photon energy were unchanged from baseline and within tolerance (1 , 3 , 2 , and 0.1 respectively) after reverting to conventional beams. Conclusion 10 MeV FLASH beams were achieved at the isocenter of the treatment room. The beam output was reproducible but requires further investigation of the ramp up time in the first 5 pulses, equivalent to 100 cGy. The eFLASH beam can irradiate both small and large subjects in minimally modified clinical settings and dose rates can be further increased by reducing the source to surface distance.


Effectus of Quantum Probability on Relational Structures. The notion of effectus from categorical logic is relevant in the emerging field of categorical probability theory. In some cases, stochastic maps are represented by maps in the Kleisli category of some probability monad. Quantum homomorphisms from combinatorics and quantum information theory are the Kleisli maps of certain sort of quantum monad. We show that the Kleisli category of this quantum monad is an effectus. This gives rise to notions of quantum probabilistic reasoning as predicates, validity, conditioning, and channels.


Distinguishing Generalized Mycielskian Graphs. A graph G is d distinguishable if there is a coloring of the vertices with d colors so that only the trivial automorphism preserves the color classes. The smallest such d is the distinguishing number, operatorname Dist (G) . The Mycielskian mu(G) of a graph G is constructed by adding a shadow vertex u i for each vertex v i of G and one additional vertex w and adding edges so that N(u i) N G(v i) cup w . The generalized Mycielskian mu t(G) is a Mycielskian graph with t layers of shadow vertices, each with edges to layers above and below. This paper examines the distinguishing number of the traditional and generalized Mycielskian graphs. Notably, if G neq K 1, K 2 and the number of isolated vertices in mu t(G) is at most operatorname Dist (G) , then operatorname Dist ( mu t(G)) le operatorname Dist (G) . This result proves and exceeds a conjecture of Alikhani and Soltani.


A Global Stochastic Maximum Principle for Forward Backward Stochastic Control Systems with Quadratic Generators. We study a stochastic optimal control problem for forward backward control systems with quadratic generators. In order to establish the first and second order variational and adjoint equations, we obtain a new estimate for one dimensional linear BSDEs with unbounded stochastic Lipschitz coefficients involving bounded mean oscillation martingales (BMO martingales for short) and prove the solvability for a class of multi dimensional BSDEs with this type. Finally, a new global stochastic maximum principle is deduced.


DREAMT Embodied Motivational Conversational Storytelling. Storytelling is fundamental to language, including culture, conversation and communication in their broadest senses. It thus emerges as an essential component of intelligent systems, including systems where natural language is not a primary focus or where we do not usually think of a story being involved. In this paper we explore the emergence of storytelling as a requirement in embodied conversational agents, including its role in educational and health interventions, as well as in a general purpose computer interface for people with disabilities or other constraints that prevent the use of traditional keyboard and speech interfaces. We further present a characterization of storytelling as an inventive fleshing out of detail according to a particular personal perspective, and propose the DREAMT model to focus attention on the different layers that need to be present in a character driven storytelling system. Most if not all aspects of the DREAMT model have arisen from or been explored in some aspect of our implemented research systems, but currently only at a primitive and relatively unintegrated level. However, this experience leads us to formalize and elaborate the DREAMT model mnemonically as follows Description Dialogue Definition Denotation Realization Representation Role Explanation Education Entertainment Actualization Activation Motivation Modelling Topicalization Transformation


Relaxation rate formula for the entropic lattice Boltzmann method. An elegant and uniform relaxation rate formula is presented for the entropic lattice Boltzmann method (ELBM). The formula not only guarantees the discrete time H theorem at numerical level but also gives full consideration to the consistency with hydrodynamics. With this novel formula, the computational cost of the ELBM is significantly reduced and the method now can be efficiently used for a broad range of hydrodynamics applications including high Renolds number flows. Moreover, we demonstrate that the grid points where flow fields change drastically are effectively marked by the formula.


Limitations on Observability of Effects in Cyber Physical Systems. Increased interconnectivity of Cyber Physical Systems, by design or otherwise, increases the cyber attack surface and attack vectors. Observing the effects of these attacks is helpful in detecting them. In this paper, we show that many attacks on such systems result in a control loop effect we term Process Model Inconsistency (PMI). Our formal approach elucidates the relationships among incompleteness, incorrectness, safety, and inconsistency of process models. We show that incomplete process models lead to inconsistency. Surprisingly, inconsistency may arise even in complete and correct models. We illustrate our approach through an Automated Teller Machine (ATM) example, and describe the practical implications of the theoretical results.


Geolocating Political Events in Text. This work introduces a general method for automatically finding the locations where political events in text occurred. Using a novel set of 8,000 labeled sentences, I create a method to link automatically extracted events and locations in text. The model achieves human level performance on the annotation task and outperforms previous event geolocation systems. It can be applied to most event extraction systems across geographic contexts. I formalize the event location linking task, describe the neural network model, describe the potential uses of such a system in political science, and demonstrate a workflow to answer an open question on the role of conventional military offensives in causing civilian casualties in the Syrian civil war.


Generalization of Jordan centralizers over matrix rings. We prove that every Jordan left centralizer over matrix ring is a left centralizer and establish that every additive mapping T on M r(R) satisfying 2T(x 2) T(x)x xT(x) , for all x in M r(R) is of a particular form T(x) alpha x , for all x in M r(R) , where alpha is in the center of the ring R . In the successive results, it is shown that additive mappings over matrix rings satisfying certain conditions becomes a two sided centralizer and has the same form as mentioned above.


Generating Empathetic Responses by Looking Ahead the User's Sentiment. An important aspect of human conversation difficult for machines is conversing with empathy, which is to understand the user's emotion and respond appropriately. Recent neural conversation models that attempted to generate empathetic responses either focused on conditioning the output to a given emotion, or incorporating the current user emotional state. However, these approaches do not factor in how the user would feel towards the generated response. Hence, in this paper, we propose Sentiment Look ahead, which is a novel perspective for empathy that models the future user emotional state. In short, Sentiment Look ahead is a reward function under a reinforcement learning framework that provides a higher reward to the generative model when the generated utterance improves the user's sentiment. We implement and evaluate three different possible implementations of sentiment look ahead and empirically show that our proposed approach can generate significantly more empathetic, relevant, and fluent responses than other competitive baselines such as multitask learning.


Colors from correlated disordered photonic systems can we outperform nature?. Living organisms have developed a wide range of appearances from iridescent to matt textures. Interestingly, angular independent structural colors, where isotropy in the scattering structure is present, only produce coloration in the blue wavelength region of the visible spectrum. One might, therefore, wonder if such observation is a limitation of the architecture of the palette of materials available in nature. Here, by exploiting numerical modeling, we discuss the origin of isotropic structural colors without restriction to a specific light scattering regime. We show that high color purity and color saturation cannot be reached in isotropic short range order structures for red hues. This conclusion holds even in the case of advanced scatterer morphologies, such as core shell particles or inverse photonic glasses explaining recent experimental findings reporting very poor performances of visual appearance for such systems.


Models for damped water waves. In this paper we derive some new weakly nonlinear asymptotic models describing viscous waves in deep water with or without surface tension effects. These asymptotic models take into account several different dissipative effects and are obtained from the free boundary problems formulated in the works of Dias, Dyachenko and Zakharov (Physics Letters A, 2008), Jiang, Ting, Perlin and Schultz (Journal of Fluid Mechanics,1996) and Wu, Liu and Yue (Journal of Fluid Mechanics, 2006).


CUNI Systems for the Unsupervised News Translation Task in WMT 2019. In this paper we describe the CUNI translation system used for the unsupervised news shared task of the ACL 2019 Fourth Conference on Machine Translation (WMT19). We follow the strategy of Artexte et al. (2018b), creating a seed phrase based system where the phrase table is initialized from cross lingual embedding mappings trained on monolingual data, followed by a neural machine translation system trained on synthetic parallel data. The synthetic corpus was produced from a monolingual corpus by a tuned PBMT model refined through iterative back translation. We further focus on the handling of named entities, i.e. the part of vocabulary where the cross lingual embedding mapping suffers most. Our system reaches a BLEU score of 15.3 on the German Czech WMT19 shared task.


Orbit counting in polarized dynamical systems. We extend recent orbit counts for finitely generated semigroups acting on mathbb P N to certain infinitely generated, polarized semigroups acting on projective varieties. We then apply these results to semigroup orbits generated by some infinite sets of unicritical polynomials.


A Survey of Resource Management for Processing in Memory and Near Memory Processing Architectures. Due to amount of data involved in emerging deep learning and big data applications, operations related to data movement have quickly become the bottleneck. Data centric computing (DCC), as enabled by processing in memory (PIM) and near memory processing (NMP) paradigms, aims to accelerate these types of applications by moving the computation closer to the data. Over the past few years, researchers have proposed various memory architectures that enable DCC systems, such as logic layers in 3D stacked memories or charge sharing based bitwise operations in DRAM. However, application specific memory access patterns, power and thermal concerns, memory technology limitations, and inconsistent performance gains complicate the offloading of computation in DCC systems. Therefore, designing intelligent resource management techniques for computation offloading is vital for leveraging the potential offered by this new paradigm. In this article, we survey the major trends in managing PIM and NMP based DCC systems and provide a review of the landscape of resource management techniques employed by system designers for such systems. Additionally, we discuss the future challenges and opportunities in DCC management.


Using observed bacteria concentration and modeled transit time under an analytical framework to estimate overall removal rate of fecal coliform in an estuary. Abundance of fecal coliform (FC) is widely used to indicate the potential presence of pathogens, the No.1 cause of water impairments in the U.S. Despite extensive monitoring efforts, assessing and modeling FC pollution still faces challenges, largely owing to the uncertainties in estimation of overall removal rate (K). This study proposes an alternative method to estimate in situ K by combining observational data, hydrodynamic simulation, and analytical solution. The method requires the observed spatial distribution of FC concentration along an estuarine channel and the numerically simulated transit time, and converts the K estimation from a temporal problem into a spatial problem, potentially reducing survey duration, effort, and cost. Application of the method gave an estimation of K 0.5 d 1 on average for the Nassawadox Creek in Chesapeake Bay. The numerical and analytical model results with the estimated K agreed well with the observation, demonstrating the credibility of the method.


Utility Driven Mining of Trend Information for Intelligent System. Useful knowledge, embedded in a database, is likely to change over time. Identifying recent changes in temporal databases can provide valuable up to date information to decision makers. Nevertheless, techniques for mining high utility patterns (HUPs) seldom consider recency as a criterion to discover patterns. Thus, the traditional utility mining framework is inadequate for obtaining up to date insights about real world data. In this paper, we address this issue by introducing a novel framework, named utility driven mining of Recent trend high Utility Patterns (RUP) in temporal databases for intelligent systems, based on user specified minimum recency and minimum utility thresholds. The utility driven RUP algorithm is based on novel global and conditional downward closure properties, and a recency utility tree. Moreover, it adopts a vertical compact recency utility list structure to store the information required by the mining process. The developed RUP algorithm recursively discovers recent HUPs. It is also fast and consumes a small amount of memory due to its pattern discovery approach that does not generate candidates. Two improved versions of the algorithm with additional pruning strategies are also designed to speed up the discovery of patterns by reducing the search space. Results of a substantial experimental evaluation show that the proposed algorithm can efficiently identify all recent high utility patterns in large scale databases, and that the improved algorithm performs best.


An Experimental Study of ILP Formulations for the Longest Induced Path Problem. Given a graph G (V,E) , the longest induced path problem asks for a maximum cardinality node subset W subseteq V such that the graph induced by W is a path. It is a long established problem with applications, e.g., in network analysis. We propose novel integer linear programming (ILP) formulations for the problem and discuss efficient implementations thereof. Comparing them with known formulations from literature, we prove that they are beneficial in theory, yielding stronger relaxations. Moreover, our experiments show their practical superiority.


Lower bounds for eigenvalues of Laplacian operator and the clamped plate problem. In this paper, we give some lower bounds for several eigenvalues. Firstly, we investigate the eigenvalues lambda i of the Laplace operator and prove a sharp lower bound. Moreover, we extent this estimate of the eigenvalues to general cases. Secondly, we study the eigenvalues Gamma i for the clamped plate problem and deliver a sharp bound for the clamped plate problem for arbitrary dimension.


Towards fully automated post event data collection and analysis pre event and post event information fusion. In post event reconnaissance missions, engineers and researchers collect perishable information about damaged buildings in the affected geographical region to learn from the consequences of the event. A typical post event reconnaissance mission is conducted by first doing a preliminary survey, followed by a detailed survey. The preliminary survey is typically conducted by driving slowly along a pre determined route, observing the damage, and noting where further detailed data should be collected. This involves several manual, time consuming steps that can be accelerated by exploiting recent advances in computer vision and artificial intelligence. The objective of this work is to develop and validate an automated technique to support post event reconnaissance teams in the rapid collection of reliable and sufficiently comprehensive data, for planning the detailed survey. The technique incorporates several methods designed to automate the process of categorizing buildings based on their key physical attributes, and rapidly assessing their post event structural condition. It is divided into pre event and post event streams, each intending to first extract all possible information about the target buildings using both pre event and post event images. Algorithms based on convolutional neural network (CNNs) are implemented for scene (image) classification. A probabilistic approach is developed to fuse the results obtained from analyzing several images to yield a robust decision regarding the attributes and condition of a target building. We validate the technique using post event images captured during reconnaissance missions that took place after hurricanes Harvey and Irma. The validation data were collected by a structural wind and coastal engineering reconnaissance team, the National Science Foundation (NSF) funded Structural Extreme Events Reconnaissance (StEER) Network.


SLOAM Semantic Lidar Odometry and Mapping for Forest Inventory. This paper describes an end to end pipeline for tree diameter estimation based on semantic segmentation and lidar odometry and mapping. Accurate mapping of this type of environment is challenging since the ground and the trees are surrounded by leaves, thorns and vines, and the sensor typically experiences extreme motion. We propose a semantic feature based pose optimization that simultaneously refines the tree models while estimating the robot pose. The pipeline utilizes a custom virtual reality tool for labeling 3D scans that is used to train a semantic segmentation network. The masked point cloud is used to compute a trellis graph that identifies individual instances and extracts relevant features that are used by the SLAM module. We show that traditional lidar and image based methods fail in the forest environment on both Unmanned Aerial Vehicle (UAV) and hand carry systems, while our method is more robust, scalable, and automatically generates tree diameter estimations.


Testing Polynomials for Vanishing on Cartesian Products of Planar Point Sets Collinearity Testing and Related Problems. We present subquadratic algorithms, in the algebraic decision tree model of computation, for detecting whether there exists a triple of points, belonging to three respective sets A , B , and C of points in the plane, that satisfy a certain polynomial equation or two equations. The best known instance of such a problem is testing for the existence of a collinear triple of points in A times B times C , a classical 3SUM hard problem that has so far defied any attempt to obtain a subquadratic solution, whether in the (uniform) real RAM model, or in the algebraic decision tree model. While we are still unable to solve this problem, in full generality, in subquadratic time, we obtain such a solution, in the algebraic decision tree model, that uses only roughly O(n 28 15 ) constant degree polynomial sign tests, for the special case where two of the sets lie on two respective one dimensional curves and the third is placed arbitrarily in the plane. Our technique is fairly general, and applies to many other problems where we seek a triple that satisfies a single polynomial equation, e.g., determining whether A times B times C contains a triple spanning a unit area triangle. This result extends recent work by Barba etal (2017) and by Chan (2018), where all three sets A , B , and C are assumed to be one dimensional. As a second application of our technique, we again have three n point sets A , B , and C in the plane, and we want to determine whether there exists a triple (a,b,c) in A times B times C that simultaneously satisfies two independent real polynomial equations. For example, this is the setup when testing for collinearity in the complex plane, when each of the sets A , B , C lies on some constant degree algebraic curve. We show that problems of this kind can be solved with roughly O(n 24 13 ) constant degree polynomial sign tests.


IN2LAAMA INertial Lidar Localisation Autocalibration And MApping. In this paper, we present INertial Lidar Localisation Autocalibration And MApping (IN2LAAMA) an offline probabilistic framework for localisation, mapping, and extrinsic calibration based on a 3D lidar and a 6 DoF IMU. Most of today's lidars collect geometric information about the surrounding environment by sweeping lasers across their field of view. Consequently, 3D points in one lidar scan are acquired at different timestamps. If the sensor trajectory is not accurately known, the scans are affected by the phenomenon known as motion distortion. The proposed method leverages preintegration with a continuous representation of the inertial measurements to characterise the system's motion at any point in time. It enables precise correction of the motion distortion without relying on any explicit motion model. The system's pose, velocity, biases, and time shift are estimated via a full batch optimisation that includes automatically generated loop closure constraints. The autocalibration and the registration of lidar data rely on planar and edge features matched across pairs of scans. The performance of the framework is validated through simulated and real data experiments.


Geometric Graph Representations and Geometric Graph Convolutions for Deep Learning on Three Dimensional (3D) Graphs. The geometry of three dimensional (3D) graphs, consisting of nodes and edges, plays a crucial role in many important applications. An excellent example is molecular graphs, whose geometry influences important properties of a molecule including its reactivity and biological activity. To facilitate the incorporation of geometry in deep learning on 3D graphs, we define three types of geometric graph representations positional, angle geometric and distance geometric. For proof of concept, we use the distance geometric graph representation for geometric graph convolutions. Further, to utilize standard graph convolution networks, we employ a simple edge weight edge distance correlation scheme, whose parameters can be fixed using reference values or determined through Bayesian hyperparameter optimization. The results of geometric graph convolutions, for the ESOL and Freesol datasets, show significant improvement over those of standard graph convolutions. Our work demonstrates the feasibility and promise of incorporating geometry, using the distance geometric graph representation, in deep learning on 3D graphs.


Learning Variable Impedance Control for Contact Sensitive Tasks. Reinforcement learning algorithms have shown great success in solving different problems ranging from playing video games to robotics. However, they struggle to solve delicate robotic problems, especially those involving contact interactions. Though in principle a policy directly outputting joint torques should be able to learn to perform these tasks, in practice we see that it has difficulty to robustly solve the problem without any given structure in the action space. In this paper, we investigate how the choice of action space can give robust performance in presence of contact uncertainties. We propose learning a policy giving as output impedance and desired position in joint space and compare the performance of that approach to torque and position control under different contact uncertainties. Furthermore, we propose an additional reward term designed to regularize these variable impedance control policies, giving them interpretability and facilitating their transfer to real systems. We present extensive experiments in simulation of both floating and fixed base systems in tasks involving contact uncertainties, as well as results for running the learned policies on a real system.


Matrix and tensor decompositions for training binary neural networks. This paper is on improving the training of binary neural networks in which both activations and weights are binary. While prior methods for neural network binarization binarize each filter independently, we propose to instead parametrize the weight tensor of each layer using matrix or tensor decomposition. The binarization process is then performed using this latent parametrization, via a quantization function (e.g. sign function) applied to the reconstructed weights. A key feature of our method is that while the reconstruction is binarized, the computation in the latent factorized space is done in the real domain. This has several advantages (i) the latent factorization enforces a coupling of the filters before binarization, which significantly improves the accuracy of the trained models. (ii) while at training time, the binary weights of each convolutional layer are parametrized using real valued matrix or tensor decomposition, during inference we simply use the reconstructed (binary) weights. As a result, our method does not sacrifice any advantage of binary networks in terms of model compression and speeding up inference. As a further contribution, instead of computing the binary weight scaling factors analytically, as in prior work, we propose to learn them discriminatively via back propagation. Finally, we show that our approach significantly outperforms existing methods when tested on the challenging tasks of (a) human pose estimation (more than 4 improvements) and (b) ImageNet classification (up to 5 performance gains).


A Survey of Knowledge based Sequential Decision Making under Uncertainty. Reasoning with declarative knowledge (RDK) and sequential decision making (SDM) are two key research areas in artificial intelligence. RDK methods reason with declarative domain knowledge, including commonsense knowledge, that is either provided a priori or acquired over time, while SDM methods (probabilistic planning and reinforcement learning) seek to compute action policies that maximize the expected cumulative utility over a time horizon both classes of methods reason in the presence of uncertainty. Despite the rich literature in these two areas, researchers have not fully explored their complementary strengths. In this paper, we survey algorithms that leverage RDK methods while making sequential decisions under uncertainty. We discuss significant developments, open problems, and directions for future work.


Residual pathologies. Several counterexamples in analysis show the existence of some special object with some sort of pathological behavior. We present three different examples where the pathological behavior is not an isolated exception, but it is the "typical" behavior of the "generic" object in a suitable class, where here generic means residual in the sense of Baire category. The first example is the revisitation of a classical result concerning approximate differentiation. The second example is the derivative loss for solutions to linear wave equations with time dependent Holder continuous propagation speed. The third result is the derivative loss for solutions to transport equations with non Lipschitz velocity field.


Multidimensional Enrichment of Spatial RDF Data for SOLAP Full Version. Large volumes of spatial data and multidimensional data are being published on the Semantic Web, which has led to new opportunities for advanced analysis, such as Spatial Online Analytical Processing (SOLAP). The RDF Data Cube (QB) and QB4OLAP vocabularies have been widely used for annotating and publishing statistical and multidimensional RDF data. Although such statistical data sets might have spatial information, such as coordinates, the lack of spatial semantics and spatial multidimensional concepts in QB4OLAP and QB prevents users from employing SOLAP queries over spatial data using SPARQL. The QB4SOLAP vocabulary, on the other hand, fully supports annotating spatial and multidimensional data on the Semantic Web and enables users to query endpoints with SOLAP operators in SPARQL. To bridge the gap between QB QB4OLAP and QB4SOLAP, we propose an RDF2SOLAP enrichment model that automatically annotates spatial multidimensional concepts with QB4SOLAP and in doing so enables SOLAP on existing QB and QB4OLAP data on the Semantic Web. Furthermore, we present and evaluate a wide range of enrichment algorithms and apply them on a non trivial real world use case involving governmental open data with complex geometry types.


Polynomially Isometric Matrices in Low Dimensions. Given two d times d matrices, say A and B , when do p(A) and p(B) have the same size'' for every polynomial p ? In this article, we provide definitive results in the cases d 2 and d 3 when the notion of size used is the spectral norm.


Exergetic Port Hamiltonian Systems Modelling Basics. Port Hamiltonian systems theory provides a structured approach to modelling, optimization and control of multiphysical systems. Yet, its relationship to thermodynamics seems to be unclear. The Hamiltonian is traditionally thought of as energy, although its meaning is exergy. This insight yields benefits 1. Links to the GENERIC structure are identified, making it relatively easy to borrow ideas from a popular framework for nonequilibrium thermodynamics. 2. The port Hamiltonian structure combined with a suitable bond graph syntax is expected to become a main ingredient in thermodynamic optimization methods akin to exergy analysis and beyond. The intuitive nature of exergy and diagrammatic language facilitates interdisciplinary communication that is necessary for implementing sustainable energy systems and processes. Port Hamiltonian systems are cyclo passive, meaning that a power balance equation immediately follows from their definition. For exergetic port Hamiltonian systems, cyclo passivity is synonymous with degradation of energy and follows from the first and the second law of thermodynamics being encoded as structural properties.


Deep learning research landscape roadmap in a nutshell past, present and future Towards deep cortical learning. The past, present and future of deep learning is presented in this work. Given this landscape roadmap, we predict that deep cortical learning will be the convergence of deep learning cortical learning which builds an artificial cortical column ultimately.


ThumbNet One Thumbnail Image Contains All You Need for Recognition. Although deep convolutional neural networks (CNNs) have achieved great success in computer vision tasks, its real world application is still impeded by its voracious demand of computational resources. Current works mostly seek to compress the network by reducing its parameters or parameter incurred computation, neglecting the influence of the input image on the system complexity. Based on the fact that input images of a CNN contain substantial redundancy, in this paper, we propose a unified framework, dubbed as ThumbNet, to simultaneously accelerate and compress CNN models by enabling them to infer on one thumbnail image. We provide three effective strategies to train ThumbNet. In doing so, ThumbNet learns an inference network that performs equally well on small images as the original input network on large images. With ThumbNet, not only do we obtain the thumbnail input inference network that can drastically reduce computation and memory requirements, but also we obtain an image downscaler that can generate thumbnail images for generic classification tasks. Extensive experiments show the effectiveness of ThumbNet, and demonstrate that the thumbnail input inference network learned by ThumbNet can adequately retain the accuracy of the original input network even when the input images are downscaled 16 times.


Prime and coprime values of polynomials. The Schinzel Hypothesis is a celebrated conjecture in number theory linking polynomial values and prime numbers. In the same vein we investigate the common divisors of values P 1(n), ldots, P s(n) of several polynomials. We deduce this coprime version of the Schinzel Hypothesis under some natural assumption, coprime polynomials assume coprime values at infinitely many integers. Consequences include a version "modulo an integer" of the original Schinzel Hypothesis, with the Goldbach conjecture, again modulo an integer, as a special case.