Domain decomposition least squares Petrov Galerkin (DD LSPG) nonlinear model reduction. A novel domain decomposition least squares Petrov Galerkin (DD LSPG) model reduction method applicable to parameterized systems of nonlinear algebraic equations (e.g., arising from discretizing a parameterized partial differential equations problem) is proposed. In contrast with previous works, we adopt an algebraically non overlapping decomposition strategy rather than a spatial decomposition strategy, which facilitates application to different spatial discretization schemes. Rather than constructing a low dimensional subspace for the entire state space in a monolithic fashion, the methodology constructs separate subspaces for the different subdomains components characterizing the original model. In the offline stage, the method constructs low dimensional bases for the interior and interface of components. In the online stage, the approach constructs an LSPG ROM for each component and enforces strong or weak compatibility on the 'ports' connecting them. We propose four different ways to construct reduced bases on the interface ports of subdomains and several ways to enforce compatibility across connecting ports. We derive a posteriori and a priori error bounds for the DD LSPG solutions. Numerical results performed on nonlinear benchmark problems in heat transfer and fluid dynamics demonstrate that the proposed method performs well in terms of both accuracy and computational cost, with different choices of basis and compatibility constraints yielding different performance profiles.


Exact electrostatic energy calculation for charged systems neutralized by uniformly distributed background charge using fast multipole method and its application to efficient free energy calculation. In the molecular dynamics calculations for the free energy of ions and ionic molecules, we often encounter wet charged molecular systems where electrical neutrality condition is broken. This causes a problem in the evaluation of electrostatic interaction under periodic boundary condition. A standard remedy for the problem is to consider a hypothetical homogeneous background charge density to neutralize the total system. Here, we present a new expression for the evaluation of electrostatic interactions for the system including the background charge by fast multipole method (FMM). Further, an efficient scheme to evaluate solute solvent interaction energy by FMM has been developed to reduce the computation of far field part. We have calculated hydration free energy of ions, Mg 2 , Na , and Cl dissolved in neutral solvent using the new expression. The calculated free energy showed a good agreement with the result using well established particle mesh Ewald method, demonstrating the validity of the present expression in the framework of FMM. An advantage of the present scheme is in an efficient free energy calculation of a large scale charged systems (particularly over million particles) based on highly parallel computations.


Polyanalytic Hardy decomposition of higher order Lipschitz functions. This paper is concerned with the problem of decomposing a higher order Lipschitz function on a closed Jordan curve Gamma into a sum of two polyanalytic functions in each open domain defined by Gamma . Our basic tools are the Hardy projections related to a singular integral operator arising in polyanalytic function theory, which, as it is proved here, represents an involution operator on the higher order Lipschitz classes. Our result generalizes the classical Hardy decomposition of Holder continuous functions on the boundary of a domain.


Scaling Systematic Literature Reviews with Machine Learning Pipelines. Systematic reviews, which entail the extraction of data from large numbers of scientific documents, are an ideal avenue for the application of machine learning. They are vital to many fields of science and philanthropy, but are very time consuming and require experts. Yet the three main stages of a systematic review are easily done automatically searching for documents can be done via APIs and scrapers, selection of relevant documents can be done via binary classification, and extraction of data can be done via sequence labelling classification. Despite the promise of automation for this field, little research exists that examines the various ways to automate each of these tasks. We construct a pipeline that automates each of these aspects, and experiment with many human time vs. system quality trade offs. We test the ability of classifiers to work well on small amounts of data and to generalise to data from countries not represented in the training data. We test different types of data extraction with varying difficulty in annotation, and five different neural architectures to do the extraction. We find that we can get surprising accuracy and generalisability of the whole pipeline system with only 2 weeks of human expert annotation, which is only 15 of the time it takes to do the whole review manually and can be repeated and extended to new data with no additional effort.


Universal to Specific Framework for Complex Action Recognition. Video based action recognition has recently attracted much attention in the field of computer vision. To solve more complex recognition tasks, it has become necessary to distinguish different levels of interclass variations. Inspired by a common flowchart based on the human decision making process that first narrows down the probable classes and then applies a "rethinking" process for finer level recognition, we propose an effective universal to specific (U2S) framework for complex action recognition. The U2S framework is composed of three subnetworks a universal network, a category specific network, and a mask network. The universal network first learns universal feature representations. The mask network then generates attention masks for confusing classes through category regularization based on the output of the universal network. The mask is further used to guide the category specific network for class specific feature representations. The entire framework is optimized in an end to end manner. Experiments on a variety of benchmark datasets, e.g., the Something Something, UCF101, and HMDB51 datasets, demonstrate the effectiveness of the U2S framework i.e., U2S can focus on discriminative spatiotemporal regions for confusing categories. We further visualize the relationship between different classes, showing that U2S indeed improves the discriminability of learned features. Moreover, the proposed U2S model is a general framework and may adopt any base recognition network.


Agent Based Simulation of Collective Cooperation From Experiment to Model. Simulation models of pedestrian dynamics have become an invaluable tool for evacuation planning. Typically crowds are assumed to stream unidirectionally towards a safe area. Simulated agents avoid collisions through mechanisms that belong to each individual, such as being repelled from each other by imaginary forces. But classic locomotion models fail when collective cooperation is called for, notably when an agent, say a first aid attendant, needs to forge a path through a densely packed group. We present a controlled experiment to observe what happens when humans pass through a dense static crowd. We formulate and test hypothesis on salient phenomena. We discuss our observations in a psychological framework. We derive a model that incorporates agents' perception and cognitive processing of a situation that needs cooperation selection from a portfolio of behaviours, such as being cooperative and a suitable action, such as swapping places. Agents' ability to successfully get through a dense crowd emerges as an effect of the psychological model.


Existence of Large deviations rate function for any S unimodal map. For an arbitrary negative Schwarzian unimodal map with non flat critical point, we establish the level 2 Large Deviation Principle (LDP) for empirical distributions. We also give an example of a multimodal map for which the level 2 LDP does not hold.


Forcing Hamiltonicity in locally finite graphs via forbidden induced subgraphs II paws. In this paper we extend a result about a sufficient condition for Hamiltonicity for finite graphs by Broersma and Veldmann to locally finite graphs. In order to do this we use topological circles within the Freudenthal compactification of a locally finite graph as infinite cycles. The condition we focus on in this paper is in terms of forbidden induced subgraphs, namely being claw free and a relaxation of being paw free.


Proxy Tasks and Subjective Measures Can Be Misleading in Evaluating Explainable AI Systems. Explainable artificially intelligent (XAI) systems form part of sociotechnical systems, e.g., human AI teams tasked with making decisions. Yet, current XAI systems are rarely evaluated by measuring the performance of human AI teams on actual decision making tasks. We conducted two online experiments and one in person think aloud study to evaluate two currently common techniques for evaluating XAI systems (1) using proxy, artificial tasks such as how well humans predict the AI's decision from the given explanations, and (2) using subjective measures of trust and preference as predictors of actual performance. The results of our experiments demonstrate that evaluations with proxy tasks did not predict the results of the evaluations with the actual decision making tasks. Further, the subjective measures on evaluations with actual decision making tasks did not predict the objective performance on those same tasks. Our results suggest that by employing misleading evaluation methods, our field may be inadvertently slowing its progress toward developing human AI teams that can reliably perform better than humans or AIs alone.


Equilibria and Systemic Risk in Saturated Networks. We undertake a fundamental study of network equilibria modeled as solutions of fixed point equations for monotone linear functions with saturation nonlinearities. The considered model extends one originally proposed to study systemic risk in networks of financial institutions interconnected by mutual obligations and is one of the simplest continuous models accounting for shock propagation phenomena and cascading failure effects. It also characterizes Nash equilibria of constrained quadratic network games with strategic complementarities. We first derive explicit expressions for network equilibria and prove necessary and sufficient conditions for their uniqueness encompassing and generalizing results available in the literature. Then, we study jump discontinuities of the network equilibria when the exogenous flows cross certain regions of measure 0 representable as graphs of continuous functions. Finally, we discuss some implications of our results in the two main motivating applications. In financial networks, this bifurcation phenomenon is responsible for how small shocks in the assets of a few nodes can trigger major aggregate losses to the system and cause the default of several agents. In constrained quadratic network games, it induces a blow up behavior of the sensitivity of Nash equilibria with respect to the individual benefits.


When HLS Meets FPGA HBM Benchmarking and Bandwidth Optimization. With the recent release of High Bandwidth Memory (HBM) based FPGA boards, developers can now exploit unprecedented external memory bandwidth. This allows more memory bounded applications to benefit from FPGA acceleration. However, we found that it is not easy to fully utilize the available bandwidth when developing some applications with high level synthesis (HLS) tools. This is due to the limitation of existing HLS tools when accessing HBM board's large number of independent external memory channels. In this paper, we measure the performance of three recent representative HBM FPGA boards (Intel's Stratix 10 MX and Xilinx's Alveo U50 U280 boards) with microbenchmarks and analyze the HLS overhead. Next, we propose HLS based optimization techniques to improve the effective bandwidth when a PE accesses multiple HBM channels or multiple PEs access an HBM channel. Our experiment demonstrates that the effective bandwidth improves by 2.4X 3.8X. We also provide a list of insights for future improvement of the HBM FPGA HLS design flow.


Fano manifolds and stability of tangent bundles. We determine the stability instability of the tangent bundles of the Fano varieties in a certain class of two orbit varieties, which are classified by Pasquier in 2009. As a consequence, we show that some of these varieties admit unstable tangent bundles, which disproves a conjecture on stability of tangent bundles of Fano manifolds.


