Polyhedral value iteration for discounted games and energy games. We present a deterministic algorithm, solving discounted games with n nodes in n O(1) cdot (2 sqrt 2 ) n time. For bipartite discounted games our algorithm runs in n O(1) cdot 2 n time. Prior to our work no deterministic algorithm running in time 2 o(n log n) regardless of the discount factor was known. We call our approach polyhedral value iteration. We rely on a well known fact that the values of a discounted game can be found from the so called optimality equations. In the algorithm we consider a polyhedron obtained by relaxing optimality equations. We iterate points on the border of this polyhedron by moving each time along a carefully chosen shift as far as possible. This continues until the current point satisfies optimality equations. Our approach is heavily inspired by a recent algorithm of Dorfman et al. (ICALP 2019) for energy games. For completeness, we present their algorithm in terms of polyhedral value iteration. Our exposition, unlike the original algorithm, does not require edge weights to be integers and works for arbitrary real weights.


Visual search over billions of aerial and satellite images. We present a system for performing visual search over billions of aerial and satellite images. The purpose of visual search is to find images that are visually similar to a query image. We define visual similarity using 512 abstract visual features generated by a convolutional neural network that has been trained on aerial and satellite imagery. The features are converted to binary values to reduce data and compute requirements. We employ a hash based search using Bigtable, a scalable database service from Google Cloud. Searching the continental United States at 1 meter pixel resolution, corresponding to approximately 2 billion images, takes approximately 0.1 seconds. This system enables real time visual search over the surface of the earth, and an interactive demo is available at https search.descarteslabs.com.


Polynomial Invariants of Singular Knots and links. We generalize the notion of the quandle polynomial to the case of singquandles. We show that the singquandle polynomial is an invariant of finite singquandles. We also construct a singular link invariant from the singquandle polynomial and show that this new singular link invariant generalizes the singquandle counting invariant. In particular, using the new polynomial invariant, we can distinguish singular links with the same singquandle counting invariant.


Analyzing the World Wide Impact of Public Health Interventions on the Transmission Dynamics of COVID 19. We analyze changes in the reproduction number, R, of COVID 19 in response to public health interventions. Our results indicate that public health measures undertaken in China reduced R from 1.5 in January to 0.4 in mid March 2020. They also suggest, however, the limitations of isolation, quarantine, and large scale attempts to limit travel. While the world wide reproduction number briefly dropped below 1 as China implemented extensive public health measures, the introduction of the virus to other nations swiftly led to an increasing world wide average value of R. In Italy, the nation hardest hit following China, social distancing measures brought the local value of R down from 3.71 to 2.51. Nonetheless, the value of R in Italy persisted at levels well above 1, allowing for ongoing transmission. By mid March 2020, as COVID 19 spread in areas without extensive public health interventions in place, the world wide value of R increased to a level similar to that of late January.


Generative models of simultaneously heavy tailed distributions of inter event times on nodes and edges. Intervals between discrete events representing human activities, as well as other types of events, often obey heavy tailed distributions, and their impacts on collective dynamics on networks such as contagion processes have been intensively studied. The literature supports that such heavy tailed distributions are present for inter event times associated with both individual nodes and individual edges in networks. However, the simultaneous presence of heavy tailed distributions of inter event times for nodes and edges is a non trivial phenomenon, and its origin has been elusive. In the present study, we propose a generative model and its variants to explain this phenomenon. We assume that each node independently transits between a high activity and low activity state according to a continuous time two state Markov process and that, for the main model, events on an edge occur at a high rate if and only if both end nodes of the edge are in the high activity state. In other words, two nodes interact frequently only when both nodes prefer to interact with others. The model produces distributions of inter event times for both individual nodes and edges that resemble heavy tailed distributions across some scales. It also produces positive correlation in consecutive inter event times, which is another stylized observation for empirical data of human activity. We expect that our modeling framework provides a useful benchmark for investigating dynamics on temporal networks driven by non Poissonian event sequences.


Optimal Stopping of Stochastic Transport Minimizing Submartingale Costs. Given a stochastic state process (X t) t and a real valued submartingale cost process (S t) t , we characterize optimal stopping times tau that minimize the expectation of S tau while realizing given initial and target distributions mu and nu , i.e., X 0 sim mu and X tau sim nu . A dual optimization problem is considered and shown to be attained under suitable conditions. The optimal solution of the dual problem then provides a contact set, which characterizes the location where optimal stopping can occur. The optimal stopping time is uniquely determined as the first hitting time of this contact set provided we assume a natural structural assumption on the pair (X t, S t) t , which generalizes the twist condition on the cost in optimal transport theory. This paper extends the Brownian motion settings studied in 15, 16 and deals with more general costs.


Recurrences and Congruences for Higher order Geometric Polynomials and Related Numbers. We obtain new recurrence relations, an explicit formula, and convolution identities for higher order geometric polynomials. These relations generalize known results for geometric polynomials, and lead to congruences for higher order geometric polynomials, particularly for p Bernoulli numbers.


Inviscid limit of the compressible Navier Stokes equations for asymptotically isothermal gases. We prove the existence of relative finite energy vanishing viscosity solutions of the one dimensional, isentropic Euler equations under the assumption of an asymptotically isothermal pressure law, that is, p( rho) rho O(1) in the limit rho to infty . This solution is obtained as the vanishing viscosity limit of classical solutions of the one dimensional, isentropic, compressible Navier Stokes equations. Our approach relies on the method of compensated compactness to pass to the limit rigorously in the nonlinear terms. Key to our strategy is the derivation of hyperbolic representation formulas for the entropy kernel and related quantities among others, a special entropy pair used to obtain higher uniform integrability estimates on the approximate solutions. Intricate bounding procedures relying on these representation formulas then yield the required compactness of the entropy dissipation measures. In turn, we prove that the Young measure generated by the classical solutions of the Navier Stokes equations reduces to a Dirac mass, from which we deduce the required convergence to a solution of the Euler equations.


Sequence space representations for spaces of entire functions with rapid decay on strips. We obtain sequence space representations for a class of Fr 'echet spaces of entire functions with rapid decay on horizontal strips. In particular, we show that the projective Gelfand Shilov spaces Sigma 1 nu and Sigma nu 1 are isomorphic to Lambda infty (n 1 ( nu 1) ) for nu 0 .


Rigidity of determinantal point processes on the unit disc with sub Bergman kernels. We give natural constructions of number rigid determinantal point processes on the unit disc mathbb D with sub Bergman kernels of the form K Lambda(z, w) sum n in Lambda (n 1) z n bar w n, quad z, w in mathbb D , with Lambda an infinite subset of the set of non negative integers. Our constructions are given both in a deterministic method and a probabilisitc method. In the deterministic method, our proofs involve the classical Bloch functions.


Analogues of Alladi's formula over global function fields. In this paper, we show an analogue of Kural, McDonald and Sah's result on Alladi's formula for global function fields. Explicitly, we show that for a global function field K , if a set S of prime divisors has a natural density delta(S) within prime divisors, then lim n to infty sum substack 1 le deg D le n D in mathfrak D (K,S) frac mu(D) D delta(S), where mu(D) is the M " o bius function on divisors and mathfrak D (K,S) is the set of all effective distinguishable divisors whose smallest prime factors are in S . As applications, we get the analogue of Dawsey's and Sweeting and Woo's results to the Chebotarev Density Theorem for function fields, and the analogue of Alladi's result to the Prime Polynomial Theorem for arithmetic progressions. We also display a connection between the M "obius function and the Fourier coefficients of modular form associated to elliptic curves. The proof of our main theorem is similar to the approach in Kural et al.'s article.


Continuous maximal covering location problems with interconnected facilities. In this paper we analyze a continuous version of the maximal covering location problem, in which the facilities are required to be interconnected by means of a graph structure in which two facilities are allowed to be linked if a given distance is not exceed. We provide a mathematical programming framework for the problem and different resolution strategies. First, we propose a Mixed Integer Non Linear Programming formulation, and derive properties of the problem that allow us to project the continuous variables out avoiding the nonlinear constraints, resulting in an equivalent pure integer programming formulation. Since the number of constraints in the integer programming formulation is large and the constraints are, in general, difficult to handle, we propose two branch cut approaches that avoid the complete enumeration of the constraints resulting in more efficient procedures. We report the results of an extensive battery of computational experiments comparing the performance of the different approaches.


Finitely generated groups acting uniformly properly on hyperbolic space. We construct an uncountable sequence of groups acting uniformly properly on hyperbolic spaces. We show that only countably many of these groups can be virtually torsion free. This gives new examples of groups acting uniformly properly on hyperbolic spaces that are not virtually torsion free and cannot be subgroups of hyperbolic groups.


Finding and counting permutations via CSPs. Permutation patterns and pattern avoidance have been intensively studied in combinatorics and computer science, going back at least to the seminal work of Knuth on stack sorting (1968). Perhaps the most natural algorithmic question in this area is deciding whether a given permutation of length n contains a given pattern of length k . In this work we give two new algorithms for this well studied problem, one whose running time is n k 4 o(k) , and a polynomial space algorithm whose running time is the better of O(1.6181 n) and O(n k 2 1 ) . These results improve the earlier best bounds of n 0.47k o(k) and O(1.79 n) due to Ahal and Rabinovich (2000) resp. Bruner and Lackner (2012) and are the fastest algorithms for the problem when k in Omega( log n ) . We show that both our new algorithms and the previous exponential time algorithms in the literature can be viewed through the unifying lens of constraint satisfaction. Our algorithms can also count, within the same running time, the number of occurrences of a pattern. We show that this result is close to optimal solving the counting problem in time f(k) cdot n o(k log k ) would contradict the exponential time hypothesis (ETH). For some special classes of patterns we obtain improved running times. We further prove that 3 increasing and 3 decreasing permutations can, in some sense, embed arbitrary permutations of almost linear length, which indicates that an algorithm with sub exponential running time is unlikely, even for patterns from these restricted classes.


