Dissipation of electron beam driven plasma wakes. Metre scale plasma wakefield accelerators have imparted energy gain approaching 10 gigaelectronvolts to single nano Coulomb electron bunches. To reach useful average currents, however, the enormous energy density that the driver deposits into the wake must be removed efficiently between shots. Yet mechanisms by which wakes dissipate their energy into surrounding plasma remain poorly understood. Here, we report ps time resolved, grazing angle optical shadowgraphic measurements and large scale particle in cell simulations of ion channels emerging from broken wakes that electron bunches from the SLAC linac generate in tenuous lithium plasma. Measurements show the channel boundary expands radially at 1 million metres per second for over a nanosecond. Simulations show that ions and electrons that the original wake propels outward, carrying 90 percent of its energy, drive this expansion by impact ionizing surrounding neutral lithium. The results provide a basis for understanding global thermodynamics of multi GeV plasma accelerators, which underlie their viability for applications demanding high average beam current.


An adaptive multiresolution ultra weak discontinuous Galerkin method for nonlinear Schrodinger equations. This paper develops a high order adaptive scheme for solving nonlinear Schrodinger equations. The solutions to such equations often exhibit solitary wave and local structures, which makes adaptivity essential in improving the simulation efficiency. Our scheme uses the ultra weak discontinuous Galerkin (DG) formulation and belongs to the framework of adaptive multiresolution schemes. Various numerical experiments are presented to demonstrate the excellent capability of capturing the soliton waves and the blow up phenomenon.


Lorenz Gauge Condition and Its Relation to Far Field Longitudinal Electric Wave. In this short paper, it is shown that Lorenz gauge condition leads to disappearance of long range longitudinal electric wave emitted by an arbitrary electrical system. In any other gauge, longitudinal electric wave would be non negligible. Implications of the obtained results are discussed.


Syn2Real Forgery Classification via Unsupervised Domain Adaptation. In recent years, image manipulation is becoming increasingly more accessible, yielding more natural looking images, owing to the modern tools in image processing and computer vision techniques. The task of the identification of forged images has become very challenging. Amongst different types of forgeries, the cases of Copy Move forgery are increasing manifold, due to the difficulties involved to detect this tampering. To tackle such problems, publicly available datasets are insufficient. In this paper, we propose to create a synthetic forged dataset using deep semantic image inpainting and copy move forgery algorithm. However, models trained on these datasets have a significant drop in performance when tested on more realistic data. To alleviate this problem, we use unsupervised domain adaptation networks to detect copy move forgery in new domains by mapping the feature space from our synthetically generated dataset. Furthermore, we improvised the F1 score on CASIA and CoMoFoD dataset to 80.3 and 78.8 , respectively. Our approach can be helpful in those cases where the classification of data is unavailable.


Higher order hesitant fuzzy Choquet integral operator and its application to multiple criteria decision making. Generally, the criteria involved in a decision making problem are interactive or inter dependent, and therefore aggregating them by the use of traditional operators which are based on additive measures is not logical. This verifies that we have to implement fuzzy measures for modelling the interaction phenomena among the criteria.On the other hand, based on the recent extension of hesitant fuzzy set, called higher order hesitant fuzzy set (HOHFS) which allows the membership of a given element to be defined in forms of several possible generalized types of fuzzy set, we encourage to propose the higher order hesitant fuzzy (HOHF) Choquet integral operator. This concept not only considers the importance of the higher order hesitant fuzzy arguments, but also it can reflect the correlations among those arguments. Then,a detailed discussion on the aggregation properties of the HOHF Choquet integral operator will be presented.To enhance the application of HOHF Choquet integral operator in decision making, we first assess the appropriate energy policy for the socio economic development. Then, the efficiency of the proposed HOHF Choquet integral operator based technique over a number of exiting techniques is further verified by employing another decision making problem associated with the technique of TODIM (an acronym in Portuguese of Interactive and Multicriteria Decision Making).


A surprising observation on the quarter plane diffraction problem. In this paper, we revisit Radlow's highly original attempt at a double Wiener Hopf solution to the canonical problem of wave diffraction by a quarter plane. Using a constructive approach, we reduce the problem to two equations, one containing his somewhat controversial ansatz, and an additional compatibility equation. We then show that despite Radlow's ansatz being erroneous, it gives surprisingly accurate results in the far field, in particular for the spherical diffraction coefficient. This unexpectedly good result is established by comparing it to results obtained by the recently established modified Smyshlyaev formulae.


Using Deep Neural Networks for Estimating Loop Unrolling Factor. Optimizing programs requires deep expertise. On one hand, it is a tedious task, because it requires a lot of tests to find out the best combination of optimizations to apply with their best factors. On the other hand, this task is critical, because it may degrade the performance of programs instead of improving it. The automatization of this task can deal with this problem and permit to obtain good results. Optimizing loops that take the most significant part of the program execution time plays a crucial role to achieve best performance. In this paper, we address Loop unrolling optimization, by proposing a deep Neural Network model to predict the optimal unrolling factor for programs written for TIRAMISU. TIRAMISU is a polyhedral framework designed to generate high performance code for multiple platforms including multicores, GPUs, and distributed machines. TIRAMISU introduces a scheduling language with novel commands to explicitly manage the complexities that arise when targeting these systems.


A Practical Tutorial on Graph Neural Networks. Graph neural networks (GNNs) have recently grown in popularity in the field of artificial intelligence (AI) due to their unique ability to ingest relatively unstructured data types as input data. Although some elements of the GNN architecture are conceptually similar in operation to traditional neural networks (and neural network variants), other elements represent a departure from traditional deep learning techniques. This tutorial exposes the power and novelty of GNNs to AI practitioners by collating and presenting details regarding the motivations, concepts, mathematics, and applications of the most common and performant variants of GNNs. Importantly, we present this tutorial concisely, alongside practical examples, thus providing a practical and accessible tutorial on the topic of GNNs.


Point of Interest Type Inference from Social Media Text. Physical places help shape how we perceive the experiences we have there. For the first time, we study the relationship between social media text and the type of the place from where it was posted, whether a park, restaurant, or someplace else. To facilitate this, we introduce a novel data set of sim 200,000 English tweets published from 2,761 different points of interest in the U.S., enriched with place type information. We train classifiers to predict the type of the location a tweet was sent from that reach a macro F1 of 43.67 across eight classes and uncover the linguistic markers associated with each type of place. The ability to predict semantic place information from a tweet has applications in recommendation systems, personalization services and cultural geography.


A computational study of positive streamers interacting with dielectrics. We use numerical simulations to study the dynamics of surface discharges, which are common in high voltage engineering. We simulate positive streamer discharges that propagate towards a dielectric surface, attach to it, and then propagate over the surface. The simulations are performed in air with a two dimensional plasma fluid model, in which a flat dielectric is placed between two plate electrodes. Electrostatic attraction is the main mechanism that causes streamers to grow towards the dielectric. Due to the net charge in the streamer head, the dielectric gets polarized, and the electric field between the streamer and the dielectric is increased. Compared to streamers in bulk gas, surface streamers have a smaller radius, a higher electric field, a higher electron density, and higher propagation velocity. A higher applied voltage leads to faster inception and faster propagation of the surface discharge. A higher dielectric permittivity leads to more rapid attachment of the streamer to the surface and a thinner surface streamer. Secondary emission coefficients are shown to play a modest role, which is due to relatively strong photoionization in air. In the simulations, a high electric field is present between the positive streamers and the dielectric surface. We show that the magnitude and decay of this field are affected by the positive ion mobility.


Approximate Numerical Integration of the Chemical Master Equation for Stochastic Reaction Networks. Numerical solution of the chemical master equation for stochastic reaction networks typically suffers from the state space explosion problem due to the curse of dimensionality and from stiffness due to multiple time scales. The dimension of the state space equals the number of molecular species involved in the reaction network and the size of the system of differential equations equals the number of states in the corresponding continuous time Markov chain, which is usually enormously huge and often even infinite. Thus, efficient numerical solution approaches must be able to handle huge, possibly infinite and stiff systems of differential equations efficiently. We present an approximate numerical integration approach that combines a dynamical state space truncation procedure with efficient numerical integration schemes for systems of ordinary differential equations including adaptive step size selection based on local error estimates. The efficiency and accuracy is demonstrated by numerical examples.


First order convergence of Milstein schemes for McKean Vlasov equations and interacting particle systems. In this paper, we derive fully implementable first order time stepping schemes for McKean Vlasov stochastic differential equations (McKean Vlasov SDEs), allowing for a drift term with super linear growth in the state component. We propose Milstein schemes for a time discretised interacting particle system associated with the McKean Vlasov equation and prove strong convergence of order 1 and moment stability, taming the drift if only a one sided Lipschitz condition holds. To derive our main results on strong convergence rates, we make use of calculus on the space of probability measures with finite second order moments. In addition, numerical examples are presented which support our theoretical findings.


A Document Skew Detection Method Using Fast Hough Transform. The majority of document image analysis systems use a document skew detection algorithm to simplify all its further processing stages. A huge amount of such algorithms based on Hough transform (HT) analysis has already been proposed. Despite this, we managed to find only one work where the Fast Hough Transform (FHT) usage was suggested to solve the indicated problem. Unfortunately, no study of that method was provided. In this work, we propose and study a skew detection algorithm for the document images which relies on FHT analysis. To measure this algorithm quality we use the dataset from the problem oriented DISEC'13 contest and its evaluation methodology. Obtained values for AED, TOP80, and CE criteria are equal to 0.086, 0.056, 68.80 respectively.


