An entrainment based model for annular wakes, with applications to airborne wind energy. Several novel wind energy systems produce wakes with annular cross sections, which are qualitatively different from the wakes with circular cross sections commonly generated by conventional horizontal axis wind turbines and by compact obstacles. Since wind farms use arrays of hundreds of turbines, good analytical wake models are essential for efficient wind farm planning. Several models already exist for circular wakes however, none have yet been proposed for annular wakes, making it impossible to estimate their array performance. We use the entrainment hypothesis to develop a reduced order model for the shape and flow velocity of an annular wake from a generic annular obstacle. Our model consists of a set of three ordinary differential equations, which we solve numerically. In addition, by assuming that the annular wake does not drift radially, we further reduce the problem to a model comprising only two differential equations, which we solve analytically. Both of our models are in good agreement with previously published large eddy simulation results.


Learning Context Graph for Person Search. Person re identification has achieved great progress with deep convolutional neural networks. However, most previous methods focus on learning individual appearance feature embedding, and it is hard for the models to handle difficult situations with different illumination, large pose variance and occlusion. In this work, we take a step further and consider employing context information for person search. For a probe gallery pair, we first propose a contextual instance expansion module, which employs a relative attention module to search and filter useful context information in the scene. We also build a graph learning framework to effectively employ context pairs to update target similarity. These two modules are built on top of a joint detection and instance feature learning framework, which improves the discriminativeness of the learned features. The proposed framework achieves state of the art performance on two widely used person search datasets.


Exact Model Reduction and Fast Forced Response Calculation in High Dimensional Nonlinear Mechanical Systems. We show how spectral submanifold (SSM) theory can be used to extract forced response curves, including isolas, without any numerical simulation in high degree of freedom, periodically forced mechanical systems. We use multivariate recurrence relations to construct the SSMs, achieving a major speed up relative to earlier autonomous SSM algorithms. The increase in computational efficiency promises to close the current gap between studying lower dimensional academic examples and analyzing larger systems obtained from finite element modeling, as we illustrate on a discretization of a damped forced beam model.


A wrinkle in time A case study in DNS poisoning. The Domain Name System (DNS) provides a translation between readable domain names and IP addresses. The DNS is a key infrastructure component of the Internet and a prime target for a variety of attacks. One of the most significant threat to the DNS's wellbeing is a DNS poisoning attack, in which the DNS responses are maliciously replaced, or poisoned, by an attacker. To identify this kind of attack, we start by an analysis of different kinds of response times. We present an analysis of typical and atypical response times, while differentiating between the different levels of DNS servers' response times, from root servers down to internal caching servers. We successfully identify empirical DNS poisoning attacks based on a novel method for DNS response timing analysis. We then present a system we developed to validate our technique that does not require any changes to the DNS protocol or any existing network equipment. Our validation system tested data from different architectures including LAN and cloud environments and real data from an Internet Service Provider (ISP). Our method and system differ from most other DNS poisoning detection methods and achieved high detection rates exceeding 99 . These findings suggest that when used in conjunction with other methods, they can considerably enhance the accuracy of these methods.


The Neumann problem for the fractional Laplacian regularity up to the boundary. We study the regularity up to the boundary of solutions to the Neumann problem for the fractional Laplacian. We prove that if u is a weak solution of ( Delta) s u f in Omega , mathcal N s u 0 in Omega c , then u is C alpha up tp the boundary for some alpha 0 . Moreover, in case s frac12 , we then show that u in C 2s 1 alpha ( overline Omega) . To prove these results we need, among other things, a delicate Moser iteration on the boundary with some logarithmic corrections. Our methods allow us to treat as well the Neumann problem for the regional fractional Laplacian, and we establish the same boundary regularity result. Prior to our results, the interior regularity for these Neumann problems was well understood, but near the boundary even the continuity of solutions was open.


Human robot Collaborative Navigation Search using Social Reward Sources. This paper proposes a Social Reward Sources (SRS) design for a Human Robot Collaborative Navigation (HRCN) task human robot collaborative search. It is a flexible approach capable of handling the collaborative task, human robot interaction and environment restrictions, all integrated on a common environment. We modelled task rewards based on unexplored area observability and isolation and evaluated the model through different levels of human robot communication. The models are validated through quantitative evaluation against both agents' individual performance and qualitative surveying of participants' perception. After that, the three proposed communication levels are compared against each other using the previous metrics.


Filiform Lie algebras with low derived length. We construct, for any integer n greater than or equal to 5, a family of complex filiform Lie algebras with derived length at most 3 and dimension n. We also give examples of n dimensional filiform Lie algebras with derived length greater than 3.


Multi scale Deep Neural Network (MscaleDNN) Methods for Oscillatory Stokes Flows in Complex Domains. In this paper, we study a multi scale deep neural network (MscaleDNN) as a meshless numerical method for computing oscillatory Stokes flows in complex domains. The MscaleDNN employs a multi scale structure in the design of its DNN using radial scalings to convert the approximation of high frequency components of the highly oscillatory Stokes solution to one of lower frequencies. The MscaleDNN solution to the Stokes problem is obtained by minimizing a loss function in terms of L2 normof the residual of the Stokes equation. Three forms of loss functions are investigated based on vorticity velocity pressure, velocity stress pressure, and velocity gradient of velocity pressure formulations of the Stokes equation. We first conduct a systematic study of the MscaleDNN methods with various loss functions on the Kovasznay flow in comparison with normal fully connected DNNs. Then, Stokes flows with highly oscillatory solutions in a 2 D domain with six randomly placed holes are simulated by the MscaleDNN. The results show that MscaleDNN has faster convergence and consistent error decays in the simulation of Kovasznay flow for all four tested loss functions. More importantly, the MscaleDNN is capable of learning highly oscillatory solutions when the normal DNNs fail to converge.


Relative phantom maps and rational homotopy. We generalize some results of Gray and McGibbon Roitberg on relations between phantom maps and rational homotopy to relative phantom maps. Since the lim 1 and the profinite completion techniques do not apply to relative phantom maps, we develop new techniques.


An Exploratory Study of COVID 19 Information on Twitter in the Greater Region. The outbreak of the COVID 19 leads to a burst of information in major online social networks (OSNs). Facing this constantly changing situation, OSNs have become an essential platform for people expressing opinions and seeking up to the minute information. Thus, discussions on OSNs may become a reflection of reality. This paper aims to figure out the distinctive characteristics of the Greater Region (GR) through conducting a data driven exploratory study of Twitter COVID 19 information in the GR and related countries using machine learning and representation learning methods. We find that tweets volume and COVID 19 cases in GR and related countries are correlated, but this correlation only exists in a particular period of the pandemic. Moreover, we plot the changing of topics in each country and region from 2020 01 22 to 2020 06 05, figuring out the main differences between GR and related countries.


Deep Reinforcement Learning for Active Flow Control around a Circular Cylinder Using Unsteady mode Plasma Actuators. Deep reinforcement learning (DRL) algorithms are rapidly making inroads into fluid mechanics, following the remarkable achievements of these techniques in a wide range of science and engineering applications. In this paper, a deep reinforcement learning (DRL) agent has been employed to train an artificial neural network (ANN) using computational fluid dynamics (CFD) data to perform active flow control (AFC) around a 2 D circular cylinder. Flow control strategies are investigated at a diameter based Reynolds number Re D 100 using advantage actor critic (A2C) algorithm by means of two symmetric plasma actuators located on the surface of the cylinder near the separation point. The DRL agent interacts with the computational fluid dynamics (CFD) environment through manipulating the non dimensional burst frequency (f ) of the two plasma actuators, and the time averaged surface pressure is used as a feedback observation to the deep neural networks (DNNs). The results show that a regular actuation using a constant non dimensional burst frequency gives a maximum drag reduction of 21.8 , while the DRL agent is able to learn a control strategy that achieves a drag reduction of 22.6 . By analyzing the flow field, it is shown that the drag reduction is accompanied with a strong flow reattachment and a significant reduction in the mean velocity magnitude and velocity fluctuations at the wake region. These outcomes prove the great capabilities of the deep reinforcement learning (DRL) paradigm in performing active flow control (AFC), and pave the way toward developing robust flow control strategies for real life applications.


Lyapunov exponents and stability properties of higher rank representations. Generalizing results of Deroin Dujardin, we introduce the notion of proximal stability for a holomorphic family ( rho lambda ) of representations Gamma to text SL (d, mathbb C ) , where Gamma is a finitely generated group, and show that it is equivalent to the pluriharmonicity of Lyapunov exponents of the family (defined using random walks).


Learning Tactile Models for Factor Graph based Estimation. We're interested in the problem of estimating object states from touch during manipulation under occlusions. In this work, we address the problem of estimating object poses from touch during planar pushing. Vision based tactile sensors provide rich, local image measurements at the point of contact. A single such measurement, however, contains limited information and multiple measurements are needed to infer latent object state. We solve this inference problem using a factor graph. In order to incorporate tactile measurements in the graph, we need local observation models that can map high dimensional tactile images onto a low dimensional state space. Prior work has used low dimensional force measurements or engineered functions to interpret tactile measurements. These methods, however, can be brittle and difficult to scale across objects and sensors. Our key insight is to directly learn tactile observation models that predict the relative pose of the sensor given a pair of tactile images. These relative poses can then be incorporated as factors within a factor graph. We propose a two stage approach first we learn local tactile observation models supervised with ground truth data, and then integrate these models along with physics and geometric factors within a factor graph optimizer. We demonstrate reliable object tracking using only tactile feedback for 150 real world planar pushing sequences with varying trajectories across three object shapes. Supplementary video https youtu.be y1kBfSmi8w0


