URIE Universal Image Enhancement for Visual Recognition in the Wild. Despite the great advances in visual recognition, it has been witnessed that recognition models trained on clean images of common datasets are not robust against distorted images in the real world. To tackle this issue, we present a Universal and Recognition friendly Image Enhancement network, dubbed URIE, which is attached in front of existing recognition models and enhances distorted input to improve their performance without retraining them. URIE is universal in that it aims to handle various factors of image degradation and to be incorporated with any arbitrary recognition models. Also, it is recognition friendly since it is optimized to improve the robustness of following recognition models, instead of perceptual quality of output image. Our experiments demonstrate that URIE can handle various and latent image distortions and improve the performance of existing models for five diverse recognition tasks when input images are degraded.


An On Line POMDP Solver for Continuous Observation Spaces. Planning under partial obervability is essential for autonomous robots. A principled way to address such planning problems is the Partially Observable Markov Decision Process (POMDP). Although solving POMDPs is computationally intractable, substantial advancements have been achieved in developing approximate POMDP solvers in the past two decades. However, computing robust solutions for problems with continuous observation spaces remains challenging. Most on line solvers rely on discretising the observation space or artificially limiting the number of observations that are considered during planning to compute tractable policies. In this paper we propose a new on line POMDP solver, called Lazy Belief Extraction for Continuous POMDPs (LABECOP), that combines methods from Monte Carlo Tree Search and particle filtering to construct a policy reprentation which doesn't require discretised observation spaces and avoids limiting the number of observations considered during planning. Experiments on three different problems involving continuous observation spaces indicate that LABECOP performs similar or better than state of the art POMDP solvers.


Let's Take This Online Adapting Scene Coordinate Regression Network Predictions for Online RGB D Camera Relocalisation. Many applications require a camera to be relocalised online, without expensive offline training on the target scene. Whilst both keyframe and sparse keypoint matching methods can be used online, the former often fail away from the training trajectory, and the latter can struggle in textureless regions. By contrast, scene coordinate regression (SCoRe) methods generalise to novel poses and can leverage dense correspondences to improve robustness, and recent work has shown how to adapt SCoRe forests between scenes, allowing their state of the art performance to be leveraged online. However, because they use features hand crafted for indoor use, they do not generalise well to harder outdoor scenes. Whilst replacing the forest with a neural network and learning suitable features for outdoor use is possible, the techniques used to adapt forests between scenes are unfortunately harder to transfer to a network context. In this paper, we address this by proposing a novel way of leveraging a network trained on one scene to predict points in another scene. Our approach replaces the appearance clustering performed by the branching structure of a regression forest with a two step process that first uses the network to predict points in the original scene, and then uses these predicted points to look up clusters of points from the new scene. We show experimentally that our online approach achieves state of the art performance on both the 7 Scenes and Cambridge Landmarks datasets, whilst running in under 300ms, making it highly effective in live scenarios.


A perspective on Microscopy Metadata data provenance and quality control. The application of microscopy in biomedical research has come a long way since Antonie van Leeuwenhoek discovered unicellular organisms. Countless innovations have positioned light microscopy as a cornerstone of modern biology and a method of choice for connecting omics datasets to their biological and clinical correlates. Still, regardless of how convincing published imaging data looks, it does not always convey meaningful information about the conditions in which it was acquired, processed, and analyzed. Adequate record keeping, reporting, and quality control are therefore essential to ensure experimental rigor and data fidelity, allow experiments to be reproducibly repeated, and promote the proper evaluation, interpretation, comparison, and re use. To this end, microscopy images should be accompanied by complete descriptions detailing experimental procedures, biological samples, microscope hardware specifications, image acquisition parameters, and image analysis procedures, as well as metrics accounting for instrument performance and calibration. However, universal, community accepted Microscopy Metadata standards and reporting specifications that would result in Findable Accessible Interoperable and Reproducible (FAIR) microscopy data have not yet been established. To understand this shortcoming and to propose a way forward, here we provide an overview of the nature of microscopy metadata and its importance for fostering data quality, reproducibility, scientific rigor, and sharing value in light microscopy. The proposal for tiered Microscopy Metadata Specifications that extend the OME Data Model put forth by the 4D Nucleome Initiative and by Bioimaging North America 1 3 as well as a suite of three complementary and interoperable tools are being developed to facilitate the process of image data documentation and are presented in related manuscripts 4 6 .


Color Constancy Convolutional Autoencoder. In this paper, we study the importance of pre training for the generalization capability in the color constancy problem. We propose two novel approaches based on convolutional autoencoders an unsupervised pre training algorithm using a fine tuned encoder and a semi supervised pre training algorithm using a novel composite loss function. This enables us to solve the data scarcity problem and achieve competitive, to the state of the art, results while requiring much fewer parameters on ColorChecker RECommended dataset. We further study the over fitting phenomenon on the recently introduced version of INTEL TUT Dataset for Camera Invariant Color Constancy Research, which has both field and non field scenes acquired by three different camera models.


Multi 2 OIE Multilingual Open Information Extraction Based on Multi Head Attention with BERT. In this paper, we propose Multi 2 OIE, which performs open information extraction (open IE) by combining BERT with multi head attention. Our model is a sequence labeling system with an efficient and effective argument extraction method. We use a query, key, and value setting inspired by the Multimodal Transformer to replace the previously used bidirectional long short term memory architecture with multi head attention. Multi 2 OIE outperforms existing sequence labeling systems with high computational efficiency on two benchmark evaluation datasets, Re OIE2016 and CaRB. Additionally, we apply the proposed method to multilingual open IE using multilingual BERT. Experimental results on new benchmark datasets introduced for two languages (Spanish and Portuguese) demonstrate that our model outperforms other multilingual systems without training data for the target languages.


Multi modal Experts Network for Autonomous Driving. End to end learning from sensory data has shown promising results in autonomous driving. While employing many sensors enhances world perception and should lead to more robust and reliable behavior of autonomous vehicles, it is challenging to train and deploy such network and at least two problems are encountered in the considered setting. The first one is the increase of computational complexity with the number of sensing devices. The other is the phenomena of network overfitting to the simplest and most informative input. We address both challenges with a novel, carefully tailored multi modal experts network architecture and propose a multi stage training procedure. The network contains a gating mechanism, which selects the most relevant input at each inference time step using a mixed discrete continuous policy. We demonstrate the plausibility of the proposed approach on our 1 6 scale truck equipped with three cameras and one LiDAR.


q-bio.BM physics.bio-ph physics.chem-ph physics.comp-ph

Estimation of Protein Ligand Unbinding Kinetics Using Non Equilibrium Targeted Molecular Dynamics Simulations. We here report on non equilibrium targeted Molecular Dynamics simulations as tool for the estimation of protein ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X ray crystallography on two small molecule compound libraries bound to the N terminal domain of the chaperone Hsp90, we show that the mean non equilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and or van der Waals interactions with the target. In the case of ligands with rigid chemical scaffold that exhibit longer residence times however, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding


Cellularity of centrosymmetric matrix algebras and Frobenius extensions. Centrosymmetric matrices of order n over an arbitrary algebra R form a subalgebra of the full n times n matrix algebra over R . It is called the centrosymmetric matrix algebra of order n over R and denoted by S n(R) . We prove (1) S n(R) is Morita equivalent to S 2(R) if n is even, and to S 3(R) if n ge 3 is odd (2) the full n times n matrix algebra over R is a separable Frobenius extension of S n(R) and (3) if R is a commutative ring, then S n(R) is a cellular R algebra in the sense of Graham Lehrer for all n ge 1 .


Universal p ary designs. We investigate p ary t designs which are simultaneously designs for all t , which we call universal p ary designs. Null universal designs are well understood due to Gordon James via the representation theory of the symmetric group. We study non null designs and determine necessary and sufficient conditions on the coefficients for such a design to exist. This allows us to classify all universal designs, up to similarity.


A codimension 2 component of the Gieseker Petri locus. We show that the Brill Noether locus M 3 18,16 is an irreducible component of the Gieseker Petri locus in genus 18 having codimension 2 in the moduli space of curves. This result disproves a conjecture predicting that the Gieseker Petri locus is always divisorial.


Cyclic framed little disks algebras, Grothendieck Verdier duality and handlebody group representations. We characterize cyclic algebras over the associative and the framed little disks operad in any symmetric monoidal bicategory. The cyclicity is appropriately treated in a homotopy coherent way. When the symmetric monoidal bicategory is specified to be a certain symmetric monoidal bicategory of linear categories subject to finiteness conditions, we prove that cyclic associative and cyclic framed little disks algebras, respectively, are equivalent to pivotal Grothendieck Verdier categories and ribbon Grothendieck Verdier categories, a type of category that was introduced by Boyarchenko Drinfeld based on Barr's notion of a autonomous category. We use these results and Costello's modular envelope construction to obtain two applications to quantum topology I) We extract a consistent system of handlebody group representations from any ribbon Grothendieck Verdier category inside a certain symmetric monoidal bicategory of linear categories and show that this generalizes the handlebody part of Lyubashenko's mapping class group representations. II) We establish a Grothendieck Verdier duality for the category extracted from a modular functor by evaluation on the circle (without any assumption on semisimplicity), thereby generalizing results of Tillmann and Bakalov Kirillov.