Cyclic Functional Mapping Self supervised correspondence between non isometric deformable shapes. We present the first utterly self supervised network for dense correspondence mapping between non isometric shapes. The task of alignment in non Euclidean domains is one of the most fundamental and crucial problems in computer vision. As 3D scanners can generate highly complex and dense models, the mission of finding dense mappings between those models is vital. The novelty of our solution is based on a cyclic mapping between metric spaces, where the distance between a pair of points should remain invariant after the full cycle. As the same learnable rules that generate the point wise descriptors apply in both directions, the network learns invariant structures without any labels while coping with non isometric deformations. We show here state of the art results by a large margin for a variety of tasks compared to known self supervised and supervised methods.


Analysis on the Pricing model for a Discrete Coupon Bond with Early redemption provision by the Structural Approach. In this paper, using the structural approach is derived a mathematical model of the discrete coupon bond with the provision that allow the holder to demand early redemption at any coupon dates prior to the maturity and based on this model is provided some analysis including min max and gradient estimates of the bond price. Using these estimates the existence and uniqueness of the default boundaries and some relationships between the design parameters of the discrete coupon bond with early redemption provision are described. Then under some assumptions the existence and uniqueness of the early redemption boundaries is proved and the analytic formula of the bond price is provided using higher binary options. Finally for our bond is provided the analysis on the duration and credit spread, which are used widely in financial reality. Our works provide a design guide of the discrete coupon bond with the early redemption provision


Mean value of real characters using a double Dirichlet series. We study the double character sum sum limits substack m leq X, m odd sum limits substack n leq Y, n odd leg mn and its smoothly weighted counterpart. An asymptotic formula with power saving error term was obtained by Conrey, Farmer and Soundararajan by applying the Poisson summation formula. The result is interesting, because the main term involves a non smooth function. In this paper, we apply the inverse Mellin transform twice and study the resulting double integral that involves a double Dirichlet series. This method has two advantages it leads to a better error term, and the surprising main term naturally arises from three residues of the double Dirichlet series.


Dynamic Inference A New Approach Toward Efficient Video Action Recognition. Though action recognition in videos has achieved great success recently, it remains a challenging task due to the massive computational cost. Designing lightweight networks is a possible solution, but it may degrade the recognition performance. In this paper, we innovatively propose a general dynamic inference idea to improve inference efficiency by leveraging the variation in the distinguishability of different videos. The dynamic inference approach can be achieved from aspects of the network depth and the number of input video frames, or even in a joint input wise and network depth wise manner. In a nutshell, we treat input frames and network depth of the computational graph as a 2 dimensional grid, and several checkpoints are placed on this grid in advance with a prediction module. The inference is carried out progressively on the grid by following some predefined route, whenever the inference process comes across a checkpoint, an early prediction can be made depending on whether the early stop criteria meets. For the proof of concept purpose, we instantiate three dynamic inference frameworks using two well known backbone CNNs. In these instances, we overcome the drawback of limited temporal coverage resulted from an early prediction by a novel frame permutation scheme, and alleviate the conflict between progressive computation and video temporal relation modeling by introducing an online temporal shift module. Extensive experiments are conducted to thoroughly analyze the effectiveness of our ideas and to inspire future research efforts. Results on various datasets also evident the superiority of our approach.


Persistence of hubs in growing random networks. We consider models of evolving networks mathcal G n n geq 0 modulated by two parameters an attachment function f mathbb N 0 to mathbb R and a (possibly random) attachment sequence m i i geq 1 . Starting with a single vertex, at each discrete step i geq 1 a new vertex v i enters the system with m i geq 1 edges which it sequentially connects to a pre existing vertex v in mathcal G i 1 with probability proportional to f( operatorname degree (v)) . We consider the problem of emergence of persistent hubs existence of a finite (a.s.) time n such that for all n geq n the identity of the maximal degree vertex (or in general the K largest degree vertices for K geq 1 ) does not change. We obtain general conditions on f and m i i geq 1 under which a persistent hub emerges, and also those under which a persistent hub fails to emerge. In the case of lack of persistence, for the specific case of trees ( m i equiv 1 for all i ), we derive asymptotics for the maximal degree and the index of the maximal degree vertex (time at which the vertex with current maximal degree entered the system) to understand the movement of the maximal degree vertex as the network evolves. A key role in the analysis is played by an inverse rate weighted martingale constructed from a continuous time embedding of the discrete time model. Asymptotics for this martingale, including concentration inequalities and moderate deviations, play a major role in the analysis of the model.


Continuous Sentences Preserved Under Reduced Products. Answering a question of Cif 'u Lopes, we give a syntactic characterization of those continuous sentences that are preserved under reduced products of metric structures. In fact, we settle this question in the wider context of general structures as introduced by the second author.


The mathbb Z p mathbb Z equivariant product isomorphism in fixed point Floer cohomology. Let p geq 2 be a prime, and mathbb F p be the field with p elements. Extending a result of Seidel for p 2, we construct an isomorphism between the Floer cohomology of an exact or Hamiltonian symplectomorphism phi, with mathbb F p coefficients, and the mathbb Z p mathbb Z equivariant Tate Floer cohomology of its p th power phi p. The construction involves a Kaledin type quasi Frobenius map, as well as a mathbb Z p mathbb Z equivariant pants product an equivariant operation with p inputs and 1 output. Our method of proof involves a spectral sequence for the action filtration, and a local mathbb Z p mathbb Z equivariant coproduct providing an inverse on the E 2 page. This strategy has the advantage of accurately describing the effect of the isomorphism on filtration levels. We describe applications to the symplectic mapping class group, as well as develop Smith theory for the persistence module of a Hamiltonian diffeomorphism phi on symplectically aspherical symplectic manifolds. We illustrate the latter by giving a new proof of the celebrated no torsion theorem of Polterovich, and by relating the growth rate of the number of periodic points of the p k th iteration of phi and its distance to the identity. Along the way, we prove a sharpening of the classical Smith inequality for actions of mathbb Z p mathbb Z .


A unified framework to generate optimized compact finite difference schemes. A unified framework to derive optimized compact schemes for a uniform grid is presented. The optimal scheme coefficients are determined analytically by solving an optimization problem to minimize the spectral error subject to equality constraints that ensure specified order of accuracy. A rigorous stability analysis for the optimized schemes is also presented. We analytically prove the relation between order of a derivative and symmetry or skew symmetry of the optimal coefficients approximating it. We also show that other types of schemes e.g., spatially explicit, and biased finite differences, can be generated as special cases of the framework.


Max plus convexity in Riesz spaces. We study max plus convexity in an Archimedean Riesz space E with an order unit un the definition of max plus convex sets is algebraic and we do not assume that E has an it a priori given topological structure. To the given unit un one can associate two equivalent norms norm cdot norm un and norm cdot norm hun on E the distance sf D hun on E associated to norm cdot norm hun is a geodesic distance for which max plus convex sets in E are geodesically closed sets. Under suitable assumptions, we establish max plus versions of some fixed points and continuous selection theorems that are well known for linear convex sets and we show that hyperspaces of compact max plus convex sets are Absolute Retracts.


Fast and Accurate Pseudoinverse with Sparse Matrix Reordering and Incremental Approach. How can we compute the pseudoinverse of a sparse feature matrix efficiently and accurately for solving optimization problems? A pseudoinverse is a generalization of a matrix inverse, which has been extensively utilized as a fundamental building block for solving linear systems in machine learning. However, an approximate computation, let alone an exact computation, of pseudoinverse is very time consuming due to its demanding time complexity, which limits it from being applied to large data. In this paper, we propose FastPI (Fast PseudoInverse), a novel incremental singular value decomposition (SVD) based pseudoinverse method for sparse matrices. Based on the observation that many real world feature matrices are sparse and highly skewed, FastPI reorders and divides the feature matrix and incrementally computes low rank SVD from the divided components. To show the efficacy of proposed FastPI, we apply them in real world multi label linear regression problems. Through extensive experiments, we demonstrate that FastPI computes the pseudoinverse faster than other approximate methods without loss of accuracy. and uses much less memory compared to full rank SVD based approach. Results imply that our method efficiently computes the low rank pseudoinverse of a large and sparse matrix that other existing methods cannot handle with limited time and space.


Dialectal Layers in West Iranian a Hierarchical Dirichlet Process Approach to Linguistic Relationships. This paper addresses a series of complex and unresolved issues in the historical phonology of West Iranian languages. The West Iranian languages (Persian, Kurdish, Balochi, and other languages) display a high degree of non Lautgesetzlich behavior. Most of this irregularity is undoubtedly due to language contact we argue, however, that an oversimplified view of the processes at work has prevailed in the literature on West Iranian dialectology, with specialists assuming that deviations from an expected outcome in a given non Persian language are due to lexical borrowing from some chronological stage of Persian. It is demonstrated that this qualitative approach yields at times problematic conclusions stemming from the lack of explicit probabilistic inferences regarding the distribution of the data Persian may not be the sole donor language additionally, borrowing at the lexical level is not always the mechanism that introduces irregularity. In many cases, the possibility that West Iranian languages show different reflexes in different conditioning environments remains under explored. We employ a novel Bayesian approach designed to overcome these problems and tease apart the different determinants of irregularity in patterns of West Iranian sound change. Our methodology allows us to provisionally resolve a number of outstanding questions in the literature on West Iranian dialectology concerning the dialectal affiliation of certain sound changes. We outline future directions for work of this sort.