Presenting Simultaneous Translation in Limited Space. Some methods of automatic simultaneous translation of a long form speech allow revisions of outputs, trading accuracy for low latency. Deploying these systems for users faces the problem of presenting subtitles in a limited space, such as two lines on a television screen. The subtitles must be shown promptly, incrementally, and with adequate time for reading. We provide an algorithm for subtitling. Furthermore, we propose a way how to estimate the overall usability of the combination of automatic translation and subtitling by measuring the quality, latency, and stability on a test set, and propose an improved measure for translation latency.


Decompositions of moduli spaces of vector bundles and graph potentials. We propose a conjectural semiorthogonal decomposition for the derived category of the moduli space of stable rank 2 bundles with fixed determinant of odd degree, independently formulated by Narasimhan. We discuss some evidence for, and furthermore propose semiorthogonal decompositions with additional structure. We also discuss two other decompositions. One is a decomposition of this moduli space in the Grothendieck ring of varieties, which relates to various known motivic decompositions. The other is the critical value decomposition of a candidate mirror Landau Ginzburg model given by graph potentials, which in turn is related under mirror symmetry to Munoz's decomposition of quantum cohomology. This corresponds to an orthogonal decomposition of the Fukaya category. We will explain how these decompositions can be seen as evidence for the conjectural semiorthogonal decomposition.


Towards Accurate Quantization and Pruning via Data free Knowledge Transfer. When large scale training data is available, one can obtain compact and accurate networks to be deployed in resource constrained environments effectively through quantization and pruning. However, training data are often protected due to privacy concerns and it is challenging to obtain compact networks without data. We study data free quantization and pruning by transferring knowledge from trained large networks to compact networks. Auxiliary generators are simultaneously and adversarially trained with the targeted compact networks to generate synthetic inputs that maximize the discrepancy between the given large network and its quantized or pruned version. We show theoretically that the alternating optimization for the underlying minimax problem converges under mild conditions for pruning and quantization. Our data free compact networks achieve competitive accuracy to networks trained and fine tuned with training data. Our quantized and pruned networks achieve good performance while being more compact and lightweight. Further, we demonstrate that the compact structure and corresponding initialization from the Lottery Ticket Hypothesis can also help in data free training.


Reverse engineering neural networks from many partial recordings. Much of neuroscience aims at reverse engineering the brain, but we only record a small number of neurons at a time. We do not currently know if reverse engineering the brain requires us to simultaneously record most neurons or if multiple recordings from smaller subsets suffice. This is made even more important by the development of novel techniques that allow recording from selected subsets of neurons, e.g. using optical techniques. To get at this question, we analyze a neural network, trained on the MNIST dataset, using only partial recordings and characterize the dependency of the quality of our reverse engineering on the number of simultaneously recorded "neurons". We find that reverse engineering of the nonlinear neural network is meaningfully possible if a sufficiently large number of neurons is simultaneously recorded but that this number can be considerably smaller than the number of neurons. Moreover, recording many times from small random subsets of neurons yields surprisingly good performance. Application in neuroscience suggests to approximate the I O function of an actual neural system, we need to record from a much larger number of neurons. The kind of scaling analysis we perform here can, and arguably should be used to calibrate approaches that can dramatically scale up the size of recorded data sets in neuroscience.


Network Optimization via Smooth Exact Penalty Functions Enabled by Distributed Gradient Computation. This paper proposes a distributed algorithm for a network of agents to solve an optimization problem with separable objective function and locally coupled constraints. Our strategy is based on reformulating the original constrained problem as the unconstrained optimization of a smooth (continuously differentiable) exact penalty function. Computing the gradient of this penalty function in a distributed way is challenging even under the separability assumptions on the original optimization problem. Our technical approach shows that the distributed computation problem for the gradient can be formulated as a system of linear algebraic equations defined by separable problem data. To solve it, we design an exponentially fast, input to state stable distributed algorithm that does not require the individual agent matrices to be invertible. We employ this strategy to compute the gradient of the penalty function at the current network state. Our distributed algorithmic solver for the original constrained optimization problem interconnects this estimation with the prescription of having the agents follow the resulting direction. Numerical simulations illustrate the convergence and robustness properties of the proposed algorithm.


A Flexible and Intelligent Framework for Remote Health Monitoring Dashboards. Developing and maintaining monitoring panels is undoubtedly the main task in the remote patient monitoring (RPM) systems. Due to the significant variations in desired functionalities, data sources, and objectives, designing an efficient dashboard that responds to the various needs in an RPM project is generally a cumbersome task to carry out. In this work, we present ViSierra, a framework for designing data monitoring dashboards in RPM projects. The abstractions and different components of this open source project are explained, and examples are provided to support our claim concerning the effectiveness of this framework in preparing fast, efficient, and accurate monitoring platforms with minimal coding. These platforms will cover all the necessary aspects in a traditional RPM project and combine them with novel functionalities such as machine learning solutions and provide better data analysis instruments for the experts to track the information.


Automated Market Makers for Decentralized Finance (DeFi). This paper compares mathematical models for automated market makers including logarithmic market scoring rule (LMSR), liquidity sensitive LMSR (LS LMSR), constant product mean sum, and others. It is shown that though LMSR may not be a good model for Decentralized Finance (DeFi) applications, LS LMSR has several advantages over constant product mean based automated market makers. However, LS LMSR requires complicated computation (i.e., logarithm and exponentiation) and the cost function curve is concave. In certain DeFi applications, it is preferred to have computationally efficient cost functions with convex curves to conform with the principle of supply and demand. This paper proposes and analyzes constant circle ellipse based cost functions for automated market makers. The proposed cost functions are computationally efficient (only requires multiplication and square root calculation) and have several advantages over widely deployed constant product cost functions. For example, the proposed market makers are more robust against front runner (slippage) attacks.


Extrapolation for multilinear compact operators and applications. This paper is devoted to studying the Rubio de Francia extrapolation for multilinear compact operators. It allows one to extrapolate the compactness of T from just one space to the full range of weighted spaces, whenever an m linear operator T is bounded on weighted Lebesgue spaces. This result is indeed established in terms of the multilinear Muckenhoupt weights A vec p , vec r , and the limited range of the L p scale. To show extrapolation theorems above, by means of a new weighted Fr ' e chet Kolmogorov theorem, we present the weighted interpolation for multilinear compact operators. To prove the latter, we also need to bulid a weighted interpolation theorem in mixed norm Lebesgue spaces. As applications, we obtain the weighted compactness of commutators of many multilinear operators, including multilinear omega Calder ' o n Zygmund operators, multilinear Fourier multipliers, bilinear rough singular integrals and bilinear Bochner Riesz means. Beyond that, we establish the weighted compactness of higher order Calder ' o n commutators, and commutators of Riesz transforms related to Schr " o dinger operators.


Unified Reinforcement Q Learning for Mean Field Game and Control Problems. We present a Reinforcement Learning (RL) algorithm to solve infinite horizon asymptotic Mean Field Game (MFG) and Mean Field Control (MFC) problems. Our approach can be described as a unified two timescale Mean Field Q learning The emph same algorithm can learn either the MFG or the MFC solution by simply tuning the ratio of two learning parameters. The algorithm is in discrete time and space where the agent not only provides an action to the environment but also a distribution of the state in order to take into account the mean field feature of the problem. Importantly, we assume that the agent can not observe the population's distribution and needs to estimate it in a model free manner. The asymptotic MFG and MFC problems are also presented in continuous time and space, and compared with classical (non asymptotic or stationary) MFG and MFC problems. They lead to explicit solutions in the linear quadratic (LQ) case that are used as benchmarks for the results of our algorithm.


Italian domination in generalized Petersen graphs. In a graph G (V,E) , each vertex v in V is labelled with 0 , 1 or 2 such that each vertex labelled with 0 is adjacent to at least one vertex labelled 2 or two vertices labelled 1 . Such kind of labelling is called an Italian dominating function (IDF) of G . The weight of an IDF f is w(f) sum v in V f(v) . The Italian domination number of G is gamma I (G) min f w(f) . Gao et al. (2019) have determined the value of gamma I(P(n,3)) . In this article, we focus on the study of the Italian domination number of generalized Petersen graphs P(n, k) , k neq3 . We determine the values of gamma I(P(n, 1)) , gamma I(P(n, 2)) and gamma I(P(n, k)) for k ge4 , k equiv2,3( bmod5) and n equiv0( bmod5) . For other P(n,k) , we present a bound of gamma I(P(n, k)) . With the obtained results, we partially solve the open problem presented by Bre v s ar et al. (2007) by giving P(n,1) is an example for which gamma I gamma r2 and characterizing P(n,2) for which gamma I(P(n,2)) gamma r2 (P(n,2)) . Moreover, our results imply P(n,1) (n equiv0( bmod 4)) is Italian, P(n,1) (n not equiv0( bmod 4)) and P(n,2) are not Italian.


Third ArchEdge Workshop Exploring the Design Space of Efficient Deep Neural Networks. This paper gives an overview of our ongoing work on the design space exploration of efficient deep neural networks (DNNs). Specifically, we cover two aspects (1) static architecture design efficiency and (2) dynamic model execution efficiency. For static architecture design, different from existing end to end hardware modeling assumptions, we conduct full stack profiling at the GPU core level to identify better accuracy latency trade offs for DNN designs. For dynamic model execution, different from prior work that tackles model redundancy at the DNN channels level, we explore a new dimension of DNN feature map redundancy to be dynamically traversed at runtime. Last, we highlight several open questions that are poised to draw research attention in the next few years.