Solution of the Unconditional Extremal Problem for a Linear Fractional Integral Functional Depending on the Parameter. The paper is devoted to the study of the unconditional extremal problem for a fractional linear integral functional defined on a set of probability distributions. In contrast to results proved earlier, the integrands of the integral expressions in the numerator and the denominator in the problem under consideration depend on a real optimization parameter vector. Thus, the optimization problem is studied on the Cartesian product of a set of probability distributions and a set of admissible values of a real parameter vector. Three statements on the extremum of a fractional linear integral functional are proved. It is established that, in all the variants, the solution of the original problem is completely determined by the extremal properties of the test function of the linear fractional integral functional this function is the ratio of the integrands of the numerator and the denominator. Possible applications of the results obtained to problems of optimal control of stochastic systems are described.


Invariant graphs and spectral type of Schr "odinger operators. In this paper we study spectral properties of Schr "odinger operators with quasi periodic potentials related to quasi periodic action minimizing trajectories for analytic twist maps. We prove that the spectrum contains a component of absolutely continuous spectrum provided that the corresponding trajectory of the twist map belongs to an analytic invariant curve.


Ricci Curvature Based Volumetric Segmentation of the Auditory Ossicles. The auditory ossicles that are located in the middle ear are the smallest bones in the human body. Their damage will result in hearing loss. It is therefore important to be able to automatically diagnose ossicles' diseases based on Computed Tomography (CT) 3D imaging. However CT images usually include the whole head area, which is much larger than the bones of interest, thus the localization of the ossicles, followed by segmentation, both play a significant role in automatic diagnosis. The commonly employed local segmentation methods require manually selected initial points, which is a highly time consuming process. We therefore propose a completely automatic method to locate the ossicles which requires neither templates, nor manual labels. It relies solely on the connective properties of the auditory ossicles themselves, and their relationship with the surrounding tissue fluid. For the segmentation task, we define a novel energy function and obtain the shape of the ossicles from the 3D CT image by minimizing this new energy. Compared to the state of the art methods which usually use the gradient operator and some normalization terms, we propose to add a Ricci curvature term to the commonly employed energy function. We compare our proposed method with the state of the art methods and show that the performance of discrete Forman Ricci curvature is superior to the others.


End to end Deep Object Tracking with Circular Loss Function for Rotated Bounding Box. The task object tracking is vital in numerous applications such as autonomous driving, intelligent surveillance, robotics, etc. This task entails the assigning of a bounding box to an object in a video stream, given only the bounding box for that object on the first frame. In 2015, a new type of video object tracking (VOT) dataset was created that introduced rotated bounding boxes as an extension of axis aligned ones. In this work, we introduce a novel end to end deep learning method based on the Transformer Multi Head Attention architecture. We also present a new type of loss function, which takes into account the bounding box overlap and orientation. Our Deep Object Tracking model with Circular Loss Function (DOTCL) shows an considerable improvement in terms of robustness over current state of the art end to end deep learning models. It also outperforms state of the art object tracking methods on VOT2018 dataset in terms of expected average overlap (EAO) metric.


DialogXL All in One XLNet for Multi Party Conversation Emotion Recognition. This paper presents our pioneering effort for emotion recognition in conversation (ERC) with pre trained language models. Unlike regular documents, conversational utterances appear alternately from different parties and are usually organized as hierarchical structures in previous work. Such structures are not conducive to the application of pre trained language models such as XLNet. To address this issue, we propose an all in one XLNet model, namely DialogXL, with enhanced memory to store longer historical context and dialog aware self attention to deal with the multi party structures. Specifically, we first modify the recurrence mechanism of XLNet from segment level to utterance level in order to better model the conversational data. Second, we introduce dialog aware self attention in replacement of the vanilla self attention in XLNet to capture useful intra and inter speaker dependencies. Extensive experiments are conducted on four ERC benchmarks with mainstream models presented for comparison. The experimental results show that the proposed model outperforms the baselines on all the datasets. Several other experiments such as ablation study and error analysis are also conducted and the results confirm the role of the critical modules of DialogXL.


Is graph based feature selection of genes better than random?. Gene interaction graphs aim to capture various relationships between genes and represent decades of biology research. When trying to make predictions from genomic data, those graphs could be used to overcome the curse of dimensionality by making machine learning models sparser and more consistent with biological common knowledge. In this work, we focus on assessing whether those graphs capture dependencies seen in gene expression data better than random. We formulate a condition that graphs should satisfy to provide a good prior knowledge and propose to test it using a Single Gene Inference' (SGI) task. We compare random graphs with seven major gene interaction graphs published by different research groups, aiming to measure the true benefit of using biologically relevant graphs in this context. Our analysis finds that dependencies can be captured almost as well at random which suggests that, in terms of gene expression levels, the relevant information about the state of the cell is spread across many genes.


Population dynamics in river networks. Natural rivers connect to each other to form networks. The geometric structure of a river network can significantly influence spatial dynamics of populations in the system. We consider a process oriented model to describe population dynamics in river networks of trees, establish the fundamental theories of the corresponding parabolic problems and elliptic problems, derive the persistence threshold by using the principal eigenvalue of the eigenvalue problem, and define the net reproductive rate to describe population persistence or extinction. By virtue of numerical simulations, we investigate the effects of hydrological, physical, and biological factors, especially the structure of the river network, on population persistence.


Learning Semantics aware Distance Map with Semantics Layering Network for Amodal Instance Segmentation. In this work, we demonstrate yet another approach to tackle the amodal segmentation problem. Specifically, we first introduce a new representation, namely a semantics aware distance map (sem dist map), to serve as our target for amodal segmentation instead of the commonly used masks and heatmaps. The sem dist map is a kind of level set representation, of which the different regions of an object are placed into different levels on the map according to their visibility. It is a natural extension of masks and heatmaps, where modal, amodal segmentation, as well as depth order information, are all well described. Then we also introduce a novel convolutional neural network (CNN) architecture, which we refer to as semantic layering network, to estimate sem dist maps layer by layer, from the global level to the instance level, for all objects in an image. Extensive experiments on the COCOA and D2SA datasets have demonstrated that our framework can predict amodal segmentation, occlusion and depth order with state of the art performance.


On the Topological Complexity of Maps. We define and develop a homotopy invariant notion for the topological complexity of a map f X to Y , denoted TC( f ), that interacts with TC( X ) and TC( Y ) in the same way cat( f ) interacts with cat( X ) and cat( Y ). Furthermore, TC( f ) and cat( f ) satisfy the same inequalities as TC( X ) and cat( X ). We compare it to other invariants defined in the papers 15,16,17,18,20 . We apply TC( f ) to studying group homomorphisms f H to G .


Strong uniqueness polynomials the complex case. The theory of strong uniqueness polynomials, satisfying the separation condition (first introduced by Fujimoto cite Fuj1 ), for complex meromorphic functions is quite complete. We construct examples of strong uniqueness polynomials which do not necessary satisfy the separation condition by constructing regular 1 forms of Wronskian type, a method introduced in cite AWW . We also use this method to produce a much easier proof in establishing the necessary and sufficient conditions for a polynomial, satisfying the separation condition, to be a strong uniqueness polynomials for meromorphic functions and rational functions.


Observation of multiramp fractional vortex beams and their properties in free space. We have experimentally investigated the evolution properties of multiramp fractional vortex beams (MFVBs) in free space, by using a fundamental Gaussian beam reflecting from a phase modulated spatial light modulator. The issue about the total vortex strength of such MFVBs is systematically addressed, and our result reveals the dependence of the total vortex strength depends on both the non integer topological charge alpha and the multiramp number m contained in initial multiramp phase structures. In the near field region, vortices contained in MFVBs are unstable and it is hard to effectively confirm the vortex strength for such fields. However, in the far field region, the evolution of vortices in fields becomes stable and the behavior of vortex strength is confirmed experimentally via measuring vortex structures by interference method. These findings give us an understanding of such complex MFVBs and may lead to potential applications in light signal process and propagation.


Efficiently Improving and Quantifying Robot Accuracy In Situ. The advancement of simulation assisted robot programming, automation of high tolerance assembly operations, and improvement of real world performance engender a need for positionally accurate robots. Despite tight machining tolerances, good mechanical design, and careful assembly, robotic arms typically exhibit average Cartesian positioning errors of several millimeters. Fortunately, the vast majority of this error can be removed in software by proper calibration of the so called "zero offsets" of a robot's joints. This research developed an automated, inexpensive, highly portable, in situ calibration method that fine tunes these kinematic parameters, thereby, improving a robot's average positioning accuracy four fold throughout its workspace. In particular, a prospective low cost motion capture system and a benchmark laser tracker were used as reference sensors for robot calibration. Bayesian inference produced optimized zero offset parameters alongside their uncertainty for data from both reference sensors. Relative and absolute accuracy metrics were proposed and applied for quantifying robot positioning accuracy. Uncertainty analysis of a validated, probabilistic robot model quantified the absolute positioning accuracy throughout its entire workspace. Altogether, three measures of accuracy conclusively revealed multi fold improvement in the positioning accuracy of the robotic arm. Bayesian inference on motion capture data yielded zero offsets and accuracy calculations comparable to those derived from laser tracker data, ultimately proving this method's viability towards robot calibration.