Human Motion Tracking by Registering an Articulated Surface to 3 D Points and Normals. We address the problem of human motion tracking by registering a surface to 3 D data. We propose a method that iteratively computes two things Maximum likelihood estimates for both the kinematic and free motion parameters of a kinematic human body representation, as well as probabilities that the data are assigned either to a body part, or to an outlier cluster. We introduce a new metric between observed points and normals on one side, and a parameterized surface on the other side, the latter being defined as a blending over a set of ellipsoids. We claim that this metric is well suited when one deals with either visual hull or visual shape observations. We illustrate the method by tracking human motions using sparse visual shape data (3 D surface points and normals) gathered from imperfect silhouettes.


On the second dual space of Banach space of vector valued little Lipschitz functions. Let (X ) be a compact metric space and (E ) be a Banach space. ( lip (X, E) ) denotes the Banach space of all (E ) valued little Lipschitz functions on (X ). We show that ( lip (X, E) ) is isometrically isomorphic to Banach space of (E ) valued Lipschitz functions ( Lip(X, E ) ) under several conditions. Moreover, we describe the isometric isomorphism from ( lip (X, E) ) to ( Lip (X, E ) ).


FRESH Interactive Reward Shaping in High Dimensional State Spaces using Human Feedback. Reinforcement learning has been successful in training autonomous agents to accomplish goals in complex environments. Although this has been adapted to multiple settings, including robotics and computer games, human players often find it easier to obtain higher rewards in some environments than reinforcement learning algorithms. This is especially true of high dimensional state spaces where the reward obtained by the agent is sparse or extremely delayed. In this paper, we seek to effectively integrate feedback signals supplied by a human operator with deep reinforcement learning algorithms in high dimensional state spaces. We call this FRESH (Feedback based REward SHaping). During training, a human operator is presented with trajectories from a replay buffer and then provides feedback on states and actions in the trajectory. In order to generalize feedback signals provided by the human operator to previously unseen states and actions at test time, we use a feedback neural network. We use an ensemble of neural networks with a shared network architecture to represent model uncertainty and the confidence of the neural network in its output. The output of the feedback neural network is converted to a shaping reward that is augmented to the reward provided by the environment. We evaluate our approach on the Bowling and Skiing Atari games in the arcade learning environment. Although human experts have been able to achieve high scores in these environments, state of the art deep learning algorithms perform poorly. We observe that FRESH is able to achieve much higher scores than state of the art deep learning algorithms in both environments. FRESH also achieves a 21.4 higher score than a human expert in Bowling and does as well as a human expert in Skiing.


Dam break in rectangular channels with different upstream downstream widths. The classic Stoker dam break problem is revisited in cases of different channel widths upstream and downstream of the dam. The channel is supposed to have a rectangular cross section and a horizontal and frictionless bottom. The system of the shallow water equations is enriched, using the width as a space dependent variable, together with the depth and the unit discharge, which conversely depend on both space and time. Such a formulation allows a quasi analytical treatment of the system, whose solution is similar to that of the classic Stoker solution when the downstream upstream depth ratio is sufficiently large, except that a further stationary contact wave exists at the dam position. When the downstream upstream depth ratio is small, the solution is richer than the Stoker solution because the critical state occurs at the dam position and the solution itself becomes resonant at the same position, where two eigenvalues are null and the strict hyperbolicity of the system is lost. The limits that identify the flow regime for channel contraction and channel expansion are discussed after showing that the nondimensional parameters governing the problem are the downstream upstream width ratio and the downstream upstream initial depth ratio. After the introduction of the previous analytical framework, a numerical analysis is also performed to evaluate a numerical method that is conceived to suitably capture rarefactions, shock waves and contact waves. A second order method is adopted, employing a Dumbser Osher Toro Riemann solver equipped with a nonlinear path. Such an original nonlinear path is shown to perform better than the classic linear path when contact waves of large amplitude must be captured, being able to obtain specific energy conservation and mass conservation at the singularity. The codes, written in MATLAB language, are made available in Mendeley Data repository.


Attractive Riesz potentials acting on hard spheres. In this paper we introduce a model for hard spheres interacting through attractive Riesz type potentials, and we study its thermodynamic limit. We show that the tail energy enforces optimal packing and round macroscopic shapes.


Learning Visually Consistent Label Embeddings for Zero Shot Learning. In this work, we propose a zero shot learning method to effectively model knowledge transfer between classes via jointly learning visually consistent word vectors and label embedding model in an end to end manner. The main idea is to project the vector space word vectors of attributes and classes into the visual space such that word representations of semantically related classes become more closer, and use the projected vectors in the proposed embedding model to identify unseen classes. We evaluate the proposed approach on two benchmark datasets and the experimental results show that our method yields significant improvements in recognition accuracy.


A note on pseudorandom Ramsey graphs. For fixed s ge 3 , we prove that if optimal K s free pseudorandom graphs exist, then the Ramsey number r(s,t) t s 1 o(1) as t rightarrow infty . Our method also improves the best lower bounds for r(C ell ,t) obtained by Bohman and Keevash from the random C ell free process by polylogarithmic factors for all odd ell geq 5 and ell in 6,10 . For ell 4 it matches their lower bound from the C 4 free process. We also prove, via a different approach, that r(C 5, t) (1 o(1))t 11 8 and r(C 7, t) (1 o(1))t 11 9 . These improve the exponent of t in the previous best results and appear to be the first examples of graphs F with cycles for which such an improvement of the exponent for r(F, t) is shown over the bounds given by the random F free process and random graphs.


Codes, Cubes, and Graphical Designs. Graphical designs are an extension of spherical designs to functions on graphs. We connect linear codes to graphical designs on cube graphs, and show that the Hamming code in particular is a highly effective graphical design. We show that even in highly structured graphs, graphical designs are distinct from the related concepts of extremal designs, maximum stable sets in distance graphs, and t designs on association schemes.


How to Present News on Social Media A Causal Analysis of Editing News Headlines for Boosting User Engagement. To reach a broader audience and optimize traffic toward news articles, media outlets commonly run social media accounts and share their content with a short text summary. Despite its importance of writing a compelling message in sharing articles, the research community does not own a sufficient understanding of what kinds of editing strategies effectively promote audience engagement. In this study, we aim to fill the gap by analyzing media outlets' current practices using a data driven approach. We first build a parallel corpus of original news articles and their corresponding tweets that eight media outlets shared. Then, we explore how those media edited tweets against original headlines and the effects of such changes. To estimate the effects of editing news headlines for social media sharing in audience engagement, we present a systematic analysis that incorporates a causal inference technique with deep learning using propensity score matching, it allows for estimating potential (dis )advantages of an editing style compared to counterfactual cases where a similar news article is shared with a different style. According to the analyses of various editing styles, we report common and differing effects of the styles across the outlets. To understand the effects of various editing styles, media outlets could apply our easy to use tool by themselves.


Relaxed parameter conditions for chemotactic collapse in logistic type parabolic elliptic Keller Segel systems. We study the finite time blow up in two variants of the parabolic elliptic Keller Segel system with nonlinear diffusion and logistic source. In n dimensional balls, we consider begin align begin cases u t nabla cdot ((u 1) m 1 nabla u u nabla v) lambda u mu u 1 kappa , 0 Delta v frac1 Omega int Omega u u end cases tag JL end align and begin align begin cases u t nabla cdot ((u 1) m 1 nabla u u nabla v) lambda u mu u 1 kappa , 0 Delta v v u, end cases tag PE end align where lambda and mu are given spatially radial nonnegative functions and m, kappa 0 are given parameters subject to further conditions. In a unified treatment, we establish a bridge between previously employed methods on blow up detection and relatively new results on pointwise upper estimates of solutions in both of the systems above and then, making use of this newly found connection, provide extended parameter ranges for m, kappa leading to the existence of finite time blow up solutions in space dimensions three and above. In particular, for constant lambda, mu 0 , we find that there are initial data which lead to blow up in (JL) if begin alignat 2 0 leq kappa min left frac 1 2 , frac n 2 n (m 1) right qquad text if m in left frac 2 n , frac 2n 2 n right) text or quad 0 leq kappa min left frac 1 2 , frac n 1 n frac m 2 right qquad text if m in left(0, frac 2 n right), end alignat and in (PE) if m in 1, frac 2n 2 n ) and begin align 0 leq kappa min left frac (m 1) n 1 2(n 1) , frac n 2 (m 1) n n(n 1) right . end align


Sparse regularization of inverse problems by operator adapted frame thresholding. We analyze sparse frame based regularization of inverse problems by means of a diagonal frame decomposition (DFD) for the forward operator, which generalizes the SVD. The DFD allows to define a non iterative (direct) operator adapted frame thresholding approach which we show to provide a convergent regularization method with linear convergence rates. These results will be compared to the well known analysis and synthesis variants of sparse ell 1 regularization which are usually implemented thorough iterative schemes. If the frame is a basis (non redundant case), the three versions of sparse regularization, namely synthesis and analysis variants of ell 1 regularization as well as the DFD thresholding are equivalent. However, in the redundant case, those three approaches are pairwise different.


Cryogenic low power CMOS analog buffer at 4.2K. A novel power efficient analog buffer at liquid helium temperature is proposed. The proposed circuit is based on an input stage consisting of two complementary differential pairs to achieve rail to rail level tracking. Results of simulation based on SMIC 0.18um CMOS technology show the high driving capability and low quiescent power consumption at cryogenic temperature. Operating at single 1.4 V supply, the circuit could achieve a slew rate of 51 V us and 93 V us for 10 pF capacitive load. The static power of the circuit is only 79uW.