Non Archimedean entire curves in projective varieties dominating an elliptic curve. Let K be an algebraically closed, complete, non Archimedean valued field of characteristic zero. We prove the non Archimedean Green Griffiths Lang conjecture for projective surfaces of irregularity one. More precisely, we prove that if X K is a groupless, projective surface that admits a dominant morphism an elliptic curve, then X is K analytically Brody hyperbolic. The main ingredient in our proof is a theorem concerning the algebraic degeneracy of non Archimedean entire curves in projective, pseudo groupless varieties admitting a dominant morphism to an elliptic curve.


Phase transition for extremes of a stochastic model with long range dependence and multiplicative noise. We consider a stochastic process with long range dependence perturbed by multiplicative noise. The marginal distributions of both the original process and the noise have regularly varying tails, with tail indices alpha, alpha' 0 , respectively. The original process is taken as the regularly varying Karlin model, a recently investigated model that has long range dependence characterized by a memory parameter beta in(0,1) . We establish limit theorems for the extremes of the model, and reveal a phase transition. In terms of the limit there are three different regimes signal dominance regime alpha alpha' beta , noise dominance regime alpha alpha' beta , and critical regime alpha alpha' beta . As for the proof, we actually establish the same phase transition phenomena for the so called Poisson Karlin model with multiplicative noise defined on generic metric spaces, and apply a Poissonization method to establish the limit theorems for the one dimensional case as a consequence.


Evolutionary Dynamic Multi objective Optimization Via Regression Transfer Learning. Dynamic multi objective optimization problems (DMOPs) remain a challenge to be settled, because of conflicting objective functions change over time. In recent years, transfer learning has been proven to be a kind of effective approach in solving DMOPs. In this paper, a novel transfer learning based dynamic multi objective optimization algorithm (DMOA) is proposed called regression transfer learning prediction based DMOA (RTLP DMOA). The algorithm aims to generate an excellent initial population to accelerate the evolutionary process and improve the evolutionary performance in solving DMOPs. When an environmental change is detected, a regression transfer learning prediction model is constructed by reusing the historical population, which can predict objective values. Then, with the assistance of this prediction model, some high quality solutions with better predicted objective values are selected as the initial population, which can improve the performance of the evolutionary process. We compare the proposed algorithm with three state of the art algorithms on benchmark functions. Experimental results indicate that the proposed algorithm can significantly enhance the performance of static multi objective optimization algorithms and is competitive in convergence and diversity.


Liouville results for fully nonlinear equations modeled on H "ormander vector fields. I. The Heisenberg group. This paper studies Liouville properties for viscosity sub and supersolutions of fully nonlinear degenerate elliptic PDEs, under the main assumption that the operator has a family of generalized subunit vector fields that satisfy the H "ormander condition. A general set of sufficient conditions is given such that all subsolutions bounded above are constant it includes the existence of a supersolution out of a big ball, that explodes at infinity. Therefore for a large class of operators the problem is reduced to finding such a Lyapunov like function. This is done here for the vector fields that generate the Heisenberg group, giving explicit conditions on the sign and size of the first and zero th order terms in the equation. The optimality of the conditions is shown via several examples. A sequel of this paper applies the methods to other Carnot groups and to Grushin geometries.


Towards Verified Stochastic Variational Inference for Probabilistic Programs. Probabilistic programming is the idea of writing models from statistics and machine learning using program notations and reasoning about these models using generic inference engines. Recently its combination with deep learning has been explored intensely, leading to the development of deep probabilistic programming languages such as Pyro. At the core of this development lie inference engines based on stochastic variational inference algorithms. When asked to find information about the posterior distribution of a model written in such a language, these algorithms convert this posterior inference query into an optimisation problem and solve it approximately by gradient ascent. In this paper, we analyse one of the most fundamental and versatile variational inference algorithms, called score estimator, using tools from denotational semantics and program analysis. We formally express what this algorithm does on models denoted by programs, and expose implicit assumptions made by the algorithm. The violation of these assumptions may lead to an undefined optimisation objective or the loss of convergence guarantee of the optimisation process. We then describe rules for proving these assumptions, which can be automated by static program analyses. Some of our rules use nontrivial facts from continuous mathematics, and let us replace requirements about integrals in the assumptions, by conditions involving differentiation or boundedness, which are much easier to prove automatically. Following our general methodology, we have developed a static program analysis for Pyro that aims at discharging the assumption about what we call model guide support match. Applied to the eight representative model guide pairs from the Pyro webpage, our analysis finds a bug in one of these cases, reveals a non standard use of an inference engine in another, and shows the assumptions are met in the remaining cases.


FedEval A Benchmark System with a Comprehensive Evaluation Model for Federated Learning. As an innovative solution for privacy preserving machine learning (ML), federated learning (FL) is attracting much attention from research and industry areas. While new technologies proposed in the past few years do evolve the FL area, unfortunately, the evaluation results presented in these works fall short in integrity and are hardly comparable because of the inconsistent evaluation metrics and the lack of a common platform. In this paper, we propose a comprehensive evaluation framework for FL systems. Specifically, we first introduce the ACTPR model, which defines five metrics that cannot be excluded in FL evaluation, including Accuracy, Communication, Time efficiency, Privacy, and Robustness. Then we design and implement a benchmarking system called FedEval, which enables the systematic evaluation and comparison of existing works under consistent experimental conditions. We then provide an in depth benchmarking study between the two most widely used FL mechanisms, FedSGD and FedAvg. The benchmarking results show that FedSGD and FedAvg both have advantages and disadvantages under the ACTPR model. For example, FedSGD is barely influenced by the none independent and identically distributed (non IID) data problem, but FedAvg suffers from a decline in accuracy of up to 9 in our experiments. On the other hand, FedAvg is more efficient than FedSGD regarding time consumption and communication. Lastly, we excavate a set of take away conclusions, which are very helpful for researchers in the FL area.


Sparsifying Transformer Models with Trainable Representation Pooling. We propose a novel method to sparsify attention in the Transformer model by learning to select the most informative token representations during the training process, thus focusing on the task specific parts of an input. A reduction of quadratic time and memory complexity to sublinear was achieved due to a robust trainable top k operator. Our experiments on a challenging long document summarization task show that even our simple baseline performs comparably to the current SOTA, and with trainable pooling, we can retain its top quality, while being 1.8 times faster during training, 4.5 times faster during inference, and up to 13 times more computationally efficient in the decoder.


Contribution of Magnetic Reconnection Events to Energy Dissipation in Magnetosheath Turbulence. By analyzing the magnetosheath measurements from MMS, we obtain the statistical results for the contribution of magnetic reconnection (MR) events at electron scales to the energy dissipation of coherent structures. The Partial Variance of Increments (PVI) method is employed to find coherent structures in the magnetic field data. The current sheet structures with reversal of magnetic field components are further selected. We consider the following criteria to identify the MR events, such as current sheet with magnetic field reversal, significant energy dissipation, and evident electron outflow velocity. Statistically, for most MR events, their PVI values are larger than that of other types of coherent structures, and their energy dissipations are also stronger than that of others. However, due to the relatively small proportion of MR events, their contribution to coherent structures' energy dissipation is relatively trivial. If taken into account the dissipation of non coherent structures, the MR's contribution to energy dissipation would be less. Hence, we suggest that MR events, though have strong dissipation locally, are not the major contributor to the energy dissipation in the magnetosheath. After analyzing the features of non MR current sheets, we propose that non MR current sheets are mainly coherent structures inherent to kinetic Alfv 'en fluctuations.


Simple dynamic word embeddings for mapping perceptions in the public sphere. Word embeddings trained on large scale historical corpora can illuminate human biases and stereotypes that perpetuate social inequalities. These embeddings are often trained in separate vector space models defined according to different attributes of interest. In this paper, we develop a unified dynamic embedding model that learns attribute specific word embeddings. We apply our model to investigate i) 20th century gender and ethnic occupation biases embedded in the Corpus of Historical American English (COHA), and ii) biases against refugees embedded in a novel corpus of talk radio transcripts containing 119 million words produced over one month across 83 stations and 64 cities. Our results shed preliminary light on scenarios when dynamic embedding models may be more suitable for representing linguistic biases than individual vector space models, and vice versa.


Faster Algorithms for Optimal Ex Ante Coordinated Collusive Strategies in Extensive Form Zero Sum Games. We focus on the problem of finding an optimal strategy for a team of two players that faces an opponent in an imperfect information zero sum extensive form game. Team members are not allowed to communicate during play but can coordinate before the game. In that setting, it is known that the best the team can do is sample a profile of potentially randomized strategies (one per player) from a joint (a.k.a. correlated) probability distribution at the beginning of the game. In this paper, we first provide new modeling results about computing such an optimal distribution by drawing a connection to a different literature on extensive form correlation. Second, we provide an algorithm that computes such an optimal distribution by only using profiles where only one of the team members gets to randomize in each profile. We can also cap the number of such profiles we allow in the solution. This begets an anytime algorithm by increasing the cap. We find that often a handful of well chosen such profiles suffices to reach optimal utility for the team. This enables team members to reach coordination through a relatively simple and understandable plan. Finally, inspired by this observation and leveraging theoretical concepts that we introduce, we develop an efficient column generation algorithm for finding an optimal distribution for the team. We evaluate it on a suite of common benchmark games. It is three orders of magnitude faster than the prior state of the art on games that the latter can solve and it can also solve several games that were previously unsolvable.