Reconstructing the thermal phonon transmission coefficient at solid interfaces in the phonon transport equation. The ab initio model for heat propagation is the phonon transport equation, a Boltzmann like kinetic equation. When two materials are put side by side, the heat that propagates from one material to the other experiences thermal boundary resistance. Mathematically, it is represented by the reflection coefficient of the phonon transport equation on the interface of the two materials. This coefficient takes different values at different phonon frequencies, between different materials. In experiments scientists measure the surface temperature of one material to infer the reflection coefficient as a function of phonon frequency. In this article, we formulate this inverse problem in an optimization framework and apply the stochastic gradient descent (SGD) method for finding the optimal solution. We furthermore prove the maximum principle and show the Lipschitz continuity of the Fr 'echet derivative. These properties allow us to justify the application of SGD in this setup.


Exponential mixing property for automorphisms of compact K "ahler manifolds. Let f be a holomorphic automorphism of a compact K "ahler manifold. Assume moreover that f admits a unique maximal dynamic degree d p with only one eigenvalue of maximal modulus. Let mu be its equilibrium measure. In this paper, we prove that mu is exponentially mixing for all d.s.h. test functions.


A Krylov Schur like method for computing the best rank (r 1,r 2,r 3) approximation of large and sparse tensors. The paper is concerned with methods for computing the best low multilinear rank approximation of large and sparse tensors. Krylov type methods have been used for this problem here block versions are introduced. For the computation of partial eigenvalue and singular value decompositions of matrices the Krylov Schur (restarted Arnoldi) method is used. We describe a generalization of this method to tensors, for computing the best low multilinear rank approximation of large and sparse tensors. In analogy to the matrix case, the large tensor is only accessed in multiplications between the tensor and blocks of vectors, thus avoiding excessive memory usage. It is proved that, if the starting approximation is good enough, then the tensor Krylov Schur method is convergent. Numerical examples are given for synthetic tensors and sparse tensors from applications, which demonstrate that for most large problems the Krylov Schur method converges faster and more robustly than higher order orthogonal iteration.


Robustness of football passing networks against continuous node and link removals. We can construct passing networks when we regard a player as a node and a pass as a link in football games. Thus, we can analyze the networks by using tools developed in network science. Among various metrics characterizing a network, centrality metrics have often been used to identify key players in a passing network. However, a tolerance to players being marked or passes being blocked in a passing network, namely the robustness of the network, has been poorly understood so far. Because the robustness of a passing network can be connected to the increase of ball possession, it would be deeply related to the outcome of a game. Here, we developed position dependent passing networks of 45 matches by 18 teams belonging to the Japan Professional Football League. Then, nodes or links were continuously removed from the passing networks by two removal methods so that we could evaluate the robustness of these networks against the removals. The results show that these passing networks commonly contain hubs (key players making passes). Then, we analyzed the most robust networks in detail and found that their full backs increase the robustness by often invoking a heavier emphasis on attack. Moreover, we showed that the robustness of the passing networks and the team performance have a positive correlation.


B 'ezoutians and injectivity of polynomial maps. We prove that an endomorphism f of affine space is injective on rational points if its B 'ezoutian is constant. Similarly, f is injective at a given rational point if its reduced B 'ezoutian is constant. We also show that if the Jacobian determinant of f is invertible, then f is injective at a given rational point if and only if its reduced B 'ezoutian is constant.


A vanishing identity on adjoint Reidemeister torsions of twist knots. For a compact oriented 3 manifold with torus boundary the adjoint Reidemeister torsion is defined as a function on the mathrm SL 2( mathbb C ) character variety depending on a choice of a boundary curve. Under reasonable assumptions, it is conjectured that the adjoint torsion satisfies a certain type of vanishing identities. In this paper, we prove that the conjecture holds for all hyperbolic twist knot exteriors by using Jacobi's residue theorem.


Design and Implementation of a Wireless SensorNetwork for Agricultural Applications. We present the design and implementation of a shortest path tree based, energy efficient data collection wireless sensor network to sense various parameters in an agricultural farm using in house developed low cost sensors. Nodes follow a synchronized, periodic sleep wake up schedule to maximize the lifetime of the network. The implemented network consists of 24 sensor nodes in a 3 acre maize farm and its performance is captured by 7 snooper nodes for different data collection intervals 10 minutes, 1 hour and 3 hours. The almost static nature of wireless links in the farm motivated us to use the same tree for a long data collection period(3 days). The imbalance in energy consumption across nodes is observed to be very small and the network architecture uses easy to implement protocols to perform different network activities including handling of node failures. We present the results and analysis of extensive tests conducted on our implementation, which provide significant insights.


Equivariant K theory and tangent spaces to Schubert varieties. Tangent spaces to Schubert varieties of type A were characterized by Lakshmibai and Seshadri. This result was extended to the other classical types by Lakshmibai. We give a uniform characterization of tangent spaces to Schubert varieties in cominuscule G P. Our results extend beyond cominuscule G P they describe the tangent space to any Schubert variety in G B at a point xB, where x is a cominuscule Weyl group element in the sense of Peterson. Our results also give partial information about the tangent space to any Schubert variety at any point. Our method is to describe the tangent spaces of Kazhdan Lusztig varieties, and then recover results for Schubert varieties. Our proof uses a relationship between weights of the tangent space of a variety with torus action, and factors of the class of the variety in torus equivariant K theory. The proof relies on a formula for Schubert classes in equivariant K theory due to Graham and Willems, as well as a theorem on subword complexes due to Knutson and Miller.


On Generalized Regularity. Szemeredi's regularity lemma is one instance in a family of regularity lemmas, replacing the definition of density of a graph by a more general coefficient. Recently, Fan Chung proved another instance, a regularity lemma for clustering graphs, and asked whether good upper bounds could be derived for the quantitative estimates it supplies. We answer this question in the negative, for every generalized regularity lemma.


Efficient Video Semantic Segmentation with Labels Propagation and Refinement. This paper tackles the problem of real time semantic segmentation of high definition videos using a hybrid GPU CPU approach. We propose an Efficient Video Segmentation(EVS) pipeline that combines (i) On the CPU, a very fast optical flow method, that is used to exploit the temporal aspect of the video and propagate semantic information from one frame to the next. It runs in parallel with the GPU. (ii) On the GPU, two Convolutional Neural Networks A main segmentation network that is used to predict dense semantic labels from scratch, and a Refiner that is designed to improve predictions from previous frames with the help of a fast Inconsistencies Attention Module (IAM). The latter can identify regions that cannot be propagated accurately. We suggest several operating points depending on the desired frame rate and accuracy. Our pipeline achieves accuracy levels competitive to the existing real time methods for semantic image segmentation(mIoU above 60 ), while achieving much higher frame rates. On the popular Cityscapes dataset with high resolution frames (2048 x 1024), the proposed operating points range from 80 to 1000 Hz on a single GPU and CPU.


Improving Investment Suggestions for Peer to Peer (P2P) Lending via Integrating Credit Scoring into Profit Scoring. In the peer to peer (P2P) lending market, lenders lend the money to the borrowers through a virtual platform and earn the possible profit generated by the interest rate. From the perspective of lenders, they want to maximize the profit while minimizing the risk. Therefore, many studies have used machine learning algorithms to help the lenders identify the "best" loans for making investments. The studies have mainly focused on two categories to guide the lenders' investments one aims at minimizing the risk of investment (i.e., the credit scoring perspective) while the other aims at maximizing the profit (i.e., the profit scoring perspective). However, they have all focused on one category only and there is seldom research trying to integrate the two categories together. Motivated by this, we propose a two stage framework that incorporates the credit information into a profit scoring modeling. We conducted the empirical experiment on a real world P2P lending data from the US P2P market and used the Light Gradient Boosting Machine (lightGBM) algorithm in the two stage framework. Results show that the proposed two stage method could identify more profitable loans and thereby provide better investment guidance to the investors compared to the existing one stage profit scoring alone approach. Therefore, the proposed framework serves as an innovative perspective for making investment decisions in P2P lending.


Effect of transcription reinitiation in stochastic gene expression. Gene expression (GE) is an inherently random or stochastic or noisy process. The randomness in different steps of GE, e.g., transcription, translation, degradation, etc., leading to cell to cell variations in mRNA and protein levels. This variation appears in organisms ranging from microbes to metazoans. Stochastic gene expression has important consequences for cellular function. The random fluctuations in protein levels produce variability in cellular behavior. It is beneficial in some contexts and harmful to others. These situations include stress response, metabolism, development, cell cycle, circadian rhythms, and aging. Different model studies e.g., constitutive, two state, etc., reveal that the fluctuations in mRNA and protein levels arise from different steps of gene expression among which the steps in transcription have the maximum effect. The pulsatile mRNA production through RNAP II based reinitiation of transcription is an important part of gene transcription. Though, the effect of that process on mRNA and protein levels is very little known. The addition of any biochemical step in the constitutive or two state process generally decreases the mean and increases the Fano factor. In this study, we have shown that the RNAP II based reinitiation process in gene transcription can have different effects on both mean and Fano factor at mRNA levels in different model systems. It decreases the mean and Fano factor both at the mRNA levels in the constitutive network whereas in other networks it can simultaneously increase or decrease both quantities or it can have mixed effect at mRNA levels. We propose that a constitutive network with reinitiation behaves like a product independent negative feedback circuit whereas other networks behave as either product independent positive or negative or mixed feedback circuit.