Hard X Ray Phase Contrast Laboratory Microscope Based on Three Block Fresnel Zone Plate Interferometer. A device based on a three block Fresnel zone plate interferometer is proposed for hard X ray phase contrast imaging. The device combines a low requirement for the coherence of the initial radiation (the interferometer operates in the amplitude division mode) with an optical magnification of the image. A numerical simulation of the image formation is carried out, taking into account the limited source interferometer distance, the size and spectral width of the X ray source. The calculations show that the proposed set up can be used as a phase contrast microscope using laboratory hard X ray sources.


physics.comp-ph cs.LG physics.data-an physics.geo-ph

Constraining subglacial processes from surface velocity observations using surrogate based Bayesian inference. Basal motion is the primary mechanism for ice flux outside Antarctica, yet a widely applicable model for predicting it in the absence of retrospective observations remains elusive. This is due to the difficulty in both observing small scale bed properties and predicting a time varying water pressure on which basal motion putatively depends. We take a Bayesian approach to these problems by coupling models of ice dynamics and subglacial hydrology and conditioning on observations of surface velocity in southwestern Greenland to infer the posterior probability distributions for eight spatially and temporally constant parameters governing the behavior of both the sliding law and hydrologic model. Because the model is computationally expensive, classical MCMC sampling is intractable. We skirt this issue by training a neural network as a surrogate that approximates the model at a sliver of the computational cost. We find that surface velocity observations establish strong constraints on model parameters relative to a prior distribution and also elucidate correlations, while the model explains 60 of observed variance. However, we also find that several distinct configurations of the hydrologic system and stress regime are consistent with observations, underscoring the need for continued data collection and model development.


Dynamic Term Modal Logic for Epistemic Social Network Dynamics (Extended Version). Logics for social networks have been studied in recent literature. This paper presents a framework based on dynamic term modal logic (DTML), a quantified variant of dynamic epistemic logic (DEL). In contrast with DEL where it is commonly known to whom agent names refer, DTML can represent dynamics with uncertainty about agent identity. We exemplify dynamics where such uncertainty and de re de dicto distinctions are key to social network epistemics. Technically, we show that DTML semantics can represent a popular class of hybrid logic epistemic social network models. We also show that DTML can encode previously discussed dynamics for which finding a complete logic was left open. As complete reduction axioms systems exist for DTML, this yields a complete system for the dynamics in question.


PASSVM A Highly Accurate Online Fast Flux Detection System. Fast Flux service networks (FFSNs) are used by adversaries to achieve a high resilient technique for their malicious servers while keeping them hidden from direct access. In this technique, a large number of botnet machines, that are known as flux agents, work as proxies to relay the traffic between end users and a malicious mothership server which is controlled by an adversary. Various mechanisms have been proposed for detecting FFSNs. Such mechanisms depend on collecting a large amount of DNS traffic traces and require a considerable amount of time to identify fast flux domains. In this paper, we propose an efficient AI based online fast flux detection system that performs highly accurate and extremely fast detection of fast flux domains. The proposed system, called PASSVM, is based on features that are associated with DNS response messages of a given domain name. The approach relies on features that are stored in two local databases, in addition to features that are extracted from the response DNS messages itself. The information in the databases are obtained from Censys search engine and IP Geolocation service. PASSVM is evaluated using three types of artificial neural networks which are Multilayer Perceptron (MLP), Radial Basis Function Network (RBF), and Support Vector Machines (SVM). Results show that SVM with RBF kernel outperformed the other two methods with an accuracy of 99.557 and a detection time of less than 18 ms.


Beyond Graph Neural Networks with Lifted Relational Neural Networks. We demonstrate a declarative differentiable programming framework based on the language of Lifted Relational Neural Networks, where small parameterized logic programs are used to encode relational learning scenarios. When presented with relational data, such as various forms of graphs, the program interpreter dynamically unfolds differentiable computational graphs to be used for the program parameter optimization by standard means. Following from the used declarative Datalog abstraction, this results into compact and elegant learning programs, in contrast with the existing procedural approaches operating directly on the computational graph level. We illustrate how this idea can be used for an efficient encoding of a diverse range of existing advanced neural architectures, with a particular focus on Graph Neural Networks (GNNs). Additionally, we show how the contemporary GNN models can be easily extended towards higher relational expressiveness. In the experiments, we demonstrate correctness and computation efficiency through comparison against specialized GNN deep learning frameworks, while shedding some light on the learning performance of existing GNN models.


Portfolio Optimization with Expectile and Omega Functions. This paper proves equivalences of portfolio optimization problems with negative expectile and omega ratio. We derive subgradients for the negative expectile as a function of the portfolio from a known dual representation of expectile and general theory about subgradients of risk measures. We also give an elementary derivation of the gradient of negative expectile under some assumptions and provide an example where negative expectile is demonstrably not differentiable. We conducted a case study and solved portfolio optimization problems with negative expectile objective and constraint.


A Non commutative Bilinear Model for Answering Path Queries in Knowledge Graphs. Bilinear diagonal models for knowledge graph embedding (KGE), such as DistMult and ComplEx, balance expressiveness and computational efficiency by representing relations as diagonal matrices. Although they perform well in predicting atomic relations, composite relations (relation paths) cannot be modeled naturally by the product of relation matrices, as the product of diagonal matrices is commutative and hence invariant with the order of relations. In this paper, we propose a new bilinear KGE model, called BlockHolE, based on block circulant matrices. In BlockHolE, relation matrices can be non commutative, allowing composite relations to be modeled by matrix product. The model is parameterized in a way that covers a spectrum ranging from diagonal to full relation matrices. A fast computation technique is developed on the basis of the duality of the Fourier transform of circulant matrices.


The Consequences of Switching Strategies in a Two Player Iterated Survival Game. We consider two player iterated survival games in which players may switch from a more cooperative behavior to a less cooperative one at some step of the game. Payoffs are survival probabilities and lone individuals have to finish the game on their own. We explore the potential of these games to support cooperation, focusing on the case in which each single step is a Prisoner's Dilemma. We find that incentives for or against cooperation depend on the number of defections at the end of the game, as opposed to the number of steps in the game. Broadly, cooperation is supported when the survival prospects of lone individuals are relatively bleak. Specifically, we find three critical values or cutoffs for the loner survival probability which, in concert with other survival parameters, determine the incentives for or against cooperation. One cutoff determines the existence of an optimal number of defections against a fully cooperative partner, one determines whether additional defections eventually become disfavored as the number of defections by the partner increases, and one determines whether additional cooperations eventually become favored as the number of defections by the partner increases. We obtain expressions for these switch points and for optimal numbers of defections against partners with various strategies. These typically involve small numbers of defections even in very long games. We show that potentially long stretches of equilibria may exist, in which there is no incentive to defect more or cooperate more. We describe how individuals find equilibria in best response walks among strategies, and establish that evolutionary stability requires there be just one such equilibrium. Otherwise, equilibria are not protected against invasion by strategies with fewer defections.


Approaching Neural Chinese Word Segmentation as a Low Resource Machine Translation Task. Supervised Chinese word segmentation has entered the deep learning era which reduces the hassle of feature engineering. Recently, some researchers attempted to treat it as character level translation which further simplified model designing and building, but there is still a performance gap between the translation based approach and other methods. In this work, we apply the best practices from low resource neural machine translation to Chinese word segmentation. We build encoder decoder models with attention, and examine a series of techniques including regularization, data augmentation, objective weighting, transfer learning and ensembling. Our method is generic for word segmentation, without the need for feature engineering or model implementation. In the closed test with constrained data, our method ties with the state of the art on the MSR dataset and is comparable to other methods on the PKU dataset.


Effective Medium Transformation the Case of Stratified Magnetic Structures. Effective medium theory replaces a given fine scale heterostructure with a homogeneous one in such a way that the physically measurable quantities, e.g. reaction fields and losses, remain approximately the same. This Letter shows that the very nature of the physical problem may change upon homogenization. A specific example is a stratified nonlinear magnetic and conducting medium, where a low frequency excitation induces eddy currents. It is shown that the appropriate coarse scale (homogeneous) model is, counter intuitively, magnetostatic, with an effective complex valued BH curve whose real and imaginary parts represent active and reactive losses in the sample. Similar situations may arise in other physical and engineering applications, notably, in diffusion problems with boundary layers.


Manipulation Trajectory Optimization with Online Grasp Synthesis and Selection. In robot manipulation, planning the motion of a robot manipulator to grasp an object is a fundamental problem. A manipulation planner needs to generate a trajectory of the manipulator arm to avoid obstacles in the environment and plan an end effector pose for grasping. While trajectory planning and grasp planning are often tackled separately, how to efficiently integrate the two planning problems remains a challenge. In this work, we present a novel method for joint motion and grasp planning. Our method integrates manipulation trajectory optimization with online grasp synthesis and selection, where we apply online learning techniques to select goal configurations for grasping, and introduce a new grasp synthesis algorithm to generate grasps online. We evaluate our planning approach and demonstrate that our method generates robust and efficient motion plans for grasping in cluttered scenes. Our video can be found at https www.youtube.com watch?v LIcACf8YkGU .