Emergence of irregular activity in networks of strongly coupled conductance based neurons. Cortical neurons are characterized by irregular firing and a broad distribution of rates. The balanced state model explains these observations with a cancellation of mean excitatory and inhibitory currents, which makes fluctuations drive firing. In networks of neurons with current based synapses, the balanced state emerges dynamically if coupling is strong, i.e. if the mean number of synapses per neuron K is large and synaptic efficacy is of order 1 sqrt K . When synapses are conductance based, current fluctuations are suppressed when coupling is strong, questioning the applicability of the balanced state idea to biological neural networks. We analyze networks of strongly coupled conductance based neurons and show that asynchronous irregular activity and broad distributions of rates emerge if synapses are of order 1 log(K) . In such networks, unlike in the standard balanced state model, current fluctuations are small and firing is maintained by a drift diffusion balance. This balance emerges dynamically, without fine tuning, if inputs are smaller than a critical value, which depends on synaptic time constants and coupling strength, and is significantly more robust to connection heterogeneities than the classical balanced state model. Our analysis makes experimentally testable predictions of how the network response properties should evolve as input increases.


ChebNet Efficient and Stable Constructions of Deep Neural Networks with Rectified Power Units using Chebyshev Approximations. In a recent paper B. Li, S. Tang and H. Yu, arXiv 1903.05858 , it was shown that deep neural networks built with rectified power units (RePU) can give better approximation for sufficient smooth functions than those with rectified linear units, by converting polynomial approximation given in power series into deep neural networks with optimal complexity and no approximation error. However, in practice, power series are not easy to compute. In this paper, we propose a new and more stable way to construct deep RePU neural networks based on Chebyshev polynomial approximations. By using a hierarchical structure of Chebyshev polynomial approximation in frequency domain, we build efficient and stable deep neural network constructions. In theory, ChebNets and the deep RePU nets based on Power series have the same upper error bounds for general function approximations. But numerically, ChebNets are much more stable. Numerical results show that the constructed ChebNets can be further trained and obtain much better results than those obtained by training deep RePU nets constructed basing on power series.


Advanced Control Methods for Particle Accelerators (ACM4PA) 2019 Workshop Report. Los Alamos is currently developing novel particle accelerator controls and diagnostics algorithms to enable higher quality beams with lower beam losses than is currently possible. The purpose of this workshop was to consider tuning and optimization challenges of a wide range of particle accelerators including linear proton accelerators such as the Los Alamos Neutron Science Center (LANSCE), rings such as the Advanced Photon Source (APS) synchrotron, free electron lasers (FEL) such as the Linac Coherent Light Source (LCLS) and LCLS II, the European X ray Free Electron Laser (EuXFEL), the Swiss FEL, and the planned MaRIE FEL, and plasma wake field accelerators such as FACET, FACET II, and AWAKE at CERN. One major challenge is an the ability to quickly create very high quality, extremely intense, custom current and energy profile beams while working with limited real time non invasive diagnostics and utilizing time varying uncertain initial beam distributions and accelerator components. Currently, a few individual accelerator labs have been developing and applying their own diagnostics tools and custom control and ML algorithms for automated machine tuning and optimization. The goal of this workshop was to bring together a group of accelerator physicists and accelerator related control and ML experts in order to define which controls and diagnostics would be most useful for existing and future accelerators and to create a plan for developing a new family of algorithms that can be shared and maintained by the community.


Reconciling vacuum laser acceleration theory and experiment. The classical theory of single electron dynamics in focused laser pulses is the foundation of both the relativistic ponderomotive force (RPF), which in turn underlies models of laser collective plasma dynamics, and the discovery of novel strong field radiation dynamics. Despite this bedrock importance, consensus eludes the community as to whether acceleration of single electrons in vacuum has been observed in experiment. We analyze the experiment of Malka et al. (1998) with respect to several features that were neglected in modeling and that can restore consistency between theory predictions and experimental data. The right or wrong pulse profile function, laser parameters, or initial electron distribution each can make or break the agreement between predictions and data. The laser phase at which the electron's interaction with the pulse begins has a large effect, explaining why much larger energies are achieved by electrons liberated in the focal region by photoionization from high Z atoms and by electrons ejected from a plasma mirror. Finally we estimate the error in a typical electron spectrum arising from fluctuating focal spot size in state of the art ultra relativistic laser facilities. Our results emphasize the importance of thoroughly characterizing laser parameters in order to achieve quantitatively accurate predictions and the precision required for discovery science.


Hexagonal Transverse Coupled Cavity VCSEL Redefining the High Speed Lasers. The vertical cavity surface emitting lasers (VCSELs) have emerged as a vital approach for realizing energy efficient, high speed optical interconnects in the data center and supercomputers. As of today, VCSEL is the most suitable for mass production in terms of cost effectiveness and reliability. However, there are still key challenges for higher speed modulation above 40 GHz. Here, a hexagonal transverse coupled cavity VCSEL adiabatically coupled through the center cavity is proposed. A 3 dB roll off modulation bandwidth of 45 GHz is demonstrated, which is five times greater than a conventional VCSEL fabricated on the same epi wafer structure. While a parity time (PT) symmetry approaches add loss to engineer the topological state of the laser system, here, a radical paradigm shift with gain introduces symmetry breaking. This idea, then enables a single mode operation with a side mode suppression ratio (SMSR) of 30 decibels and signal to noise ratio (SNR) of 45 decibels. The energy distribution inside the coupled cavity system is also redistributed to provide a coherent gain in a spatially separated system. Consequently, throughput power is three times higher than that of the conventional VCSEL.


Access Control Mechanisms in Named Data Networks A Comprehensive Survey. Information Centric Networking (ICN) has recently emerged as a prominent candidate for the Future Internet Architecture (FIA) that addresses existing issues with the host centric communication model of the current TCP IP based Internet. Named Data Networking (NDN) is one of the most recent and active ICN architectures that provides a clean slate approach for Internet communication. NDN provides intrinsic content security where security is directly provided to the content instead of communication channel. Among other security aspects, Access Control (AC) rules specify the privileges for the entities that can access the content. In TCP IP based AC systems, due to the client server communication model, the servers control which client can access a particular content. In contrast, ICN based networks use content names to drive communication and decouple the content from its original location. This phenomenon leads to the loss of control over the content causing different challenges for the realization of efficient AC mechanisms. To date, considerable efforts have been made to develop various AC mechanisms in NDN. In this paper, we provide a detailed and comprehensive survey of the AC mechanisms in NDN. We follow a holistic approach towards AC in NDN where we first summarize the ICN paradigm, describe the changes from channel based security to content based security and highlight different cryptographic algorithms and security protocols in NDN. We then classify the existing AC mechanisms into two main categories Encryption based AC and Encryption independent AC. Each category has different classes based on the working principle of AC (e.g., Attribute based AC, Name based AC, Identity based AC, etc). Finally, we present the lessons learned from the existing AC mechanisms and identify the challenges of NDN based AC at large, highlighting future research directions for the community.


A note on singular Hermitian Yang Mills connections. We give an example of a homogeneous reflexive sheaf over mathbb C 3 which admits a non conical Hermitian Yang Mills connection. This is expected to model bubbling phenomenon along complex codimension 2 submanifolds when the Fueter section takes zero value.


An Ethical Highlighter for People Centric Dataset Creation. Important ethical concerns arising from computer vision datasets of people have been receiving significant attention, and a number of datasets have been withdrawn as a result. To meet the academic need for people centric datasets, we propose an analytical framework to guide ethical evaluation of existing datasets and to serve future dataset creators in avoiding missteps. Our work is informed by a review and analysis of prior works and highlights where such ethical challenges arise.


A Natural Quadratic Approach to the Generalized Graph Layering Problem. We propose a new exact approach to the generalized graph layering problem that is based on a particular quadratic assignment formulation. It expresses, in a natural way, the associated layout restrictions and several possible objectives, such as a minimum total arc length, minimum number of reversed arcs, and minimum width, or the adaptation to a specific drawing area. Our computational experiments show a competitive performance compared to prior exact models.


Multiple weak solutions for a kind of time dependent equation involving singularity. The existence of at least three weak solutions for a kind of nonlinear time dependent equation is studied. In fact, we consider the case that the source function has singularity at origin. To this aim, the variational methods and the well known critical points theorem are main tools.


Improvement on the crossing number of crossing critical graphs. The crossing number of a graph G is the minimum number of edge crossings over all drawings of G in the plane. A graph G is k crossing critical if its crossing number is at least k , but if we remove any edge of G , its crossing number drops below k . There are examples of k crossing critical graphs that do not have drawings with exactly k crossings. Richter and Thomassen proved in 1993 that if G is k crossing critical, then its crossing number is at most 2.5k 16 . We improve this bound to 2k 6 sqrt k 44 .


Liouville property and non negative Ollivier curvature on graphs. For graphs with non negative Ollivier curvature, we prove the Liouville property, i.e., every bounded harmonic function is constant. Moreover, we improve Ollivier's results on concentration of the measure under positive Ollivier curvature.


Hochschild cohomology of generalised Grassmannians. We compute the Hochschild Kostant Rosenberg decomposition of the Hochschild cohomology of generalised Grassmannians, i.e. partial flag varieties associated to maximal parabolic subgroups in a simple algebraic group. We explain how the decomposition is concentrated in global sections for so called (co)minuscule and (co)adjoint generalised Grassmannians, and we suggest that for (almost) all other cases the same vanishing of the higher cohomology does not hold. Our methods give an explicit partial description of the Gerstenhaber algebra structure for the Hochschild cohomology of cominuscule and adjoint generalised Grassmannians. We also consider the case of adjoint partial flag varieties in type A, which are associated to certain submaximal parabolic subgroups.