Computational prediction of RNA tertiary structures using machine learning methods. RNAs play crucial and versatile roles in biological processes. Computational prediction approaches can help to understand RNA structures and their stabilizing factors, thus providing information on their functions, and facilitating the design of new RNAs. Machine learning (ML) techniques have made tremendous progress in many fields in the past few years. Although their usage in protein related fields has a long history, the use of ML methods in predicting RNA tertiary structures is new and rare. Here, we review the recent advances of using ML methods on RNA structure predictions and discuss the advantages and limitation, the difficulties and potentials of these approaches when applied in the field.


Reflection positivity on spheres. In this article we specialize a construction of a reflection positive Hilbert space due to Dimock and Jaffe Ritter to the sphere mathbb S n . We determine the resulting Osterwalder Schrader Hilbert space, a construction that can be viewed as the step from euclidean to relativistic quantum field theory. We show that this process gives rise to an irreducible unitary spherical representation of the orthochronous Lorentz group G c mathrm O 1,n ( mathbb R ) uparrow and that the representations thus obtained are the irreducible unitary spherical representations of this group. A key tool is a certain complex domain Xi , known as the crown of the hyperboloid, containing a half sphere mathbb S n and the hyperboloid mathbb H n as totally real submanifolds. This domain provides a bridge between those two manifolds when we study unitary representations of G c in spaces of holomorphic functions on Xi . We connect this analysis with the boundary components which are the de Sitter space and the Lorentz cone of future pointing light like vectors.


On the roots of the Poupard and Kreweras polynomials. The Poupard polynomials are polynomials in one variable with integer coefficients, with some close relationship to Bernoulli and tangent numbers. They also have a combinatorial interpretation. We prove that every Poupard polynomial has all its roots on the unit circle. We also obtain the same property for another sequence of polynomials introduced by Kreweras and related to Genocchi numbers. This is obtained through a general statement about some linear operators acting on palindromic polynomials.


physics.med-ph cs.LG physics.bio-ph q-bio.GN

Next Generation Radiogenomics Sequencing for Prediction of EGFR and KRAS Mutation Status in NSCLC Patients Using Multimodal Imaging and Machine Learning Approaches. Aim In the present work, we aimed to evaluate a comprehensive radiomics framework that enabled prediction of EGFR and KRAS mutation status in NSCLC cancer patients based on PET and CT multi modalities radiomic features and machine learning (ML) algorithms. Methods Our study involved 211 NSCLC cancer patient with PET and CTD images. More than twenty thousand radiomic features from different image feature sets were extracted Feature value was normalized to obtain Z scores, followed by student t test students for comparison, high correlated features were eliminated and the False discovery rate (FDR) correction were performed Six feature selection methods and twelve classifiers were used to predict gene status in patient and model evaluation was reported on independent validation sets (68 patients). Results The best predictive power of conventional PET parameters was achieved by SUVpeak (AUC 0.69, P value 0.0002) and MTV (AUC 0.55, P value 0.0011) for EGFR and KRAS, respectively. Univariate analysis of radiomics features improved prediction power up to AUC 75 (q value 0.003, Short Run Emphasis feature of GLRLM from LOG preprocessed image of PET with sigma value 1.5) and AUC 0.71 (q value 0.00005, The Large Dependence Low Gray Level Emphasis from GLDM in LOG preprocessed image of CTD sigma value 5) for EGFR and KRAS, respectively. Furthermore, the machine learning algorithm improved the perdition power up to AUC 0.82 for EGFR (LOG preprocessed of PET image set with sigma 3 with VT feature selector and SGD classifier) and AUC 0.83 for KRAS (CT image set with sigma 3.5 with SM feature selector and SGD classifier). Conclusion We demonstrated that radiomic features extracted from different image feature sets could be used for EGFR and KRAS mutation status prediction in NSCLC patients, and showed that they have more predictive power than conventional imaging parameters.


Non trivial informational closure of a Bayesian hyperparameter. We investigate the non trivial informational closure (NTIC) of a Bayesian hyperparameter inferring the underlying distribution of an identically and independently distributed finite random variable. For this we embed both the Bayesian hyper parameter updating process and the random data process into a Markov chain. The original publication by Bertschinger et al. (2006) mentioned that NTIC may be able to capture an abstract notion of modeling that is agnostic to the specific internal structure of and existence of explicit representations within the modeling process. The Bayesian hyperparameter is of interest since it has a well defined interpretation as a model of the data process and at the same time its dynamics can be specified without reference to this interpretation. On the one hand we show explicitly that the NTIC of the hyperparameter increases indefinitely over time. On the other hand we attempt to establish a connection between a quantity that is a feature of the interpretation of the hyperparameter as a model, namely the information gain, and the one step pointwise NTIC which is a quantity that does not depend on this interpretation. We find that in general we cannot use the one step pointwise NTIC as an indicator for information gain. We hope this exploratory work can lead to further rigorous studies of the relation between NTIC and modeling.


A Knotted Meta molecule with 2 D Isotropic Optical Activity Rotating the Incident Polarization by 90 deg . Optical activity is the ability of chiral materials to rotate linearly polarized (LP) electromagnetic waves. Because of their intrinsic asymmetry, traditional chiral molecules usually lack isotropic performance, or at best only possess a weak form of chirality. Here we introduce a knotted chiral meta molecule that exhibits optical activity corresponding to a 90 deg polarization rotation of the incident waves. More importantly, arising from the continuous multi fold rotational symmetry of the chiral torus knot structure, the observed polarization rotation behavior is found to be independent of how the incident wave is polarized. In other words, the proposed chiral knot structure possesses two dimensional (2 D) isotropic optical activity as illustrated in Fig. 1, which has been experimentally validated in the microwave spectrum. The proposed chiral torus knot represents the most optically active meta molecule reported to date that is intrinsically isotropic to the incident polarization.


Surface bundles and the section conjecture. We formulate a tropical analogue of Grothendieck's section conjecture that for every stable graph G of genus g 2, and every field k, the generic curve with reduction type G over k satisfies the section conjecture. We prove many cases of this conjecture. In so doing we produce many examples of curves satisfying the section conjecture over fields of geometric interest, and then over p adic fields and number fields via a Chebotarev argument. We construct two Galois cohomology classes o 1 and o 2, which obstruct the existence of pi 1 sections and hence of rational points. The first is an abelian obstruction, closely related to the period of a curve and to a cohomology class on the moduli space of curves M g studied by Morita. The second is a 2 nilpotent obstruction and appears to be new. We study the degeneration of these classes via topological techniques, and we produce examples of surface bundles over surfaces where these classes obstruct sections. We then use these constructions to produce curves over p adic fields and number fields where each class obstructs pi 1 sections and hence rational points. Among our geometric results are a new proof of the section conjecture for the generic curve of genus g 2, and a proof of the section conjecture for the generic curve of even genus with a rational divisor class of degree one (where the obstruction to the existence of a section is genuinely non abelian).


Introductory E M Lab Manual for Virtual Teaching. Introductory electricity and magnetism lab manual was designed to use with virtual Physics II class. The lab manual consists of experiments on electrostatics, electric potential and energy, current and resistance, DC circuits, electromagnetism, and AC circuits. Virtual experiments were based on simulations. Open educational resources (OER) were used for all experiments. Virtual experiments were designed to simulate in person physical lab experiments. Special emphasis was given to computational data analysis with excel. Formatted excel sheets per each lab were given to students and step by step calculation in excel were explained during the synchronous class. Learning management system (LMS) was used to fully web enhance the lab class. Virtual labs were delivered by using live video conference technology and recorded lab sessions were added to LMS. Lab class were tested with both virtual delivery methods (synchronous and asynchronous). Student learning outcomes (understand, apply, analyze and evaluate) were studied with detailed lab reports and end of the semester lab based written exam which confirmed the virtual lab class was as effective as the in person physical lab class.


Nonuniform Multiresolution Analysis Associated with Linear Canonical Transform. The linear canonical transform (LCT) has attained respectable status within a short span and is being broadly employed across several disciplines of science and engineering including signal processing, optical and radar systems, electrical and communication systems, quantum physics etc, mainly due to the extra degrees of freedom and simple geometrical manifestation. In this article, we introduce a novel multiresolution analysis (LCT NUMRA) on the spectrum Omega big 0, r N big 2 mathbb Z , where N geqq 1 is an integer and r is an odd integer with 1 leqq r leqq 2N 1, such that r and N are relatively prime, by intertwining the ideas of nonuniform MRA and linear canonical transforms. We first develop nonuniform multiresolution analysis associated with the LCT and then drive an algorithm to construct an LCT NUMRA starting from a linear canonical low pass filter Lambda 0 M ( omega) with suitable conditions. Nevertheless, to extend the scope of the present study, we construct the associated wavelet packets for such an MRA and investigate their properties by means of the linear canonical transforms.


On List Decoding of 5G NR Polar Codes. The 5th generation wireless systems (5G) standardization process of the 3rd generation partnership project (3GPP) chose polar codes as a channel coding scheme for the control channel. In case of downlink control information, polar codes are concatenated with distributed distributed cyclic redundancy check (CRC). Whereas CRC bits allow to improve the performance of successive cancellation list (SCL) decoders by improving distance properties, distributed CRC bits allow for path pruning and decoding early termination. In this paper, we show how to take advantage of the distributed CRC to improve SCL decoding, analyzing various schemes having different earlytermination and error correction properties. Simulation results compare the proposed decoding schemes, showing different tradeoffs between error correction performance and early termination with different decoder parameters.