Conversational Transfer Learning for Emotion Recognition. Recognizing emotions in conversations is a challenging task due to the presence of contextual dependencies governed by self and inter personal influences. Recent approaches have focused on modeling these dependencies primarily via supervised learning. However, purely supervised strategies demand large amounts of annotated data, which is lacking in most of the available corpora in this task. To tackle this challenge, we look at transfer learning approaches as a viable alternative. Given the large amount of available conversational data, we investigate whether generative conversational models can be leveraged to transfer affective knowledge for detecting emotions in context. We propose an approach, TL ERC, where we pre train a hierarchical dialogue model on multi turn conversations (source) and then transfer its parameters to a conversational emotion classifier (target). In addition to the popular practice of using pre trained sentence encoders, our approach also incorporates recurrent parameters that model inter sentential context across the whole conversation. Based on this idea, we perform several experiments across multiple datasets and find improvement in performance and robustness against limited training data. TL ERC also achieves better validation performances in significantly fewer epochs. Overall, we infer that knowledge acquired from dialogue generators can indeed help recognize emotions in conversations.


An Eulerian model for sea spray transport and evaporation. Reliable estimates of the fluxes of momentum, heat, and moisture at the air sea interface are essential for accurate long term climate projections, as well as the prediction of short term weather events such as tropical cyclones. In recent years, it has been suggested that these estimates need to incorporate an accurate description of the transport of sea spray within the atmospheric boundary layer and the drop induced fluxes of moisture, momentum, and heat, so that the resulting effects on the atmospheric flow can be evaluated. In this paper we propose a model based on a theoretical and mathematical framework inspired from kinetic gas theory. This approach reconciles the Lagrangian nature of the spray transport with the Eulerian description of the atmosphere. In turn, this enables a relatively straightforward inclusion of the spray fluxes and the resulting spray effects on the atmospheric flow. A comprehensive dimensional analysis has led us to identify the spray effects that are most likely to influence the temperature, moisture, and speed of the airflow. We also provide an example application to illustrate the capabilities of the model in specific environmental conditions. Finally, suggestions for future work are offered.


Hybrid Score and Rank level Fusion for Person Identification using Face and ECG Data. Uni modal identification systems are vulnerable to errors in sensor data collection and are therefore more likely to misidentify subjects. For instance, relying on data solely from an RGB face camera can cause problems in poorly lit environments or if subjects do not face the camera. Other identification methods such as electrocardiograms (ECG) have issues with improper lead connections to the skin. Errors in identification are minimized through the fusion of information gathered from both of these models. This paper proposes a methodology for combining the identification results of face and ECG data using Part A of the BioVid Heat Pain Database containing synchronized RGB video and ECG data on 87 subjects. Using 10 fold cross validation, face identification was 98.8 accurate, while the ECG identification was 96.1 accurate. By using a fusion approach the identification accuracy improved to 99.8 . Our proposed methodology allows for identification accuracies to be significantly improved by using disparate face and ECG models that have non overlapping modalities.


Group Re Identification with Multi grained Matching and Integration. The task of re identifying groups of people underdifferent camera views is an important yet less studied problem.Group re identification (Re ID) is a very challenging task sinceit is not only adversely affected by common issues in traditionalsingle object Re ID problems such as viewpoint and human posevariations, but it also suffers from changes in group layout andgroup membership. In this paper, we propose a novel conceptof group granularity by characterizing a group image by multi grained objects individual persons and sub groups of two andthree people within a group. To achieve robust group Re ID,we first introduce multi grained representations which can beextracted via the development of two separate schemes, i.e. onewith hand crafted descriptors and another with deep neuralnetworks. The proposed representation seeks to characterize bothappearance and spatial relations of multi grained objects, and isfurther equipped with importance weights which capture varia tions in intra group dynamics. Optimal group wise matching isfacilitated by a multi order matching process which in turn,dynamically updates the importance weights in iterative fashion.We evaluated on three multi camera group datasets containingcomplex scenarios and large dynamics, with experimental resultsdemonstrating the effectiveness of our approach. The published dataset can be found in url http min.sjtu.edu.cn lwydemo GroupReID.html


On the geometry of the Pappas Rapoport models for PEL Shimura varieties. We study integral models, so called Pappas Rapoport or splitting models, of some PEL Shimura Varieties whose data are ramified at a prime p. We show that except in a specific case, these models are smooth when there is no level at p, and we study their special fiber modulo p. We show that interesting stratifications appear, and we show that the mu ordinary locus is open and dense, following a strategy of Wedhorn in the unramified case.


Predicting Engagement with the Internet Research Agency's Facebook and Instagram Campaigns around the 2016 U.S. Presidential Election. The Russian Internet Research Agency's (IRA) online interference campaign in the 2016 U.S. presidential election represents a turning point in the trajectory of democratic elections in the digital age. What can we learn about how the IRA engages U.S. audiences, ahead of the 2020 U.S. presidential election? We provide the first in depth analysis of the relationships between IRA content characteristics and user engagement on Facebook and Instagram around the 2016 election. We find that content targeting right wing and non Black marginalised groups had the strongest positive association with engagement on both Facebook and Instagram, in contrast to findings from the IRA campaign on Twitter and to some previous commentary in the media. Higher engagement was associated with posting later in the 2015 2017 period and using less text on both platforms, using negative wording and not including links on Facebook, and using fewer hashtags on Instagram. The sub audiences and sub issues associated with most engagement differed across the platforms.


RM CVaR Regularized Multiple beta CVaR Portfolio. The problem of finding the optimal portfolio for investors is called the portfolio optimization problem. Such problem mainly concerns the expectation and variability of return (i.e., mean and variance). Although the variance would be the most fundamental risk measure to be minimized, it has several drawbacks. Conditional Value at Risk (CVaR) is a relatively new risk measure that addresses some of the shortcomings of well known variance related risk measures, and because of its computational efficiencies, it has gained popularity. CVaR is defined as the expected value of the loss that occurs beyond a certain probability level ( beta ). However, portfolio optimization problems that use CVaR as a risk measure are formulated with a single beta and may output significantly different portfolios depending on how the beta is selected. We confirm even small changes in beta can result in huge changes in the whole portfolio structure. In order to improve this problem, we propose RM CVaR Regularized Multiple beta CVaR Portfolio. We perform experiments on well known benchmarks to evaluate the proposed portfolio. Compared with various portfolios, RM CVaR demonstrates a superior performance of having both higher risk adjusted returns and lower maximum drawdown.


Distributed Multi agent Video Fast forwarding. In many intelligent systems, a network of agents collaboratively perceives the environment for better and more efficient situation awareness. As these agents often have limited resources, it could be greatly beneficial to identify the content overlapping among camera views from different agents and leverage it for reducing the processing, transmission and storage of redundant unimportant video frames. This paper presents a consensus based distributed multi agent video fast forwarding framework, named DMVF, that fast forwards multi view video streams collaboratively and adaptively. In our framework, each camera view is addressed by a reinforcement learning based fast forwarding agent, which periodically chooses from multiple strategies to selectively process video frames and transmits the selected frames at adjustable paces. During every adaptation period, each agent communicates with a number of neighboring agents, evaluates the importance of the selected frames from itself and those from its neighbors, refines such evaluation together with other agents via a system wide consensus algorithm, and uses such evaluation to decide their strategy for the next period. Compared with approaches in the literature on a real world surveillance video dataset VideoWeb, our method significantly improves the coverage of important frames and also reduces the number of frames processed in the system.


Finding Needles in a Moving Haystack Prioritizing Alerts with Adversarial Reinforcement Learning. Detection of malicious behavior is a fundamental problem in security. One of the major challenges in using detection systems in practice is in dealing with an overwhelming number of alerts that are triggered by normal behavior (the so called false positives), obscuring alerts resulting from actual malicious activity. While numerous methods for reducing the scope of this issue have been proposed, ultimately one must still decide how to prioritize which alerts to investigate, and most existing prioritization methods are heuristic, for example, based on suspiciousness or priority scores. We introduce a novel approach for computing a policy for prioritizing alerts using adversarial reinforcement learning. Our approach assumes that the attackers know the full state of the detection system and dynamically choose an optimal attack as a function of this state, as well as of the alert prioritization policy. The first step of our approach is to capture the interaction between the defender and attacker in a game theoretic model. To tackle the computational complexity of solving this game to obtain a dynamic stochastic alert prioritization policy, we propose an adversarial reinforcement learning framework. In this framework, we use neural reinforcement learning to compute best response policies for both the defender and the adversary to an arbitrary stochastic policy of the other. We then use these in a double oracle framework to obtain an approximate equilibrium of the game, which in turn yields a robust stochastic policy for the defender. Extensive experiments using case studies in fraud and intrusion detection demonstrate that our approach is effective in creating robust alert prioritization policies.


Qualitative properties of spatial epidemiological models. We study the qualitative properties of a spatial diffusive heterogeneous SIR model, that appears in mathematical epidemiology to describe the spread of an infectious disease in a population. The model we consider consists in a system of parabolic PDEs. In the first part of the paper, we give a criterion that ensures whether or not an epidemic propagates in a given population. We show how the features of the disease and of the population (rates of infection and of recovery, localisation and diffusivity of individuals) influence the propagation of the epidemic. In particular, we prove that there are situations where "slowing down" the individuals can trigger an epidemic that would not propagate otherwise. In the second part of the paper, we show how the spatial diffusive SIR model qualitatively differs from the usual, purely temporal, SIR model.


Toward a 6G AI Native Air Interface. Each generation of cellular communication systems is marked by a defining disruptive technology of its time, such as orthogonal frequency division multiplexing (OFDM) for 4G or Massive multiple input multiple output (MIMO) for 5G. Since artificial intelligence (AI) is the defining technology of our time, it is natural to ask what role it could play for 6G. While it is clear that 6G must cater to the needs of large distributed learning systems, it is less certain if AI will play a defining role in the design of 6G itself. The goal of this article is to paint a vision of a new air interface which is partially designed by AI to enable optimized communication schemes for any hardware, radio environment, and application.