A Subjective Interestingness measure for Business Intelligence explorations. This paper addresses the problem of defining a subjective interestingness measure for BI exploration. Such a measure involves prior modeling of the belief of the user. The complexity of this problem lies in the impossibility to ask the user about the degree of belief in each element composing their knowledge prior to the writing of a query. To this aim, we propose to automatically infer this user belief based on the user's past interactions over a data cube, the cube schema and other users past activities. We express the belief under the form of a probability distribution over all the query parts potentially accessible to the user, and use a random walk to learn this distribution. This belief is then used to define a first Subjective Interestingness measure over multidimensional queries. Experiments conducted on simulated and real explorations show how this new subjective interestingness measure relates to prototypical and real user behaviors, and that query parts offer a reasonable proxy to infer user belief.


1 Laplacian type problems with strongly singular nonlinearities and gradient terms. We show optimal existence, nonexistence and regularity results for nonnegative solutions to Dirichlet problems as begin cases displaystyle Delta 1 u g(u) D u h(u)f text in Omega, u 0 text on partial Omega, end cases where Omega is an open bounded subset of mathbb R N , f geq 0 belongs to L N( Omega) , and g and h are continuous functions that may blow up at zero. As a noteworthy fact we show how a non trivial interaction mechanism between the two nonlinearities g and h produces remarkable regularizing effects on the solutions. The sharpness of our main results is discussed through the use of appropriate explicit examples.


Disjunctive cuts in Mixed Integer Conic Optimization. This paper studies disjunctive cutting planes in Mixed Integer Conic Programming. Building on conic duality, we formulate a cut generating conic program for separating disjunctive cuts, and investigate the impact of the normalization condition on its resolution. In particular, we show that a careful selection of normalization guarantees its solvability and conic strong duality. Then, we highlight the shortcomings of separating conic infeasible points in an outer approximation context, and propose conic extensions to the classical lifting and monoidal strengthening procedures. Finally, we assess the computational behavior of various normalization conditions in terms of gap closed, computing time and cut sparsity. In the process, we show that our approach is competitive with the internal lift and project cuts of a state of the art solver.


Off diagonal estimates for commutators of bi parameter singular integrals. We study the boundedness of commutators of bi parameter singular integrals between mixed spaces b,T L p 1 L p 2 to L q 1 L q 2 in the off diagonal situation q i,p i in(1, infty) where we also allow q i not p i. Boundedness is fully characterized for several arrangements of the integrability exponents with some open problems presented.


Dynamics and Allocation of Transaction Cost in Multiagent Industrial Symbiosis. This paper discusses the dynamics of Transaction Cost (TC) in Industrial Symbiosis Institutions (ISI) and provides a fair and stable mechanism for TC allocation among the involved firms in a given ISI. In principle, industrial symbiosis, as an implementation of the circular economy paradigm in the context of industrial relation, is a practice aiming at reducing the material energy footprint of the firm. The well engineered form of this practice is proved to decrease the transaction costs at a collective level. This can be achieved using information systems for identifying potential synergies, evaluating mutually beneficial ones, implementing the contracts, and governing the behavior of the established relations. Then the question is "how to distribute the costs for maintaining such an information system in a fair and stable manner?" We see such a cost as a collective transaction cost and employ an integrated method rooted in cooperative game theory and multiagent systems research to develop a fair and stable allocation mechanism for it. The novelty is twofold in developing analytical multiagent methods for capturing the dynamics of transaction costs in industrial symbiosis and in presenting a novel game theoretic mechanism for its allocation in industrial symbiosis institutions. While the former contributes to the theories of industrial symbiosis (methodological contribution), the latter supports decision makers aiming to specify fair and stable industrial symbiosis contracts (practical contribution).


KLDivNet An unsupervised neural network for multi modality image registration. Multi modality image registration is one of the most underlined processes in medical image analysis. Recently, convolutional neural networks (CNNs) have shown significant potential in deformable registration. However, the lack of voxel wise ground truth challenges the training of CNNs for an accurate registration. In this work, we propose a cross modality similarity metric, based on the KL divergence of image variables, and implement an efficient estimation method using a CNN. This estimation network, referred to as KLDivNet, can be trained unsupervisedly. We then embed the KLDivNet into a registration network to achieve the unsupervised deformable registration for multi modality images. We employed three datasets, i.e., AAL Brain, LiTS Liver and Hospital Liver, with both the intra and inter modality image registration tasks for validation. Results showed that our similarity metric was effective, and the proposed registration network delivered superior performance compared to the state of the art methods.


Homological aspects of oriented hypergraphs. We define an algebraic setup of homology for hypergraphs, which defaults to simplicial homology in the case of graphs, and study its basic properties. As part of our study we define algebraic spanning trees of hypergraphs, along with fundamental cuts and cycles playing their usual roles.


LSM Learning Subspace Minimization for Low level Vision. We study the energy minimization problem in low level vision tasks from a novel perspective. We replace the heuristic regularization term with a learnable subspace constraint, and preserve the data term to exploit domain knowledge derived from the first principle of a task. This learning subspace minimization (LSM) framework unifies the network structures and the parameters for many low level vision tasks, which allows us to train a single network for multiple tasks simultaneously with completely shared parameters, and even generalizes the trained network to an unseen task as long as its data term can be formulated. We demonstrate our LSM framework on four low level tasks including interactive image segmentation, video segmentation, stereo matching, and optical flow, and validate the network on various datasets. The experiments show that the proposed LSM generates state of the art results with smaller model size, faster training convergence, and real time inference.


GeoBlocks A Query Cache Accelerated Data Structure for Spatial Aggregation over Polygons. As individual traffic and public transport in cities are changing, city authorities need to analyze urban geospatial data to improve transportation and infrastructure. To that end, they highly rely on spatial aggregation queries that extract summarized information from point data (e.g., Uber rides) contained in a given polygonal region (e.g., a city neighborhood). To support such queries, current analysis tools either allow only predefined aggregates on predefined regions and are thus unsuitable for exploratory analyses, or access the raw data to compute aggregate results on the fly, which severely limits the interactivity. At the same time, existing pre aggregation techniques are inadequate since they maintain aggregates over rectangular regions. As a result, when applied over arbitrary polygonal regions, they induce an approximation error that cannot be bounded. In this paper, we introduce GeoBlocks, a novel pre aggregating data structure that supports spatial aggregation over arbitrary polygons. GeoBlocks closely approximate polygons using a set of fine grained grid cells and, in contrast to prior work, allow to bound the approximation error by adjusting the cell size. Furthermore, GeoBlocks employ a trie like cache that caches aggregate results of frequently queried regions, thereby dynamically adapting to the skew inherently present in query workloads and improving performance over time. In summary, GeoBlocks outperform on the fly aggregation by up to three orders of magnitude, achieving the sub second query latencies required for interactive exploratory analytics.


The Picard groups of inclusions of C algebras induced by equivalence bimodules. Let A and B be sigma unital C algebras and X and Y an A A equivalence bimodule and a B B equivalence bimodule, respectively. Also, let A rtimes X mathbb Z and B rtimes Y mathbb Z be the crossed products of A and B by X and Y , respectively. Furthermore, let A subset A rtimes X mathbb Z and B subset B rtimes Y mathbb Z be the inclusions of C algebras induced by X and Y , respectively. We suppose that A' cap M(A rtimes X mathbb Z ) mathbb C 1 . In this paper we shall show that the inclusions A subset A rtimes X mathbb Z and B subset B rtimes Y mathbb Z are strongly Morita equivalent if and only if there is an A B equivalence bimodule M such that Y cong widetilde M otimes A X otimes A M or widetilde Y cong widetilde M otimes A X otimes A M as B B equivalence bimodules, where widetilde M and widetilde Y are the dual B A equivalence bimodule and the dual B B equivalence bimodule of M and Y , respectively. Applying this result, we shall compute the Picard group of the inclusion A subset A rtimes X mathbb Z under the assumption that A' cap M(A rtimes X mathbb Z ) mathbb C 1 .


Combinatorial Decision Dags A Natural Computational Model for General Intelligence. A novel computational model (CoDD) utilizing combinatory logic to create higher order decision trees is presented. A theoretical analysis of general intelligence in terms of the formal theory of pattern recognition and pattern formation is outlined, and shown to take especially natural form in the case where patterns are expressed in CoDD language. Relationships between logical entropy and algorithmic information, and Shannon entropy and runtime complexity, are shown to be elucidated by this approach. Extension to the quantum computing case is also briefly discussed.


Crouzeix's Conjecture and related problems. In this paper, we establish several results related to Crouzeix's conjecture. We show that the conjecture holds for contractions with eigenvalues that are sufficiently well separated. This separation is measured by the so called separation constant, which is defined in terms of the pseudohyperbolic metric. Moreover, we study general properties of related extremal functions and associated vectors. Throughout, compressions of the shift serve as illustrating examples which also allow for refined results.


To Be Announced. In this survey we review dynamic epistemic logics with modalities for quantification over information change. Of such logics we present complete axiomatizations, focussing on axioms involving the interaction between knowledge and such quantifiers, we report on their relative expressivity, on decidability and on the complexity of model checking and satisfiability, and on applications. We focus on open problems and new directions for research.


Learning Representations for Predicting Future Activities. Foreseeing the future is one of the key factors of intelligence. It involves understanding of the past and current environment as well as decent experience of its possible dynamics. In this work, we address future prediction at the abstract level of activities. We propose a network module for learning embeddings of the environment's dynamics in a self supervised way. To take the ambiguities and high variances in the future activities into account, we use a multi hypotheses scheme that can represent multiple futures. We demonstrate the approach by classifying future activities on the Epic Kitchens and Breakfast datasets. Moreover, we generate captions that describe the future activities