Learning a metacognition for object perception. Beyond representing the external world, humans also represent their own cognitive processes. In the context of perception, this metacognition helps us identify unreliable percepts, such as when we recognize that we are seeing an illusion. Here we propose MetaGen, a model for the unsupervised learning of metacognition. In MetaGen, metacognition is expressed as a generative model of how a perceptual system produces noisy percepts. Using basic principles of how the world works (such as object permanence, part of infants' core knowledge), MetaGen jointly infers the objects in the world causing the percepts and a representation of its own perceptual system. MetaGen can then use this metacognition to infer which objects are actually present in the world. On simulated data, we find that MetaGen quickly learns a metacognition and improves overall accuracy, outperforming models that lack a metacognition.


The Sobolev embedding constant on Lie groups. In this paper we estimate the Sobolev embedding constant on general noncompact Lie groups, for sub Riemannian inhomogeneous Sobolev spaces endowed with a left invariant measure. The bound that we obtain, up to a constant depending only on the group and its sub Riemannian structure, reduces to the best known bound for the classical inhomogeneous Sobolev embedding constant on mathbb R d . As an application, we prove local and global Moser Trudinger inequalities.


On modelling bicycle power for velodromes Part I Formulation for individual pursuits. For a moving bicycle, the power can be modelled as a response to the propulsion of the centre of mass of the bicycle cyclist system. On a velodrome, an accurate modelling of power requires a distinction between the trajectory of the wheels and the trajectory of the centre of mass. We formulate and examine an individual pursuit model that takes into account the aforementioned distinction. In doing so, we provide details of the invoked physical principles and mathematical derivations, with an emphasis on their limitations. We assume that a velodrome consists of two parallel straights and two semicircular arcs. We neglect the effects of the track inclination along the straights and assume the track inclination along the curves to be constant. For either segment, we consider two distinct black line speeds. For the latter, the lean angle expression is derived based on a noninertial frame of the cyclist. Among conclusions quantified by this model is the fact that a constant cadence approach to an individual pursuit does not minimize the required power.


Domain Generalization Using a Mixture of Multiple Latent Domains. When domains, which represent underlying data distributions, vary during training and testing processes, deep neural networks suffer a drop in their performance. Domain generalization allows improvements in the generalization performance for unseen target domains by using multiple source domains. Conventional methods assume that the domain to which each sample belongs is known in training. However, many datasets, such as those collected via web crawling, contain a mixture of multiple latent domains, in which the domain of each sample is unknown. This paper introduces domain generalization using a mixture of multiple latent domains as a novel and more realistic scenario, where we try to train a domain generalized model without using domain labels. To address this scenario, we propose a method that iteratively divides samples into latent domains via clustering, and which trains the domain invariant feature extractor shared among the divided latent domains via adversarial learning. We assume that the latent domain of images is reflected in their style, and thus, utilize style features for clustering. By using these features, our proposed method successfully discovers latent domains and achieves domain generalization even if the domain labels are not given. Experiments show that our proposed method can train a domain generalized model without using domain labels. Moreover, it outperforms conventional domain generalization methods, including those that utilize domain labels.


Detecting Botnet Attacks in IoT Environments An Optimized Machine Learning Approach. The increased reliance on the Internet and the corresponding surge in connectivity demand has led to a significant growth in Internet of Things (IoT) devices. The continued deployment of IoT devices has in turn led to an increase in network attacks due to the larger number of potential attack surfaces as illustrated by the recent reports that IoT malware attacks increased by 215.7 from 10.3 million in 2017 to 32.7 million in 2018. This illustrates the increased vulnerability and susceptibility of IoT devices and networks. Therefore, there is a need for proper effective and efficient attack detection and mitigation techniques in such environments. Machine learning (ML) has emerged as one potential solution due to the abundance of data generated and available for IoT devices and networks. Hence, they have significant potential to be adopted for intrusion detection for IoT environments. To that end, this paper proposes an optimized ML based framework consisting of a combination of Bayesian optimization Gaussian Process (BO GP) algorithm and decision tree (DT) classification model to detect attacks on IoT devices in an effective and efficient manner. The performance of the proposed framework is evaluated using the Bot IoT 2018 dataset. Experimental results show that the proposed optimized framework has a high detection accuracy, precision, recall, and F score, highlighting its effectiveness and robustness for the detection of botnet attacks in IoT environments.


Top DB Net Top DropBlock for Activation Enhancement in Person Re Identification. Person Re Identification is a challenging task that aims to retrieve all instances of a query image across a system of non overlapping cameras. Due to the various extreme changes of view, it is common that local regions that could be used to match people are suppressed, which leads to a scenario where approaches have to evaluate the similarity of images based on less informative regions. In this work, we introduce the Top DB Net, a method based on Top DropBlock that pushes the network to learn to focus on the scene foreground, with special emphasis on the most task relevant regions and, at the same time, encodes low informative regions to provide high discriminability. The Top DB Net is composed of three streams (i) a global stream encodes rich image information from a backbone, (ii) the Top DropBlock stream encourages the backbone to encode low informative regions with high discriminative features, and (iii) a regularization stream helps to deal with the noise created by the dropping process of the second stream, when testing the first two streams are used. Vast experiments on three challenging datasets show the capabilities of our approach against state of the art methods. Qualitative results demonstrate that our method exhibits better activation maps focusing on reliable parts of the input images.


Some identities on r central factorial numbers and r central Bell polynomials. In this paper, we introduce the extended r central factorial numbers of the second and first kinds and the extended r central Bell polynomials, as extended versions and central analogues of some previously introduced numbers and polynomials. Then we study various properties and identities related to these numbers and polynomials and also their connections.


Data Diversification A Simple Strategy For Neural Machine Translation. We introduce Data Diversification a simple but effective strategy to boost neural machine translation (NMT) performance. It diversifies the training data by using the predictions of multiple forward and backward models and then merging them with the original dataset on which the final NMT model is trained. Our method is applicable to all NMT models. It does not require extra monolingual data like back translation, nor does it add more computations and parameters like ensembles of models. Our method achieves state of the art BLEU scores of 30.7 and 43.7 in the WMT'14 English German and English French translation tasks, respectively. It also substantially improves on 8 other translation tasks 4 IWSLT tasks (English German and English French) and 4 low resource translation tasks (English Nepali and English Sinhala). We demonstrate that our method is more effective than knowledge distillation and dual learning, it exhibits strong correlation with ensembles of models, and it trades perplexity off for better BLEU score. We have released our source code at https github.com nxphi47 data diversification


Proceedings of the SREC (Social Robots in Therapy and Care) Workshop at HRI 2019. Robot Assisted Therapy (RAT) has successfully been used in Human Robot Interaction (HRI) research by including social robots in health care interventions by virtue of their ability to engage human users in both social and emotional dimensions. Robots used for these tasks must be designed with several user groups in mind, including both individuals receiving therapy and care professionals responsible for the treatment. These robots must also be able to perceive their context of use, recognize human actions and intentions, and follow the therapeutic goals to perform meaningful and personalized treatment. Effective interactions require for robots to be capable of coordinated, timely behavior in response to social cues. This means being able to estimate and predict levels of engagement, attention, intentionality and emotional state during human robot interactions. An additional challenge for social robots in therapy and care is the wide range of needs and conditions the different users can have during their interventions, even if they may share the same pathologies their current requirements and the objectives of their therapies can varied extensively. Therefore, it becomes crucial for robots to adapt their behaviors and interaction scenario to the specific needs, preferences and requirements of the patients they interact with. This personalization should be considered in terms of the robot behavior and the intervention scenario and must reflect the needs, preferences and requirements of the user.


FEQA A Question Answering Evaluation Framework for Faithfulness Assessment in Abstractive Summarization. Neural abstractive summarization models are prone to generate content inconsistent with the source document, i.e. unfaithful. Existing automatic metrics do not capture such mistakes effectively. We tackle the problem of evaluating faithfulness of a generated summary given its source document. We first collected human annotations of faithfulness for outputs from numerous models on two datasets. We find that current models exhibit a trade off between abstractiveness and faithfulness outputs with less word overlap with the source document are more likely to be unfaithful. Next, we propose an automatic question answering (QA) based metric for faithfulness, FEQA, which leverages recent advances in reading comprehension. Given question answer pairs generated from the summary, a QA model extracts answers from the document non matched answers indicate unfaithful information in the summary. Among metrics based on word overlap, embedding similarity, and learned language understanding models, our QA based metric has significantly higher correlation with human faithfulness scores, especially on highly abstractive summaries.


A Note on Generalized Algebraic Theories and Categories with Families. We give a new syntax independent definition of the notion of a generalized algebraic theory as an initial object in a category of categories with families (cwfs) with extra structure. To this end we define inductively how to build a valid signature Sigma for a generalized algebraic theory and the associated category of cwfs with a Sigma structure and cwf morphisms that preserve this structure on the nose. Our definition refers to uniform families of contexts, types, and terms, a purely semantic notion. Furthermore, we show how to syntactically construct initial cwfs with Sigma structures. This result can be viewed as a generalization of Birkhoff's completeness theorem for equational logic. It is obtained by extending Castellan, Clairambault, and Dybjer's construction of an initial cwf. We provide examples of generalized algebraic theories for monoids, categories, categories with families, and categories with families with extra structure for some type formers of dependent type theory. The models of these are internal monoids, internal categories, and internal categories with families (with extra structure) in a category with families.


Goal Recognition over Imperfect Domain Models. Goal recognition is the problem of recognizing the intended goal of autonomous agents or humans by observing their behavior in an environment. Over the past years, most existing approaches to goal and plan recognition have been ignoring the need to deal with imperfections regarding the domain model that formalizes the environment where autonomous agents behave. In this thesis, we introduce the problem of goal recognition over imperfect domain models, and develop solution approaches that explicitly deal with two distinct types of imperfect domains models (1) incomplete discrete domain models that have possible, rather than known, preconditions and effects in action descriptions and (2) approximate continuous domain models, where the transition function is approximated from past observations and not well defined. We develop novel goal recognition approaches over imperfect domains models by leveraging and adapting existing recognition approaches from the literature. Experiments and evaluation over these two types of imperfect domains models show that our novel goal recognition approaches are accurate in comparison to baseline approaches from the literature, at several levels of observability and imperfections.


Vortex dynamics for 2D Euler flows with unbounded vorticity. It is well known that the dynamics of vortices in an ideal incompressible two dimensional fluid contained in a bounded not necessarily simply connected smooth domain is described by the Kirchhoff Routh point vortex system. In this paper, we revisit the classical problem of how well solutions to the Euler equations approximate these vortex dynamics and extend previous rigorous results to the case where the vorticity field is unbounded. More precisely, we establish estimates for the 2 Wasserstein distance between the vorticity and the empirical measure associated with the point vortex dynamics. In particular, we derive an estimate on the order of weak convergence of the Euler solutions to the solutions of the point vortex system.


The Importance and the Limitations of Sim2Real for Robotic Manipulation in Precision Agriculture. In recent years Sim2Real approaches have brought great results to robotics. Techniques such as model based learning or domain randomization can help overcome the gap between simulation and reality, but in some situations simulation accuracy is still needed. An example is agricultural robotics, which needs detailed simulations, both in terms of dynamics and visuals. However, simulation software is still not capable of such quality and accuracy. Current Sim2Real techniques are helpful in mitigating the problem, but for these specific tasks they are not enough.


Making Adversarial Examples More Transferable and Indistinguishable. Fast gradient sign attack series are popular methods that are used to generate adversarial examples. However, most of the approaches based on fast gradient sign attack series cannot balance the indistinguishability and transferability due to the limitations of the basic sign structure. To address this problem, we propose a method, called Adam Iterative Fast Gradient Tanh Method (AI FGTM), to generate indistinguishable adversarial examples with high transferability. Besides, smaller kernels and dynamic step size are also applied to generate adversarial examples for further increasing the attack success rates. Extensive experiments on an ImageNet compatible dataset show that our method generates more indistinguishable adversarial examples and achieves higher attack success rates without extra running time and resource. Our best transfer based attack NI TI DI AITM can fool six classic defense models with an average success rate of 89.3 and three advanced defense models with an average success rate of 82.7 , which are higher than the state of the art gradient based attacks. Additionally, our method can also reduce nearly 20 mean perturbation. We expect that our method will serve as a new baseline for generating adversarial examples with better transferability and indistinguishability.


Smooth digraphs modulo primitive positive constructability and cyclic loop conditions. Finite smooth digraphs, that is, finite directed graphs without sources and sinks, can be partially ordered via pp constructability. We give a complete description of this poset and, in particular, we prove that it is a distributive lattice. Moreover, we show that in order to separate two smooth digraphs in our poset it suffices to show that the polymorphism clone of one of the digraphs satisfies a prime cyclic loop condition that is not satisfied by the polymorphism clone of the other. Furthermore, we prove that the poset of cyclic loop ordered by their strength for clones is a distributive lattice, too.


An Alternating Rank K Nonnegative Least Squares Framework (ARkNLS) for Nonnegative Matrix Factorization. Nonnegative matrix factorization (NMF) is a prominent technique for data dimensionality reduction that has been widely used for text mining, computer vision, pattern discovery, and bioinformatics. In this paper, a framework called ARkNLS (Alternating Rank k Nonnegativity constrained Least Squares) is proposed for computing NMF. First, a recursive formula for the solution of the rank k nonnegativity constrained least squares (NLS) is established. This recursive formula can be used to derive the closed form solution for the Rank k NLS problem for any integer k ge 1. As a result, each subproblem for an alternating rank k nonnegative least squares framework can be obtained based on this closed form solution. Assuming that all matrices involved in rank k NLS in the context of NMF computation are of full rank, two of the currently best NMF algorithms HALS (hierarchical alternating least squares) and ANLS BPP (Alternating NLS based on Block Principal Pivoting) can be considered as special cases of ARkNLS with k 1 and k r for rank r NMF, respectively. This paper is then focused on the framework with k 3, which leads to a new algorithm for NMF via the closed form solution of the rank 3 NLS problem. Furthermore, a new strategy that efficiently overcomes the potential singularity problem in rank 3 NLS within the context of NMF computation is also presented. Extensive numerical comparisons using real and synthetic data sets demonstrate that the proposed algorithm provides state of the art performance in terms of computational accuracy and cpu time.


A Recursive Algebraic Coloring Technique for Hardware Efficient Symmetric Sparse Matrix Vector Multiplication. The symmetric sparse matrix vector multiplication (SymmSpMV) is an important building block for many numerical linear algebra kernel operations or graph traversal applications. Parallelizing SymmSpMV on today's multicore platforms with up to 100 cores is difficult due to the need to manage conflicting updates on the result vector. Coloring approaches can be used to solve this problem without data duplication, but existing coloring algorithms do not take load balancing and deep memory hierarchies into account, hampering scalability and full chip performance. In this work, we propose the recursive algebraic coloring engine (RACE), a novel coloring algorithm and open source library implementation, which eliminates the shortcomings of previous coloring methods in terms of hardware efficiency and parallelization overhead. We describe the level construction, distance k coloring, and load balancing steps in RACE, use it to parallelize SymmSpMV, and compare its performance on 31 sparse matrices with other state of the art coloring techniques and Intel MKL on two modern multicore processors. RACE outperforms all other approaches substantially and behaves in accordance with the Roofline model. Outliers are discussed and analyzed in detail. While we focus on SymmSpMV in this paper, our algorithm and software is applicable to any sparse matrix operation with data dependencies that can be resolved by distance k coloring.


Predicting the Critical Number of Layers for Hierarchical Support Vector Regression. Hierarchical support vector regression (HSVR) models a function from data as a linear combination of SVR models at a range of scales, starting at a coarse scale and moving to finer scales as the hierarchy continues. In the original formulation of HSVR, there were no rules for choosing the depth of the model. In this paper, we observe in a number of models a phase transition in the training error the error remains relatively constant as layers are added, until a critical scale is passed, at which point the training error drops close to zero and remains nearly constant for added layers. We introduce a method to predict this critical scale a priori with the prediction based on the support of either a Fourier transform of the data or the Dynamic Mode Decomposition (DMD) spectrum. This allows us to determine the required number of layers prior to training any models.


Implementing geometric complexity theory On the separation of orbit closures via symmetries. Understanding the difference between group orbits and their closures is a key difficulty in geometric complexity theory (GCT) While the GCT program is set up to separate certain orbit closures, many beautiful mathematical properties are only known for the group orbits, in particular close relations with symmetry groups and invariant spaces, while the orbit closures seem much more difficult to understand. However, in order to prove lower bounds in algebraic complexity theory, considering group orbits is not enough. In this paper we tighten the relationship between the orbit of the power sum polynomial and its closure, so that we can separate this orbit closure from the orbit closure of the product of variables by just considering the symmetry groups of both polynomials and their representation theoretic decomposition coefficients. In a natural way our construction yields a multiplicity obstruction that is neither an occurrence obstruction, nor a so called vanishing ideal occurrence obstruction. All multiplicity obstructions so far have been of one of these two types. Our paper is the first implementation of the ambitious approach that was originally suggested in the first papers on geometric complexity theory by Mulmuley and Sohoni (SIAM J Comput 2001, 2008) Before our paper, all existence proofs of obstructions only took into account the symmetry group of one of the two polynomials (or tensors) that were to be separated. In our paper the multiplicity obstruction is obtained by comparing the representation theoretic decomposition coefficients of both symmetry groups. Our proof uses a semi explicit description of the coordinate ring of the orbit closure of the power sum polynomial in terms of Young tableaux, which enables its comparison to the coordinate ring of the orbit.


Modeling and simulation of heat source trajectories through phase change materials. The modeling and simulation of heat source trajectories through phase change materials is a relevant problem both for space exploration and for terrestrial climate research, among other fields. In space, the DLR and NASA are both interested in exploring beneath the surfaces of icy moons, primarily Enceladus and Europa, where conditions may alloy for extraterrestrial life. On Earth, unique sub glacial aquatic ecosystems offer potential for geo biological discoveries. Unfortunately, existing ice drilling technology is dirty and cumbersome. Melting probes are a clean and compact alternative technology which use heaters to melt through the ice. A melting probe's trajectory can be controlled with differential heating. Successful trajectory control requires advancements not only in the modeling and simulation of the ambient dynamics, but also of the probe's coupled rigid body dynamics. Fundamentally, the rigid body dynamics can be modeled by the equations of motion but this approach is prohibitively complex. This work proposes an approach which exploits that the motion of the probe is driven by contact with the evolving liquid solid interface. From this perspective, an energy minimization problem is formulated. The general mathematical problem is formulated as two split operators, respectively for the rigid body dynamics and the ambient dynamics. These operators are coupled with feasibility constraints which ensure that the probe does not penetrate the solid. Concrete examples are shown both for the energy minimization problem and for the unsteady ambient dynamics. Finally, an algorithm is presented for the temporal coupling of the split operators, which is implemented using Python and C . Example trajectories are shown, including the dynamic response of the probe velocity to a rapid change in the heat flux.


Light driven interactions in spatial predator prey model with toxicant chemotaxis. We develop and analyze a spatial temporal model of light driven ecotoxicological processes, motivated by an aquatic predator prey system of algae and textsl Daphnia subject to a contaminant. Population dynamics are driven by light, which is periodic in time and varies with spatial depth. The existence and uniqueness of spatial and temporal dependent periodic solutions are shown and analytical functions of the solutions under parameter constraints are presented. We conduct Turing stability analyses of solutions with respect to perturbations of initial conditions. Given a perturbation to a periodic equilibrium state, we show the system will return to this equilibrium state as long as motility is fast enough and or the reservoir depth is shallow enough. Analytical results assume some Dirichlet boundary conditions that match the periodic equilibrium state, however numerical simulations with more relaxed boundary conditions capture similar periodic solutions. The work sheds light onto spatially dependent population dynamics that are driven by periodic forces, such as light levels.


On the structure of the Sally module and the second normal Hilbert coefficient. The Hilbert coefficients of the normal filtration give important geometric information on the base ring like the pseudo rationality. The Sally module was introduced by W.V. Vasconcelos and it is useful to connect the Hilbert coefficients to the homological properties of the associated graded module of a Noetherian filtration. In this paper we give a complete structure of the Sally module in the case the second normal Hilbert coefficient attains almost minimal value in an analytically unramified Cohen Macaulay local ring. As a consequence, in this case we present a complete description of the Hilbert function of the associated graded ring of the normal filtration. A deep analysis of the vanishing of the third Hilbert coefficient has been necessary. This study is related to a long standing conjecture stated by S. Itoh.


On higher rank Khovanskii Teissier inequalities. We shall discuss a higher rank Khovanskii Teissier inequality, generalizing a theorem of Li. In the course of the proof, we develop new Hodge Riemann bilinear relations in certain mixed settings, which in themselves slightly extend the existing results and imply new Khovanskii Teissier type inequalities and log concavity results.


e positivity of vertical strip LLT polynomials. In this article we prove the e positivity of G mathbf nu X q 1 when G mathbf nu X q is a vertical strip LLT polynomial. This property has been conjectured by Alexandersson and Panova, and by Garsia, Haglund, Qiu and Romero, and it implies several e positivities conjectured by them and also by Bergeron. We make use of a result of Carlsson and Mellit that shows that a vertical strip LLT polynomial can be obtained by applying certain compositions of operators of the Dyck path algebra to the constant 1 . Our proof gives in fact an algorithm to expand these symmetric functions in the elementary basis, and it shows, as a byproduct, that these compositions of operators are actually multiplication operators.


On conditional connectivity of the Cartesian product of cycles. The conditional h vertex( h edge) connectivity of a connected graph H of minimum degree k h is the size of a smallest vertex(edge) set F of H such that H F is a disconnected graph of minimum degree at least h. Let G be the Cartesian product of r geq 1 cycles, each of length at least four and let h be an integer such that 0 leq h leq 2r 2 . In this paper, we determine the conditional h vertex connectivity and the conditional h edge connectivity of the graph G. We prove that both these connectivities are equal to (2r h)a h r , where a h r is the number of vertices of a smallest h regular subgraph of G.


A Thurston compactification of the space of stability conditions. We propose a compactification of the moduli space of Bridgeland stability conditions of a triangulated category. Under mild conditions on the triangulated category, we conjecture that this compactification is a real manifold with boundary, on which the action of the auto equivalence group extends continuously. The key ingredient in the compactification is an embedding of the stability space into an infinite dimensional projective space. We study this embedding in detail in the case of 2 Calabi Yau categories associated to quivers, and prove our conjectures in the A 2 and widehat A 1 cases. Central to our analysis is a detailed understanding of Harder Narasimhan multiplicities and how they transform under auto equivalences. We achieve this by introducing a structure called a Harder Narasimhan automaton and constructing examples for A 2 and widehat A 1 .


On derivatives of Kato's Euler system for elliptic curves. In this paper we study a new conjecture concerning Kato's Euler system of zeta elements for elliptic curves E over mathbb Q . This conjecture, which we refer to as the Generalized Perrin Riou Conjecture', predicts a precise congruence relation between a Darmon type derivative' of the zeta element of E over an arbitrary real abelian field and the critical value of an appropriate higher derivative of the L function of E over mathbb Q . We prove that the conjecture specializes in the relevant case of analytic rank one to recover Perrin Riou's conjecture on the logarithm of Kato's zeta element. Under mild hypotheses we also prove that the order of vanishing' part of the conjecture is valid in arbitrary rank. An Iwasawa theoretic analysis of our approach leads to the formulation and proof of a natural higher rank generalization of Rubin's formula concerning derivatives of p adic L functions. In addition, we establish a concrete and apparently new connection between the p part of the classical Birch and Swinnerton Dyer Formula and the Iwasawa Main Conjecture in arbitrary rank and for arbitrary reduction at p . In a forthcoming paper we will show that the Generalized Perrin Riou Conjecture implies (in arbitrary rank) the conjecture of Mazur and Tate concerning congruences for modular elements and, by using this approach, we are able to give a proof, under certain mild and natural hypotheses, that the Mazur Tate Conjecture is valid in analytic rank one.


Semi automatic System for Title Construction. In this paper, we propose a semi automatic system for title construction from scientific abstracts. The system extracts and recommends impactful words from the text, which the author can creatively use to construct an appropriate title for the manuscript. The work is based on the hypothesis that keywords are good candidates for title construction. We extract important words from the document by inducing a supervised keyword extraction model. The model is trained on novel features extracted from graph of text representation of the document. We empirically show that these graph based features are capable of discriminating keywords from non keywords. We further establish empirically that the proposed approach can be applied to any text irrespective of the training domain and corpus. We evaluate the proposed system by computing the overlap between extracted keywords and the list of title words for documents, and we observe a macro averaged precision of 82 .


II 20 Intelligent and pragmatic analytic categorization of image collections. We introduce II 20 (Image Insight 2020), a multimedia analytics approach for analytic categorization of image collections. Advanced visualizations for image collections exist, but they need tight integration with a machine model to support analytic categorization. Directly employing computer vision and interactive learning techniques gravitates towards search. Analytic categorization, however, is not machine classification (the difference between the two is called the pragmatic gap) a human adds redefines deletes categories of relevance on the fly to build insight, whereas the machine classifier is rigid and non adaptive. Analytic categorization that brings the user to insight requires a flexible machine model that allows dynamic sliding on the exploration search axis, as well as semantic interactions. II 20 brings 3 major contributions to multimedia analytics on image collections and towards closing the pragmatic gap. Firstly, a machine model that closely follows the user's interactions and dynamically models her categories of relevance. II 20's model, in addition to matching and exceeding the state of the art w. r. t. relevance, allows the user to dynamically slide on the exploration search axis without additional input from her side. Secondly, the dynamic, 1 image at a time Tetris metaphor that synergizes with the model. It allows the model to analyze the collection by itself with minimal interaction from the user and complements the classic grid metaphor. Thirdly, the fast forward interaction, allowing the user to harness the model to quickly expand ("fast forward") the categories of relevance, expands the multimedia analytics semantic interaction dictionary. Automated experiments show that II 20's model outperforms the state of the art and also demonstrate Tetris's analytic quality. User studies confirm that II 20 is an intuitive, efficient, and effective multimedia analytics tool.


Synthesis of Boolean Networks from Biological Dynamical Constraints using Answer Set Programming. Boolean networks model finite discrete dynamical systems with complex behaviours. The state of each component is determined by a Boolean function of the state of (a subset of) the components of the network. This paper addresses the synthesis of these Boolean functions from constraints on their domain and emerging dynamical properties of the resulting network. The dynamical properties relate to the existence and absence of trajectories between partially observed configurations, and to the stable behaviours (fixpoints and cyclic attractors). The synthesis is expressed as a Boolean satisfiability problem relying on Answer Set Programming with a parametrized complexity, and leads to a complete non redundant characterization of the set of solutions. Considered constraints are particularly suited to address the synthesis of models of cellular differentiation processes, as illustrated on a case study. The scalability of the approach is demonstrated on random networks with scale free structures up to 100 to 1,000 nodes depending on the type of constraints.


More Hierarchy in Route Planning Using Edge Hierarchies. A highly successful approach to route planning in networks (particularly road networks) is to identify a hierarchy in the network that allows faster queries after some preprocessing that basically inserts additional "shortcut" edges into a graph. In the past there has been a succession of techniques that infer a more and more fine grained hierarchy enabling increasingly more efficient queries. This appeared to culminate in contraction hierarchies that assign one hierarchy level to each vertex. In this paper we show how to identify an even more fine grained hierarchy that assigns one level to each edge of the network. Our findings indicate that this can lead to considerably smaller search spaces in terms of visited edges. Currently, this rarely implies improved query times so that it remains an open question whether edge hierarchies can lead to consistently improved performance. However, we believe that the technique as such is a noteworthy enrichment of the portfolio of available techniques that might prove useful in the future.


On abstract F systems. A graph theoretic model for paradoxes involving a falsity predicate and its application to argumentation frameworks. F systems are digraphs that enable to model sentences that predicate the falsity of other sentences. Paradoxes like the Liar and Yablo's can be analyzed with that tool to find graph theoretic patterns. In this paper we present the F systems model abstracting from all the features of the language in which the represented sentences are expressed. All that is assumed is the existence of sentences and the binary relation '... affirms the falsity of ...' among them. The possible existence of non referential sentences is also considered. To model the sets of all the sentences that can jointly be valued as true we introduce the notion of conglomerate, the existence of which guarantees the absence of paradox. Conglomerates also enable to characterize referential contradictions, i.e. sentences that can only be false under a classical valuation due to the interactions with other sentences in the model. A Kripke's style fixed point characterization of groundedness is offered and fixed points which are complete (meaning that every sentence is deemed either true or false) and consistent (meaning that no sentence is deemed true and false) are put in correspondence with conglomerates. Furthermore, argumentation frameworks are special cases of F systems. We show the relation between local conglomerates and admissible sets of arguments and argue about the usefulness of the concept for argumentation theory.


Vaporization frequency response to pressure oscillations an approximate analytical solution for mixed injection regimes. This work is devoted to a theoretical analysis of mass frequency response to pressure oscillations of a spray of repetitively injected drops into a combustion chamber. A single stationary spherical droplet continuously fed with the same liquid fuel so that its volume remains constant despite the evaporation, the so called 'mean droplet' in the Heidmann analogy, represents this vaporizing spray. The feeding is realized with a liquid liquid heat transfer coefficient by using a source point placed at the mean droplet centre, in such a way that only radial thermal convection and conduction effects are allowed inside the droplet during the process. This feeding procedure is now viewed as a proper boundary condition that is a mixed or a generalized feeding regime controlling the liquid fuel injection process into the combustion chamber. Drawing upon a linear analysis based on the Rayleigh criterion, the evaporating mass response factor is evaluated. Effects due to the variation of the heat transfer coefficient and that of the process characteristic times are analysed. An abrupt increase appears in the response when a fuel thermodynamic coefficient approaches a particular value.


A fast memoryless predictive algorithm in a chain of recurrent neural networks. In the recent publication (arxiv 2007.08063v2 cs.LG ) a fast prediction algorithm for a single recurrent network (RN) was suggested. In this manuscript we generalize this approach to a chain of RNs and show that it can be implemented in natural neural systems. When the network is used recursively to predict sequence of values the proposed algorithm does not require to store the original input sequence. It increases robustness of the new approach compared to the standard moving expanding window predictive procedure. We consider requirements on trained networks that allow to implement the proposed algorithm and discuss them in the neuroscience context.


Extension of order bounded operators. Assume that a normed lattice E is order dense majorizing of a vector lattice E t . There is an extension norm Vert. Vert t for E t and we extend some lattice and topological properties of normed lattice (E, Vert. Vert) to new normed lattice (E t, Vert. Vert t ). For a Dedekind complete Banach lattice F , T t is an extension of T from E t into F whenever T is an order bounded operator from E into F . For each positive operator T , we have Vert T Vert Vert T t Vert and we show that T t is a lattice homomorphism from E t into F and moreover T t in mathcal L n(E t,F) whenever 0 leq T in mathcal L n(E,F) and T(x wedge y) Tx wedge Ty for each 0 leq x,y in E . We also extend some lattice and topological properties of T in mathcal L b(E,F) to the extension operator T t in mathcal L b(E t,F) .


On unitriangular basic sets for symmetric and alternating groups. We study the modular representation theory of the symmetric and alternating groups. One of the most natural ways to label the irreducible representations of a given group or algebra in the modular case is to show the unitriangularity of the decomposition matrices, that is, the existence of a unitriangular basic set. We study several ways to obtain such sets in the general case of a symmetric algebra. We apply our results to the symmetric groups and to their Hecke algebras and thus obtain new ways to label the simple modules for these objects. Finally, we show that these sets do not always exist in the case of the alternating groups by studying two explicit cases in characteristic 3.


Deep Differentiable Grasp Planner for High DOF Grippers. We present an end to end algorithm for training deep neural networks to grasp novel objects. Our algorithm builds all the essential components of a grasping system using a forward backward automatic differentiation approach, including the forward kinematics of the gripper, the collision between the gripper and the target object, and the metric for grasp poses. In particular, we show that a generalized Q1 grasp metric is defined and differentiable for inexact grasps generated by a neural network, and the derivatives of our generalized Q1 metric can be computed from a sensitivity analysis of the induced optimization problem. We show that the derivatives of the (self )collision terms can be efficiently computed from a watertight triangle mesh of low quality. Altogether, our algorithm allows for the computation of grasp poses for high DOF grippers in an unsupervised mode with no ground truth data, or it improves the results in a supervised mode using a small dataset. Our new learning algorithm significantly simplifies the data preparation for learning based grasping systems and leads to higher qualities of learned grasps on common 3D shape datasets 7, 49, 26, 25 , achieving a 22 higher success rate on physical hardware and a 0.12 higher value on the Q1 grasp quality metric.


Robust Policy Synthesis for Uncertain POMDPs via Convex Optimization. We study the problem of policy synthesis for uncertain partially observable Markov decision processes (uPOMDPs). The transition probability function of uPOMDPs is only known to belong to a so called uncertainty set, for instance in the form of probability intervals. Such a model arises when, for example, an agent operates under information limitation due to imperfect knowledge about the accuracy of its sensors. The goal is to compute a policy for the agent that is robust against all possible probability distributions within the uncertainty set. In particular, we are interested in a policy that robustly ensures the satisfaction of temporal logic and expected reward specifications. We state the underlying optimization problem as a semi infinite quadratically constrained quadratic program (QCQP), which has finitely many variables and infinitely many constraints. Since QCQPs are non convex in general and practically infeasible to solve, we resort to the so called convex concave procedure to convexify the QCQP. Even though convex, the resulting optimization problem still has infinitely many constraints and is NP hard. For uncertainty sets that form convex polytopes, we provide a transformation of the problem to a convex QCQP with finitely many constraints. We demonstrate the feasibility of our approach by means of several case studies that highlight typical bottlenecks for our problem. In particular, we show that we are able to solve benchmarks with hundreds of thousands of states, hundreds of different observations, and we investigate the effect of different levels of uncertainty in the models.


Computer Aided Generation of N shift RWS. Generating schedules for shift workers is essential for many employers, whether the employer is a small or a large industrial complex, research laboratory, or other businesses involving shift works. Previous methods for creating rotational workforce schedules included interactions between the schedule maker and the algorithm, including defining the length of sequences of consecutive days of working shifts. In this method, an algorithm takes into account inputs (or constraints) from the schedule maker and then presents the possible solutions (incl. that all shifts must be filled, working hours per week, minimal resting time, etc.) in a first phase. The schedule maker can then select which solutions are most feasible to proceed with in the second phase, where the final schedules are then constructed and exported.


Revisiting Reynolds and Nusselt numbers in turbulent thermal convection. In this paper, we extend Grossmann and Lohse's (GL) model Phys. Rev. Lett. bf 86 , 3316 (2001) for the predictions of Reynolds number (Re) and Nusselt number (Nu) in turbulent Rayleigh B ' e nard convection (RBC). Towards this objective, we use functional forms for the prefactors of the dissipation rates in the bulk and the boundary layers. The functional forms arise due to inhibition of nonlinear interactions in the presence of walls and buoyancy compared to free turbulence, along with a deviation of viscous boundary layer profile from Prandtl Blasius theory. We perform 60 numerical runs on a three dimensional unit box for a range of Rayleigh numbers (Ra) and Prandtl numbers (Pr) and determine the aforementioned functional forms using machine learning. The revised predictions are in better agreement with the past numerical and experimental results than those of the GL model, especially for extreme Prandtl numbers.


Examples of reflective projective billiards and outer ghost billiards. In the class of projective billiards, which contains the usual billiards, we exhibit counter examples to Ivrii's conjecture, which states that in any planar billiard with smooth boundary the set of periodic orbits has zero measure. The counter examples are polygons admitting a 2 parameters family of n periodic orbits, with n being either 3 or any even number grower than 4 .


Refining 6D Object Pose Predictions using Abstract Render and Compare. Robotic systems often require precise scene analysis capabilities, especially in unstructured, cluttered situations, as occurring in human made environments. While current deep learning based methods yield good estimates of object poses, they often struggle with large amounts of occlusion and do not take inter object effects into account. Vision as inverse graphics is a promising concept for detailed scene analysis. A key element for this idea is a method for inferring scene parameter updates from the rasterized 2D scene. However, the rasterization process is notoriously difficult to invert, both due to the projection and occlusion process, but also due to secondary effects such as lighting or reflections. We propose to remove the latter from the process by mapping the rasterized image into an abstract feature space learned in a self supervised way from pixel correspondences. Using only a light weight inverse rendering module, this allows us to refine 6D object pose estimations in highly cluttered scenes by optimizing a simple pixel wise difference in the abstract image representation. We evaluate our approach on the challenging YCB Video dataset, where it yields large improvements and demonstrates a large basin of attraction towards the correct object poses.


Classification of finite Morse index solutions of higher order Gelfand Liouville equation. We classify finite Morse index solutions of the following Gelfand Liouville equation begin equation ( Delta) s u e u text in mathbb R n, end equation for 1 s 2 and s 2 via a novel monotonicity formula and technical blow down analysis. We show that the above equation does not admit any finite Morse index solution with ( Delta) s 2 u vanishes at infinity provided n 2s and begin equation label 1.condition frac Gamma( frac n 2s 4 ) 2 Gamma( frac n 2s 4 ) 2 frac Gamma( frac n 2 ) Gamma(1 s) Gamma( frac n 2s 2 ) , end equation where Gamma is the classical Gamma function. The cases of s 1 and s 2 are settled by Dancer and Farina cite df,d and Dupaigne et al. cite dggw , respectively, using Moser iteration arguments established by Crandall and Rabinowitz cite CR . The case of 0 s 1 is established by Hyder Yang in cite hy applying arguments provided in cite ddw,fw .


Multi Task Reinforcement Learning based Mobile Manipulation Control for Dynamic Object Tracking and Grasping. Agile control of mobile manipulator is challenging because of the high complexity coupled by the robotic system and the unstructured working environment. Tracking and grasping a dynamic object with a random trajectory is even harder. In this paper, a multi task reinforcement learning based mobile manipulation control framework is proposed to achieve general dynamic object tracking and grasping. Several basic types of dynamic trajectories are chosen as the task training set. To improve the policy generalization in practice, random noise and dynamics randomization are introduced during the training process. Extensive experiments show that our policy trained can adapt to unseen random dynamic trajectories with about 0.1m tracking error and 75 grasping success rate of dynamic objects. The trained policy can also be successfully deployed on a real mobile manipulator.


An institutional study on plan quality and variation of manual forward planning for Gamma Knife radiosurgery for vestibular schwannoma. Due to the complexity and cumbersomeness of Gamma Knife (GK) manual forward planning, the quality of the resulting treatment plans heavily depends on the planners skill, experience and the amount of effort devoted to plan development. Hence, GK plan quality may vary significantly among institutions and planners, and even for a same planner at different cases. This is particularly a concern for challenging cases with complicated geometry, such as vestibular schwannoma cases. The purpose of this retrospective study is to investigate the plan quality and variation in the manually forward planned, clinically acceptable GK treatment plans of 22 previous vestibular schwannoma cases. Considering the impacts of different patient geometry and different trade offs among the planning objectives in GK planning, it is difficult to objectively assess the plan quality across different cases. To reduce these confounding factors on plan quality assessment, we employed our recently developed multiresolution level inverse planning algorithm to generate a golden plan for each case, which is expected to be on or close to the pareto surface with a similar trade off as used in the manual plan. The plan quality of the manual plan is then quantified in terms of its deviation from the golden plan. A scoring criterion between 0 100 was designed to calculate a final score for each manual plan to simplify our analysis. Large quality variation was observed in these 22 cases, with two cases having a score lower than 75, three cases scoring between 80 and 85, two cases between 85 and 90, eight cases between 90 and 95, and seven cases higher than 95. Inter and intra planner variability was also observed in our study. This large variation in GK manual planning deserves high attention, and merits further investigation on how to reduce the variation in GK treatment plan quality.


Capacity and Optimal Resource Allocation for IRS assisted Multi user Communication Systems. The fundamental capacity limits of intelligent reflecting surface (IRS) assisted multi user wireless communication systems are investigated in this paper. Specifically, the capacity and rate regions for both capacity achieving non orthogonal multiple access (NOMA) and orthogonal multiple access (OMA) transmission schemes are characterized under the constraints of a maximum number of IRS reconfiguration times. For NOMA, the ideal case with asymptotically large number of IRS reconfiguration times is considered, where the optimal solution is obtained by employing the Lagrange duality method. Inspired by this result, an inner bound of the capacity region for the general case with a finite number of IRS reconfiguration times is derived. For OMA, the optimal transmission strategy for the ideal case is to serve each individual user alternatingly with its effective channel power gain maximized. Based on this result, a rate region inner bound for the general case is derived. Finally, numerical results are provided to show that i) a significant capacity and rate region improvement can be achieved by using IRS ii) the capacity gain can be further improved by dynamically configuring the IRS reflection matrix.