BAGS An automatic homework grading system using the pictures taken by smart phones. Homework grading is critical to evaluate teaching quality and effect. However, it is usually time consuming to grade the homework manually. In automatic homework grading scenario, many optical mark reader (OMR) based solutions which require specific equipments have been proposed. Although many of them can achieve relatively high accuracy, they are less convenient for users. In contrast, with the popularity of smart phones, the automatic grading system which depends on the image photographed by phones becomes more available. In practice, due to different photographing angles or uneven papers, images may be distorted. Moreover, most of images are photographed under complex backgrounds, making answer areas detection more difficult. To solve these problems, we propose BAGS, an automatic homework grading system which can effectively locate and recognize handwritten answers. In BAGS, all the answers would be written above the answer area underlines (AAU), and we use two segmentation networks based on DeepLabv3 to locate the answer areas. Then, we use the characters recognition part to recognize students' answers. Finally, the grading part is designed for the comparison between the recognized answers and the standard ones. In our test, BAGS correctly locates and recognizes the handwritten answers in 91 of total answer areas.


Internet of Underwater Things and Big Marine Data Analytics A Comprehensive Survey. The Internet of Underwater Things (IoUT) is an emerging communication ecosystem developed for connecting underwater objects in maritime and underwater environments. The IoUT technology is intricately linked with intelligent boats and ships, smart shores and oceans, automatic marine transportations, positioning and navigation, underwater exploration, disaster prediction and prevention, as well as with intelligent monitoring and security. The IoUT has an influence at various scales ranging from a small scientific observatory, to a midsized harbor, and to covering global oceanic trade. The network architecture of IoUT is intrinsically heterogeneous and should be sufficiently resilient to operate in harsh environments. This creates major challenges in terms of underwater communications, whilst relying on limited energy resources. Additionally, the volume, velocity, and variety of data produced by sensors, hydrophones, and cameras in IoUT is enormous, giving rise to the concept of Big Marine Data (BMD), which has its own processing challenges. Hence, conventional data processing techniques will falter, and bespoke Machine Learning (ML) solutions have to be employed for automatically learning the specific BMD behavior and features facilitating knowledge extraction and decision support. The motivation of this paper is to comprehensively survey the IoUT, BMD, and their synthesis. It also aims for exploring the nexus of BMD with ML. We set out from underwater data collection and then discuss the family of IoUT data communication techniques with an emphasis on the state of the art research challenges. We then review the suite of ML solutions suitable for BMD handling and analytics. We treat the subject deductively from an educational perspective, critically appraising the material surveyed.


Maximum sparse induced subgraphs of the binomial random graph with given number of edges. We prove that a.a.s. the maximum size of an induced subtree of the binomial random graph G(n,p) is concentrated in 2 consecutive points. We also prove that, given a non negative integer valued function t(k) varepsilon k 2 , under a certain smoothness condition on this function, a.a.s. the maximum size k of an induced subgraph with exactly t(k) edges of G(n,p) is concentrated in 2 consecutive points as well.


High current ion beam formation with strongly inhomogeneous electrostatic field. A new approach to the development of extraction systems capable of forming ion beams with previously inaccessible intensity is proposed. The use of inhomogeneous accelerating field allows to improve the ion beam formation efficiency significantly. The increase of electric field magnitude is achieved by changing the shape of the electrodes only, without increasing the accelerating voltage and decreasing the interelectrode distance. The comparison is made between a new extraction system and a flat traditional one, which is the most common. The use of a new electrode geometry allows to increase the lifetime of the electrodes in sources of intense beams operating in a continuous wave mode. For electron cyclotron resonance ion sources, results demonstrate the possibility to form high quality ion beams with a current density of more than 1 A cm 2 .


Multi IF An Approach to Anomaly Detection in Self Driving Systems. Autonomous driving vehicles (ADVs) are implemented with rich software functions and equipped with many sensors, which in turn brings broad attack surface. Moreover, the execution environment of ADVs is often open and complex. Hence, ADVs are always at risk of safety and security threats. This paper proposes a fast method called Multi IF, using multiple invocation features of system calls to detect anomalies in self driving systems. Since self driving functions take most of the computation resources and upgrade frequently, Multi IF is designed to work under such resource constraints and support frequent updates. Given the collected sequences of system calls, the combination of different syntax patterns is used to analyze and construct feature vectors of those sequences. By taking the feature vectors as inputs, one class support vector machine is adopted to determine whether the current sequence of system calls is abnormal, which is trained with the feature vectors from the normal sequences. The evaluations on both simulated and real data prove that the proposed method is effective in identifying the abnormal behavior after minutes of feature extraction and training. Further comparisons with the existing methods on the ADFA LD data set also validate that the proposed approach achieves a higher accuracy with less time overhead.


Measurement of magnetic fields using the voltage generated by a vibrating wire. A vibrating wire may be used as an instrument with a variety of applications, one of which is the measurement of magnetic fields. Often, the magnetic fields are determined by measuring the amplitude of the wire vibration under the action of a Lorentz force. Though generally adequate, this approach may be inconvenient in certain circumstances. One of these occurs when it is necessary to measure the amplitude of high frequency vibration, as the amplitude is expected to decrease linearly with frequency, and thus becomes harder to measure. Another example may be found in situations where the sensor must operate over a wide range of vibration frequencies. In this case the sensor will be unresponsive to specific frequencies of wire vibration, which are determined by the placement of the sensor. This means that for the instrument to be robust, the sensor must be precisely mobile, or multiple sensors must be used. Here a technique which may be used to supplement the displacement sensor is described. This technique makes use of the voltage generated by the motion of the wire in the magnetic field under measurement. It is predicted that the technique may be more suitable for measurements requiring high frequency vibration, and is sensitive to all frequencies of vibration. Measurements of a magnetic field obtained using this technique are compared to those found using only a displacement sensor, and the benefits and drawbacks of the technique are discussed.


Named entity recognition in chemical patents using ensemble of contextual language models. Chemical patent documents describe a broad range of applications holding key reaction and compound information, such as chemical structure, reaction formulas, and molecular properties. These informational entities should be first identified in text passages to be utilized in downstream tasks. Text mining provides means to extract relevant information from chemical patents through information extraction techniques. As part of the Information Extraction task of the Cheminformatics Elsevier Melbourne University challenge, in this work we study the effectiveness of contextualized language models to extract reaction information in chemical patents. We assess transformer architectures trained on a generic and specialised corpora to propose a new ensemble model. Our best model, based on a majority ensemble approach, achieves an exact F1 score of 92.30 and a relaxed F1 score of 96.24 . The results show that ensemble of contextualized language models can provide an effective method to extract information from chemical patents.


A Dual Path Model With Adaptive Attention For Vehicle Re Identification. In recent years, attention models have been extensively used for person and vehicle re identification. Most re identification methods are designed to focus attention on key point locations. However, depending on the orientation, the contribution of each key point varies. In this paper, we present a novel dual path adaptive attention model for vehicle re identification (AAVER). The global appearance path captures macroscopic vehicle features while the orientation conditioned part appearance path learns to capture localized discriminative features by focusing attention on the most informative key points. Through extensive experimentation, we show that the proposed AAVER method is able to accurately re identify vehicles in unconstrained scenarios, yielding state of the art results on the challenging dataset VeRi 776. As a byproduct, the proposed system is also able to accurately predict vehicle key points and shows an improvement of more than 7 over state of the art. The code for key point estimation model is available at https github.com Pirazh Vehicle Key Point Orientation Estimation.


O( log log n) Worst Case Local Decoding and Update Efficiency for Data Compression. This paper addresses the problem of data compression with local decoding and local update. A compression scheme has worst case local decoding d wc if any bit of the raw file can be recovered by probing at most d wc bits of the compressed sequence, and has update efficiency of u wc if a single bit of the raw file can be updated by modifying at most u wc bits of the compressed sequence. This article provides an entropy achieving compression scheme for memoryless sources that simultaneously achieves O( log log n) local decoding and update efficiency. Key to this achievability result is a novel succinct data structure for sparse sequences which allows efficient local decoding and local update. Under general assumptions on the local decoder and update algorithms, a converse result shows that d wc and u wc must grow as Omega( log log n) .


New methodology for facilitating the food wastage quantification. Identifying gaps and data inconsistencies. The work aims at providing a new methodology to facilitate the process of quantifying the food waste according to European standards all along the agrifood chain combining information that is becoming available at local level. This new methodology generates straightforward and easy to interpret results for the decision making process in the framework of the quantification of the food waste at local and supralocal scale and it provides adequate procedures which are easy adaptable to the specific circumstances in each municipality. Moreover, this method could have applications for larger territorial contexts, as the national scale, detecting possible points for improvement of the current official figures at this respect.


Lexicographic Groebner bases of bivariate polynomials modulo a univariate one. Let T(x) in k x be a monic non constant polynomial and write R k x (T) the quotient ring. Consider two bivariate polynomials a(x, y), b(x, y) in R y . In a first part, T p e is assumed to be the power of an irreducible polynomial p. A new algorithm that computes a minimal lexicographic Groebner basis of the ideal ( a, b, p e), is introduced. A second part extends this algorithm when T is general through the "local global" principle realized by a generalization of "dynamic evaluation", restricted so far to a polynomial T that is squarefree. The algorithm produces splittings according to the case distinction "invertible nilpotent", extending the usual "invertible zero" in classic dynamic evaluation. This algorithm belongs to the Euclidean family, the core being a subresultant sequence of a and b modulo T. In particular no factorization or Groebner basis computations are necessary. The theoretical background relies on Lazard's structural theorem for lexicographic Groebner bases in two variables. An implementation is realized in Magma. Benchmarks show clearly the benefit, sometimes important, of this approach compared to the Groebner bases approach.