MIMA MAPPER Induced Manifold Alignment for Semi Supervised Fusion of Optical Image and Polarimetric SAR Data. Multi modal data fusion has recently been shown promise in classification tasks in remote sensing. Optical data and radar data, two important yet intrinsically different data sources, are attracting more and more attention for potential data fusion. It is already widely known that, a machine learning based methodology often yields excellent performance. However, the methodology relies on a large training set, which is very expensive to achieve in remote sensing. The semi supervised manifold alignment (SSMA), a multi modal data fusion algorithm, has been designed to amplify the impact of an existing training set by linking labeled data to unlabeled data via unsupervised techniques. In this paper, we explore the potential of SSMA in fusing optical data and polarimetric SAR data, which are multi sensory data sources. Furthermore, we propose a MAPPER induced manifold alignment (MIMA) for semi supervised fusion of multi sensory data sources. Our proposed method unites SSMA with MAPPER, which is developed from the emerging topological data analysis (TDA) field. To our best knowledge, this is the first time that SSMA has been applied on fusing optical data and SAR data, and also the first time that TDA has been applied in remote sensing. The conventional SSMA derives a topological structure using k nearest neighbor (kNN), while MIMA employs MAPPER, which considers the field knowledge and derives a novel topological structure through the spectral clustering in a data driven fashion. Experiment results on data fusion with respect to land cover land use classification and local climate zone classification suggest superior performance of MIMA.


Heatmap Regression via Randomized Rounding. Heatmap regression has become the mainstream methodology for deep learning based semantic landmark localization, including in facial landmark localization and human pose estimation. Though heatmap regression is robust to large variations in pose, illumination, and occlusion in unconstrained settings, it usually suffers from a sub pixel localization problem. Specifically, considering that the activation point indices in heatmaps are always integers, quantization error thus appears when using heatmaps as the representation of numerical coordinates. Previous methods to overcome the sub pixel localization problem usually rely on high resolution heatmaps. As a result, there is always a trade off between achieving localization accuracy and computational cost, where the computational complexity of heatmap regression depends on the heatmap resolution in a quadratic manner. In this paper, we formally analyze the quantization error of vanilla heatmap regression and propose a simple yet effective quantization system to address the sub pixel localization problem. The proposed quantization system induced by the randomized rounding operation 1) encodes the fractional part of numerical coordinates into the ground truth heatmap using a probabilistic approach during training and 2) decodes the predicted numerical coordinates from a set of activation points during testing. We prove that the proposed quantization system for heatmap regression is unbiased and lossless. Experimental results on popular facial landmark localization datasets (WFLW, 300W, COFW, and AFLW) and human pose estimation datasets (MPII and COCO) demonstrate the effectiveness of the proposed method for efficient and accurate semantic landmark localization. Code is available at http github.com baoshengyu H3R.


81 Kr dating at the Guliya ice cap, Tibetan Plateau. We present radiometric 81 Kr dating results for ice samples collected at the outlets of the Guliya ice cap in the western Kunlun Mountains of the Tibetan Plateau. This first application of 81 Kr dating on mid latitude glacier ice was made possible by recent advances in Atom Trap Trace Analysis, particularly a reduction in the required sample size down to 1 mu L STP of krypton. Eight ice blocks were sampled from the bottom of the glacier at three different sites along the southern edges. The 81 Kr data yield upper age limits in the range of 15 74 ka (90 confidence level). This is an order of magnitude lower than the ages exceeding 500 ka which the previous 36 Cl data suggest for the bottom of the Guliya ice core. It is also significantly lower than the widely used chronology up to 110 ka established for the upper part of the core based on delta 18 O in the ice.


Planning a Sequence of Base Positions for a Mobile Manipulator to Perform Multiple Pick and Place Tasks. In this paper, we present a planner that plans a sequence of base positions for a mobile manipulator to efficiently and robustly collect objects stored in distinct trays. We achieve high efficiency by exploring the common areas where a mobile manipulator can grasp objects stored in multiple trays simultaneously and move the mobile manipulator to the common areas to reduce the time needed for moving the mobile base. We ensure robustness by optimizing the base position with the best clearance to positioning uncertainty so that a mobile manipulator can complete the task even if there is a certain deviation from the planned base positions. Besides, considering different styles of object placement in the tray, we analyze feasible schemes for dynamically updating the base positions based on either the remaining objects or the target objects to be picked in one round of the tasks. In the experiment part, we examine our planner on various scenarios, including different object placement (1) Regularly placed toy objects (2) Randomly placed industrial parts and different schemes for online execution (1) Apply globally static base positions (2) Dynamically update the base positions. The experiment results demonstrate the efficiency, robustness and feasibility of the proposed method.


Complex Far Field Geometries Determine the Stability of Solid Tumor Growth with Chemotaxis. In this paper, we develop a sharp interface tumor growth model to study the effect of the tumor microenvironment using a complex far field geometry that mimics a heterogeneous distribution of vasculature. Together with different nutrient uptake rates inside and outside the tumor, this introduces variability in spatial diffusion gradients. Linear stability analysis suggests that the uptake rate in the tumor microenvironment, together with chemotaxis, may induce unstable growth, especially when the nutrient gradients are large. We investigate the fully nonlinear dynamics using a spectrally accurate boundary integral method. Our nonlinear simulations reveal that vascular heterogeneity plays an important role in the development of morphological instabilities that range from fingering and chain like morphologies to compact, plate like shapes in two dimensions.


Engineering Faster Sorters for Small Sets of Items. Sorting a set of items is a task that can be useful by itself or as a building block for more complex operations. That is why a lot of effort has been put into finding sorting algorithms that sort large sets as fast as possible. But the more sophisticated and complex the algorithms become, the less efficient they are for small sets of items due to large constant factors. We aim to determine if there is a faster way than insertion sort to sort small sets of items to provide a more efficient base case sorter. We looked at sorting networks, at how they can improve the speed of sorting few elements, and how to implement them in an efficient manner by using conditional moves. Since sorting networks need to be implemented explicitly for each set size, providing networks for larger sizes becomes less efficient due to increased code sizes. To also enable the sorting of slightly larger base cases, we adapted sample sort to Register Sample Sort, to break down those larger sets into sizes that can in turn be sorted by sorting networks. From our experiments we found that when sorting only small sets, the sorting networks outperform insertion sort by a factor of at least 1.76 for any array size between six and sixteen, and by a factor of 2.72 on average across all machines and array sizes. When integrating sorting networks as a base case sorter into Quicksort, we achieved far less performance improvements, which is probably due to the networks having a larger code size and cluttering the L1 instruction cache. But for x86 machines with a larger L1 instruction cache of 64 KiB or more, we obtained speedups of 12.7 when using sorting networks as a base case sorter in std sort. In conclusion, the desired improvement in speed could only be achieved under special circumstances, but the results clearly show the potential of using conditional moves in the field of sorting algorithms.


Supporting On Stack Replacement in Unstructured Languages by Loop Reconstruction and Extraction. On stack replacement (OSR) is a common technique employed by dynamic compilers to reduce program warm up time. OSR allows switching from interpreted to compiled code during the execution of this code. The main targets are long running loops, which need to be represented explicitly, with dedicated information about condition and body, to be optimized at run time. Bytecode interpreters, however, represent control flow implicitly via unstructured jumps and thus do not exhibit the required high level loop representation. To enable OSR also for jump based often called unstructured languages, we propose the partial reconstruction of loops in order to explicitly represent them in a bytecode interpreter. Besides an outline of the general idea, we implemented our approach in Sulong, a bytecode interpreter for LLVM bitcode, which allows the execution of C C . We conducted an evaluation with a set of C benchmarks, which showed speed ups in warm up of up to 9x for certain benchmarks. This facilitates execution of programs with long running loops in rarely called functions, which would yield significant slowdown without OSR. While shown with a prototype implementation, the overall idea of our approach is generalizable for all bytecode interpreters.


Hierarchical Affordance Discovery using Intrinsic Motivation. To be capable of lifelong learning in a real life environment, robots have to tackle multiple challenges. Being able to relate physical properties they may observe in their environment to possible interactions they may have is one of them. This skill, named affordance learning, is strongly related to embodiment and is mastered through each person's development each individual learns affordances differently through their own interactions with their surroundings. Current methods for affordance learning usually use either fixed actions to learn these affordances or focus on static setups involving a robotic arm to be operated. In this article, we propose an algorithm using intrinsic motivation to guide the learning of affordances for a mobile robot. This algorithm is capable to autonomously discover, learn and adapt interrelated affordances without pre programmed actions. Once learned, these affordances may be used by the algorithm to plan sequences of actions in order to perform tasks of various difficulties. We then present one experiment and analyse our system before comparing it with other approaches from reinforcement learning and affordance learning.


False Detection (Positives and Negatives) in Object Detection. Object detection is a very important function of visual perception systems. Since the early days of classical object detection based on HOG to modern deep learning based detectors, object detection has improved in accuracy. Two stage detectors usually have higher accuracy than single stage ones. Both types of detectors use some form of quantization of the search space of rectangular regions of image. There are far more of the quantized elements than true objects. The way these bounding boxes are filtered out possibly results in the false positive and false negatives. This empirical experimental study explores ways of reducing false positives and negatives with labelled data.. In the process also discovered insufficient labelling in Openimage 2019 Object Detection dataset.