Generalization to Novel Objects using Prior Relational Knowledge. To solve tasks in new environments involving objects unseen during training, agents must reason over prior information about those objects and their relations. We introduce the Prior Knowledge Graph network, an architecture for combining prior information, structured as a knowledge graph, with a symbolic parsing of the visual scene, and demonstrate that this approach is able to apply learned relations to novel objects whereas the baseline algorithms fail. Ablation experiments show that the agents ground the knowledge graph relations to semantically relevant behaviors. In both a Sokoban game and the more complex Pacman environment, our network is also more sample efficient than the baselines, reaching the same performance in 5 10x fewer episodes. Once the agents are trained with our approach, we can manipulate agent behavior by modifying the knowledge graph in semantically meaningful ways. These results suggest that our network provides a framework for agents to reason over structured knowledge graphs while still leveraging gradient based learning approaches.


Hankel determinants of linear combinations of moments of orthogonal polynomials. We prove evaluations of Hankel determinants of linear combinations of moments of orthogonal polynomials (or, equivalently, of generating functions for Motzkin paths), thus generalising known results for Catalan numbers.


Indicators of Open Access for universities. This paper presents a first attempt to analyse Open Access integration at the institutional level. For this, we combine information from Unpaywall and the Leiden Ranking to offer basic OA indicators for universities. We calculate the overall number of Open Access publications for 930 universities worldwide. OA indicators are also disaggregated by green, gold and hybrid Open Access. We then explore differences between and within countries and offer a general ranking of universities based on the proportion of their output which is openly accessible.


Mobile Robot Path Planning in Dynamic Environments A Survey. There are many challenges for robot navigation in densely populated dynamic environments. This paper presents a survey of the path planning methods for robot navigation in dense environments. Particularly, the path planning in the navigation framework of mobile robots is composed of global path planning and local path planning, with regard to the planning scope and the executability. Within this framework, the recent progress of the path planning methods is presented in the paper, while examining their strengths and weaknesses. Notably, the recently developed Velocity Obstacle method and its variants that serve as the local planner are analyzed comprehensively. Moreover, as a model free method that is widely used in current robot applications, the reinforcement learning based path planning algorithms are detailed in this paper.


Powering Electronic Devices from Salt Gradients in AA Battery Sized Stacks of Hydrogel Infused Paper. Strongly electric fish use gradients of ions within their bodies to generate stunning external electrical discharges the most powerful of these organisms, the Atlantic torpedo ray, can produce pulses of over 1 kW from its electric organs. Despite extensive study of this phenomenon in nature, the development of artificial power generation schemes based on ion gradients for portable, wearable, or implantable human use has remained out of reach. Previously, inspired by the electric eel, we developed an artificial electric organ that generated electricity from ion gradients within stacked hydrogels and, like the eel, was optimized to deliver large voltages that exceeded 100 V. Due to its high internal resistance, the current of this power source was, however, too low to power standard electronics. Here we introduce an artificial electric organ that takes inspiration from the unique morphologies of torpedo rays for maximal current output. This power source uses a hybrid material of hydrogel infused paper to create, organize, and reconfigure stacks of thin, arbitrarily large gel films both in series and in parallel. The resulting increase in electrical power by almost two orders of magnitude compared to the original eel inspired design makes it possible to power electronic devices and establishes that biology's mechanism of generating significant electrical power can now be realized from benign and soft materials in a portable size.


Marked point processes and intensity ratios for limit order book modeling. This paper extends the analysis of Muni Toke and Yoshida (2020) to the case of marked point processes. We consider multiple marked point processes with intensities defined by three multiplicative components, namely a common baseline intensity, a state dependent component specific to each process, and a state dependent component specific to each mark within each process. We show that for specific mark distributions, this model is a combination of the ratio models defined in Muni Toke and Yoshida (2020). We prove convergence results for the quasi maximum and quasi Bayesian likelihood estimators of this model and provide numerical illustrations of the asymptotic variances. We use these ratio processes in order to model transactions occuring in a limit order book. Model flexibility allows us to investigate both state dependency (emphasizing the role of imbalance and spread as significant signals) and clustering. Calibration, model selection and prediction results are reported for high frequency trading data on multiple stocks traded on Euronext Paris. We show that the marked ratio model outperforms other intensity based methods (such as "pure" Hawkes based methods) in predicting the sign and aggressiveness of market orders on financial markets.


A General Stabilization Bound for Influence Propagation in Graphs. We study the stabilization time of a wide class of processes on graphs, in which each node can only switch its state if it is motivated to do so by at least a frac 1 lambda 2 fraction of its neighbors, for some 0 lambda 1 . Two examples of such processes are well studied dynamically changing colorings in graphs in majority processes, nodes switch to the most frequent color in their neighborhood, while in minority processes, nodes switch to the least frequent color in their neighborhood. We describe a non elementary function f( lambda) , and we show that in the sequential model, the worst case stabilization time of these processes can completely be characterized by f( lambda) . More precisely, we prove that for any epsilon 0 , O(n 1 f( lambda) epsilon ) is an upper bound on the stabilization time of any proportional majority minority process, and we also show that there are graph constructions where stabilization indeed takes Omega(n 1 f( lambda) epsilon ) steps.


New structurally unstable families of planar vector fields. We study global bifurcations in generic 3 parameter families of vector fields on S 2 . In the recent article arXiv 1506.06797 , Ilyashenko, Kudryashov, and Schurov show that 3 parameter unfoldings of vector fields with the polycycle "tears of the heart" are structurally unstable. We consider 3 parameter unfoldings of vector fields with separatrix graphs "ears" and "glasses", and prove that these families are structurally unstable as well. We also study in more details the classical bifurcation of a saddle loop, and use it as a building block in our main example.


Regulation of interferon production as a potential strategy for COVID 19 treatment. Regulating the upstream of the cytokines production could be a promising strategy to the treatment of COVID 19. We suggest to pay more attention to the dysregulated IFN I production in COVID 19 and to considerate cGAS, ALK and STING as potential therapeutic targets preventing cytokine storm. Approved drugs like suramin and ALK inhibitors are worthy of clinical trials.


On Khovanov Homology and Related Invariants. This paper begins with a survey of some applications of Khovanov homology to low dimensional topology, with an eye toward extending these results to mathfrak sl (n) homologies. We extend Levine and Zemke's ribbon concordance obstruction from Khovanov homology to mathfrak sl (n) homology for n geq 2 , including the universal mathfrak sl (2) and mathfrak sl (3) homology theories. Inspired by Alishahi and Dowlin's bounds for the unknotting number coming from Khovanov homology and relying on spectral sequence arguments, we produce bounds on the alternation number of a knot. Lee and Bar Natan spectral sequences also provide lower bounds on Turaev genus.


Multi scale Dynamic Feature Encoding Network for Image Demoireing. The prevalence of digital sensors, such as digital cameras and mobile phones, simplifies the acquisition of photos. Digital sensors, however, suffer from producing Moire when photographing objects having complex textures, which deteriorates the quality of photos. Moire spreads across various frequency bands of images and is a dynamic texture with varying colors and shapes, which pose two main challenges in demoireing an important task in image restoration. In this paper, towards addressing the first challenge, we design a multi scale network to process images at different spatial resolutions, obtaining features in different frequency bands, and thus our method can jointly remove moire in different frequency bands. Towards solving the second challenge, we propose a dynamic feature encoding module (DFE), embedded in each scale, for dynamic texture. Moire pattern can be eliminated more effectively via DFE.Our proposed method, termed Multi scale convolutional network with Dynamic feature encoding for image DeMoireing (MDDM), can outperform the state of the arts in fidelity as well as perceptual on benchmarks.


Unsupervised Learning from Video with Deep Neural Embeddings. Because of the rich dynamical structure of videos and their ubiquity in everyday life, it is a natural idea that video data could serve as a powerful unsupervised learning signal for training visual representations in deep neural networks. However, instantiating this idea, especially at large scale, has remained a significant artificial intelligence challenge. Here we present the Video Instance Embedding (VIE) framework, which extends powerful recent unsupervised loss functions for learning deep nonlinear embeddings to multi stream temporal processing architectures on large scale video datasets. We show that VIE trained networks substantially advance the state of the art in unsupervised learning from video datastreams, both for action recognition in the Kinetics dataset, and object recognition in the ImageNet dataset. We show that a hybrid model with both static and dynamic processing pathways is optimal for both transfer tasks, and provide analyses indicating how the pathways differ. Taken in context, our results suggest that deep neural embeddings are a promising approach to unsupervised visual learning across a wide variety of domains.


Least Squares Optimization from Theory to Practice. Nowadays, Non Linear Least Squares embodies the foundation of many Robotics and Computer Vision systems. The research community deeply investigated this topic in the last years, and this resulted in the development of several open source solvers to approach constantly increasing classes of problems. In this work, we propose a unified methodology to design and develop efficient Least Squares Optimization algorithms, focusing on the structures and patterns of each specific domain. Furthermore, we present a novel open source optimization system, that addresses transparently problems with a different structure and designed to be easy to extend. The system is written in modern C and can run efficiently on embedded systems. Source code https srrg.gitlab.io srrg2 solver.html. We validated our approach by conducting comparative experiments on several problems using standard datasets. The results show that our system achieves state of the art performances in all tested scenarios.