Invariant subspaces and exact solutions for some types of scalar and coupled time space fractional diffusion equations. We explain how the invariant subspace method can be extended to a scalar and coupled system of time space fractional partial differential equations. The effectiveness and applicability of the method have been illustrated through time space (i) fractional diffusion convection equation, (ii) fractional reaction diffusion equation, (iii) fractional diffusion equation with source term, (iv) two coupled system of the fractional diffusion equation, (v) two coupled system of fractional stationary transonic plane parallel gas flow equation and (vi) three coupled system of fractional Hirota Satsuma KdV equation. Also, we explicitly presented how to derive more than one exact solution of the equations as mentioned above using the invariant subspace method.


Implied and Realized Volatility A Study of Distributions and the Distribution of Difference. We study distributions of realized variance (squared realized volatility) and squared implied volatility, as represented by VIX and VXO indices. We find that Generalized Beta distribution provide the best fits. These fits are much more accurate for realized variance than for squared VIX and VXO possibly another indicator that the latter have deficiencies in predicting the former. We also show that there are noticeable differences between the distributions of the 1970 2017 realized variance and its 1990 2017 portion, for which VIX and VXO became available. This may be indicative of a feedback effect that implied volatility has on realized volatility. We also discuss the distribution of the difference between squared implied volatility and realized variance and show that, at the basic level, it is consistent with Pearson's correlations obtained from linear regression.


Mean field limit of Age and Leaky Memory dependent Hawkes processes. We propose a mean field model of interacting point processes where each process has a memory of the time elapsed since its last event (age) and its recent past (leaky memory), generalizing Age dependent Hawkes processes. The model is motivated by interacting nonlinear Hawkes processes with Markovian self interaction and networks of spiking neurons with adaptation and short term synaptic plasticity. By proving propagation of chaos and using a path integral representation for the law of the limit process, we show that, in the mean field limit, the empirical measure of the system follows a multidimensional nonlocal transport equation.


Periodic Higgs bundles over curves. In this article, we study periodic Higgs bundles and their applications. We obtain the following results i). an elliptic curve has infinitely many primes of supersingular reduction if and only if any periodic Higgs bundle over it is a direct sum of torsion line bundles ii). the uniformizing de Rham bundle attached to a generic projective hyperbolic curve is not one periodic, and it is motivic iff it admits a modular embedding (e.g. Shimura curves, triangle curves) iii). there is an explicit Deuring Eichler mass formula for the Newton jumping locus a Shimura curve of Hodge type. We propose the periodic Higgs conjecture, which would imply an arithmetic Simpson correspondence. The conjecture holds in rank one case.


Towards a Metric for Automated Conversational Dialogue System Evaluation and Improvement. We present "AutoJudge", an automated evaluation method for conversational dialogue systems. The method works by first generating dialogues based on self talk, i.e. dialogue systems talking to itself. Then, it uses human ratings on these dialogues to train an automated judgement model. Our experiments show that AutoJudge correlates well with the human ratings and can be used to automatically evaluate dialogue systems, even in deployed systems. In a second part, we attempt to apply AutoJudge to improve existing systems. This works well for re ranking a set of candidate utterances. However, our experiments show that AutoJudge cannot be applied as reward for reinforcement learning, although the metric can distinguish good from bad dialogues. We discuss potential reasons, but state here already that this is still an open question for further research.


Efficient near field to far field transformations for quasinormal modes of optical cavities and plasmonic resonators. We describe an efficient near field to far field transformation for optical quasinormal modes, which are the dissipative modes of open cavities and plasmonic resonators with complex eigenfrequencies. As an application of the theory, we show how one can compute the reservoir modes (or regularized quasinormal modes) outside the resonator, which are essential to use in both classical and quantum optics. We subsequently demonstrate how to efficiently compute the quantum optical parameters necessary in the theory of quantized quasinormal modes Franke et al., Phys. Rev. Lett. 122, 213901 (2019) . To confirm the accuracy of our technique, we directly compare with a Dyson equation approach currently used in the literature (in regimes where this is possible), and demonstrate several order of magnitude improvement for the calculation run times. We also introduce an efficient pole approximation for computing the quantized quasinormal mode parameters, since they require an integration over a range of frequencies. Using this approach, we show how to compute regularized quasinormal modes and quantum optical parameters for a full 3D metal dimer in under one minute on a standard desktop computer. Our technique is exemplified by studying the quasinormal modes of metal dimers and a hybrid structure consisting of a gold dimer on top of a photonic crystal beam. In the latter example, we show how to compute the quantum optical parameters that describe a pronounced Fano resonance, using structural geometries that cannot practically be solved using a Dyson equation approach. All calculations for the spontaneous emission rates are confirmed with full dipole calculations in Maxwell's equations and are shown to be in excellent agreement.


The explicit formula for solution of anomalous diffusion equation in the multi dimensional space. This paper intends on obtaining the explicit solution of n dimensional anomalous diffusion equation in the infinite domain with non zero initial condition and vanishing condition at infinity. It is shown that this equation can be derived from the parabolic integro differential equation with memory in which the kernel is t alpha E 1 alpha, 1 alpha ( t 1 alpha ), alpha in(0, 1), where E alpha, beta is the Mittag Liffler function. Based on Laplace and Fourier transforms the properties of the Fox H function and convolution theorem, explicit solution for anomalous diffusion equation is obtained.


On the multiplier problem for the ball on graded Lie groups. In this note, we consider a non commutative analogy of the classical Fefferman multiplier problem for the ball. More precisely, if chi is the characteristic function of the unit interval I 0,1 , we investigate a family of differential operators mathcal R on a graded Lie group G, for which the multipliers chi( mathcal R ) are bounded on L p(G), if and only if p 2.


Taxonomical hierarchy of canonicalized relations from multiple Knowledge Bases. This work addresses two important questions pertinent to Relation Extraction (RE). First, what are all possible relations that could exist between any two given entity types? Second, how do we define an unambiguous taxonomical (is a) hierarchy among the identified relations? To address the first question, we use three resources Wikipedia Infobox, Wikidata, and DBpedia. This study focuses on relations between person, organization and location entity types. We exploit Wikidata and DBpedia in a data driven manner, and Wikipedia Infobox templates manually to generate lists of relations. Further, to address the second question, we canonicalize, filter, and combine the identified relations from the three resources to construct a taxonomical hierarchy. This hierarchy contains 623 canonical relations with highest contribution from Wikipedia Infobox followed by DBpedia and Wikidata. The generated relation list subsumes an average of 85 of relations from RE datasets when entity types are restricted.


Axiomatisability of the Class of Monolithic Groups in a Variety of Nilpotent Groups. The class of all subdirectly irreducible groups belonging to a variety generated by a finite nilpotent group can be axiomatised by a finite set of elementary sentences.


The Value of Information in Selfish Routing. Path selection by selfish agents has traditionally been studied by comparing social optima and equilibria in the Wardrop model, i.e., by investigating the Price of Anarchy in selfish routing. In this work, we refine and extend the traditional selfish routing model in order to answer questions that arise in emerging path aware Internet architectures. The model enables us to characterize the impact of different degrees of congestion information that users possess. Furthermore, it allows us to analytically quantify the impact of selfish routing, not only on users, but also on network operators. Based on our model, we show that the cost of selfish routing depends on the network topology, the perspective (users versus network operators), and the information that users have. Surprisingly, we show analytically and empirically that less information tends to lower the Price of Anarchy, almost to the optimum. Our results hence suggest that selfish routing has modest social cost even without the dissemination of path load information.


NutCracker at WNUT 2020 Task 2 Robustly Identifying Informative COVID 19 Tweets using Ensembling and Adversarial Training. We experiment with COVID Twitter BERT and RoBERTa models to identify informative COVID 19 tweets. We further experiment with adversarial training to make our models robust. The ensemble of COVID Twitter BERT and RoBERTa obtains a F1 score of 0.9096 (on the positive class) on the test data of WNUT 2020 Task 2 and ranks 1st on the leaderboard. The ensemble of the models trained using adversarial training also produces similar result.


Utilizing the redundant constraints for the uplift payment elimination. A power market with non convexities may not have an equilibrium price for power that provides economic stability of the centralized dispatch outcome. In this case, the market players are entitled to receive the uplift payments that compensate the economic profit lost when following the centralized dispatch. We consider a special class of the (possibly non linear) redundant constraints that are redundant not only on the feasible set of the centralized dispatch optimization problem (and, therefore, do not change the centralized dispatch outcome) but also on the larger set obtained when the power balance constraint is relaxed. We show that the Lagrangian relaxation of these redundant constraints may reduce the uplift payments without changing the duality gap. For any given market price (or a pricing algorithm that sets the producer revenue as a function of its output volume) in a uninode multi period power market with fixed load, we explicitly construct a family of the redundant constraints that do not change the maximum profit of the producer and result in zero uplift payment. We show that the introduction and subsequent Lagrangian relaxation of just one redundant constraint in the centralized dispatch problem suffice to eliminate the uplift payments for all the producers. In the case of the convex hull pricing method, the introduction of these redundant constraints affects neither the duality gap nor the market price for power. The results can be straightforwardly generalized to a power market with the price sensitive load.