On a limit theorem for a non linear scaling. In this note, we proved that weak limits, of sums of independent positive identically distributed random variables which are re normalized by a non linear shrinking transform max(0, x r) , are either degenerate or (some) compound Poisson distributions.


A new splitting algorithm for dynamical low rank approximation motivated by the fibre bundle structure of matrix manifolds. In this paper, we propose a new splitting algorithm for dynamical low rank approximation motivated by the fibre bundle structure of the set of fixed rank matrices. We first introduce a geometric description of the set of fixed rank matrices which relies on a natural parametrization of matrices. More precisely, it is endowed with the structure of analytic principal bundle, with an explicit description of local charts. For matrix differential equations, we introduce a first order numerical integrator working in local coordinates. The resulting algorithm can be interpreted as a particular splitting of the projection operator onto the tangent space of the low rank matrix manifold. It is proven to be exact in some particular case. Numerical experiments confirm this result and illustrate the robustness of the proposed algorithm.


Fragility of Layout Based and Visual GUI Test Scripts An Assessment Study on a Hybrid Mobile Application. Context Albeit different approaches exist for automated GUI testing of hybrid mobile applications, the practice appears to be not so commonly adopted by developers. A possible reason for such a low diffusion can be the fragility of the techniques, i.e. the frequent need for maintaining test cases when the GUI of the app is changed. Goal In this paper, we perform an assessment of the maintenance needed by test cases for a hybrid mobile app, and the related fragility causes. Methods We evaluated a small test suite with a Layout based testing tool (Appium) and a Visual one (EyeAutomate) and observed the changes needed by tests during the co evolution with the GUI of the app. Results We found that 20 Layout based test methods and 30 Visual test methods had to be modified at least once, and that each release induced fragilities in 3 4 of the test methods. Conclusion Fragility of GUI tests can induce relevant maintenance efforts in test suites of large applications. Several principal causes for fragilities have been identified for the tested hybrid application, and guidelines for developers are deduced from them.