Robust Zero Shot Cross Domain Slot Filling with Example Values. Task oriented dialog systems increasingly rely on deep learning based slot filling models, usually needing extensive labeled training data for target domains. Often, however, little to no target domain training data may be available, or the training and target domain schemas may be misaligned, as is common for web forms on similar websites. Prior zero shot slot filling models use slot descriptions to learn concepts, but are not robust to misaligned schemas. We propose utilizing both the slot description and a small number of examples of slot values, which may be easily available, to learn semantic representations of slots which are transferable across domains and robust to misaligned schemas. Our approach outperforms state of the art models on two multi domain datasets, especially in the low data setting.


Error analysis for a fractional derivative parabolic problem on quasi graded meshes using barrier functions. An initial boundary value problem with a Caputo time derivative of fractional order alpha in(0,1) is considered, solutions of which typically exhibit a singular behaviour at an initial time. For this problem, we give a simple and general numerical stability analysis using barrier functions, which yields sharp pointwise in time error bounds on quasi graded temporal meshes with arbitrary degree of grading. L1 type and Alikhanov type discretization in time are considered. In particular, those results imply that milder (compared to the optimal) grading yields optimal convergence rates in positive time. Semi discretizations in time and full discretizations are addressed. The theoretical findings are illustrated by numerical experiments.


Mixed multiplicities and projective degrees of rational maps. We consider the notion of mixed multiplicities for multigraded modules by using Hilbert series, and this is later applied to study the projective degrees of rational maps. We use a general framework to determine the projective degrees of a rational map via a computation of the multiplicity of the saturated special fiber ring. As specific applications, we provide explicit formulas for all the projective degrees of rational maps determined by perfect ideals of height two or by Gorenstein ideals of height three.


Extensional proofs in a propositional logic modulo isomorphisms. System I is a proof language for a fragment of propositional logic where isomorphic propositions, such as A wedge B and B wedge A , or A Rightarrow(B wedge C) and (A Rightarrow B) wedge(A Rightarrow C) are made equal. System I enjoys the strong normalisation property. This is sufficient to prove the existence of empty types, but not to prove the introduction property (every closed term in normal form is an introduction). Moreover, a severe restriction had to be made on the types of the variables in order to obtain the existence of empty types. We show here that adding eta expansion rules to System I permits to drop this restriction, and yields a strongly normalizing calculus with enjoying the introduction property.


Space Time Medium Functions as a Perfect Antenna Mixer Amplifier Transceiver. We show that a space time varying medium can function as a front end transceiver, i.e., an antenna mixer amplifier. Such a unique functionality is endowed by space time surface waves associated with complex space time wave vectors in a subluminal space time medium. The proposed structure introduces pure frequency up and down conversions and with very weak undesired time harmonics. In contrast to other recently proposed space time mixers, a large frequency up down conversion ratio, associated with gain is achievable. Furthermore, as the structure does not operate based on progressive energy transition between the space time modulation and the incident wave, it possesses a subwavelength thickness (metasurface). Such a multi functional, highly efficient and compact medium is expected to find various applications in modern wireless telecommunication systems.


Homogeneous Bases for Demazure Modules. We study the PBW filtration on the Demazure modules V r gamma ( lambda) associated to reflections r gamma at positive roots in type A n , long roots in type C n , short roots in type B n and positive roots not involving the simple root alpha n 1 in type D n . In each type, we construct a normal polytope labelling a basis for the associated graded space operatorname gr V r gamma ( lambda) and prove that the annihilating ideal of operatorname gr V r gamma ( lambda) is a monomial ideal.


Triple Junction and Grain Boundary Influences on Climate Signals in Polar Ice. The Climate Change Institute W. M. Keck Laser Ice Facility laser ablation inductively coupled plasma spectrometer (LA ICP MS) yields a sample every 121 micrometers, a resolution on the scale of ice crystal triple junctions and grain boundaries in ice cores. Recent publications suggest that these features can allow amplification of impurity concentrations, and allow migration through veins potentially obscuring climate signals preserved in polar ice. LA ICP MS data reveal that such features modify these signals by less than 6 in the case of Na, Ca, and Fe based on the examination of GISP2 ice deposited at the beginning of the Holocene.


The conjugate gradient method with various viewpoints. Connections of the conjugate gradient (CG) method with other methods in computational mathematics are surveyed, including the connections with the conjugate direction method, the subspace optimization method and the quasi Newton method BFGS in numrical optimization, and the Lanczos method in numerical linear algebra. Two sequences of polynomials related to residual vectors and conjugate vectors are reviewed, where the residual polynomials are similar to orthogonal polynomials in the approximation theory and the roots of the polynomials reveal certain information of the coefficient matrix. The convergence rates of the steepest descent and CG are reconsidered in a viewpoint different from textbooks. The connection of infinite dimensional CG with finite dimensional preconditioned CG is also reviewed via numerical solution of an elliptic equation.


Modal model theory. We introduce the subject of modal model theory, where one studies a mathematical structure within a class of similar structures under an extension concept that gives rise to mathematically natural notions of possibility and necessity. A statement varphi is possible in a structure (written Diamond varphi ) if varphi is true in some extension of that structure, and varphi is necessary (written Box varphi ) if it is true in all extensions of the structure. A principal case for us will be the class Mod(T) of all models of a given theory T all graphs, all groups, all fields, or what have you considered under the substructure relation. In this article, we aim to develop the resulting modal model theory. The class of all graphs is a particularly insightful case illustrating the remarkable power of the modal vocabulary, for the modal language of graph theory can express connectedness, k colorability, finiteness, countability, size continuum, size aleph 1 , aleph 2 , aleph omega , beth omega , first beth fixed point, first beth hyper fixed point and much more. A graph obeys the maximality principle Diamond Box varphi(a) to varphi(a) with parameters if and only if it satisfies the theory of the countable random graph, and it satisfies the maximality principle for sentences if and only if it is universal for finite graphs.


Mesoscopic architecture enhances communication across the Macaque connectome revealing structure function correspondence in the brain. Analyzing the brain in terms of organizational structures at intermediate scales provides an approach to negotiate the complexity arising from interactions between its large number of components. Focusing on a wiring diagram that spans the cortex, basal ganglia and thalamus of the Macaque brain, we provide a mesoscopic level description of the topological architecture of one of the most well studied mammalian connectomes. The robust modules we identify each comprise densely inter connected cortical and sub cortical areas that play complementary roles in executing specific cognitive functions. We find that physical proximity between areas is insufficient to explain the modular organization, as similar mesoscopic structures can be obtained even after factoring out the effect of distance constraints on the connectivity. We observe that the distribution profile of brain areas, classified in terms of their intra and inter modular connectivity, is conserved across the principal cortical subdivisions, as well as, sub cortical structures. In particular provincial hubs, which have significantly higher number of connections with members of their module, but relatively less well connected to other modules, are the only class that exhibits homophily, i.e., a discernible preference to connect to each other. By considering a process of diffusive propagation we demonstrate that this architecture, instead of localizing the activity, facilitates rapid communication across the connectome. By supplementing the topological information about the Macaque connectome with physical locations, volumes and functions of the constituent areas and analyzing this augmented dataset, we reveal a counter intuitive role played by the modular architecture of the brain in promoting global interaction.


On the Length of Monotone Paths in Polyhedra. Motivated by the problem of bounding the number of iterations of the Simplex algorithm we investigate the possible lengths of monotone paths followed by the Simplex method inside the oriented graphs of polyhedra (oriented by the objective function). We consider both the shortest and the longest monotone paths and estimate the monotone diameter and height of polyhedra. Our analysis applies to transportation polytopes, matroid polytopes, matching polytopes, shortest path polytopes, and the TSP, among others. We begin by showing that combinatorial cubes have monotone and Bland pivot height bounded by their dimension and that in fact all monotone paths of zonotopes are no larger than the number of edge directions of the zonotope. We later use this to show that several polytopes have polynomial size pivot height, for all pivot rules. In contrast, we show that many well known combinatorial polytopes have exponentially long monotone paths. Surprisingly, for some famous pivot rules, e.g., greatest improvement and steepest edge, these same polytopes have polynomial size simplex paths.


Interactions between Large Molecules Puzzle for Reference Quantum Mechanical Methods. Quantum mechanical methods are widely used for understanding molecular interactions throughout biology, chemistry, and materials science. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations CCSD(T) are two state of the art and trusted wavefunction methods that have been categorically shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely used semi empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. Approaching such well converged predictive power in larger molecules has motivated major developments in CCSD(T) as well as DMC algorithms in the past years, resulting in orders of magnitude time to solution reductions. Here, we show that CCSD(T) and DMC interaction energies are not in consistent agreement for a set of polarizable supramolecules. Whilst agreement is found for some of the complexes, in a few key systems disagreements of up to 8 kcal mol remained. This leads to differences of up to 6 orders of magnitude in the corresponding binding association constant at room temperature for systems which are well within the accustomed domain of applicability for both methods. These findings thus indicate that more caution is required when aiming at reproducible non covalent interactions between extended molecules. Our data contradicts the expectation that the most comprehensive and robust wavefunction methods predict identical non covalent interactions and indicate an unsolved challenge for benchmark approaches.