On COVID 19 Modelling. This contribution analyzes the COVID 19 outbreak by comparably simple mathematical and numerical methods. The final goal is to predict the peak of the epidemic outbreak per country with a reliable technique. This is done by an algorithm motivated by standard SIR models and aligned with the standard data provided by the Johns Hopkins University. To reconstruct data for the unregistered Infected, the algorithm uses current values of the infection fatality rate and a data driven estimation of a specific form of the recovery rate. All other ingredients are data driven as well. Various examples of predictions are provided for illustration.


An Air Suspension to Demonstrate the Properties of Torsion Balances with Fibres of Zero Length. We report on the design and characterisation of an air bearing suspension that has been constructed to highlight the properties of torsion balances with fibres of zero length. A float is levitated on this suspension and its rotational and translational motion in the horizontal plane of the laboratory is controlled using magnetic actuators. We demonstrate the in situ electromagnetic tuning of the float's centre of buoyancy to an accuracy of pm 0.3 mm, which was limited by the noise in the air bearing. The rotational stiffness of the float, which is approximately zero by design, was also measured. We compare the observed behaviour of the float with the predictions of a detailed model of the statics of the float actuator system. Finally, we briefly discuss the application of these ideas and results to the construction of sensitive devices for the measurement of weak forces with short ranges.


Spatio temporal Tubelet Feature Aggregation and Object Linking in Videos. This paper addresses the problem of how to exploit spatio temporal information available in videos to improve the object detection precision. We propose a two stage object detector called FANet based on short term spatio temporal feature aggregation to give a first detection set, and long term object linking to refine these detections. Firstly, we generate a set of short tubelet proposals containing the object in N consecutive frames. Then, we aggregate RoI pooled deep features through the tubelet using a temporal pooling operator that summarizes the information with a fixed size output independent of the number of input frames. On top of that, we define a double head implementation that we feed with spatio temporal aggregated information for spatio temporal object classification, and with spatial information extracted from the current frame for object localization and spatial classification. Furthermore, we also specialize each head branch architecture to better perform in each task taking into account the input data. Finally, a long term linking method builds long tubes using the previously calculated short tubelets to overcome detection errors. We have evaluated our model in the widely used ImageNet VID dataset achieving a 80.9 mAP, which is the new state of the art result for single models. Also, in the challenging small object detection dataset USC GRAD STDdb, our proposal outperforms the single frame baseline by 5.4 mAP.


Algebraic k systems of curves. A collection Delta of simple closed curves on an orientable surface is an algebraic k system if the algebraic intersection number langle alpha, beta rangle is equal to k in absolute value for every alpha , beta in Delta distinct. Generalizing a theorem of MRT14 we compute that the maximum size of an algebraic k system of curves on a surface of genus g is 2g 1 when g ge 3 or k is odd, and 2g otherwise. To illustrate the tightness in our assumptions, we present a construction of curves pairwise geometrically intersecting twice whose size grows as g 2 .


Variational Inference with Mixture Model Approximation Robotic Applications. We propose a method to approximate the distribution of robot configurations satisfying multiple objectives. Our approach uses variational inference, a popular method in Bayesian computation, which has several advantages over sampling based techniques. To be able to represent the complex and multimodal distribution of configurations, we propose to use a mixture model as approximate distribution, an approach that has gained popularity recently. In this work, we show the interesting properties of this approach and how it can be applied to a wide range of problems in robotics.


Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics. A comprehensive microscopic understanding of ambient liquid water is a major challenge for ab initio simulations as it simultaneously requires an accurate quantum mechanical description of the underlying potential energy surface (PES) as well as extensive sampling of configuration space. Due to the presence of light atoms (e.g., H or D), nuclear quantum fluctuations lead to observable changes in the structural properties of liquid water (e.g., isotope effects), and therefore provide yet another challenge for ab initio approaches. In this work, we demonstrate that the combination of dispersion inclusive hybrid density functional theory (DFT), the Feynman discretized path integral (PI) approach, and machine learning (ML) constitutes a versatile ab initio based framework that enables extensive sampling of both thermal and nuclear quantum fluctuations on a quite accurate underlying PES. In particular, we employ the recently developed deep potential molecular dynamics (DPMD) model a neural network representation of the ab initio PES in conjunction with a PI approach based on the generalized Langevin equation (PIGLET) to investigate how isotope effects influence the structural properties of ambient liquid H 2 O and D 2 O. Through a detailed analysis of the interference differential cross sections as well as several radial and angular distribution functions, we demonstrate that this approach can furnish a semi quantitative prediction of these subtle isotope effects.


The Cesium Source Calibration and Monitoring System of the ATLAS Tile Calorimeter Design, Construction and Results. This article describes the design, construction and use of a calibration and monitoring system, based on movable 137Cs gamma ray sources, for the ATLAS Tile Calorimeter (TileCal). The sources, propelled by a water based liquid through tubes that traverse all the calorimeter's cells, produce signals that precisely characterise the response of each tile, thereby providing very granular and accurate data on the response of TileCal to particles. The system has been used to guide and control the quality of the optical instrumentation of all TileCal modules, to set and equalise the dynamic range of the response to physics data, and to set the energy scale of the readout system. In the ATLAS cavern, periodic measurements of the whole detector's response to 137Cs sources allow monitoring the uniformity and stability of all the calorimeter's cells as well as maintaining precise knowledge of its energy calibration. The design of the source hydraulic drive system's hardware and software, the data acquisition system and the data processing algorithms are described. Finally, the results of this two decade program are shown.


Global phase amplitude description of oscillatory dynamics via the parameterization method. In this paper we use the parameterization method to provide a complete description of the dynamics of an n dimensional oscillator beyond the classical phase reduction. The parameterization method allows, via efficient algorithms, to obtain a parameterization of the attracting invariant manifold of the limit cycle in terms of the phase amplitude variables. The method has several advantages. It provides analytically a Fourier Taylor expansion of the parameterization up to any order, as well as a simplification of the dynamics that allows for a numerical globalization of the manifolds. Thus, one can obtain the local and global isochrons and isostables, including the slow attracting manifold, up to high accuracy, which offer a geometrical portrait of the oscillatory dynamics. Furthermore, it provides straightforwardly the infinitesimal Phase and Amplitude Response Functions, that is, the extended infinitesimal Phase and Amplitude Response Curves, which monitor the phase and amplitude shifts beyond the asymptotic state. Thus, the methodology presented yields an accurate description of the phase dynamics for perturbations not restricted to the limit cycle but to its attracting invariant manifold. Finally, we explore some strategies to reduce the dimension of the dynamics, including the reduction of the dynamics to the slow stable submanifold. We illustrate our methods by applying them to different three dimensional single neuron and neural population models in neuroscience.


Spatio Temporal Graph for Video Captioning with Knowledge Distillation. Video captioning is a challenging task that requires a deep understanding of visual scenes. State of the art methods generate captions using either scene level or object level information but without explicitly modeling object interactions. Thus, they often fail to make visually grounded predictions, and are sensitive to spurious correlations. In this paper, we propose a novel spatio temporal graph model for video captioning that exploits object interactions in space and time. Our model builds interpretable links and is able to provide explicit visual grounding. To avoid unstable performance caused by the variable number of objects, we further propose an object aware knowledge distillation mechanism, in which local object information is used to regularize global scene features. We demonstrate the efficacy of our approach through extensive experiments on two benchmarks, showing our approach yields competitive performance with interpretable predictions.


Scalable Derivative Free Optimization for Nonlinear Least Squares Problems. Derivative free or zeroth order optimization (DFO) has gained recent attention for its ability to solve problems in a variety of application areas, including machine learning, particularly involving objectives which are stochastic and or expensive to compute. In this work, we develop a novel model based DFO method for solving nonlinear least squares problems. We improve on state of the art DFO by performing dimensionality reduction in the observational space using sketching methods, avoiding the construction of a full local model. Our approach has a per iteration computational cost which is linear in problem dimension in a big data regime, and numerical evidence demonstrates that, compared to existing software, it has dramatically improved runtime performance on overdetermined least squares problems.


Formula Transformers and Combinatorial Test Generators for Propositional Intuitionistic Theorem Provers. We develop combinatorial test generation algorithms for progressively more powerful theorem provers, covering formula languages ranging from the implicational fragment of intuitionistic logic to full intuitionistic propositional logic. Our algorithms support exhaustive and random generators for formulas of these logics. To provide known to be provable formulas, via the Curry Howard formulas as types correspondence, we use generators for typable lambda terms and combinator expressions. Besides generators for several classes of formulas, we design algorithms that restrict formula generation to canonical representatives among equiprovable formulas and introduce program transformations that reduce formulas to equivalent formulas of a simpler structure. The same transformations, when applied in reverse, create harder formulas that can catch soundness or incompleteness bugs. To test the effectiveness of the testing framework itself, we describe use cases for deriving lightweight theorem provers for several of these logics and for finding bugs in known theorem provers. Our Prolog implementation available at https github.com ptarau TypesAndProofs and a subset of formula generators and theorem provers, implemented in Python is available at https github.com ptarau PythonProvers. Keywords term and formula generation algorithms, Prolog based theorem provers, formulas as types, type inference and type inhabitation, combinatorial testing, finding bugs in theorem provers.


Lyubeznik numbers, F modules and modules of generalized fractions. This paper presents an algorithm for calculation of the Lyubeznik numbers of a local ring which is a homomorphic image of a regular local ring R of prime characteristic. The methods used employ Lyubeznik's F modules over R , particularly his F finite F modules, and also the modules of generalized fractions of Sharp and Zakeri. It is shown that many modules of generalized fractions over R have natural structures as F modules these lead to F module structures on certain local cohomology modules over R , which are exploited, in conjunction with F module structures on injective R modules that result from work of Huneke and Sharp, to compute Lyubeznik numbers. The resulting algorithm has been implemented in Macaulay2.