Scalable Reinforcement Learning for Multi Agent Networked Systems. We study reinforcement learning (RL) in a setting with a network of agents whose states and actions interact in a local manner where the objective is to find localized policies such that the (discounted) global reward is maximized. A fundamental challenge in this setting is that the state action space size scales exponentially in the number of agents, rendering the problem intractable for large networks. In this paper, we propose a Scalable Actor Critic (SAC) framework that exploits the network structure and finds a localized policy that is an O( rho kappa ) approximation of a stationary point of the objective for some rho in(0,1) , with complexity that scales with the local state action space size of the largest kappa hop neighborhood of the network. We illustrate our model and approach using examples from wireless communication, epidemics and traffic.


Singularity of pm 1 matrices and asymptotics of the number of threshold functions. Two results concerning the number of threshold functions P(2, n) and the probability mathbb P n that a random n times n Bernoulli matrix is singular are established. We introduce a supermodular function eta bigstar n 2 bf RP n fin to mathbb Z geq 0 , defined on finite subsets of bf RP n, that allows us to obtain a lower bound for P(2, n) in terms of mathbb P n 1 . This, together with L.Schl "afli's famous upper bound, give us asymptotics P(2, n) thicksim 2 2 n 1 choose n , quad n to infty. Also, the validity of the long standing conjecture concerning mathbb P n is proved mathbb P n thicksim (n 1) 22 1 n , quad n to infty .


Stationary Reflection and the failure of SCH. In this paper we prove that from large cardinals it is consistent that there is a singular strong limit cardinal nu such that the singular cardinal hypothesis fails at nu and every collection of fewer than mathrm cf ( nu) stationary subsets of nu reflects simultaneously. For uncountable cofinality, this situation was not previously known to be consistent. Using different methods, we reduce the upper bound on the consistency strength of this situation for mathrm cf ( nu) omega to below a single partially supercompact cardinal. The previous upper bound of infinitely many supercompact cardinals was due to Sharon.


A P k 2 polynomial lifting operator on polygons and polyhedrons. A P k 2 polynomial lifting operator is defined on polygons and polyhedrons. It lifts discontinuous polynomials inside the polygon polyhedron and on the faces to a one piece P k 2 polynomial. With this lifting operator, we prove that the weak Galerkin finite element solution, after this lifting, converges at two orders higher than the optimal order, in both L 2 and H 1 norms. The theory is confirmed by numerical solutions of 2D and 3D Poisson equations.


Deep Reinforcement Learning of Transition States. Combining reinforcement learning (RL) and molecular dynamics (MD) simulations, we propose a machine learning approach (RL ddag ) to automatically unravel chemical reaction mechanisms. In RL ddag , locating the transition state of a chemical reaction is formulated as a game, where a virtual player is trained to shoot simulation trajectories connecting the reactant and product. The player utilizes two functions, one for value estimation and the other for policy making, to iteratively improve the chance of winning this game. We can directly interpret the reaction mechanism according to the value function. Meanwhile, the policy function enables efficient sampling of the transition paths, which can be further used to analyze the reaction dynamics and kinetics. Through multiple experiments, we show that RL ddag can be trained tabula rasa hence allows us to reveal chemical reaction mechanisms with minimal subjective biases.


Learning and Evaluating Representations for Deep One class Classification. We present a two stage framework for deep one class classification. We first learn self supervised representations from one class data, and then build one class classifiers on learned representations. The framework not only allows to learn better representations, but also permits building one class classifiers that are faithful to the target task. We argue that classifiers inspired by the statistical perspective in generative or discriminative models are more effective than existing approaches, such as a normality score from a surrogate classifier. We thoroughly evaluate different self supervised representation learning algorithms under the proposed framework for one class classification. Moreover, we present a novel distribution augmented contrastive learning that extends training distributions via data augmentation to obstruct the uniformity of contrastive representations. In experiments, we demonstrate state of the art performance on visual domain one class classification benchmarks, including novelty and anomaly detection. Finally, we present visual explanations, confirming that the decision making process of deep one class classifiers is intuitive to humans. The code is available at https github.com google research deep representation one class.


Improved FPT Algorithms for Deletion to Forest like Structures. The Feedback Vertex Set problem is undoubtedly one of the most well studied problems in Parameterized Complexity. In this problem, given an undirected graph G and a non negative integer k , the objective is to test whether there exists a subset S subseteq V(G) of size at most k such that G S is a forest. After a long line of improvement, recently, Li and Nederlof SODA, 2020 designed a randomized algorithm for the problem running in time mathcal O star (2.7 k) . In the Parameterized Complexity literature, several problems around Feedback Vertex Set have been studied. Some of these include Independent Feedback Vertex Set (where the set S should be an independent set in G ), Almost Forest Deletion and Pseudoforest Deletion. In Pseudoforest Deletion, each connected component in G S has at most one cycle in it. However, in Almost Forest Deletion, the input is a graph G and non negative integers k, ell in mathbb N , and the objective is to test whether there exists a vertex subset S of size at most k , such that G S is ell edges away from a forest. In this paper, using the methodology of Li and Nederlof SODA, 2020 , we obtain the current fastest algorithms for all these problems. In particular we obtain following randomized algorithms. 1) Independent Feedback Vertex Set can be solved in time mathcal O star (2.7 k) . 2) Pseudo Forest Deletion can be solved in time mathcal O star (2.85 k) . 3) Almost Forest Deletion can be solved in mathcal O star ( min 2.85 k cdot 8.54 ell,2.7 k cdot 36.61 ell,3 k cdot 1.78 ell ) .


A Hierarchical Reasoning Graph Neural Network for The Automatic Scoring of Answer Transcriptions in Video Job Interviews. We address the task of automatically scoring the competency of candidates based on textual features, from the automatic speech recognition (ASR) transcriptions in the asynchronous video job interview (AVI). The key challenge is how to construct the dependency relation between questions and answers, and conduct the semantic level interaction for each question answer (QA) pair. However, most of the recent studies in AVI focus on how to represent questions and answers better, but ignore the dependency information and interaction between them, which is critical for QA evaluation. In this work, we propose a Hierarchical Reasoning Graph Neural Network (HRGNN) for the automatic assessment of question answer pairs. Specifically, we construct a sentence level relational graph neural network to capture the dependency information of sentences in or between the question and the answer. Based on these graphs, we employ a semantic level reasoning graph attention network to model the interaction states of the current QA session. Finally, we propose a gated recurrent unit encoder to represent the temporal question answer pairs for the final prediction. Empirical results conducted on CHNAT (a real world dataset) validate that our proposed model significantly outperforms text matching based benchmark models. Ablation studies and experimental results with 10 random seeds also show the effectiveness and stability of our models.


Aspect Based Sentiment Analysis in Education Domain. Analysis of a large amount of data has always brought value to institutions and organizations. Lately, people's opinions expressed through text have become a very important aspect of this analysis. In response to this challenge, a natural language processing technique known as Aspect Based Sentiment Analysis (ABSA) has emerged. Having the ability to extract the polarity for each aspect of opinions separately, ABSA has found itself useful in a wide range of domains. Education is one of the domains in which ABSA can be successfully utilized. Being able to understand and find out what students like and don't like most about a course, professor, or teaching methodology can be of great importance for the respective institutions. While this task represents a unique NLP challenge, many studies have proposed different approaches to tackle the problem. In this work, we present a comprehensive review of the existing work in ABSA with a focus in the education domain. A wide range of methodologies are discussed and conclusions are drawn.


Stability of the linear complementarity problem properties under interval uncertainty. We consider the linear complementarity problem with uncertain data modeled by intervals, representing the range of possible values. Many properties of the linear complementarity problem (such as solvability, uniqueness, convexity, finite number of solutions etc.) are reflected by the properties of the constraint matrix. In order that the problem has desired properties even in the uncertain environment, we have to be able to check them for all possible realizations of interval data. This leads us to the robust properties of interval matrices. In particular, we will discuss S matrix, Z matrix, copositivity, semimonotonicity, column sufficiency, principal nondegeneracy, R 0 matrix and R matrix. We characterize the robust properties and also suggest efficiently recognizable subclasses.


Comonadic semantics for guarded fragments. In previous work, Abramsky, Dawar and Wang (LiCS 2017) and Abramsky and Shah (CSL 2018) have shown how a range of model comparison games which play a central role in finite model theory, including Ehrenfeucht Fraisse, pebbling, and bisimulation games, can be captured in terms of resource indexed comonads on the category of relational structures. Moreover, the coalgebras for these comonads capture important combinatorial parameters such as tree width and tree depth. The present paper extends this analysis to quantifier guarded fragments of first order logic. We give a systematic account, covering atomic, loose and clique guards. In each case, we show that coKleisli morphisms capture winning strategies for Duplicator in the existential guarded bisimulation game, while back and forth bisimulation, and hence equivalence in the full guarded fragment, is captured by spans of open morphisms. We study the coalgebras for these comonads, and show that they correspond to guarded tree decompositions. We relate these constructions to a syntax free setting, with a comonad on the category of hypergraphs.


Characterization of passivity in Mueller matrices. Except for very particular and artificial experimental configurations, linear transformations of the state of polarization of an electromagnetic wave result in a reduction of the intensity of the exiting wave with respect to the incoming one. This natural passive behavior imposes certain mathematical restrictions on the corresponding Mueller matrices associated to the said transformations. Although the general conditions for passivity in Mueller matrices were presented in a previous paper J. J. Gil, J. Opt. Soc. Am. A 17, 328 334 (2000) , the demonstration was incomplete. In this paper, the set of two necessary and sufficient conditions for a Mueller matrix to represent a passive medium are determined and demonstrated on the basis of its arbitrary decomposition as a convex combination of nondepolarizing and passive pure Mueller matrices. The procedure followed to solve the problem provides also an appropriate framework to identify the Mueller matrix that, among the family of proportional passive Mueller matrices, exhibits the maximal physically achievable intensity transmittance. Beyond the theoretical interest on the rigorous characterization of passivity, the results obtained, when applied to absolute Mueller polarimetry, also provide a criterion to discard those experimentally measured Mueller matrices that do not satisfy the passivity criterion.