Dynamic model of HIV infection with immune system response of T lymphocytes, B cells and dendritic cells a review. A dynamic model of non lineal time dependent ordinary differential equations (ODE) has been applied to the interactions of a HIV infection with the immune system cells. This model has been simplified into two compartments lymph node and peripheral blood. The model includes CD4 T lymphocytes in several states (quiescent Q, naive N and activated T), cytotoxic CD8 T cells, B cells and dendritic cells. Cytokines and immunoglobulins specific for each antigen (i.e. gp41 or p24) have been also included in the model, modelling the atraction effect of CD4 T cells to the infected area and the reduction of virus concentration by immunoglobulins. HIV virus infection of CD4 T lymphocytes is modelled in several stages before fusion as HIV attached (H) and after fusion as non permissive abortively infected (M), and permissive latently infected (L) and permissive actively infected (I). These equations have been implemented in a C Python interface application, called Immune System app, which runs Open Modelica software to solve the ODE system through a 4th order Runge Kutta numerical approximation. Results of the simulation show that although HIV virus concentration in both compartments is lower than 10 10 virus mu L after t 2 years, quiescent lymphocytes reach an equilibrium with a concentration lower than the initial conditions, due to the latency state, which serves as a reservoir in time of virus production. As a conclusion, this model can provide reliable results in other conditions, such as antiviral therapies.


Autonomics In Search of a Foundation for Next Generation Autonomous Systems. The potential benefits of autonomous systems have been driving intensive development of such systems, and of supporting tools and methodologies. However, there are still major issues to be dealt with before such development becomes commonplace engineering practice, with accepted and trustworthy deliverables. We argue that a solid, evolving, publicly available, community controlled foundation for developing next generation autonomous systems is a must. We discuss what is needed for such a foundation, identify a central aspect thereof, namely, decision making, and focus on three main challenges (i) how to specify autonomous system behavior and the associated decisions in the face of unpredictability of future events and conditions and the inadequacy of current languages for describing these (ii) how to carry out faithful simulation and analysis of system behavior with respect to rich environments that include humans, physical artifacts, and other systems, and (iii) how to engineer systems that combine executable model driven techniques and data driven machine learning techniques. We argue that autonomics, i.e., the study of unique challenges presented by next generation autonomous systems, and research towards resolving them, can introduce substantial contributions and innovations in system engineering and computer science.


Efficient implementation of immersed boundary lattice Boltzmann method for massive particle laden flows Part I Serial computing. Immersed boundary lattice Boltzmann method (IB LBM) has been widely used for simulation of particle laden flows recently. However, it was limited to small scale simulations with no more than O(103) particles. Here, we expand IB LBM for massive particle laden flows with more than O(104) particles by two sequential works. First is the Part I serial computing on a single CPU core and following the Part II parallel computing on many CPU cores. In this Part I paper, a highly efficient and localized implementation of IB LBM is proposed for serial computing. We optimize in three main aspects swap algorithm for incompressible LBM, local grid to point algorithm for IBM and improved grid search algorithm for particle pair short range interaction. In addition, symmetry algorithm is proposed for the half calculation of LB collision and external force term. The computational performance on a single CPU core is analyzed. Different scales of two dimensional (2D) and three dimensional (3D) particles settling in closed cavities are used for testing. The solid volume fraction is varied from 0 to 0.40. Simulation results demonstrate that all calculation parts are dramatically decreased by the improved algorithm. For the particle free flows, the Mega Lattice Site Update per Second (MLUPS) can be achieved up to 36 (2D) and 12 (3D) using the improved algorithm. For the particle laden flows, MLUPS can be achieved no lower than 15 (2D) and 7 (3D) in the simulations of dense flows. At last, we discuss the potential of the new algorithms for the high performance computation of the large scale systems of particle laden flows with MPI parallel technique.


An Introduction to Mordell Weil Theorem. This article is an introduction to Mordell Weil theorem. In this article, I introduced some basic properties about ellptic curves and proved the theorem in two different ways.


A note on rank two stable bundles over surfaces. Let pi X longrightarrow C be a fibration with reduced fibers over a curve C and consider a polarization H on the surface X . Let E be a stable vector bundle of rank 2 on C . We prove that the pullback pi E is a H stable bundle over X . This result allows us to relate the corresponding moduli spaces of stable bundles mathcal M C(2,d) and mathcal M X,H (2,df,0) through an injective morphism. We study the induced morphism at the level of Brill Noether loci to construct examples of Brill Noether loci on fibered surfaces. Results concerning the emptiness of Brill Noether loci follow as a consequence of a generalization of Clifford's Theorem for rank two bundles on surfaces.


On Lp viscosity solutions of bilateral obstacle problems with unbounded ingredients. The global equi continuity estimate on Lp viscosity solutions of bilateral obstacle problems with unbounded ingredients is established when obstacles are merely continuous. The existence of Lp viscosity solutions is established via an approximation of given data. The local H "older continuity estimate on the first derivative of Lp viscosity solutions is shown when the obstacles belong to C 1, beta , and p n.


Privacy Preserving Distributed Machine Learning with Federated Learning. Edge computing and distributed machine learning have advanced to a level that can revolutionize a particular organization. Distributed devices such as the Internet of Things (IoT) often produce a large amount of data, eventually resulting in big data that can be vital in uncovering hidden patterns, and other insights in numerous fields such as healthcare, banking, and policing. Data related to areas such as healthcare and banking can contain potentially sensitive data that can become public if they are not appropriately sanitized. Federated learning (FedML) is a recently developed distributed machine learning (DML) approach that tries to preserve privacy by bringing the learning of an ML model to data owners'. However, literature shows different attack methods such as membership inference that exploit the vulnerabilities of ML models as well as the coordinating servers to retrieve private data. Hence, FedML needs additional measures to guarantee data privacy. Furthermore, big data often requires more resources than available in a standard computer. This paper addresses these issues by proposing a distributed perturbation algorithm named as DISTPAB, for privacy preservation of horizontally partitioned data. DISTPAB alleviates computational bottlenecks by distributing the task of privacy preservation utilizing the asymmetry of resources of a distributed environment, which can have resource constrained devices as well as high performance computers. Experiments show that DISTPAB provides high accuracy, high efficiency, high scalability, and high attack resistance. Further experiments on privacy preserving FedML show that DISTPAB is an excellent solution to stop privacy leaks in DML while preserving high data utility.


Seeing Things from a Different Angle Discovering Diverse Perspectives about Claims. One key consequence of the information revolution is a significant increase and a contamination of our information supply. The practice of fact checking won't suffice to eliminate the biases in text data we observe, as the degree of factuality alone does not determine whether biases exist in the spectrum of opinions visible to us. To better understand controversial issues, one needs to view them from a diverse yet comprehensive set of perspectives. For example, there are many ways to respond to a claim such as "animals should have lawful rights", and these responses form a spectrum of perspectives, each with a stance relative to this claim and, ideally, with evidence supporting it. Inherently, this is a natural language understanding task, and we propose to address it as such. Specifically, we propose the task of substantiated perspective discovery where, given a claim, a system is expected to discover a diverse set of well corroborated perspectives that take a stance with respect to the claim. Each perspective should be substantiated by evidence paragraphs which summarize pertinent results and facts. We construct PERSPECTRUM, a dataset of claims, perspectives and evidence, making use of online debate websites to create the initial data collection, and augmenting it using search engines in order to expand and diversify our dataset. We use crowd sourcing to filter out noise and ensure high quality data. Our dataset contains 1k claims, accompanied with pools of 10k and 8k perspective sentences and evidence paragraphs, respectively. We provide a thorough analysis of the dataset to highlight key underlying language understanding challenges, and show that human baselines across multiple subtasks far outperform ma chine baselines built upon state of the art NLP techniques. This poses a challenge and opportunity for the NLP community to address.


Generic symmetry forced infinitesimal rigidity translations and rotations. We characterize the combinatorial types of symmetric frameworks in the plane that are minimally generically symmetry forced infinitesimally rigid when the symmetry group consists of rotations and translations. Along the way, we use tropical geometry to show how a construction of Edmonds that associates a matroid to a submodular function can be used to give a description of the algebraic matroid of a Hadamard product of two linear spaces in terms of the matroids of each linear space. This leads to new, short, proofs of Laman's theorem, and a theorem of Jord 'a n, Kaszanitzky, and Tanigawa, and Malestein and Theran characterizing the minimally generically symmetry forced rigid graphs in the plane when the symmetry group contains only rotations.


Contraction Clustering (RASTER) A Very Fast Big Data Algorithm for Sequential and Parallel Density Based Clustering in Linear Time, Constant Memory, and a Single Pass. Clustering is an essential data mining tool for analyzing and grouping similar objects. In big data applications, however, many clustering algorithms are infeasible due to their high memory requirements and or unfavorable runtime complexity. In contrast, Contraction Clustering (RASTER) is a single pass algorithm for identifying density based clusters with linear time complexity. Due to its favorable runtime and the fact that its memory requirements are constant, this algorithm is highly suitable for big data applications where the amount of data to be processed is huge. It consists of two steps (1) a contraction step which projects objects onto tiles and (2) an agglomeration step which groups tiles into clusters. This algorithm is extremely fast in both sequential and parallel execution. Our quantitative evaluation shows that a sequential implementation of RASTER performs significantly better than various standard clustering algorithms. Furthermore, the parallel speedup is significant on a contemporary workstation, an implementation in Rust processes a batch of 500 million points with 1 million clusters in less than 50 seconds on one core. With 8 cores, the algorithm is about four times faster.