On line and on board planning and perception for quadrupedal locomotion. We present a legged motion planning approach for quadrupedal locomotion over challenging terrain. We decompose the problem into body action planning and footstep planning. We use a lattice representation together with a set of defined body movement primitives for computing a body action plan. The lattice representation allows us to plan versatile movements that ensure feasibility for every possible plan. To this end, we propose a set of rules that define the footstep search regions and footstep sequence given a body action. We use Anytime Repairing A (ARA ) search that guarantees bounded suboptimal plans. Our main contribution is a planning approach that generates on line versatile movements. Experimental trials demonstrate the performance of our planning approach in a set of challenging terrain conditions. The terrain information and plans are computed on line and on board.


Radiation instability of a relativistic electron beam into a split cavity resonator. The radiation instability in a split cavity asymmetric resonator is considered for the relativistic case. The space charge of an electron beam is taken into account. In the small signal approximation, the energy loss by particles passing through the resonator and the modulation of beam current are investigated. Based on analytical and numerical calculations, it is shown that the symmetric configuration of a split cavity resonator provides the highest rate of instability growth. It is shown that the beam modulation and the efficiency of energy transfer from particles to electromagnetic field decrease with the increase in the initial electron energy. The increase in beam density has a positive effect on the radiation instability growth. It is important to take into account the obtained results when developing generators of electromagnetic radiation and systems for modulating the beam current based on a split cavity resonator.


Crack initiation of printed lines predicted with digital image correlation. Printing of metallic films has been preferred for roll to roll processes over vacuum technologies due to faster processing times and lower processing costs. Films can be produced by depositing inks containing suspended metallic particles within a solvent and then heating the films to both remove the solvent and sinter the particles. The resulting printed structure, electrical and mechanical behavior of the printed films has been studied to better understand their electro mechanical response to loading and eventual brittle fracture. This study evaluated the electro mechanical behavior of 1.25 mu m printed Ag films using in situ resistance and in situ imaging methods. Digital image correlation was utilized with confocal laser scanning microscope images to better visualize crack initiation during tensile straining. This technique showed that cracks initiated earlier in the thicker areas of the film (crests) than in lower areas (troughs) because of a higher density of printing defects and the increased thickness.


Short simplex paths in lattice polytopes. The goal of this paper is to design a simplex algorithm for linear programs on lattice polytopes that traces short' simplex paths from any given vertex to an optimal one. We consider a lattice polytope P contained in 0,k n and defined via m linear inequalities. Our first contribution is a simplex algorithm that reaches an optimal vertex by tracing a path along the edges of P of length in O(n 4 k log(nk) . The length of this path is independent from m and it is the best possible up to a polynomial function. In fact, it is only polynomially far from the worst case diameter, which roughly grows as a linear function in n and k . Motivated by the fact that most known lattice polytopes are defined via 0, pm 1 constraint matrices, our second contribution is an iterative algorithm which exploits the largest absolute value alpha of the entries in the constraint matrix. We show that the length of the simplex path generated by the iterative algorithm is in O(n 2k log(nk alpha)) . In particular, if alpha is bounded by a polynomial in n, k , then the length of the simplex path is in O(n 2k log(nk)) . For both algorithms, the number of arithmetic operations needed to compute the next vertex in the path is polynomial in n , m and log k . If k is polynomially bounded by n and m , the algorithm runs in strongly polynomial time.


Hardware and Interference Limited Cooperative CR NOMA Networks under Imperfect SIC and CSI. The conflation of cognitive radio (CR) and nonorthogonal multiple access (NOMA) concepts is a promising approach to fulfil the massive connectivity goals of future networks given the spectrum scarcity. Accordingly, this letter investigates the outage performance of imperfect cooperative CR NOMA networks under hardware impairments and interference. Our analysis is involved with the derivation of the end to end outage probability (OP) for secondary NOMA users by accounting for imperfect channel state information (CSI), as well as the residual interference caused by successive interference cancellation (SIC) errors and coexisting primary secondary users. The numerical results validated by Monte Carlo simulations show that CR NOMA network provides a superior outage performance over orthogonal multiple access. As imperfections become more significant, CR NOMA is observed to deliver relatively poor outage performance.


Benchmarking Tether based UAV Motion Primitives. This paper proposes and benchmarks two tether based motion primitives for tethered UAVs to execute autonomous flight with proprioception only. Tethered UAVs have been studied mainly due to power and safety considerations. Tether is either not included in the UAV motion (treated same as free flying UAV) or only in terms of station keeping and high speed steady flight. However, feedback from and control over the tether configuration could be utilized as a set of navigational tools for autonomous flight, especially in GPS denied environments and without vision based exteroception. In this work, two tether based motion primitives are proposed, which can enable autonomous flight of a tethered UAV. The proposed motion primitives are implemented on a physical tethered UAV for autonomous path execution with motion capture ground truth. The navigational performance is quantified and compared. The proposed motion primitives make tethered UAV a mobile and safe autonomous robot platform. The benchmarking results suggest appropriate usage of the two motion primitives for tethered UAVs with different path plans.


On Multiple Access in Queue Length Sensitive Systems. We consider transmission of packets over queue length sensitive unreliable links, where packets are randomly corrupted through a noisy channel whose transition probabilities are modulated by the queue length. The goal is to characterize the capacity of this channel. We particularly consider multiple access systems, where transmitters dispatch encoded symbols over a system that is a superposition of continuous time GI k GI 1 queues. A server receives and processes symbols in order of arrivals with queue length dependent noise. We first determine the capacity of single user queue length dependent channels. Further, we characterize the best and worst dispatch processes for GI M 1 queues and the best and worst service processes for M GI 1 queues. Then, the multiple access channel capacity is obtained using point processes. When the number of transmitters is large and each arrival process is sparse, the superposition of arrivals approaches a Poisson point process. In characterizing the Poisson approximation, we show that the capacity of the multiple access system converges to that of a single user M GI 1 queue length dependent system, and an upper bound on the convergence rate is obtained. This implies that the best and worst server behaviors of single user M GI 1 queues are preserved in the sparse multiple access case.


Hairy graphs to ribbon graphs via a fixed source graph complex. We show that the hairy graph complex (HGC n,n ,d) appears as an associated graded complex of the oriented graph complex (OGC n 1 ,d) , subject to the filtration on the number of targets, or equivalently sources, called the fixed source graph complex. The fixed source graph complex (OGC 1,d 0) maps into the ribbon graph complex RGC , which models the moduli space of Riemann surfaces with marked points. The full differential d on the oriented graph complex OGC n 1 corresponds to the deformed differential d h on the hairy graph complex HGC n,n , where h adds a hair. This deformed complex (HGC n,n ,d h) is already known to be quasi isomorphic to standard Kontsevich's graph complex GC 2 n . This gives a new connection between the standard and the oriented version of Kontsevich's graph complex.


A two stage approach for beam hardening artifact reduction in low dose dental CBCT. This paper presents a two stage method for beam hardening artifact correction of dental cone beam computerized tomography (CBCT). The proposed artifact reduction method is designed to improve the quality of maxillofacial imaging, where soft tissue details are not required. Compared to standard CT, the additional difficulty of dental CBCT comes from the problems caused by offset detector, FOV truncation, and low signal to noise ratio due to low X ray irradiation. To address these problems, the proposed method primarily performs a sinogram adjustment in the direction of enhancing data consistency, considering the situation according to the FOV truncation and offset detector. This sinogram correction algorithm significantly reduces beam hardening artifacts caused by high density materials such as teeth, bones, and metal implants, while tending to amplify special types of noise. To suppress such noise, a deep convolutional neural network is complementarily used, where CT images adjusted by the sinogram correction are used as the input of the neural network. Numerous experiments validate that the proposed method successfully reduces beam hardening artifacts and, in particular, has the advantage of improving the image quality of teeth, associated with maxillofacial CBCT imaging.


Deterministic particle approximation of aggregation diffusion equations on unbounded domains. We consider a one dimensional aggregation diffusion equation, which is the gradient flow in the Wasserstein space of a functional with competing attractive repulsive interactions. We prove that the fully deterministic particle approximations with piecewise constant densities introduced in cite Di Francesco Rosini starting from general bounded initial densities converge strongly in L 1 to bounded weak solutions of the PDE. In particular, the result is achieved in unbounded domains and for arbitrary nonnegative bounded initial densities, thus extending the results in cite Gosse Toscani, Matthes Osberger, Mathes Soellner (in which a no vacuum condition is required) and giving an alternative approach to cite Carrillo Craig Patacchini in the one dimensional case, including also subquadratic and superquadratic diffusions.


Multi Decoder Attention Model with Embedding Glimpse for Solving Vehicle Routing Problems. We present a novel deep reinforcement learning method to learn construction heuristics for vehicle routing problems. In specific, we propose a Multi Decoder Attention Model (MDAM) to train multiple diverse policies, which effectively increases the chance of finding good solutions compared with existing methods that train only one policy. A customized beam search strategy is designed to fully exploit the diversity of MDAM. In addition, we propose an Embedding Glimpse layer in MDAM based on the recursive nature of construction, which can improve the quality of each policy by providing more informative embeddings. Extensive experiments on six different routing problems show that our method significantly outperforms the state of the art deep learning based models.