Transition Frequencies and Hyperfine Structure in 113,115 In Application of a Liquid Metal Ion Source for Collinear Laser Spectroscopy. We demonstrate the first application of a liquid metal ion source for collinear laser spectroscopy in proof of principle measurements on naturally abundant In . The superior beam quality, i.e., the actively stabilized current and energy of a beam with very low transverse emittance, allowed us to perform precision spectroscopy on the 5s 2 1 mathrm S 0 rightarrow 5s5p 3 mathrm P 1 intercombination transition in 115 In , which is to our knowledge the slowest transition measured with collinear fluorescence laser spectroscopy so far. By applying collinear and anticollinear spectroscopy, we improved the center of gravity frequency nu mathrm cg 1 ,299 ,617 ,759. ,3 ,(1.2) and the hyperfine constants A 6957.19 ,(28) ,MHz and B 443.7 ,(2.4) ,MHz by more than two orders of magnitude. A similar accuracy was reached for 113 In in combination with literature data and the isotope shift between both naturally abundant isotopes was deduced to nu( 113 mathrm In ) nu( 115 mathrm In ) 696.3 ,(3.1) ,MHz. Nuclear alignment induced by optical pumping in a preparation section of the ion beamline was demonstrated as a pump and probe approach to provide sharp features on top of the Doppler broadened resonance profile.


Endpoint estimates for the maximal function over prime numbers. Given an ergodic dynamical system (X, mathcal B , mu, T) , we prove that for each function f belonging to the Orlicz space L( log L) 2( log log L)(X, mu) , the ergodic averages frac 1 pi(N) sum p in mathbb P N f big(T p x big), converge for mu almost all x in X , where mathbb P N is the set of prime numbers not larger that N and pi(N) mathbb P N .


PATSQL Efficient Synthesis of SQL Queries from Example Tables with Quick Inference of Projected Columns. SQL is one of the most popular tools for data analysis, and it is now used by an increasing number of users without having expertise in databases. Several studies have proposed programming by example approaches to help such non experts to write correct SQL queries. While existing methods support a variety of SQL features such as aggregation and nested query, they suffer a significant increase in computational cost as the scale of example tables increases. In this paper, we propose an efficient algorithm utilizing properties known in relational algebra to synthesize SQL queries from input and output tables. Our key insight is that a projection operator in a program sketch can be lifted above other operators by applying transformation rules in relational algebra, while preserving the semantics of the program. This enables a quick inference of appropriate columns in the projection operator, which is an essential component in synthesis but causes combinatorial explosions in prior work. We also introduce a novel form of constraints and its top down propagation mechanism for efficient sketch completion. We implemented this algorithm in our tool PATSQL and evaluated it on 226 queries from prior benchmarks and Kaggle's tutorials. As a result, PATSQL solved 68 of the benchmarks and found 89 of the solutions within a second. Our tool is available at https naist se.github.io patsql .


Multicommodity Flows in Planar Graphs with Demands on Faces. We consider the problem of multicommodity flows in planar graphs. Seymour showed that if the union of supply and demand graphs is planar, then the cut condition is sufficient for routing demands. Okamura Seymour showed that if all demands are incident on one face, then again cut condition is sufficient for routing demands. We consider a common generalization of these settings where the end points of each demand are on the same face of the planar graph. We show that if the source sink pairs on each face of the graph are such that sources and sinks appear contiguously on the cycle bounding the face, then the flow cut gap is at most 3. We come up with a notion of approximating demands on a face by convex combination of laminar demands to prove this result.


A Summary of Formal Specification and Verification of Autonomous Robotic Systems. Autonomous robotic systems are complex, hybrid, and often safety critical this makes their formal specification and verification uniquely challenging. Though commonly used, testing and simulation alone are insufficient to ensure the correctness of, or provide sufficient evidence for the certification of, autonomous robotics. Formal methods for autonomous robotics have received some attention in the literature, but no resource provides a current overview. This short paper summarises the contributions of Luckcuck 2019, which surveys the state of the art in formal specification and verification for autonomous robotics.


Coupled Cluster Externally Corrected by Adaptive Configuration Interaction. An externally corrected coupled cluster (CC) method, where an adaptive configuration interaction (ACI) wave function provides the external cluster amplitudes, named ACI CC, is presented. By exploiting the connection between configuration interaction and coupled cluster through cluster analysis, the higher order T3 and T4 terms obtained from ACI are used to augment the T1 and T2 amplitude equations from traditional coupled cluster. These higher order contributions are kept frozen during the coupled cluster iterations and do not contribute to an increased cost with respect to CCSD. We have benchmarked this method on three closed shell systems beryllium dimer, carbonyl oxide, and cyclobutadiene, with good results compared to other corrected coupled cluster methods. In all cases, the inclusion of these external corrections improved upon the "gold standard" CCSD(T) results, indicating that ACI CCSD(T) can be used to assess strong correlation effects in a system and as an inexpensive starting point for more complex external corrections.


The generalised principal eigenvalue of time periodic nonlocal dispersal operators and applications. This paper is mainly concerned with the generalised principal eigenvalue for time periodic nonlocal dispersal operators. We first establish the equivalence between two different characterisations of the generalised principal eigenvalue. We further investigate the dependence of the generalised principal eigenvalue on the frequency, the dispersal rate and the dispersal spread. Finally, these qualitative results for time periodic linear operators are applied to time periodic nonlinear KPP equations with nonlocal dispersal, focusing on the effects of the frequency, the dispersal rate and the dispersal spread on the existence and stability of positive time periodic solutions to nonlinear equations.


The relative trace formula in analytic number theory. We discuss a variety of applications of the relative trace formula of Kuznetsov type in analytic number theory and the theory of automorphic forms.


Spectral properties of anti heptadiagonal persymmetric Hankel matrices. In this paper we express the eigenvalues of anti heptadiagonal persymmetric Hankel matrices as the zeros of explicit polynomials giving also a representation of its eigenvectors. We present also an expression depending on localizable parameters to compute its integer powers. In particular, an explicit formula not depending on any unknown parameter for the inverse of anti heptadiagonal persymmetric Hankel matrices is provided.


On simultaneous rational approximation to a real number and its integral powers, II. For a positive integer n and a real number xi , let lambda n ( xi) denote the supremum of the real numbers lambda for which there are arbitrarily large positive integers q such that q xi , q xi 2 , ldots , q xi n are all less than q lambda . Here, cdot denotes the distance to the nearest integer. We establish new results on the Hausdorff dimension of the set of real numbers xi such that lambda n ( xi) is equal (or greater than or equal) to a given value.


Nota T 'ecnica dos Modelos Implementados pelo Coletivo Covid19br para Proje c c oes de Cen 'arios Futuros da Pandemia COVID 19 no Brasil. This technical note aims to provide a brief introduction to the projection models used by the group to project future scenarios for states and municipalities in real time, according to the disease's behavior in previous days. However, the parameters can be modified by the user to design customized scenarios. The proposed model begins with the calculation of the basic reproduction number for the state or municipality based on the incidence of cases in the last 12 days. Once this is done, the epidemiological curve is projected using the SEIR compartmentalized epidemic model, in possession of this curve, part of the newly projected infected enter a simulation model for health systems in queuing theory, aiming to project future occupations and collapses.


Dimensions of the popcorn graph. The 'popcorn function' isThe popcorn function' is a well known and important example in real analysis with many interesting features. We prove that the box dimension of the graph of the popcorn function is 4 3, as well as computing the Assouad dimension and Assouad spectrum. The main ingredients include Duffin Schaeffer type estimates from Diophantine approximation and the Chung Erd H o s inequality from probability theory.


Centrally Essential Endomorphism Rings of Abelian Groups. We study Abelian groups A with centrally essential endomorphism ring text End ,A . If A is a such group which is either a torsion group or a non reduced group, then the ring text End ,A is commutative. We give examples of Abelian torsion free groups of finite rank with non commutative centrally essential endomorphism rings. The paper will appear in Communications in Algebra.


Multidimensional topological Galois theory. In this preprint we present an outline of the multidimensional version of topological Galois theory. The theory studies topological obstruction to solvability of equations "in finite terms" (i.e. to their solvability by radicals, by elementary functions, by quadratures and so on). This preprint is based on the author's book on topological Galois theory. It contains definitions, statements of results and comments to them. Basically no proofs are presented. This preprint was written as a part of the comments to a new edition (in preparation) of the classical book "Integration in finite terms" by J.F. Ritt.


From Welcome New Gabbers to the Pittsburgh Synagogue Shooting The Evolution of Gab. Gab, an online social media platform with very little content moderation, has recently come to prominence as an alt right community and a haven for hate speech. We document the evolution of Gab since its inception until a Gab user carried out the most deadly attack on the Jewish community in US history. We investigate Gab language use, study how topics evolved over time, and find that the shooters' posts were among the most consistently anti Semitic on Gab, but that hundreds of other users were even more extreme.