A Novel DNN Training Framework via Data Sampling and Multi Task Optimization. Conventional DNN training paradigms typically rely on one training set and one validation set, obtained by partitioning an annotated dataset used for training, namely gross training set, in a certain way. The training set is used for training the model while the validation set is used to estimate the generalization performance of the trained model as the training proceeds to avoid over fitting. There exist two major issues in this paradigm. Firstly, the validation set may hardly guarantee an unbiased estimate of generalization performance due to potential mismatching with test data. Secondly, training a DNN corresponds to solve a complex optimization problem, which is prone to getting trapped into inferior local optima and thus leads to undesired training results. To address these issues, we propose a novel DNN training framework. It generates multiple pairs of training and validation sets from the gross training set via random splitting, trains a DNN model of a pre specified structure on each pair while making the useful knowledge (e.g., promising network parameters) obtained from one model training process to be transferred to other model training processes via multi task optimization, and outputs the best, among all trained models, which has the overall best performance across the validation sets from all pairs. The knowledge transfer mechanism featured in this new framework can not only enhance training effectiveness by helping the model training process to escape from local optima but also improve on generalization performance via implicit regularization imposed on one model training process from other model training processes. We implement the proposed framework, parallelize the implementation on a GPU cluster, and apply it to train several widely used DNN models. Experimental results demonstrate the superiority of the proposed framework over the conventional training paradigm.