Volatility Estimation of General Gaussian Ornstein Uhlenbeck Process. In this article we study the asymptotic behaviour of the realized quadratic variation of a process int 0 t u s dG H s , where u is a beta H "older continuous process with beta 1 H and G H is a self similar Gaussian process with parameters H in(0,3 4) . We prove almost sure convergence uniformly in time, and a stable weak convergence for the realized quadratic variation. As an application, we construct strongly consistent estimator for the integrated volatility parameter in a model driven by G H .


Spectral analysis of mixing in 2D high Reynolds flows. We use spectral analysis of Eulerian and Lagrangian dynamics to study the advective mixing in an incompressible 2D bounded cavity flow. A significant property of such a rotational flow at high Reynolds numbers is that mixing in its core is slower than wall adjacent areas and corner eddies. We explain this property by considering the resonance between frequencies of unsteady motion captured by the Koopman spectral analysis of the velocity field and the circulation frequency of Lagrangian tracers in the mean flow. In high Reynolds rotational 2D flows, the vorticity in the rotational core is uniformly distributed, which leads to uniform distribution of circulation periods in the mean flow, i.e., the kinematics in the core of mean flow is like rigid body rotation. When this rigid" core is exposed to velocity fluctuations arising from bifurcations at high Reynolds, it shows more resilience toward resonance in Lagrangian motion and hence mixes more slowly compared to other areas. We also show how our qualitative resonance argument extends to chaotic flows where the classical tools of dynamical systems are not applicable.


Determining geopotential difference via relativistic precise point positioning time comparison A case study using simulated observations. According to general relativity theory (GRT), the geopotential difference (GD) can be determined by comparing the change in time difference between precise clocks using the precise point positioning (PPP) time transfer technique, referred to as the relativistic PPP time comparison approach. We focused on high precision time comparison between two precise clocks for determining the GD using the relativistic PPP time transfer,and conducted simulation experiments to validate the approach. In the experiments, we consider three cases to evaluate the performance of the approach using clocks with different stabilities, namely, the frequency stabilities of the clocks equipped at three selected ground stations are respectively (Case 1), (Case 2), and (Case 3) at time period. Conclusions are drawn from the experimental results. First, high precision clocks can significantly improve the accuracy for PPP time transfer, but the improvement is limited by measurement noises. Compared to Case 1, the long term stabilities of OPMT BRUX as well as PTBB BRUX are improved in Cases 2 and 3. The frequency stabilities of Cases 1 3 are approximately 4.28 10 16, 4.00 10 17, and 3.22 10 17 at 10 day averaging time for OPMT BRUX, respectively, and for PTBB BRUX, these values are approximately 3.73 10 16, 8.17 10 17, and 4.64 10 17. Second, the geopotential difference between any two stations can be determined at the decimeter level, with its accuracy being consistent with the stabilities of the time links in Cases 1 3. In Case 3, the determined geopotential differences between OPMT and BRUX deviate from the EIGEN 6C4 model values by 0.64 m2 s2 with an uncertainty of 1.11 m2 s2, whereas the deviation error between PTBB and BRUX is 0.76 m2 s2 with an uncertainty of 1.79 m2 s2. The approach proposed in this study can be also applied to testing GRT.


Moves on k graphs preserving Morita equivalence. We initiate the program of extending to higher rank graphs ( k graphs) the geometric classification of directed graph C algebras, as completed in the 2016 paper of Eilers, Restorff, Ruiz, and Sorensen ERRS16 . To be precise, we identify four "moves," or modifications, one can perform on a k graph Lambda , which leave invariant the Morita equivalence class of its C algebra C ( Lambda) . These moves insplitting, delay, sink deletion, and reduction are inspired by the moves for directed graphs described by Sorensen S o13 and Bates Pask BP04 . Because of this, our perspective on k graphs focuses on the underlying directed graph. We consequently include two new results, Theorem 2.3 and Lemma 2.9, about the relationship between a k graph and its underlying directed graph.


A hybrid discrete continuum approach to model Turing pattern formation. Since its introduction in 1952, Turing's (pre )pattern theory ("the chemical basis of morphogenesis") has been widely applied to a number of areas in developmental biology. The related pattern formation models normally comprise a system of reaction diffusion equations for interacting chemical species ("morphogens"), whose heterogeneous distribution in some spatial domain acts as a template for cells to form some kind of pattern or structure through, for example, differentiation or proliferation induced by the chemical pre pattern. Here we develop a hybrid discrete continuum modelling framework for the formation of cellular patterns via the Turing mechanism. In this framework, a stochastic individual based model of cell movement and proliferation is combined with a reaction diffusion system for the concentrations of some morphogens. As an illustrative example, we focus on a model in which the dynamics of the morphogens are governed by an activator inhibitor system that gives rise to Turing pre patterns. The cells then interact with morphogens in their local area through either of two forms of chemically dependent cell action chemotaxis and chemically controlled proliferation. We begin by considering such a hybrid model posed on static spatial domains, and then turn to the case of growing domains. In both cases, we formally derive the corresponding deterministic continuum limit and show that that there is an excellent quantitative match between the spatial patterns produced by the stochastic individual based model and its deterministic continuum counterpart, when sufficiently large numbers of cells are considered. This paper is intended to present a proof of concept for the ideas underlying the modelling framework, with the aim to then apply the related methods to the study of specific patterning and morphogenetic processes in the future.


Do Truly Unidirectional Surface Plasmon Polaritons Exist?. In this work, we revisit the topic of surface waves on nonreciprocal plasmonic structures, and clarify whether strictly unidirectional surface plasmon polaritons are allowed to exist in this material platform. By investigating different three dimensional configurations and frequency regimes, we theoretically show that, while conventional surface magneto plasmons are not strictly unidirectional due to nonlocal effects, consistent with recent predictions made in the literature, another important class of one way surface plasmon polaritons, existing at an interface with an opaque isotropic material, robustly preserve their unidirectionality even in the presence of nonlocality, and for arbitrarily small levels of dissipation. We also investigate the extreme behavior of terminated unidirectional wave guiding structures, for both classes of surface waves, and discuss their counter intuitive implications.


Self sweeping ytterbium doped fiber laser based on a fiber saturable absorber. Generally speaking, the self sweeping effect relies on the dynamical grating formed in a gain fiber. Here, the normal self sweeping was generated in a pump free ytterbium doped fiber which serves as a fiber saturable absorber and is introduced to the laser cavity by a circulator in this experiment. The sweeping rate and the sweeping range alter as usual, both of which can be controlled by the pump power. Further, a new self pulse signal is observed and discussed in this work, which shows the difference of the self sweeping effects between active fiber and fiber saturable absorber.


Regularisation, optimisation, subregularity. Regularisation theory in Banach spaces, and non norm squared regularisation even in finite dimensions, generally relies upon Bregman divergences to replace norm convergence. This is comparable to the extension of first order optimisation methods to Banach spaces. Bregman divergences can, however, be somewhat suboptimal in terms of descriptiveness. Using the concept of (strong) metric subregularity, previously used to prove the fast local convergence of optimisation methods, we show norm convergence in Banach spaces and for non norm squared regularisation. For problems such as total variation regularised image reconstruction, the metric subregularity reduces to a geometric condition on the ground truth flat areas in the ground truth have to compensate for the fidelity term not having second order growth within the kernel of the forward operator. Our approach to proving such regularisation results is based on optimisation formulations of inverse problems. As a side result of the regularisation theory that we develop, we provide regularisation complexity results for optimisation methods how many steps N delta of the algorithm do we have to take for the approximate solutions to converge as the corruption level delta searrow 0 ?


Non uniqueness in law for two dimensional Navier Stokes equations with diffusion weaker than a full Laplacian. We study the two dimensional Navier Stokes equations forced by random noise with a diffusive term generalized via a fractional Laplacian that has a positive exponent strictly less than one. Because intermittent jets are inherently three dimensional, we instead adapt the theory of intermittent form of the two dimensional stationary flows to the stochastic approach presented by Hofmanov acute mathrm a , Zhu Zhu (2019, arXiv 1912.11841 math.PR ) and prove its non uniqueness in law.


De(con)struction of the lazy F loop improving performance of Smith Waterman alignment. Striped variation of the Smith Waterman algorithm is known as extremely efficient and easily adaptable for the SIMD architectures. However, the potential for improvement has not been exhausted yet. The popular Lazy F loop heuristic requires additional memory access operations, and the worst case performance of the loop could be as bad as the nonvectorized version. We demonstrate the progression of the lazy F loop transformations that improve the loop performance, and ultimately eliminate the loop completely. Our algorithm achieves the best asymptotic performance of all scan based SW algorithms O(n p log(p)), and is very efficient in practice.


Basis Prediction Networks for Effective Burst Denoising with Large Kernels. Bursts of images exhibit significant self similarity across both time and space. This motivates a representation of the kernels as linear combinations of a small set of basis elements. To this end, we introduce a novel basis prediction network that, given an input burst, predicts a set of global basis kernels shared within the image and the corresponding mixing coefficients which are specific to individual pixels. Compared to state of the art techniques that output a large tensor of per pixel spatiotemporal kernels, our formulation substantially reduces the dimensionality of the network output. This allows us to effectively exploit comparatively larger denoising kernels, achieving both significant quality improvements (over 1dB PSNR) and faster run times over state of the art methods.