Herriott cavity assisted all optical atomic vector magnetometer. We report an all optical atomic vector magnetometer using dual Bell Bloom optical pumping beams in a Rb vapor cell. This vector magnetometer consists of two orthogonal optical pumping beams, with amplitude modulations at 85 Rb and 87 Rb Larmor frequencies respectively. We simultaneously detect atomic signals excited by these two pumping beams using a single probe beam in the third direction, and extract the field orientation information using the phase delays between the modulated atomic signals and the driving beams. By adding a Herriott cavity inside the vapor cell, we improve the magnetometer sensitivity. We study the performance of this vector magnetometer in a magnetic field ranging from 100 mG to 500 mG, and demonstrate a field angle sensitivity better than 10 mu rad Hz 1 2 above 10 Hz.


Unexpectedly strong diamagnetism and superparamagnetism of aromatic peptides due to self assembling and cations. There is a considerable amount of work that shows the biomagnetism of organic components without ferromagnetic components at the molecular level, but it is of great challenge to cover the giant gap of biomagnetism between their experimental and theoretical results. Here, we show that the diamagnetism of an aromatic peptide, the AYFFF, is greatly enhanced for about 11 times by self assembling, reaching two orders of magnitude higher than the mass susceptibility of pure water. Moreover, the AYFFF self assemblies further mixed with ZnCl2 solution of sufficiently high concentrations display superparamagnetism, with the mass susceptibility reaching more than two orders of magnitude higher than the absolute value of pure water, which may approach the mass susceptibility of ferromagnetism. The aromatic rings in the peptide molecules and the cations are the keys to such a strong diamagnetism and superparamagnetism of aromatic peptides.


Ranking the spreading influence of nodes in complex networks based on mixing degree centrality and local structure. The safety and robustness of the network have attracted the attention of people from all walks of life, and the damage of several key nodes will lead to extremely serious consequences. In this paper, we proposed the clustering H index mixing (CHM) centrality based on the H index of the node itself and the relative distance of its neighbors. Starting from the node itself and combining with the topology around the node, the importance of the node and its spreading capability were determined. In order to evaluate the performance of the proposed method, we use Susceptible Infected Recovered (SIR) model, monotonicity and resolution as the evaluation standard of experiment. Experimental results in artificial networks and real world networks show that CHM centrality has excellent performance in identifying node importance and its spreading capability.


Tree limits and limits of random trees. We explore the tree limits recently defined by Elek and Tardos. In particular, we find tree limits for many classes of random trees. We give general theorems for three classes of conditional Galton Watson trees and simply generated trees, for split trees and generalized split trees (as defined here), and for trees defined by a continuous time branching process. These general results include, for example, random labelled trees, ordered trees, random recursive trees, preferential attachment trees, and binary search trees.


physics.class-ph cs.NA math.DS math.NA math.OC

Finding the strongest stable weightless column with a follower load and relocatable concentrated masses. We consider the problem of optimal placement of concentrated masses along a massless elastic column that is clamped at one end and loaded by a nonconservative follower force at the free end. The goal is to find the largest possible interval such that the variation in the loading parameter within this interval preserves stability of the structure. The stability constraint is nonconvex and nonsmooth, making the optimization problem quite challenging. We give a detailed analytical treatment for the case of two masses, arguing that the optimal parameter configuration approaches the flutter and divergence boundaries of the stability region simultaneously. Furthermore, we conjecture that this property holds for any number of masses, which in turn suggests a simple formula for the maximal load interval for n masses. This conjecture is strongly supported by extensive computational results, obtained using the recently developed open source software package GRANSO (GRadient based Algorithm for Non Smooth Optimization) to maximize the load interval subject to an appropriate formulation of the nonsmooth stability constraint. We hope that our work will provide a foundation for new approaches to classical long standing problems of stability optimization for nonconservative elastic systems arising in civil and mechanical engineering.


Construction of Const Dimension Codes from Serval Parallel Lift MRD Code. In this paper, we generalize the method of using two parallel versions of the lifted MRD code from the existing work 1 . The Delsarte theorem of the rank distribution of MRD codes is an important part to count codewords in our construction. We give a new generalize construction to the following bounds if n k d, then Aq(n k,k,d) q n(k frac d 2 1) sum r frac d 2 k frac d 2 A r(Q q(n,k, frac d 2 )). On this basis, we also give a construction of constant dimension subspace codes from several parallel versions of lifted MRD codes. This construction contributes to a new lower bounds for Aq((s 1)k n,d,k).


Construction of local A packets. In this paper, we reformulate Moeglin's explicit construction of local A packets of split odd special orthogonal groups and symplectic groups. By this reformulation together with results of the previous paper with Minguez, we can compute the A packets explicitly. Also, we give a non vanishing criterion of our parametrization, and an algorithm to compute certain derivatives. Finally, we prove a formula for the Aubert duals of irreducible representations of Arthur type.


Top Fourier coefficients of residual Eisenstein series on symplectic or metaplectic groups induced from Speh representations. We consider the residues at the poles in the right half plane of Eisenstein series, on symplectic groups, or their double covers, induced from Speh representations. We show that for each such pole, there is a unique maximal nilpotent orbit, attached to Fourier coefficients admitted by the corresponding residual representation. We find this orbit in each case.


Microalgae based bioremediation of water contaminated by pesticides in peri urban agricultural areas. The present study evaluated the capacity of a semi closed, tubular horizontal photobioreactor (PBR) to remove pesticides from agricultural run off. The study was carried out in July to study its efficiency under the best conditions (highest solar irradiation). A total of 51 pesticides, including 10 transformation products, were selected and investigated based on their consumption rate and environmental relevance. Sixteen of them were detected in the agricultural run off, and the estimated removal efficiencies ranged from negative values, obtained for 3 compounds, namely terbutryn, diuron, and imidacloprid, to 100 , achieved for 10 compounds. The acidic herbicide MCPA was removed by 88 on average, and the insecticides 2,4 D and diazinon showed variable removals, between 100 and negative values. The environmental risk associated with the compounds still present in the effluent of the PBR was evaluated using hazard quotients (HQs), calculated using the average and highest measured concentrations of the compounds. HQ values 10 (meaning high risk) were obtained for imidacloprid (21), between 1 and 10 (meaning moderate risk) for 2,4 D (2.8), diazinon (4.6) and terbutryn (1.5), and 1 (meaning low risk) for the remaining compounds diuron, linuron, and MCPA. The PBR treatment yielded variable removals depending on the compound, similar to conventional wastewater treatment plants. This study provides new data on the capacity of icroalgae based treatment systems to eliminate a wide range of priority pesticides under real environmental conditions.


Stein Variational Model Predictive Control. Decision making under uncertainty is critical to real world, autonomous systems. Model Predictive Control (MPC) methods have demonstrated favorable performance in practice, but remain limited when dealing with complex probability distributions. In this paper, we propose a generalization of MPC that represents a multitude of solutions as posterior distributions. By casting MPC as a Bayesian inference problem, we employ variational methods for posterior computation, naturally encoding the complexity and multi modality of the decision making problem. We present a Stein variational gradient descent method to estimate the posterior directly over control parameters, given a cost function and observed state trajectories. We show that this framework leads to successful planning in challenging, non convex optimal control problems.


Three dimensional tracking of finite size spheres in a turbulent boundary layer. The motion of individual magnetic wax spheres with specific gravities of 1.003, 1.050 and 1.150 was investigated in turbulent boundary layers with Re tau 700 and 1300 ( d 60 and 120). The spheres were marked with dots all over the surface to monitor their translation and rotation via high speed stereoscopic imaging. Upon release from rest on a smooth wall, each sphere typically accelerated strongly over a streamwise distance of one boundary layer thickness before approaching an approximate terminal velocity. Spheres with sufficient net upward force lifted off of the wall once released before descending back towards the wall. These spheres mostly translated with the fluid above the wall, undergoing saltation or resuspension, with minimal rotation about all axes. By contrast, spheres that did not lift off after release mainly slid along the wall. As they propagated downstream, they began to roll forward with occasional lift off events of smaller magnitude. All of the lift off activities observed were limited to the buffer and logarithmic layers. Both translation and rotation of the spheres were significantly affected by the wall turbulence.


Global F splitting of surfaces admitting an int amplified endomorphim. In this paper, we study the global F splitting of varieties admitting an int amplified endomoprhism. We prove that surfaces admitting an int amplified endomorphism are of dense globally F split type and, in particular, of Calabi Yau type.


Behavioural Equivalence via Modalities for Algebraic Effects. The paper investigates behavioural equivalence between programs in a call by value functional language extended with a signature of (algebraic) effect triggering operations. Two programs are considered as being behaviourally equivalent if they enjoy the same behavioural properties. To formulate this, we define a logic whose formulas specify behavioural properties. A crucial ingredient is a collection of modalities expressing effect specific aspects of behaviour. We give a general theory of such modalities. If two conditions, openness and decomposability, are satisfied by the modalities then the logically specified behavioural equivalence coincides with a modality defined notion of applicative bisimilarity, which can be proven to be a congruence by a generalisation of Howe's method. We show that the openness and decomposability conditions hold for several examples of algebraic effects nondeterminism, probabilistic choice, global store and input output.


Compact reduction in Lipschitz free spaces. We prove a general principle satisfied by weakly precompact sets of Lipschitz free spaces. By this principle, certain infinite dimensional phenomena in Lipschitz free spaces over general metric spaces may be reduced to the same phenomena in free spaces over their compact subsets. As easy consequences we derive several new and some known results. The main new results are mathcal F(X) is weakly sequentially complete for every superreflexive Banach space X , and mathcal F(M) has the Schur property and the approximation property for every scattered complete metric space M .