Analysis of simultaneous inpainting and geometric separation based on sparse decomposition. Natural images are often the superposition of various parts of different geometric characteristics. For instance, an image might be a mixture of cartoon and texture structures. In addition, images are often given with missing data. In this paper, we develop a method for simultaneously decomposing an image to its two underlying parts and inpainting the missing data. Our separation inpainting method is based on and l 1 minimization approach, using two dictionaries, each sparsifying one of the image parts but not the other. We introduce a comprehensive convergence analysis of our method, in a general setting, utilizing the concepts of joint concentration, clustered sparsity, and cluster coherence. As the main application of our theory, we consider the problem of separating and inpainting an image to a cartoon and texture parts.


Optimal Questionnaires for Screening of Strategic Agents. During the COVID 19 pandemic the health authorities at airports and train stations try to screen and identify the travellers possibly exposed to the virus. However, many individuals avoid getting tested and hence may misreport their travel history. This is a challenge for the health authorities who wish to ascertain the truly susceptible cases in spite of this strategic misreporting. We investigate the problem of questioning travellers to classify them for further testing when the travellers are strategic or are unwilling to reveal their travel histories. We show there are fundamental limits to how many travel histories the health authorities can recover. can be correctly classified by any probing mechanism.


Perfect graphs for domination games. Let gamma g(G) and gamma tg (G) be the game domination number and the total game domination number of a graph G , respectively. Then G is gamma g perfect (resp. gamma tg perfect), if every induced subgraph F of G satisfies gamma g(F) gamma(F) (resp. gamma tg (F) gamma t(F) ). A recursive characterization of gamma g perfect graphs is derived. The characterization yields a polynomial recognition algorithm for gamma g perfect graphs. It is proved that every minimally gamma g imperfect graph has domination number 2 . All minimally gamma g imperfect triangle free graphs are determined. It is also proved that gamma tg perfect graphs are precisely overline 2P 3 free cographs.


Monte Carlo Tree Search Based Tactical Maneuvering. In this paper we explore the application of simultaneous move Monte Carlo Tree Search (MCTS) based online framework for tactical maneuvering between two unmanned aircrafts. Compared to other techniques, MCTS enables efficient search over long horizons and uses self play to select best maneuver in the current state while accounting for the opponent aircraft tactics. We explore different algorithmic choices in MCTS and demonstrate the framework numerically in a simulated 2D tactical maneuvering application.


Image Inpainting with Learnable Bidirectional Attention Maps. Most convolutional network (CNN) based inpainting methods adopt standard convolution to indistinguishably treat valid pixels and holes, making them limited in handling irregular holes and more likely to generate inpainting results with color discrepancy and blurriness. Partial convolution has been suggested to address this issue, but it adopts handcrafted feature re normalization, and only considers forward mask updating. In this paper, we present a learnable attention map module for learning feature renormalization and mask updating in an end to end manner, which is effective in adapting to irregular holes and propagation of convolution layers. Furthermore, learnable reverse attention maps are introduced to allow the decoder of U Net to concentrate on filling in irregular holes instead of reconstructing both holes and known regions, resulting in our learnable bidirectional attention maps. Qualitative and quantitative experiments show that our method performs favorably against state of the arts in generating sharper, more coherent and visually plausible inpainting results. The source code and pre trained models will be available.


Product Formula of Multiple Integrals of Levy Process. We derive a product formula for the multiple stochastic integrals with respect to Levy process. The idea is to use exponential vectors and the polarization technique which greatly simplify the argument.


Targeted Physical World Attention Attack on Deep Learning Models in Road Sign Recognition. Real world traffic sign recognition is an important step towards building autonomous vehicles, most of which highly dependent on Deep Neural Networks (DNNs). Recent studies demonstrated that DNNs are surprisingly susceptible to adversarial examples. Many attack methods have been proposed to understand and generate adversarial examples, such as gradient based attack, score based attack, decision based attack, and transfer based attacks. However, most of these algorithms are ineffective in real world road sign attack, because (1) iteratively learning perturbations for each frame is not realistic for a fast moving car and (2) most optimization algorithms traverse all pixels equally without considering their diverse contribution. To alleviate these problems, this paper proposes the targeted attention attack (TAA) method for real world road sign attack. Specifically, we have made the following contributions (1) we leverage the soft attention map to highlight those important pixels and skip those zero contributed areas this also helps to generate natural perturbations, (2) we design an efficient universal attack that optimizes a single perturbation noise based on a set of training images under the guidance of the pre trained attention map, (3) we design a simple objective function that can be easily optimized, (4) we evaluate the effectiveness of TAA on real world data sets. Experimental results validate that the TAA method improves the attack successful rate (nearly 10 ) and reduces the perturbation loss (about a quarter) compared with the popular RP2 method. Additionally, our TAA also provides good properties, e.g., transferability and generalization capability. We provide code and data to ensure the reproducibility https github.com AdvAttack RoadSignAttack.


Regularity results for a nonlinear elliptic parabolic system with oscillating coefficients. In this paper we study the initial boundary value problem for the system mbox div ( sigma(u) nabla varphi) 0, u t Delta u sigma(u) nabla varphi 2 . This problem is known as the thermistor problem which models the electrical heating of conductors. Our assumptions on sigma(u) leave open the possibility that liminf u rightarrow infty sigma(u) 0 , while limsup u rightarrow infty sigma(u) is large. This means that sigma(u) can oscillate wildly between 0 and a large positive number as u rightarrow infty . Thus our degeneracy is fundamentally different from the one that is present in porous medium type of equations. We obtain a weak solution (u, varphi) with nabla varphi , nabla u in L infty by first establishing a uniform upper bound for e varepsilon u for some small varepsilon . This leads to an inequality in nabla varphi , from whence follows the regularity result. This approach enables us to avoid first proving the H " o lder continuity of varphi in the space variables, which would have required that the elliptic coefficient sigma(u) be an A 2 weight. As it is known, the latter implies that ln sigma(u) is "nearly bounded".


Classical Theta Lifts for Higher Metaplectic Covering Group. The classical theta correspondence establishes a relationship between automorphic representations on special orthogonal groups and automorphic representations on symplectic groups or their double covers. This correspondence is achieved by using as integral kernel a theta series on the metaplectic double cover of a symplectic group that is constructed from the Weil representation. There is also an analogous local correspondence. In this work we present an extension of the classical theta correspondence to higher degree metaplectic covers of symplectic and special orthogonal groups. The key issue here is that for higher degree covers there is no analogue of the Weil representation, and additional ingredients are needed. Our work reflects a broader paradigm constructions in automorphic forms that work for algebraic groups or their double covers should often extend to higher degree metaplectic covers.


The Fourier Mukai transform made easy. We propose a slightly modified definition for the Fourier Mukai transform (on abelian varieties) that makes it much easier to remember various formulas. As an application, we give relatively short proofs for two important theorems the characterization of GV sheaves in terms of vanishing, due to Hacon and fact that M regularity implies (continuous) global generation, due to Pareschi and Popa.


Spray pyrolysis deposited La1 xSrxCoO3 thin films for potential non volatile memory applications. In this work thin films of the La1 xSrxCoO3 (0.05 x 0.26) compound were grown, employing the so called spray pyrolysis process. The as grown thin films exhibit polycrystalline microstructure, with uniform grain size distribution, and observable porosity. Regarding their electrical transport properties, the produced thin films show semiconducting like behavior, regardless the Sr doping level, which is most likely due to both the oxygen deficiencies and the grainy nature of the films. Furthermore, room temperature current voltage (I V) measurements reveal stable resistance switching behavior, which is well explained in terms of space charge limited conduction mechanism. The presented experimental results provide essential evidence regarding the engagement of low cost, industrial scale methods of growing perovskite transition metal oxide thin films, for potential applications in random access memory devices.


Lefschetz Bott fibrations on line bundles over symplectic manifolds. We describe Lefschetz Bott fibrations on complex line bundles over symplectic manifolds explicitly. As an application, we construct more than one strong symplectic filling of the link of the A k type singularity. In the appendix, we show that the total space of a Lefschetz Bott fibration over the unit disk serves as a strong symplectic filling of a contact manifold compatible with an open book induced by the fibration.


Mirostat A Neural Text Decoding Algorithm that Directly Controls Perplexity. Neural text decoding is important for generating high quality texts using language models. To generate high quality text, popular decoding algorithms like top k, top p (nucleus), and temperature based sampling truncate or distort the unreliable low probability tail of the language model. Though these methods generate high quality text after parameter tuning, they are ad hoc. Not much is known about the control they provide over the statistics of the output, which is important since recent reports show text quality is highest for a specific range of likelihoods. Here, first we provide a theoretical analysis of perplexity in top k, top p, and temperature sampling, finding that cross entropy behaves approximately linearly as a function of p in top p sampling whereas it is a nonlinear function of k in top k sampling, under Zipfian statistics. We use this analysis to design a feedback based adaptive top k text decoding algorithm called mirostat that generates text (of any length) with a predetermined value of perplexity, and thereby high quality text without any tuning. Experiments show that for low values of k and p in top k and top p sampling, perplexity drops significantly with generated text length, which is also correlated with excessive repetitions in the text (the boredom trap). On the other hand, for large values of k and p, we find that perplexity increases with generated text length, which is correlated with incoherence in the text (confusion trap). Mirostat avoids both traps experiments show that cross entropy has a near linear relation with repetition in generated text. This relation is almost independent of the sampling method but slightly dependent on the model used. Hence, for a given language model, control over perplexity also gives control over repetitions. Experiments with human raters for fluency, coherence, and quality further verify our findings.