Inclusion among commutators of elementary subgroups. In the present paper we continue the study of the elementary commutator subgroups E(n,A),E(n,B) , where A and B are two sided ideals of an associative ring R , n ge 3 . First, we refine and expand a number of the auxiliary results, both classical ones, due to Bass, Stein, Mason, Stothers, Tits, Vaserstein, van der Kallen, Stepanov, as also some of the intermediate results in our joint works with Hazrat, and our own recent papers 40,41 . The gimmick of the present paper is an explicit triple congruence for elementary commutators t ij (ab),t ji (c) , where a,b,c belong to three ideals A,B,C of R . In particular, it provides a sharper counterpart of the three subgroups lemma at the level of ideals. We derive some further striking corollaries thereof, such as a complete description of generic lattice of commutator subgroups E(n,I r),E(n,I s) , new inclusions among multiple elementary commutator subgroups, etc.


A new variant of the Optimised Schwarz Method for arbitrary non overlapping subdomain partitions. We consider a scalar wave propagation in harmonic regime modelled by Helmholtz equation with heterogeneous coefficients. Using the Multi Trace Formalism (MTF), we propose a new variant of the Optimized Schwarz Method (OSM) that can accomodate the presence of cross points in the subdomain partition. This leads to the derivation of a strongly coercive formulation of our Helmholtz problem posed on the union of all interfaces. The corresponding operator takes the form "identity contraction".


A Synopses Data Engine for Interactive Extreme Scale Analytics. In this work, we detail the design and structure of a Synopses Data Engine (SDE) which combines the virtues of parallel processing and stream summarization towards delivering interactive analytics at extreme scale. Our SDE is built on top of Apache Flink and implements a synopsis as a service paradigm. In that it achieves (a) concurrently maintaining thousands of synopses of various types for thousands of streams on demand, (b) reusing maintained synopses among various concurrent workflows, (c) providing data summarization facilities even for cross (Big Data) platform workflows, (d) pluggability of new synopses on the fly, (e) increased potential for workflow execution optimization. The proposed SDE is useful for interactive analytics at extreme scales because it enables (i) enhanced horizontal scalability, i.e., not only scaling out the computation to a number of processing units available in a computer cluster, but also harnessing the processing load assigned to each by operating on carefully crafted data summaries, (ii) vertical scalability, i.e., scaling the computation to very high numbers of processed streams and (iii) federated scalability i.e., scaling the computation beyond single clusters and clouds by controlling the communication required to answer global queries posed over a number of potentially geo dispersed clusters.


Heterogeneous Gaussian Mechanism Preserving Differential Privacy in Deep Learning with Provable Robustness. In this paper, we propose a novel Heterogeneous Gaussian Mechanism (HGM) to preserve differential privacy in deep neural networks, with provable robustness against adversarial examples. We first relax the constraint of the privacy budget in the traditional Gaussian Mechanism from (0, 1 to (0, infty), with a new bound of the noise scale to preserve differential privacy. The noise in our mechanism can be arbitrarily redistributed, offering a distinctive ability to address the trade off between model utility and privacy loss. To derive provable robustness, our HGM is applied to inject Gaussian noise into the first hidden layer. Then, a tighter robustness bound is proposed. Theoretical analysis and thorough evaluations show that our mechanism notably improves the robustness of differentially private deep neural networks, compared with baseline approaches, under a variety of model attacks.


LIP Local Importance based Pooling. Spatial downsampling layers are favored in convolutional neural networks (CNNs) to downscale feature maps for larger receptive fields and less memory consumption. However, for discriminative tasks, there is a possibility that these layers lose the discriminative details due to improper pooling strategies, which could hinder the learning process and eventually result in suboptimal models. In this paper, we present a unified framework over the existing downsampling layers (e.g., average pooling, max pooling, and strided convolution) from a local importance view. In this framework, we analyze the issues of these widely used pooling layers and figure out the criteria for designing an effective downsampling layer. According to this analysis, we propose a conceptually simple, general, and effective pooling layer based on local importance modeling, termed as em Local Importance based Pooling (LIP). LIP can automatically enhance discriminative features during the downsampling procedure by learning adaptive importance weights based on inputs. Experiment results show that LIP consistently yields notable gains with different depths and different architectures on ImageNet classification. In the challenging MS COCO dataset, detectors with our LIP ResNets as backbones obtain a consistent improvement ( ge 1.4 ) over the vanilla ResNets, and especially achieve the current state of the art performance in detecting small objects under the single scale testing scheme.


Optically levitated rotor at its thermal limit of frequency stability. Optically levitated rotors are prime candidates for torque sensors whose precision is limited by the fluctuationsof the rotation frequency. In this work we investigate an optically levitated rotor at its fundamental thermal limitof frequency stability, where rotation frequency fluctuations arise solely due to coupling to the thermal bath.


Multimodal 3D Object Detection from Simulated Pretraining. The need for simulated data in autonomous driving applications has become increasingly important, both for validation of pretrained models and for training new models. In order for these models to generalize to real world applications, it is critical that the underlying dataset contains a variety of driving scenarios and that simulated sensor readings closely mimics real world sensors. We present the Carla Automated Dataset Extraction Tool (CADET), a novel tool for generating training data from the CARLA simulator to be used in autonomous driving research. The tool is able to export high quality, synchronized LIDAR and camera data with object annotations, and offers configuration to accurately reflect a real life sensor array. Furthermore, we use this tool to generate a dataset consisting of 10 000 samples and use this dataset in order to train the 3D object detection network AVOD FPN, with finetuning on the KITTI dataset in order to evaluate the potential for effective pretraining. We also present two novel LIDAR feature map configurations in Bird's Eye View for use with AVOD FPN that can be easily modified. These configurations are tested on the KITTI and CADET datasets in order to evaluate their performance as well as the usability of the simulated dataset for pretraining. Although insufficient to fully replace the use of real world data, and generally not able to exceed the performance of systems fully trained on real data, our results indicate that simulated data can considerably reduce the amount of training on real data required to achieve satisfactory levels of accuracy.


Approximate Evaluation of First Order Counting Queries. Kuske and Schweikardt introduced the very expressive first order counting logic FOC(P) to model database queries with counting operations. They showed that there is an efficient model checking algorithm on graphs with bounded degree, while Grohe and Schweikardt showed that probably no such algorithm exists for trees of bounded depth. We analyze the fragment FO( 0 ) of this logic. While we remove for example subtraction and comparison between two non atomic counting terms, this logic remains quite expressive We allow nested counting and comparison between counting terms and arbitrarily large numbers. Our main result is an approximation scheme of the model checking problem for FO( 0 ) that runs in linear fpt time on structures with bounded expansion. This scheme either gives the correct answer or says "I do not know." The latter answer may only be given if small perturbations in the number symbols of the formula could make it both satisfied and unsatisfied. This is complemented by showing that exactly solving the model checking problem for FO( 0 ) is already hard on trees of bounded depth and just slightly increasing the expressiveness of FO( 0 ) makes even approximation hard on trees.


On weak metric dimension of digraphs. Using the two way distance, we introduce the concepts of weak metric dimension of a strongly connected digraph Gamma . We first establish lower and upper bounds for the number of arcs in Gamma by using the diameter and weak metric dimension of Gamma , and characterize all digraphs attaining the lower or upper bound. Then we study a digraph with weak metric dimension 1 and classify all vertex transitive digraphs having weak metric dimension 1 . Finally, all digraphs of order n with weak metric dimension n 1 or n 2 are determined.


An experiment to address conceptual difficulties in slipping and rolling problems. A bicycle wheel that was initially spinning freely was placed in contact with a rough surface and a digital film was made of its motion. Using Tracker software for video analysis, we obtained the velocity vectors for several points on the wheel, in the frame of reference of the laboratory as well as in a relative frame of reference having as its origin the wheel s center of mass. The velocity of the wheel s point of contact with the floor was also determined obtaining then a complete picture of the kinematic state of the wheel in both frames of reference. An empirical approach of this sort to problems in mechanics can contribute to overcoming the considerable difficulties they entail.


QoE driven Coupled Uplink and Downlink Rate Adaptation for 360 degree Video Live Streaming. 360 degree video provides an immersive 360 degree viewing experience and has been widely used in many areas. The 360 degree video live streaming systems involve capturing, compression, uplink (camera to video server) and downlink (video server to user) transmissions. However, few studies have jointly investigated such complex systems, especially the rate adaptation for the coupled uplink and downlink in the 360 degree video streaming under limited bandwidth constraints. In this letter, we propose a quality of experience (QoE) driven 360 degree video live streaming system, in which a video server performs rate adaptation based on the uplink and downlink bandwidths and information concerning each user's real time field of view (FOV). We formulate it as a nonlinear integer programming problem and propose an algorithm, which combines the Karush Kuhn Tucker (KKT) condition and branch and bound method, to solve it. The numerical results show that the proposed optimization model can improve users' QoE significantly in comparison with other baseline schemes.


Adapting Workflow Management Systems to BFT Blockchains The YAWL Example. Blockchain technology provides an auditable and tamper proof distributed storage infrastructure for information records. This can be leveraged to support distributed workflow management. Compared to proof of work consensus, popularized by Bitcoin and Ethereum, blockchains based on BFT (byzantine fault tolerance) ordering consensus trade scalability for immediacy and finality of consensus. This makes them easier to use as distribution infrastructure, as applications need not be adapted to deal with eventual consistency and delayed consensus of proof of work blockchains. Hence, applications such as workflow engines can be easily ported to such a blockchain infrastructure to take advantage of their decentralized integrity assurance and information distribution model. In this paper we describe how the YAWL workflow engine can be used on a BFT based blockchain infrastructure to enable collaborative workflows across different organizations.