Privacy, Altruism, and Experience Estimating the Perceived Value of Internet Data for Medical Uses. People increasingly turn to the Internet when they have a medical condition. The data they create during this process is a valuable source for medical research and for future health services. However, utilizing these data could come at a cost to user privacy. Thus, it is important to balance the perceived value that users assign to these data with the value of the services derived from them. Here we describe experiments where methods from Mechanism Design were used to elicit a truthful valuation from users for their Internet data and for services to screen people for medical conditions. In these experiments, 880 people from around the world were asked to participate in an auction to provide their data for uses differing in their contribution to the participant, to society, and in the disease they addressed. Some users were offered monetary compensation for their participation, while others were asked to pay to participate. Our findings show that 99 of people were willing to contribute their data in exchange for monetary compensation and an analysis of their data, while 53 were willing to pay to have their data analyzed. The average perceived value users assigned to their data was estimated at US 49. Their value to screen them for a specific cancer was US 22 while the value of this service offered to the general public was US 22. Participants requested higher compensation when notified that their data would be used to analyze a more severe condition. They were willing to pay more to have their data analyzed when the condition was more severe, when they had higher education or if they had recently experienced a serious medical condition.


ALTER Auxiliary Text Rewriting Tool for Natural Language Generation. In this paper, we describe ALTER, an auxiliary text rewriting tool that facilitates the rewriting process for natural language generation tasks, such as paraphrasing, text simplification, fairness aware text rewriting, and text style transfer. Our tool is characterized by two features, i) recording of word level revision histories and ii) flexible auxiliary edit support and feedback to annotators. The text rewriting assist and traceable rewriting history are potentially beneficial to the future research of natural language generation.


On Frobenius Perron Dimension. We propose a notion of Frobenius Perron dimension for certain free mathbb Z modules of infinite rank and compute it for the mathbb Z modules of finite dimensional complex representations of unitary groups with nonnegative dominant weights. We also provide a lower bound for the Frobenius Perron dimension of Schubert classes in the quantum cohomology of complex Grassmannians.


On products of primes and square free integers in arithmetic progressions. We obtain an asymptotic formula for the number of ways to represent every reduced residue class as a product of a prime and square free integer. This may be considered as a relaxed version of a conjecture of Erd "os, Odlyzko, and S 'ark "ozy.


VarGNet Variable Group Convolutional Neural Network for Efficient Embedded Computing. In this paper, we propose a novel network design mechanism for efficient embedded computing. Inspired by the limited computing patterns, we propose to fix the number of channels in a group convolution, instead of the existing practice that fixing the total group numbers. Our solution based network, named Variable Group Convolutional Network (VarGNet), can be optimized easier on hardware side, due to the more unified computing schemes among the layers. Extensive experiments on various vision tasks, including classification, detection, pixel wise parsing and face recognition, have demonstrated the practical value of our VarGNet.


On the structure of weak solutions to scalar conservation laws with finite entropy production. We consider weak solutions with finite entropy production to the scalar conservation law begin equation partial t u mathrm div x F(u) 0 quad mbox in (0,T) times mathbb R d. end equation Building on the kinetic formulation we prove under suitable nonlinearity assumption on f that the set of non Lebesgue points of u has Hausdorff dimension at most d . A notion of Lagrangian representation for this class of solutions is introduced and this allows for a new interpretation of the entropy dissipation measure.


Mapping in a cycle Sinkhorn regularized unsupervised learning for point cloud shapes. We propose an unsupervised learning framework with the pretext task of finding dense correspondences between point cloud shapes from the same category based on the cycle consistency formulation. In order to learn discriminative pointwise features from point cloud data, we incorporate in the formulation a regularization term based on Sinkhorn normalization to enhance the learned pointwise mappings to be as bijective as possible. Besides, a random rigid transform of the source shape is introduced to form a triplet cycle to improve the model's robustness against perturbations. Comprehensive experiments demonstrate that the learned pointwise features through our framework benefits various point cloud analysis tasks, e.g. partial shape registration and keypoint transfer. We also show that the learned pointwise features can be leveraged by supervised methods to improve the part segmentation performance with either the full training dataset or just a small portion of it.


Irreducibility of x n a . A. Capelli gave a necessary and sufficient condition for the reducibility of x n a over Q . In this article, we are providing an alternate elementary proof for the same.


Packing of mixed hyperarborescences with flexible roots via matroid intersection. Given a mixed hypergraph mathcal F (V, mathcal A cup mathcal E ) , functions f,g V rightarrow mathbb Z and an integer k , a packing of k spanning mixed hyperarborescences is called (k,f,g) flexible if every v in V is the root of at least f(v) and at most g(v) of the mixed hyperarborescences. We give a characterization of the mixed hypergraphs admitting such packings. This generalizes results of Frank and, more recently, Gao and Yang. Our approach is based on matroid intersection, generalizing a construction of Edmonds. We also obtain an algorithm for finding a minimum weight solution to the above mentioned problem.


Can Unconditional Language Models Recover Arbitrary Sentences?. Neural network based generative language models like ELMo and BERT can work effectively as general purpose sentence encoders in text classification without further fine tuning. Is it possible to adapt them in a similar way for use as general purpose decoders? For this to be possible, it would need to be the case that for any target sentence of interest, there is some continuous representation that can be passed to the language model to cause it to reproduce that sentence. We set aside the difficult problem of designing an encoder that can produce such representations and, instead, ask directly whether such representations exist at all. To do this, we introduce a pair of effective, complementary methods for feeding representations into pretrained unconditional language models and a corresponding set of methods to map sentences into and out of this representation space, the reparametrized sentence space. We then investigate the conditions under which a language model can be made to generate a sentence through the identification of a point in such a space and find that it is possible to recover arbitrary sentences nearly perfectly with language models and representations of moderate size without modifying any model parameters.


Bootstrap Aggregation for Point based Generalized Membership Inference Attacks. An efficient scheme is introduced that extends the generalized membership inference attack to every point in a model's training data set. Our approach leverages data partitioning to create variable sized training sets for the reference models. We then train an attack model for every single training example for a reference model configuration based upon output for each individual point. This allows us to quantify the membership inference attack vulnerability of each training data point. Using this approach, we discovered that smaller amounts of reference model training data led to a stronger attack. Furthermore, the reference models do not need to be of the same architecture as the target model, providing additional attack efficiencies. The attack may also be performed by an adversary even when they do not have the complete original data set.


Dedekind sums arising from newform Eisenstein series. For primitive non trivial Dirichlet characters chi 1 and chi 2 , we study the weight zero newform Eisenstein series E chi 1, chi 2 (z,s) at s 1 . The holomorphic part of this function has a transformation rule that we express in finite terms as a generalized Dedekind sum. This gives rise to the explicit construction (in finite terms) of elements of H 1( Gamma 0(N), mathbb C ) . We also give a short proof of the reciprocity formula for this Dedekind sum.


Managing aquatic invasions optimal locations and operating times for watercraft inspection stations. Aquatic invasive species (AIS) cause significant ecological and economic damages around the world. A major spread mechanism for AIS is traffic of boaters transporting their watercraft from invaded to uninvaded waterbodies. To inhibit the spread of AIS, several Canadian provinces and American states set up watercraft inspection stations at roadsides, where potentially infested boats are screened for AIS and, if necessary, decontaminated. However, since budgets for AIS control are limited, watercraft inspection stations can only be operated at specific locations and daytimes. Though theoretical studies provide managers with general guidelines for AIS management, more specific results are needed to determine when and where watercraft inspections would be most effective. This is the subject of this paper. We show how linear integer programming techniques can be used to optimize watercraft inspection policies under budget constraints. We introduce our approach as a general framework and apply it to the prevention of the spread of zebra and quagga mussels (Dreissena spp.) to the Canadian province British Columbia. We consider a variety of scenarios and show how variations in budget constraints, propagule sources, and model uncertainty affect the optimal policy. Based on these results, we identify simple, generally applicable principles for optimal AIS management.


Dense 3D Face Decoding over 2500FPS Joint Texture Shape Convolutional Mesh Decoders. 3D Morphable Models (3DMMs) are statistical models that represent facial texture and shape variations using a set of linear bases and more particular Principal Component Analysis (PCA). 3DMMs were used as statistical priors for reconstructing 3D faces from images by solving non linear least square optimization problems. Recently, 3DMMs were used as generative models for training non linear mappings ( ie, regressors) from image to the parameters of the models via Deep Convolutional Neural Networks (DCNNs). Nevertheless, all of the above methods use either fully connected layers or 2D convolutions on parametric unwrapped UV spaces leading to large networks with many parameters. In this paper, we present the first, to the best of our knowledge, non linear 3DMMs by learning joint texture and shape auto encoders using direct mesh convolutions. We demonstrate how these auto encoders can be used to train very light weight models that perform Coloured Mesh Decoding (CMD) in the wild at a speed of over 2500 FPS.


Manifolds of Triangulations, braid groups of manifolds and the groups Gamma n k . The spaces of triangulations of a given manifold have been widely studied. The celebrated theorem of Pachner cite Pachner says that any two triangulations of a given manifold can be connected by a sequence of bistellar moves, or Pachner moves, see also cite GKZ,Nabutovsky . In the present paper we consider groups which naturally appear when considering the set of triangulations with fixed number of simplices of maximal dimension. There are three ways of introducing this groups the geometrical one, which depends on the metric, the topological one, and the combinatorial one. The second one can be thought of as a braid group'' of the manifold and, by definition, is an invariant of the topological type of manifold in a similar way, one can construct the smooth version. We construct a series of groups Gamma n k corresponding to Pachner moves of (k 2) dimensional manifolds and construct a canonical map from the braid group of any k dimensional manifold to Gamma n k thus getting topological smooth invariants of these manifolds.