Subterranean MAV Navigation based on Nonlinear MPC with Collision Avoidance Constraints. Micro Aerial Vehicles (MAVs) navigation in subterranean environments is gaining attention in the field of aerial robotics, however there are still multiple challenges for collision free navigation in such harsh environments. This article proposes a novel baseline solution for collision free navigation with Nonlinear Model Predictive Control (NMPC). In the proposed method, the MAV is considered as a floating object, where the velocities on the x , y axes and the position on altitude are the references for the NMPC to navigate along the tunnel, while the NMPC avoids the collision by considering kinematics of the obstacles based on measurements from a 2D lidar. Moreover, a novel approach for correcting the heading of the MAV towards the center of the mine tunnel is proposed, while the efficacy of the suggested framework has been evaluated in multiple field trials in an underground mine in Sweden.


Stable connectivity over a base. Morel's stable connectivity theorems state that for any connective S 1 spectrum F of motivic spaces (Nisnevich simplicial sheaves) over an arbitrary field, the spectrum L mathbb A 1 (F) is connective, and the same property for mathbb P 1 spectra of motivic spaces. Here L mathbb A 1 denotes the mathbb A 1 localisation in the category of motivic spectra over a field k . Originally the same property was conjectured for the case of motivic S 1 spectra over a base scheme S .In view of Ayoub's conterexamples the modified version of conjecture states that L mathbb A 1 (F) is ( d) connective for any connective F , where d mathrm dim S is the Krull dimension. The conjecture is proven under the infiniteness assumption on the residue fields for the cases of Dedekind schemes by J. Schmidt and F. Strunk and noetherian domains of arbitrary dimension by N. Deshmukh, A. Hogadi, G. Kulkarni and S. Yadavand. In the article we prove the result or general base with out the assumption on the residue fields. So by the result for any smooth scheme X over a base scheme S of Krull dimension d the Nisnevich sheaves of S 1 stable motivic homotopy groups pi i S 1 (X) and mathbb P 1 stable motivic homotopy groups pi i j,j mathbb P 1 (X) vanishes for all i d .