Matrix Deviation Inequality for ell p Norm, 1 leq p leq 2 . We prove the matrix deviation inequality for ell p norm, 1 leq p leq 2 , and i.i.d ensemble sub gaussian random matrix, which is a random matrix with i.i.d mean zero, unit variance, sub gaussian entries. As an application of this, we extend the Johnson Lindenstrauss lemma from ell 2 norm to ell p norm.


Flow dynamics in the closure region of an internal ship air cavity. This work is dedicated to providing a detailed account of the flow dynamics in the closure region of an internal ship air cavity. A geometrically simple multiwave test cavity is considered, and a simulation of the flow is conducted using large eddy simulation coupled with an algebraic Volume of Fluid interface capturing method. Results reveal that the flow in the closure region is highly unsteady and turbulent. The main cause of this is established to be the pressure gradient occurring due to the stagnation of the flow on the beach wall of the cavity. The pressure gradient leads to a steep incline in the mean location of the air water interface, which, in turn, leads to the flow separating from it and forming a recirculation zone, in which air and water are mixed. The separated flow becomes turbulent, which further facilitates the mixing and entrainment of air. Swarms of air bubbles leak periodically. Upstream of the closure region, the phase and length of the wave are found to be well predicted using existing approximations based on linear flow theory. However, for the corresponding prediction of the amplitude of the wave the agreement is worse, with the estimates under predicting the simulation results.


Anubhuti An annotated dataset for emotional analysis of Bengali short stories. Thousands of short stories and articles are being written in many different languages all around the world today. Bengali, or Bangla, is the second highest spoken language in India after Hindi and is the national language of the country of Bangladesh. This work reports in detail the creation of Anubhuti the first and largest text corpus for analyzing emotions expressed by writers of Bengali short stories. We explain the data collection methods, the manual annotation process and the resulting high inter annotator agreement of the dataset due to the linguistic expertise of the annotators and the clear methodology of labelling followed. We also address some of the challenges faced in the collection of raw data and annotation process of a low resource language like Bengali. We have verified the performance of our dataset with baseline Machine Learning as well as a Deep Learning model for emotion classification and have found that these standard models have a high accuracy and relevant feature selection on Anubhuti. In addition, we also explain how this dataset can be of interest to linguists and data analysts to study the flow of emotions as expressed by writers of Bengali literature.


Harmonic Forms on the Kodaira Thurston Manifold. We introduce an effective method to solve the bar partial harmonic forms on the Kodaira Thurston manifold endowed with an almost complex structure and an Hermitian metric. Using the Weil Brezin transform, we reduce the elliptic PDE system to countably many linear ODE systems. By solving a fundamental problem on linear ODE systems, the problem of finding bar partial harmonic forms is equivalent to a generalised Gauss circle problem. We demonstrate two remarkable applications. First, the dimension of the almost complex bar partial Hodge numbers on the Kodaira Thurston manifold could be arbitrarily large. Second, Hodge numbers vary with different choices of Hermitian metrics. This answers a question of Kodaira and Spencer in Hirzebruch's 1954 problem list.


Internal Wave Attractors in 3D Geometries a dynamical systems approach. We study the propagation in three dimensions of internal waves using ray tracing methods and traditional dynamical systems theory. The wave propagation on a cone that generalizes the Saint Andrew's cross justifies the introduction of an angle of propagation that allows to describe the position of the wave ray in the horizontal plane. Considering the evolution of this reflection angle for waves that repeatedly reflect off an inclined slope, a new trapping mechanism emerges that displays the tendency to align this angle with the upslope gradient. In the rather simple geometry of a translationally invariant canal, we show first that this configuration leads to trapezium shaped attractors, very similar to what has been extensively studied in two dimensions. However, we also establish a direct link between the trapping and the existence of two dimensional attractors. In a second stage, considering a geometry that is not translationally invariant, closer to realistic configurations, we prove that although there are no two dimensional attractors, one can find a structure in three dimensional space with properties similar to internal wave attractors a one dimensional attracting manifold. Moreover, as this structure is unique, it should be easy to visualize in laboratory experiments since energy injected in the domain would eventually be confined to a very thin region in three dimensional space, for which reason it is called a super attractor.


Towards Good Practices for Multi Person Pose Estimation. Multi Person Pose Estimation is an interesting yet challenging task in computer vision. In this paper, we conduct a series of refinements with the MSPN and PoseFix Networks, and empirically evaluate their impact on the final model performance through ablation studies. By taking all the refinements, we achieve 78.7 on the COCO test dev dataset and 76.3 on the COCO test challenge dataset.


Linear Schr "odinger equation with an almost periodic potential. We study the reducibility of a Linear Schr "odinger equation subject to a small unbounded almost periodic perturbation which is analytic in time and space. Under appropriate assumptions on the smallness, analiticity and on the frequency of the almost periodic perturbation, we prove that such an equation is reducible to constant coefficients via an analytic almost periodic change of variables. This implies control of both Sobolev and Analytic norms for the solution of the corresponding Schr "odinger equation for all times.


PodSumm Podcast Audio Summarization. The diverse nature, scale, and specificity of podcasts present a unique challenge to content discovery systems. Listeners often rely on text descriptions of episodes provided by the podcast creators to discover new content. Some factors like the presentation style of the narrator and production quality are significant indicators of subjective user preference but are difficult to quantify and not reflected in the text descriptions provided by the podcast creators. We propose the automated creation of podcast audio summaries to aid in content discovery and help listeners to quickly preview podcast content before investing time in listening to an entire episode. In this paper, we present a method to automatically construct a podcast summary via guidance from the text domain. Our method performs two key steps, namely, audio to text transcription and text summary generation. Motivated by a lack of datasets for this task, we curate an internal dataset, find an effective scheme for data augmentation, and design a protocol to gather summaries from annotators. We fine tune a PreSumm 10 model with our augmented dataset and perform an ablation study. Our method achieves ROUGE F(1 2 L) scores of 0.63 0.53 0.63 on our dataset. We hope these results may inspire future research in this direction.


Regularity criteria via one directional derivative of the velocity in anisotropic Lebesgue spaces to the 3D Navier Stokes equations. In this paper, we consider the regularity criterion for 3D incompressible Navier Stokes equations in terms of one directional derivative of the velocity in anisotropic Lebesgue spaces. More precisely, it is proved that u becomes a regular solution if the partial 3u satisfies int T 0 frac left left left partial 3 u(t) right L p x 1 right L q x 2 right beta L r x 3 1 ln left( partial 3u left(t right) L 2 e right) dt infty, text where frac 2 beta frac 1 p frac 1 q frac 1 r 1 text and 2 p, q, r leq infty, 1 left( frac 1 p frac 1 q frac 1 r right) geq 0 .


Minmax Hierarchies, Minimal Fibrations and a PDE based Proof of the Willmore Conjecture. We introduce a general scheme that permits to generate successive min max problems for producing critical points of higher and higher indices to Palais Smale Functionals in normal Banach manifolds equipped with complete Finsler structures. We call the resulting tree of minmax problems a minmax hierarchy. We give several examples and in particular we explain how to implement this scheme in the framework of the viscosity method introduced by the author some years ago in order to give a new proof of the Willmore conjecture after the famous result by Marques and Neves.


Sums of Cubes in Quaternion Rings. We investigate a version of Waring's Problem over quaternion rings, focusing on cubes in quaternion rings with integer coefficients. We determine the global upper and lower bounds for the number of cubes necessary to represent all such quaternions.


Solution to the Gliding Tone Problem. The solution is given to the classical problem of an oscillator driven by a sinusoid of steadily varying frequency. A closed analytical expression is obtained in the case where the Q factor of the oscillator is high, equivalent to the rotating wave approximation of atomic physics. In this case all independent variables of the system combine into a single parameter. The results are compared with previous work series and other approximations, numerical calculations, graphical solutions and analogue simulations. Attention is paid to the distortion of the resonance specifically frequency shift and amplitude attenuation consequent upon the finite sweep rate. The frequency shift is interpreted as a delay in the appearance of the resonance peak to leading order this time delay is twice the oscillator's ring down time. Measurements on a high Q oscillator are consistent with the mathematical solution.


A short proof that the L p diameter of Diff 0(S,area) is infinite. We give a short proof that the L p diameter of the group of area preserving diffeomorphisms isotopic to the identity of a compact surface is infinite.


Modeling stripe formation on growing zebrafish tailfins. As zebrafish develop, black and gold stripes form across their skin due to the interactions of brightly colored pigment cells. These characteristic patterns emerge on the growing fish body, as well as on the anal and caudal fins. While wild type stripes form parallel to a horizontal marker on the body, patterns on the tailfin gradually extend distally outward. Interestingly, several mutations lead to altered body patterns without affecting fin stripes. Through an exploratory modeling approach, our goal is to help better understand these differences between body and fin patterns. By adapting a prior agent based model of cell interactions on the fish body, we present an in silico study of stripe development on tailfins. Our main result is a demonstration that two cell types can produce stripes on the caudal fin. We highlight several ways that bone rays, growth, and the body fin interface may be involved in patterning, and we raise questions for future work related to pattern robustness.