Genus one cobordisms between torus knots. We determine the pairs of torus knots that have a genus one cobordism between them, with one notable exception. This is done by combining obstructions using nu from the Heegaard Floer knot complex and explicit constructions of cobordisms. As an application, we determine the pairs of torus knots related by a single crossing change. Also, we determine the pairs of Thurston Bennequin number maximizing Legendrian torus knots that have a genus one exact Lagrangian cobordism, with one exception.


L 2 times L 2 times L 2 to L 2 3 boundedness for trilinear multiplier operator. This paper discusses the boundedness of the trilinear multiplier operator T m (f 1,f 2,f 3) , when the multiplier satisfies a certain degree of smoothness but with no decaying condition and is L q integrable with an admissible range of q . The boundedness is stated in the terms of m L q . In particular, begin equation T m L 2 times L 2 times L 2 to L 2 3 lesssim m L q q 3 . end equation


A categorical duality for algebras of partial functions. We prove a categorical duality between a class of abstract algebras of partial functions and a class of (small) topological categories. The algebras are the isomorphs of collections of partial functions closed under the operations of composition, antidomain, range, and preferential union (or 'override'). The topological categories are those whose space of objects is a Stone space, source map is a local homeomorphism, target map is open, and all of whose arrows are epimorphisms.


Theta divisors and Ulrich bundles on Geometrically ruled surfaces. We consider the following question for which invariants g and e is there a geometrically ruled surface S rightarrow C over a curve C of genus g with invariant e such that S is the support of an Ulrich line bundle with respect to a very ample line bundle? A surprising relation between the existence of certain proper Theta divisors on some moduli spaces of vector bundles on C with the existence of Ulrich line bundles on S will be the key to completely solve the above question. The relation is realized by translating the vanishing conditions characterizing Ulrich line bundles to specific geometric conditions on the symmetric powers of the defining vector bundle of a given ruled surface. This general principle leads to some finer existence results of Ulrich line bundles in particular cases. Another focus is on the rank two case where, with very few exceptions, we show the existence of large families of special Ulrich bundles on arbitrary polarized ruled surfaces.


Tur 'an and Ramsey type results for unavoidable subgraphs. We study Tur 'an and Ramsey type problems on edge colored graphs. An edge colored graph is called em varepsilon balanced if each color class contains at least an varepsilon proportion of its edges. Given a family mathcal F of edge colored graphs, the Ramsey function R( varepsilon, mathcal F ) is the smallest n for which any varepsilon balanced K n must contain a copy of an F in mathcal F , and the Tur 'an function mathrm ex ( varepsilon, n, mathcal F ) is the maximum number of edges in an n vertex varepsilon balanced graph which avoids all of mathcal F . In this paper, we consider this Tur 'an function for several classes of edge colored graphs, we show that the Ramsey function is linear for bounded degree graphs, and we prove a theorem that gives a relationship between the two parameters.


Rigidity for perimeter inequality under spherical symmetrisation. Necessary and sufficient conditions for rigidity of the perimeter inequality under spherical symmetrisation are given. That is, a characterisation for the uniqueness (up to orthogonal transformations) of the extremals is provided. This is obtained through a careful analysis of the equality cases, and studying fine properties of the circular symmetrisation, which was firstly introduced by P 'olya in 1950.


Graph Spanners by Sketching in Dynamic Streams and the Simultaneous Communication Model. Graph sketching is a powerful technique introduced by the seminal work of Ahn, Guha and McGregor'12 on connectivity in dynamic graph streams that has enjoyed considerable attention in the literature since then, and has led to near optimal dynamic streaming algorithms for many fundamental problems such as connectivity, cut and spectral sparsifiers and matchings. Interestingly, however, the sketching and dynamic streaming complexity of approximating the shortest path metric of a graph is still far from well understood. Besides a direct k pass implementation of classical spanner constructions (recently improved to lfloor frac k2 rfloor 1 passes by Fernandez, Woodruff and Yasuda'20) the state of the art amounts to a O( log k) pass algorithm of Ahn, Guha and McGregor'12, and a 2 pass algorithm of Kapralov and Woodruff'14. In particular, no single pass algorithm is known, and the optimal tradeoff between the number of passes, stretch and space complexity is open. In this paper we introduce several new graph sketching techniques for approximating the shortest path metric of the input graph. We give the first em single pass sketching algorithm for constructing graph spanners we show how to obtain a widetilde O (n frac23 ) spanner using widetilde O (n) space, and in general a widetilde O (n frac23(1 alpha) ) spanner using widetilde O (n 1 alpha ) space for every alpha in 0, 1 , a tradeoff that we think may be close optimal. We also give new spanner construction algorithms for any number of passes, simultaneously improving upon all prior work on this problem. Finally, we study the simultaneous communication model and propose the first protocols with low per player information.


Optimal Randomized First Order Methods for Least Squares Problems. We provide an exact analysis of a class of randomized algorithms for solving overdetermined least squares problems. We consider first order methods, where the gradients are pre conditioned by an approximation of the Hessian, based on a subspace embedding of the data matrix. This class of algorithms encompasses several randomized methods among the fastest solvers for least squares problems. We focus on two classical embeddings, namely, Gaussian projections and subsampled randomized Hadamard transforms (SRHT). Our key technical innovation is the derivation of the limiting spectral density of SRHT embeddings. Leveraging this novel result, we derive the family of normalized orthogonal polynomials of the SRHT density and we find the optimal pre conditioned first order method along with its rate of convergence. Our analysis of Gaussian embeddings proceeds similarly, and leverages classical random matrix theory results. In particular, we show that for a given sketch size, SRHT embeddings exhibits a faster rate of convergence than Gaussian embeddings. Then, we propose a new algorithm by optimizing the computational complexity over the choice of the sketching dimension. To our knowledge, our resulting algorithm yields the best known complexity for solving least squares problems with no condition number dependence.


A framework for step wise explaining how to solve constraint satisfaction problems. We explore the problem of step wise explaining how to solve constraint satisfaction problems, with a use case on logic grid puzzles. More specifically, we study the problem of explaining the inference steps that one can take during propagation, in a way that is easy to interpret for a person. Thereby, we aim to give the constraint solver explainable agency, which can help in building trust in the solver by being able to understand and even learn from the explanations. The main challenge is that of finding a sequence of simple explanations, where each explanation should aim to be as cognitively easy as possible for a human to verify and understand. This contrasts with the arbitrary combination of facts and constraints that the solver may use when propagating. We propose the use of a cost function to quantify how simple an individual explanation of an inference step is, and identify the explanation production problem of finding the best sequence of explanations of a CSP. Our approach is agnostic of the underlying constraint propagation mechanisms, and can provide explanations even for inference steps resulting from combinations of constraints. In case multiple constraints are involved, we also develop a mechanism that allows to break the most difficult steps up and thus gives the user the ability to zoom in on specific parts of the explanation. Our proposed algorithm iteratively constructs the explanation sequence by using an optimistic estimate of the cost function to guide the search for the best explanation at each step. Our experiments on logic grid puzzles show the feasibility of the approach in terms of the quality of the individual explanations and the resulting explanation sequences obtained.


Compact quasi Einstein manifolds with boundary. The goal of this article is to study compact quasi Einstein manifolds with boundary. We provide boundary estimates for compact quasi Einstein manifolds simi lar to previous results obtained for static and V static spaces. In addition, we show that compact quasi Einstein manifolds with connected boundary and satisfying a suitable pinching condition must be isometric, up to scaling, to the standard hemisphere mathbb S n .


A Faithful Discretization of the Augmented Persistent Homology Transform. The persistent homology transform (PHT) represents a shape with a multiset of persistence diagrams parameterized by the sphere of directions in the ambient space. In this work, we describe a finite set of diagrams that discretize the PHT such that it faithfully represents the underlying shape. We provide a discretization that is exponential in the dimension of the shape (making it Furthermore, we provide an output sensitive algorithm that is, the algorithm reports the discretization in time proportional to the size of the discretization. Finally, our approach relies only on knowing the heights and dimensions of topological events, meaning that it can be adapted to provide discretizations of other dimension returning topological transforms, including the Betti curve transform.


Entropy of Toeplitz systems over residually finite groups. The purpose of this work is to bound sofic topological entropy of Toeplitz systems over residually finite groups and to prove the Krieger Theorem about attaining arbitrary entropy by the Toeplitz systems. To achieve these results, we discuss certain properties of the sofic topological entropy in the context of finitely indexed normal subgroups of the group. It will help us to formulate results almost independently of the natural sofic approximation sequences of residually finite groups.


Exponential Convergence Rate for the Asymptotic Optimality of Whittle Index Policy. We evaluate the performance of Whittle index policy for restless Markovian bandits, when the number of bandits grows. It is proven in 30 that this performance is asymptotically optimal if the bandits are indexable and the associated deterministic system has a global attractor fixed point. In this paper we show that, under the same conditions, the convergence rate is exponential in the number of bandits, unless the fixed point is singular (to be defined later). Our proof is based on the nature of the deterministic equation governing the stochastic system We show that it is a piecewise affine continuous dynamical system inside the simplex of the empirical measure of the bandits. Using simulations and numerical solvers, we also investigate the cases where the conditions for the exponential rate theorem are violated, notably when attracting limit cycles appear, or when the fixed point is singular. We illustrate our theorem on a Markovian fading channel model, which has been well studied in the literature. Finally, we extend our synchronous model results to the asynchronous model.