Statistical properties of periodic points for infinitely renormalizable unimodal maps. For an infinitely renormalizable negative Schwarzian unimodal map f with non flat critical point, we analyze statistical properties of periodic points as the periods tend to infinity. Introducing a weight function varphi which is a continuous or a geometric potential varphi beta log f' ( beta in mathbb R ), we establish the level 2 Large Deviation Principle for weighted periodic points. From this, we deduce that all weighted periodic points equidistribute with respect to equilibrium states for the potential varphi . In particular, it follows that all periodic points are equidistributed with respect to measures of maximal entropy, and all periodic points weighted with their Lyapunov exponents are equidistributed with respect to the post critical measure supported on the attracting Cantor set.


Model order reduction of flow based on a modular geometrical approximation of blood vessels. We are interested in a reduced order method for the efficient simulation of blood flow in arteries. The blood dynamics is modeled by means of the incompressible Navier Stokes equations. Our algorithm is based on an approximated domain decomposition of the target geometry into a number of subdomains obtained from the parametrized deformation of geometrical building blocks (e.g. straight tubes and model bifurcations). On each of these building blocks, we build a set of spectral functions by proper orthogonal decomposition of a large number of snapshots of finite element solutions (offline phase). The global solution of the Navier Stokes equations on a target geometry is then found by coupling linear combinations of these local basis functions by means of spectral Lagrange multipliers (online phase). Being that the number of reduced degrees of freedom is considerably smaller than their finite element counterpart, this approach allows us to significantly decrease the size of the linear system to be solved in each iteration of the Newton Raphson algorithm. We achieve large speedups with respect to the full order simulation (in our numerical experiments, the gain is at least of one order of magnitude and grows inversely with respect to the reduced basis size), whilst still retaining satisfactory accuracy for most cardiovascular simulations.