Numerical study on the effect of geometric approximation error in the numerical solution of PDEs using a high order curvilinear mesh. When time dependent partial differential equations (PDEs) are solved numerically in a domain with curved boundary or on a curved surface, mesh error and geometric approximation error caused by the inaccurate location of vertices and other interior grid points, respectively, could be the main source of the inaccuracy and instability of the numerical solutions of PDEs. The role of these geometric errors in deteriorating the stability and particularly the conservation properties are largely unknown, which seems to necessitate very fine meshes especially to remove geometric approximation error. This paper aims to investigate the effect of geometric approximation error by using a high order mesh with negligible geometric approximation error, even for high order polynomial of order p. To achieve this goal, the high order mesh generator from CAD geometry called NekMesh is adapted for surface mesh generation in comparison to traditional meshes with non negligible geometric approximation error. Two types of numerical tests are considered. Firstly, the accuracy of differential operators is compared for various p on a curved element of the sphere. Secondly, by applying the method of moving frames, four different time dependent PDEs on the sphere are numerically solved to investigate the impact of geometric approximation error on the accuracy and conservation properties of high order numerical schemes for PDEs on the sphere.


q Supercongruences from the q Saalsch " u tz identity. In terms of the q Saalsch " u tz identity and the Chinese remainder theorem for coprime polynomials, we establish some q supercongruences modulo the third power of a cyclotomic polynomial. In particular, we give a q analogue of a formula due to Long and Ramakrishna Adv. Math. 290 (2016), 773 808 .


On Products of Closed Subsets in Free Groups. We present examples of closed subsets of a free group such that their product is not closed in the profinite topology. We discuss how to characterize a subset of a free group which is closed in the profinite topology and its product with any finitely generated subgroup of a free group is also closed in the profinite topology.


Time resolved XUV Opacity Measurements of Warm Dense Aluminium. The free free opacity in plasmas is fundamental to our understanding of energy transport in stellar interiors and for inertial confinement fusion research. However, theoretical predictions in the challenging dense plasma regime are conflicting and there is a dearth of accurate experimental data to allow for direct model validation. Here we present time resolved transmission measurements in solid density Al heated by an XUV free electron laser. We use a novel functional optimization approach to extract the temperature dependent absorption coefficient directly from an oversampled pool of single shot measurements, and find a pronounced enhancement of the opacity as the plasma is heated to temperatures of order the Fermi energy. Plasma heating and opacity enhancement is observed on ultrafast time scales, within the duration of the femtosecond XUV pulse. We attribute further rises in the opacity on ps timescales to melt and the formation of warm dense matter.


People Still Care About Facts Twitter Users Engage More with Factual Discourse than Misinformation A Comparison Between COVID and General Narratives on Twitter. Misinformation entails the dissemination of falsehoods that leads to the slow fracturing of society via decreased trust in democratic processes, institutions, and science. The public has grown aware of the role of social media as a superspreader of untrustworthy information, where even pandemics have not been immune. In this paper, we focus on COVID 19 misinformation and examine a subset of 2.1M tweets to understand misinformation as a function of engagement, tweet content (COVID 19 vs. non COVID 19 related), and veracity (misleading or factual). Using correlation analysis, we show the most relevant feature subsets among over 126 features that most heavily correlate with misinformation or facts. We found that (i) factual tweets, regardless of whether COVID related, were more engaging than misinformation tweets and (ii) features that most heavily correlated with engagement varied depending on the veracity and content of the tweet.


Generating Question Titles for Stack Overflow from Mined Code Snippets. Stack Overflow has been heavily used by software developers as a popular way to seek programming related information from peers via the internet. The Stack Overflow community recommends users to provide the related code snippet when they are creating a question to help others better understand it and offer their help. Previous studies have shown that a significant number of these questions are of low quality and not attractive to other potential experts in Stack Overflow. These poorly asked questions are less likely to receive useful answers and hinder the overall knowledge generation and sharing process. Considering one of the reasons for introducing low quality questions in SO is that many developers may not be able to clarify and summarize the key problems behind their presented code snippets due to their lack of knowledge and terminology related to the problem, and or their poor writing skills, in this study we propose an approach to assist developers in writing high quality questions by automatically generating question titles for a code snippet using a deep sequence to sequence learning approach. Our approach is fully data driven and uses an attention mechanism to perform better content selection, a copy mechanism to handle the rare words problem and a coverage mechanism to eliminate word repetition problem. We evaluate our approach on Stack Overflow datasets over a variety of programming languages (e.g., Python, Java, Javascript, C and SQL) and our experimental results show that our approach significantly outperforms several state of the art baselines in both automatic and human evaluation. We have released our code and datasets to facilitate other researchers to verify their ideas and inspire the follow up work.


Learning Non Monotonic Automatic Post Editing of Translations from Human Orderings. Recent research in neural machine translation has explored flexible generation orders, as an alternative to left to right generation. However, training non monotonic models brings a new complication how to search for a good ordering when there is a combinatorial explosion of orderings arriving at the same final result? Also, how do these automatic orderings compare with the actual behaviour of human translators? Current models rely on manually built biases or are left to explore all possibilities on their own. In this paper, we analyze the orderings produced by human post editors and use them to train an automatic post editing system. We compare the resulting system with those trained with left to right and random post editing orderings. We observe that humans tend to follow a nearly left to right order, but with interesting deviations, such as preferring to start by correcting punctuation or verbs.


On highly accurate computations of the Coulomb two center systems with unit charges. Results of our recent highly accurate computations of the Coulomb two center systems with the unit electrical charges X X e and X Y e are discussed. In particular, we have determined (to very high accuracy) the total energies of the ground 1 s sigma states in the two center adiabatic (or molecular) H 2 , D 2 , HD , HT , T 2 and DT ions. In these computations we applied the new masses of hydrogen isotopes, which have been measured in the recent high energy experiments. We also derived (and tested) our accurate mass interpolation formula for the total energies of the model two center X X e ions with very heavy masses of the point X particles ( M X ge 100,000 m e ). By using this formula we analyze the overall accuracy and validity of the adiabatic two center approximation in application to the three body atomic (or Coulomb) systems.


On iterations of rational functions over perfect fields. Let mathbb K be a perfect field of characterstic p ge 0 and let R in mathbb K(x) be a rational function. This paper studies the number Delta alpha, R (n) of distinct solutions of R (n) (x) alpha over the algebraic closure overline mathbb K of mathbb K , where alpha in overline mathbb K and R (n) is the n fold composition of R with itself. With the exception of some pairs ( alpha, R) , we prove that Delta alpha, R (n) c alpha, R cdot d n O alpha, R (1) for some 0 c alpha, R le 1 d . The number d is readily obtained from R and we provide estimates on c alpha, R . Moreover we prove that the exceptional pairs ( alpha, R) satisfy Delta alpha, R (n) le 2 for every n ge 0 , and we fully describe them. We also discuss further questions and propose some problems in the case where mathbb K is finite.


An Improved Morse Index Bound of Min Max Minimal Hypersurfaces. In this paper, we give an improved Morse index bound of minimal hypersurfaces from Almgren Pitts min max construction in any closed Riemannian manifold M n 1 (n 1 geq 3 ), which generalizes a result by X. Zhou cite zhou multiplicity 2019 for 3 leq n 1 leq 7 . The novel techniques are the construction of hierarchical deformations and the restrictive min max theory. These techniques do not rely on bumpy metrics, and thus could be adapted to many other min max settings.


GoodPoint unsupervised learning of keypoint detection and description. This paper introduces a new algorithm for unsupervised learning of keypoint detectors and descriptors, which demonstrates fast convergence and good performance across different datasets. The training procedure uses homographic transformation of images. The proposed model learns to detect points and generate descriptors on pairs of transformed images, which are easy for it to distinguish and repeatedly detect. The trained model follows SuperPoint architecture for ease of comparison, and demonstrates similar performance on natural images from HPatches dataset, and better performance on retina images from Fundus Image Registration Dataset, which contain low number of corner like features. For HPatches and other datasets, coverage was also computed to provide better estimation of model quality.


The automorphism groups of the profinite braid groups. In this paper we determine the automorphism groups of the profinite braid groups with four or more strings in terms of the profinite Grothendieck Teichm "uller group.


Evidence for increasing frequency of extreme coastal sea levels. Projections of extreme sea levels (ESLs) are critical for managing coastal risks, but are made complicated by deep uncertainties. One key uncertainty is the choice of model structure used to estimate coastal hazards. Differences in model structural choices contribute to uncertainty in estimated coastal hazard, so it is important to characterize how model structural choice affects estimates of ESL. Here, we present a collection of 36 ESL data sets, from tide gauge stations along the United States East and Gulf Coasts. The data are processed using both annual block maxima and peaks over thresholds approaches for modeling distributions of extremes. We use these data sets to fit a suite of potentially nonstationary extreme value models by covarying the ESL statistics with multiple climate variables. We demonstrate how this data set enables inquiry into deep uncertainty surrounding coastal hazards. For all of the models and sites considered here, we find that accounting for changes in the frequency of coastal extreme sea levels provides a better fit than using a stationary extreme value model.