The Generalized Independent and Dominating Set Problems on Unit Disk Graphs. In this article, we study a generalized version of the maximum independent set and minimum dominating set problems, namely, the maximum d distance independent set problem and the minimum d distance dominating set problem on unit disk graphs for a positive integer d 0 . We first show that the maximum d distance independent set problem and the minimum d distance dominating set problem belongs to NP hard class. Next, we propose a simple polynomial time constant factor approximation algorithms and PTAS for both the problems.


3D Pose Detection in Videos Focusing on Occlusion. In this work, we build upon existing methods for occlusion aware 3D pose detection in videos. We implement a two stage architecture that consists of the stacked hourglass network to produce 2D pose predictions, which are then inputted into a temporal convolutional network to produce 3D pose predictions. To facilitate prediction on poses with occluded joints, we introduce an intuitive generalization of the cylinder man model used to generate occlusion labels. We find that the occlusion aware network is able to achieve a mean per joint position error 5 mm less than our linear baseline model on the Human3.6M dataset. Compared to our temporal convolutional network baseline, we achieve a comparable mean per joint position error of 0.1 mm less at reduced computational cost.


Filtering Rules for Flow Time Minimization in a Parallel Machine Scheduling Problem. This paper studies the scheduling of jobs of different families on parallel machines with qualification constraints. Originating from semiconductor manufacturing, this constraint imposes a time threshold between the execution of two jobs of the same family. Otherwise, the machine becomes disqualified for this family. The goal is to minimize both the flow time and the number of disqualifications. Recently, an efficient constraint programming model has been proposed. However, when priority is given to the flow time objective, the efficiency of the model can be improved. This paper uses a polynomial time algorithm which minimize the flow time for a single machine relaxation where disqualifications are not considered. Using this algorithm one can derived filtering rules on different variables of the model. Experimental results are presented showing the effectiveness of these rules. They improve the competitiveness with the mixed integer linear program of the literature.


The complexity of 3 colouring H colourable graphs. We study the complexity of approximation on satisfiable instances for graph homomorphism problems. For a fixed graph H , the H colouring problem is to decide whether a given graph has a homomorphism to H . By a result of Hell and Ne v s et v r il, this problem is NP hard for any non bipartite graph H . In the context of promise constraint satisfaction problems, Brakensiek and Guruswami conjectured that this hardness result extends to promise graph homomorphism as follows fix any non bipartite graph H and another graph G with a homomorphism from H to G , it is NP hard to find a homomorphism to G from a given H colourable graph. Arguably, the two most important special cases of this conjecture are when H is fixed to be the complete graph on 3 vertices (and G is any graph with a triangle) and when G is the complete graph on 3 vertices (and H is any 3 colourable graph). The former case is equivalent to the notoriously difficult approximate graph colouring problem. In this paper, we confirm the Brakensiek Guruswami conjecture for the latter case. Our proofs rely on a novel combination of the universal algebraic approach to promise constraint satisfaction, that was recently developed by Barto, Bul 'in and the authors, with some ideas from algebraic topology.


AutoEmb Automated Embedding Dimensionality Search in Streaming Recommendations. Deep learning based recommender systems (DLRSs) often have embedding layers, which are utilized to lessen the dimensionality of categorical variables (e.g. user item identifiers) and meaningfully transform them in the low dimensional space. The majority of existing DLRSs empirically pre define a fixed and unified dimension for all user item embeddings. It is evident from recent researches that different embedding sizes are highly desired for different users items according to their popularity. However, manually selecting embedding sizes in recommender systems can be very challenging due to the large number of users items and the dynamic nature of their popularity. Thus, in this paper, we propose an AutoML based end to end framework (AutoEmb), which can enable various embedding dimensions according to the popularity in an automated and dynamic manner. To be specific, we first enhance a typical DLRS to allow various embedding dimensions then we propose an end to end differentiable framework that can automatically select different embedding dimensions according to user item popularity finally we propose an AutoML based optimization algorithm in a streaming recommendation setting. The experimental results based on widely used benchmark datasets demonstrate the effectiveness of the AutoEmb framework.


Euclidean Geometry and Elliptic Curves. In this paper, we demonstrate the intimate relationships among some geometric figures and the families of elliptic curves with positive ranks. These geometric figures include textit textbf Heron triangles , textit textbf Brahmagupta quadrilaterals and textit textbf Bicentric quadrilaterals . Firstly, we investigate the important properties of these figures and then utilizing these properties, we show that how to construct various families of elliptic curves with different positive ranks having different torsion subgroups.


An efficient Kinetic Monte Carlo to study analyte capture by a nanopore Transients, boundary conditions and time dependent fields. To better understand the capture process by a nanopore, we introduce an efficient Kinetic Monte Carlo (KMC) algorithm that can simulate long times and large system sizes by mapping the dynamic of a point like particle in a 3D spherically symmetric system onto the 1D biased random walk. Our algorithm recovers the steady state analytical solution and allows us to study time dependent processes such as transients. Simulation results show that the steady state depletion zone near pore is barely larger than the pore radius and narrows at higher field intensities as a result, the time to reach steady state is much smaller than the time required to empty a zone of the size of the capture radius lambda e . When the sample reservoir has a finite size, a second depletion region propagates inward from the outer wall, and the capture rate starts decreasing when it reaches the capture radius lambda e . We also note that the flatness of the electric field near the pore, which is often neglected, induces a traffic jam that can increase the transient time by several orders of magnitude. Finally, we propose a new proof of concept scheme to separate two analytes of the same mobility but different diffusion coefficients using time varying fields.


CM Method and Expansion of Numbers. We show that an iterative method for computing the center of mass (CM) of q units of mass, placed on a unit interval 0, 1 along the x axis, give rise to a simple procedure for expanding rational numbers less than unity in powers of r s 1, with r, s, integers larger than 0. The method is then extended to all numbers, real or complex, though the procedure for none rational numbers is more time consuming. We also show how our method provides a natural way to generalize Jacobsthal numbers. Moreover, the method provides a way to generate infinitely many sequences of numbers, of which many play an important rule in mathematical sciences and engineering, to name few, Jacobsthal sequence, Fibonacci sequence, and Pell sequence.


MIMICS A Large Scale Data Collection for Search Clarification. Search clarification has recently attracted much attention due to its applications in search engines. It has also been recognized as a major component in conversational information seeking systems. Despite its importance, the research community still feels the lack of a large scale data for studying different aspects of search clarification. In this paper, we introduce MIMICS, a collection of search clarification datasets for real web search queries sampled from the Bing query logs. Each clarification in MIMICS is generated by a Bing production algorithm and consists of a clarifying question and up to five candidate answers. MIMICS contains three datasets (1) MIMICS Click includes over 400k unique queries, their associated clarification panes, and the corresponding aggregated user interaction signals (i.e., clicks). (2) MIMICS ClickExplore is an exploration data that includes aggregated user interaction signals for over 60k unique queries, each with multiple clarification panes. (3) MIMICS Manual includes over 2k unique real search queries. Each query clarification pair in this dataset has been manually labeled by at least three trained annotators. It contains graded quality labels for the clarifying question, the candidate answer set, and the landing result page for each candidate answer. MIMICS is publicly available for research purposes, thus enables researchers to study a number of tasks related to search clarification, including clarification generation and selection, user engagement prediction for clarification, click models for clarification, and analyzing user interactions with search clarification.


Uniform estimates for the planning problem with potential. In this paper, we study a priori estimates for a first order mean field planning problem with a potential. In the theory of mean field games (MFGs), a priori estimates play a crucial role to prove the existence of classical solutions. In particular, uniform bounds for the density of players' distribution and its inverse are of utmost importance. Here, we investigate a priori bounds for those quantities for a planning problem with a non vanishing potential. The presence of a potential raises non trivial difficulties, which we overcome by exploring a displacement convexity property for the mean field planning problem with a potential together with Moser's iteration method. We show that if the potential satisfies a certain smallness condition, then a displacement convexity property holds. This property enables L q bounds for the density. In the one dimensional case, the displacement convexity property also gives L q bounds for the inverse of the density. Finally, using these L q estimates and Moser's iteration method, we obtain L infty estimates for the density of the distribution of the players and its inverse.