Improvement on corrosion resistance and biocompability of ZK60 magnesium alloy by carboxyl ion implantation. Magnesium alloys have been considered to be potential biocompatible metallic materials. Further improvement on the anti corrosion is expected to make this type of materials more suitable for biomedical applications in the fields of orthopedics, cardiovascular surgery and others. In this paper, we introduce a method of carboxyl ion (COOH ) implantation to reduce the degradation of ZK60 Mg alloy and improve its functionality in physiological environment. X ray photoelectron spectroscopy (XPS) and atomic force microscopy (AFM) experiments show the formation of a smooth layer containing carbaxylic group, carbonate, metal oxides and hydroxides on the ion implanted alloy surface. Corrosion experiments and in vitro cytotoxicity tests demonstrate that the ion implantation treatment can both reduce the corrosion rate and improve the biocompatibility of the alloy. The promising results indicate that organic functional group ion implantation may be a practical method of improving the biological and corrosion properties of magnesium alloys.


Nonpositive Curvature of the quantomorphism group and quasigeostrophic motion. In this paper, we compute the sectional curvature of the quantomorphism group mathcal D q(M) whose geodesic equation is the quasi geostrophic (QG) equation in geophysics and oceanography, for flows with a stream function depending on only one variable. Using this explicit formula, we will also derive a criterion for the curvature operator to be nonpositive and discuss the role of the Froude number and the Rossby number on curvature. The main technique to obtain a usable formula is a simplification of Arnold's general formula in the case where a vector field is close to a Killing field, and then use the Green's function explicitly. We show that nonzero Froude number and Rossby numbers both tend to stabilize flows in the Lagrangian sense.


Data Driven Game Development Ethical Considerations. In recent years, the games industry has made a major move towards data driven development, using data analytics and player modeling to inform design decisions. Data driven techniques are beneficial as they allow for the study of player behavior at scale, making them very applicable to modern digital game development. However, with this move towards data driven decision making comes a number of ethical concerns. Previous work in player modeling as well as work in the fields of AI and machine learning have demonstrated several ways in which algorithmic decision making can be flawed due to data or algorithmic bias or lack of data from specific groups. Further, black box algorithms create a trust problem due to lack of interpretability and transparency of the results or models developed based on the data, requiring blind faith in the results. In this position paper, we discuss several factors affecting the use of game data in the development cycle. In addition to issues raised by previous work, we also raise issues with algorithms marginalizing certain player groups and flaws in the resulting models due to their inability to reason about situational factors affecting players' decisions. Further, we outline some work that seeks to address these problems and identify some open problems concerning ethics and game data science.


Entropy Estimation of Physically Unclonable Functions via Chow Parameters. A physically unclonable function (PUF) is an electronic circuit that produces an intrinsic identifier in response to a challenge. These identifiers depend on uncontrollable variations of the manufacturing process, which make them hard to predict or to replicate. Various security protocols leverage on such intrinsic randomness for authentification, cryptographic key generation, anti counterfeiting, etc. Evaluating the entropy of PUFs (for all possible challenges) allows one to assess the security properties of such protocols. In this paper, we estimate the probability distribution of certain kinds of PUFs composed of n delay elements. This is used to evaluate relevant R 'enyi entropies and determine how they increase with n . Such a problem was known to have extremely high complexity (in the order of 2 2 n ) and previous entropy estimations were carried out up to n 7 . Making the link with the theory of Boolean threshold functions, we leverage on the representation by Chow parameters to estimate probability distributions up to n 10 . The resulting Shannon entropy of the PUF is close to the max entropy, which is asymptotically quadratic in n .


Development of Adaptive Frame Reservation Scheme and Naive Persistent State Co Located Coexistence Controller. Future broadband networks need to provide high capacity at low cost with increased revenue through enhanced services. WiMAX came up as one of the leading technologies, however, the 2.3 GHz and 2.5 GHz frequency bands allocated create two serious coexistence issues with the adjacent 2.4 GHz ISM band. First problem is to address radio interfaces that are located on two independent platforms and still possess the potential for mutual interference owing to close proximity to each other. The Adaptive Frame Reservation Scheme presented here extends the CTS frame reservation signaling defined in 802.11 specifications to a demand based and adaptive scheme. Second issue is to address the coexistence problem in multi radio platforms where two or more radios are co located, creating an even worse interference scenario. This can be managed by hardware signaling that can be made available between radio interfaces through OS control. The development of a smart Co located Coexistence Controller is explored which continuously receives transmission, reception and sleep requests from attached interfaces and in return grant permissions.


A short proof of the blow up lemma for approximate decompositions. Kim, K "uhn, Osthus and Tyomkyn (Trans. Amer. Math. Soc. 371 (2019), 4655 4742) greatly extended the well known blow up lemma of Koml 'os, S 'ark "ozy and Szemer 'edi by proving a blow up lemma for approximate decompositions' which states that multipartite quasirandom graphs can be almost decomposed into any collection of bounded degree graphs with the same multipartite structure and slightly fewer edges. This result has already been used by Joos, Kim, K "uhn and Osthus to prove the tree packing conjecture due to Gy 'arf 'as and Lehel from 1976 and Ringel's conjecture from 1963 for bounded degree trees as well as implicitly in the recent resolution of the Oberwolfach problem (asked by Ringel in 1967) by Glock, Joos, Kim, K "uhn and Osthus. Here we present a new and significantly shorter proof of the blow up lemma for approximate decompositions. In fact, we prove a more general theorem that yields packings with stronger quasirandom properties so that it can be combined with Keevash's results on designs to obtain results of the following form. For all varepsilon 0 , r in mathbb N and all large n (such that r divides n 1 ), there is a decomposition of K n into any collection of r regular graphs H 1, ldots,H (n 1) r on n vertices provided that H 1, ldots,H varepsilon n contain each at least varepsilon n vertices in components of size at most varepsilon 1 .


Controlling acoustic waves using magnetoelastic Fano resonances. We propose and analyze theoretically a class of energy efficient magneto elastic devices for analogue signal processing. The signals are carried by transverse acoustic waves while the bias magnetic field controls their scattering from a magneto elastic slab. By tuning the bias field, one can alter the resonant frequency at which the propagating acoustic waves hybridize with the magnetic modes, and thereby control transmission and reflection coefficients of the acoustic waves. The scattering coefficients exhibit Breit Wigner Fano resonant behaviour akin to inelastic scattering in atomic and nuclear physics. Employing oblique incidence geometry, one can effectively enhance the strength of magnetoelastic coupling, and thus countermand the magnetic losses due to the Gilbert damping. We apply our theory to discuss potential benefits and issues in realistic systems and suggest further routes to enhance performance of the proposed devices.


Four generated direct powers of partition lattices and authentication. For an integer n geq 5 , H. Strietz (1975) and L. Z 'adori (1986) proved that the lattice Part (n) of all partitions of 1,2, dots,n is four generated. Developing L. Z 'adori's particularly elegant construction further, we prove that even the k th direct power Part (n) k of Part (n) is four generated for many but only finitely many exponents k . E.g., Part (n) k is four generated for every k leq 3 cdot 10 89 , and it has a four element generating set that is not an antichain for every k leq 1.4 cdot 10 34 . In connection with these results, we outline a protocol how to use these lattices in authentication and secret key cryptography.


Introspective Learning by Distilling Knowledge from Online Self explanation. In recent years, many explanation methods have been proposed to explain individual classifications of deep neural networks. However, how to leverage the created explanations to improve the learning process has been less explored. As the privileged information, the explanations of a model can be used to guide the learning process of the model itself. In the community, another intensively investigated privileged information used to guide the training of a model is the knowledge from a powerful teacher model. The goal of this work is to leverage the self explanation to improve the learning process by borrowing ideas from knowledge distillation. We start by investigating the effective components of the knowledge transferred from the teacher network to the student network. Our investigation reveals that both the responses in non ground truth classes and class similarity information in teacher's outputs contribute to the success of the knowledge distillation. Motivated by the conclusion, we propose an implementation of introspective learning by distilling knowledge from online self explanations. The models trained with the introspective learning procedure outperform the ones trained with the standard learning procedure, as well as the ones trained with different regularization methods. When compared to the models learned from peer networks or teacher networks, our models also show competitive performance and requires neither peers nor teachers.


Conformal Perturbation of Heat Kernels with applications. Let (M, g) be a smooth n dimensional Riemannian manifold for n ge 2 . Consider the conformal perturbation tilde g h g where h is a smooth bounded positive function on M . Denote by tilde p t(x,y) the heat kernel of manifolds (M, tilde g ) . In this paper, we derive the upper bounds and gradient estimates of tilde p t(x,y) .


Experimental study on wave attenuation by vegetation by means of particle imaging velocimetry. Coastal vegetation serves as an ecological defense and plays an important role when facing marine disaster such as storm surge. In this study, we investigated experimentally the wave attenuation, including wave height and wave force decay, and flow field within rigid vegetation by means of particle imaging velocimetry (PIV). Drag coefficient of vegetation is closely related the wave attenuation by vegetation. Based on Morison equation, a new modified direct measurement method to derive the drag coefficient of vegetation was proposed in this study. This method can eliminate the potential errors and obtain more accurate drag coefficient to have a well fitted relationship with Reynolds number compared with the previous methods. More importantly, we showed the flow field including instantaneous velocity filed in different wave phase within the vegetation and turbulent kinetic energy (TKE) field to describe the mechanisms of wave energy dissipation. In addition, a prediction model to reproduce the instantaneous TKE was proposed by the work input of vegetation. The results illustrated that all these characteristics of parameters can help us understand the effect of the coastal vegetation on wave attenuation better.


Multi level Cross modal Interaction Network for RGB D Salient Object Detection. Depth cues with affluent spatial information have been proven beneficial in boosting salient object detection (SOD), while the depth quality directly affects the subsequent SOD performance. However, it is inevitable to obtain some low quality depth cues due to limitations of its acquisition devices, which can inhibit the SOD performance. Besides, existing methods tend to combine RGB images and depth cues in a direct fusion or a simple fusion module, which makes they can not effectively exploit the complex correlations between the two sources. Moreover, few methods design an appropriate module to fully fuse multi level features, resulting in cross level feature interaction insufficient. To address these issues, we propose a novel Multi level Cross modal Interaction Network (MCINet) for RGB D based SOD. Our MCI Net includes two key components 1) a cross modal feature learning network, which is used to learn the high level features for the RGB images and depth cues, effectively enabling the correlations between the two sources to be exploited and 2) a multi level interactive integration network, which integrates multi level cross modal features to boost the SOD performance. Extensive experiments on six benchmark datasets demonstrate the superiority of our MCI Net over 14 state of the art methods, and validate the effectiveness of different components in our MCI Net. More important, our MCI Net significantly improves the SOD performance as well as has a higher FPS.


Symmetry constrained Rectification Network for Scene Text Recognition. Reading text in the wild is a very challenging task due to the diversity of text instances and the complexity of natural scenes. Recently, the community has paid increasing attention to the problem of recognizing text instances with irregular shapes. One intuitive and effective way to handle this problem is to rectify irregular text to a canonical form before recognition. However, these methods might struggle when dealing with highly curved or distorted text instances. To tackle this issue, we propose in this paper a Symmetry constrained Rectification Network (ScRN) based on local attributes of text instances, such as center line, scale and orientation. Such constraints with an accurate description of text shape enable ScRN to generate better rectification results than existing methods and thus lead to higher recognition accuracy. Our method achieves state of the art performance on text with both regular and irregular shapes. Specifically, the system outperforms existing algorithms by a large margin on datasets that contain quite a proportion of irregular text instances, e.g., ICDAR 2015, SVT Perspective and CUTE80.


A Self Matched Leaky Wave Antenna for Ultrahigh Field MRI with Low SAR. The technology of magnetic resonance imaging is developing towards higher magnetic fields to improve resolution and contrast. However, whole body imaging at 7 T or even higher fields remains challenging due to wave interference, tissue inhomogneities and high RF power deposition. Nowadays, proper RF excitation of a human body in prostate and cardiac MRI is only possible to achieve by using phased arrays of antennas attached to the body (so called surface coils). Due to safety concerns, the design of such coils aims to minimize the local specific absorption rate (SAR) keeping the highest possible RF signal in the region of interest. All previously demonstrated approaches were based on resonant structures such as e. g.dipoles, capacitively loaded loops, TEM line sections. In this study, we show that there is a better compromise between the transmit signal and the local SAR using non resonant surface coils due to weaker RF near fields in the close proximity of their conductors. With this aim, we propose and experimentally demonstrate a first leaky wave surface coil implemented as a periodically slotted microstrip transmission line. Due to its non resonant radiation, the proposed coil induces only slightly over half the peak local SAR compared to a state of the art dipole coil, but has the same transmit efficiency in prostate imaging at 7 T. Unlike other coils, the leaky wave coil intrinsically matches its input impedance to the averaged wave impedance of body tissues in a broad frequency range, which makes it very attractive for future clinical applications of 7 T MRI.


Simultaneous Stabilization of Traffic Flow on Two Connected Roads. In this paper we develop a boundary state feedback control law for a traffic flow network system in its most fundamental form one incoming and one outgoing road connected by a junction. The macroscopic traffic dynamics on each road segment are governed by Aw Rascle Zhang (ARZ) model, consisting of second order nonlinear partial differential equations (PDEs) for traffic density and velocity. Different equilibrium road conditions are considered for the connected segments. For stabilization of the stop and go traffic congestion on the two roads, we consider a ramp metering located at the connecting junction. The traffic flow rate entering from the on ramp to the mainline junction is actuated. The objective is to simultaneously stabilize the upstream and downstream traffic to a given spatially uniform constant steady state. We design a full state feedback control law for this under actuated network of two systems of two hetero directional linear first order hyperbolic PDEs interconnected through the boundary condition (junction). The exponential stability is validated by numerical simulation.


Unsupervised Domain Adaptation for Mobile Semantic Segmentation based on Cycle Consistency and Feature Alignment. The supervised training of deep networks for semantic segmentation requires a huge amount of labeled real world data. To solve this issue, a commonly exploited workaround is to use synthetic data for training, but deep networks show a critical performance drop when analyzing data with slightly different statistical properties with respect to the training set. In this work, we propose a novel Unsupervised Domain Adaptation (UDA) strategy to address the domain shift issue between real world and synthetic representations. An adversarial model, based on the cycle consistency framework, performs the mapping between the synthetic and real domain. The data is then fed to a MobileNet v2 architecture that performs the semantic segmentation task. An additional couple of discriminators, working at the feature level of the MobileNet v2, allows to better align the features of the two domain distributions and to further improve the performance. Finally, the consistency of the semantic maps is exploited. After an initial supervised training on synthetic data, the whole UDA architecture is trained end to end considering all its components at once. Experimental results show how the proposed strategy is able to obtain impressive performance in adapting a segmentation network trained on synthetic data to real world scenarios. The usage of the lightweight MobileNet v2 architecture allows its deployment on devices with limited computational resources as the ones employed in autonomous vehicles.


Mitigating Leakage in Federated Learning with Trusted Hardware. In federated learning, multiple parties collaborate in order to train a global model over their respective datasets. Even though cryptographic primitives (e.g., homomorphic encryption) can help achieve data privacy in this setting, some partial information may still be leaked across parties if this is done non judiciously. In this work, we study the federated learning framework of SecureBoost Cheng et al., FL IJCAI'19 as a specific such example, demonstrate a leakage abuse attack based on its leakage profile, and experimentally evaluate the effectiveness of our attack. We then propose two secure versions relying on trusted execution environments. We implement and benchmark our protocols to demonstrate that they are 1.2 5.4X faster in computation and need 5 49X less communication than SecureBoost.


Real zeros of algebraic polynomials with nonidentical dependent random coefficients. The expected number of real zeros of an algebraic polynomial a 0 a 1x a 2x 2 a 3x 3 .... a n 1 x n 1 depends on the types of random coefficients, with large n. In this article, we show that when the random coefficients a i i 1 n 1 are assumed to be negatively dependent with var(a i) sigma 2i and correlation between any two coefficients for i neq j, assumed to be rho ij rho i j , where 0 rho frac 1 3 , then the expected number of real zeros is asymptotically equal to frac 2 pi sigma logn.


An orthogonality relation for GL(4,R). Orthogonality is a fundamental theme in representation theory and Fourier analysis. An orthogonality relation for characters of finite abelian groups (now recognized as an orthogonality relation on mathrm GL (1) ) was used by Dirichlet to prove infinitely many primes in arithmetic progressions. Orthogonality relations for mathrm GL (2) and mathrm GL (3) have been worked on by many researchers with a broad range of applications to number theory. We present here, for the first time, very explicit orthogonality relations for the real group mathrm GL (4, mathbb R ) with a power savings error term. The proof requires novel techniques in the computation of the geometric side of the Kuznetsov trace formula. An appendix by Bingrong Huang gives new bounds for the relevant Kloosterman sums.


Inverse Problems of Determining Coefficients of the Fractional Partial Differential Equations. When considering fractional diffusion equation as model equation in analyzing anomalous diffusion processes, some important parameters in the model, for example, the orders of the fractional derivative or the source term, are often unknown, which requires one to discuss inverse problems to identify these physical quantities from some additional information that can be observed or measured practically. This chapter investigates several kinds of inverse coefficient problems for the fractional diffusion equation.


Im2Pencil Controllable Pencil Illustration from Photographs. We propose a high quality photo to pencil translation method with fine grained control over the drawing style. This is a challenging task due to multiple stroke types (e.g., outline and shading), structural complexity of pencil shading (e.g., hatching), and the lack of aligned training data pairs. To address these challenges, we develop a two branch model that learns separate filters for generating sketchy outlines and tonal shading from a collection of pencil drawings. We create training data pairs by extracting clean outlines and tonal illustrations from original pencil drawings using image filtering techniques, and we manually label the drawing styles. In addition, our model creates different pencil styles (e.g., line sketchiness and shading style) in a user controllable manner. Experimental results on different types of pencil drawings show that the proposed algorithm performs favorably against existing methods in terms of quality, diversity and user evaluations.


CONICOPF Conic relaxations for AC optimal power flow computations. Computational speed and global optimality are key needs for practical algorithms for the optimal power flow problem. Two convex relaxations offer a favorable trade off between the standard second order cone and the standard semidefinite relaxations for large scale meshed networks in terms of optimality gap and computation time the tight and cheap relaxation (TCR) and the quadratic convex relaxation (QCR). We compare these relaxations on 60 PGLib OPF test cases with up to 1,354 buses under three operating conditions and show that TCR dominates QCR on all 20 typical (TYP) test cases, on 18 out of 20 active power increase (API), and 12 out of 20 small angle difference (SAD). Selected state of the art conic relaxations are implemented in the new MATLAB based package CONICOPF available on GitHub.


Algebraic Language Theory for Eilenberg Moore Algebras. We develop an algebraic language theory based on the notion of an Eilenberg Moore algebra. In comparison to previous such frameworks the main contribution is the support for algebras with infinitely many sorts and the connection to logic in form of so called definable algebras'.


Facial Expressions Analysis Under Occlusions Based on Specificities of Facial Motion Propagation. Although much progress has been made in the facial expression analysis field, facial occlusions are still challenging. The main innovation brought by this contribution consists in exploiting the specificities of facial movement propagation for recognizing expressions in presence of important occlusions. The movement induced by an expression extends beyond the movement epicenter. Thus, the movement occurring in an occluded region propagates towards neighboring visible regions. In presence of occlusions, per expression, we compute the importance of each unoccluded facial region and we construct adapted facial frameworks that boost the performance of per expression binary classifier. The output of each expression dependant binary classifier is then aggregated and fed into a fusion process that aims constructing, per occlusion, a unique model that recognizes all the facial expressions considered. The evaluations highlight the robustness of this approach in presence of significant facial occlusions.


Higher Nash blowups of normal toric varieties in prime characteristic. We prove that the higher Nash blowup of a normal toric variety defined over a field of positive characteristic is an isomorphism if and only if it is non singular. We also extend a result by R. Toh Yama which shows that higher Nash blowups do not give a one step resolution of the A 3 singularity. These results were previously known only in characteristic zero.


Learning to Detect and Retrieve Objects from Unlabeled Videos. Learning an object detector or retrieval requires a large data set with manual annotations. Such data sets are expensive and time consuming to create and therefore difficult to obtain on a large scale. In this work, we propose to exploit the natural correlation in narrations and the visual presence of objects in video, to learn an object detector and retrieval without any manual labeling involved. We pose the problem as weakly supervised learning with noisy labels, and propose a novel object detection paradigm under these constraints. We handle the background rejection by using contrastive samples and confront the high level of label noise with a new clustering score. Our evaluation is based on a set of 11 manually annotated objects in over 5000 frames. We show comparison to a weakly supervised approach as baseline and provide a strongly labeled upper bound.


Nash equilibria in games over graphs equipped with a communication mechanism. We study pure Nash equilibria in infinite duration games on graphs, with partial visibility of actions but communication (based on a graph) among the players. We show that a simple communication mechanism consisting in reporting the deviator when seeing it and propagating this information is sufficient for characterizing Nash equilibria. We propose an epistemic game construction, which conveniently records important information about the knowledge of the players. With this abstraction, we are able to characterize Nash equilibria which follow the simple communication pattern via winning strategies. We finally discuss the size of the construction, which would allow efficient algorithmic solutions to compute Nash equilibria in the original game.


Polytope duality for families of K3 surfaces and coupling. We study a relation between coupling introduced by Ebeling and the polytope duality among families of K3 surfaces.


Data Selection for training Semantic Segmentation CNNs with cross dataset weak supervision. Training convolutional networks for semantic segmentation with strong (per pixel) and weak (per bounding box) supervision requires a large amount of weakly labeled data. We propose two methods for selecting the most relevant data with weak supervision. The first method is designed for finding visually similar images without the need of labels and is based on modeling image representations with a Gaussian Mixture Model (GMM). As a byproduct of GMM modeling, we present useful insights on characterizing the data generating distribution. The second method aims at finding images with high object diversity and requires only the bounding box labels. Both methods are developed in the context of automated driving and experimentation is conducted on Cityscapes and Open Images datasets. We demonstrate performance gains by reducing the amount of employed weakly labeled images up to 100 times for Open Images and up to 20 times for Cityscapes.


Max Min Lyapunov Functions for Switched Systems and Related Differential Inclusions. Starting from a finite family of continuously differentiable positive definite functions, we study conditions under which a function obtained by max min combinations is a Lyapunov function, establishing stability for two kinds of nonlinear dynamical systems a) Differential inclusions where the set valued right hand side comprises the convex hull of a finite number of vector fields, and b) Autonomous switched systems with a state dependent switching signal. We investigate generalized notions of directional derivatives for these max min functions, and use them in deriving stability conditions with various degrees of conservatism, where more conservative conditions are numerically more tractable. The proposed constructions also provide nonconvex Lyapunov functions, which are shown to be useful for systems with state dependent switching that do not admit a convex Lyapunov function. Several examples are included to illustrate the results.


MAT Multi Fingered Adaptive Tactile Grasping via Deep Reinforcement Learning. Vision based grasping systems typically adopt an open loop execution of a planned grasp. This policy can fail due to many reasons, including ubiquitous calibration error. Recovery from a failed grasp is further complicated by visual occlusion, as the hand is usually occluding the vision sensor as it attempts another open loop regrasp. This work presents MAT, a tactile closed loop method capable of realizing grasps provided by a coarse initial positioning of the hand above an object. Our algorithm is a deep reinforcement learning (RL) policy optimized through the clipped surrogate objective within a maximum entropy RL framework to balance exploitation and exploration. The method utilizes tactile and proprioceptive information to act through both fine finger motions and larger regrasp movements to execute stable grasps. A novel curriculum of action motion magnitude makes learning more tractable and helps turn common failure cases into successes. Careful selection of features that exhibit small sim to real gaps enables this tactile grasping policy, trained purely in simulation, to transfer well to real world environments without the need for additional learning. Experimentally, this methodology improves over a vision only grasp success rate substantially on a multi fingered robot hand. When this methodology is used to realize grasps from coarse initial positions provided by a vision only planner, the system is made dramatically more robust to calibration errors in the camera robot transform.


Tropical degenerations and stable rationality. We use the motivic obstruction to stable rationality introduced by Shinder and the first named author to establish several new classes of stably irrational hypersurfaces and complete intersections. In particular, we show that very general quartic fivefolds and complete intersections of a quadric and a cubic in mathbb P 6 are stably irrational. An important new ingredient is the use of tropical degeneration techniques.


Derandomized compressed sensing with nonuniform guarantees for ell 1 recovery. We extend the techniques of H " u gel, Rauhut and Strohmer (Found. Comput. Math., 2014) to show that for every delta in(0,1 , there exists an explicit random m times N partial Fourier matrix A with m s operatorname polylog (N epsilon) and entropy s delta operatorname polylog (N epsilon) such that for every s sparse signal x in mathbb C N , there exists an event of probability at least 1 epsilon over which x is the unique minimizer of z 1 subject to Az Ax . The bulk of our analysis uses tools from decoupling to estimate the extreme singular values of the submatrix of A whose columns correspond to the support of x .


Coloring of Graphs Avoiding Bicolored Paths of a Fixed Length. The problem of finding the minimum number of colors to color a graph properly without containing any bicolored copy of a fixed family of subgraphs has been widely studied. Most well known examples are star coloring and acyclic coloring of graphs (Gr "unbaum, 1973) where bicolored copies of P 4 and cycles are not allowed, respectively. In this paper, we introduce a variation of these problems and study proper coloring of graphs not containing a bicolored path of a fixed length and provide general bounds for all graphs. A P k coloring of an undirected graph G is a proper vertex coloring of G such that there is no bicolored copy of P k in G, and the minimum number of colors needed for a P k coloring of G is called the P k chromatic number of G, denoted by s k(G). We provide bounds on s k(G) for all graphs, in particular, proving that for any graph G with maximum degree d geq 2, and k geq4, s k(G) O(d frac k 1 k 2 ). Moreover, we find the exact values for the P k chromatic number of the products of some cycles and paths for k 5,6.


Why and How to Balance Alignment and Diversity of Requirements Engineering Practices in Automotive. In large scale automotive companies, various requirements engineering (RE) practices are used across teams. RE practices manifest in Requirements Information Models (RIM) that define what concepts and information should be captured for requirements. Collaboration of practitioners from different parts of an organization is required to define a suitable RIM that balances support for diverse practices in individual teams with the alignment needed for a shared view and team support on system level. There exists no guidance for this challenging task. This paper presents a mixed methods study to examine the role of RIMs in balancing alignment and diversity of RE practices in four automotive companies. Our analysis is based on data from systems engineering tools, 11 semi structured interviews, and a survey to validate findings and suggestions. We found that balancing alignment and diversity of RE practices is important to consider when defining RIMs. We further investigated enablers for this balance and actions that practitioners take to achieve it. From these factors, we derived and evaluated recommendations for managing RIMs in practice that take into account the lifecycle of requirements and allow for diverse practices across sub disciplines in early development, while enforcing alignment of requirements that are close to release.


Precise Morrey regularity of the weak solutions to a kind of quasilinear systems with discontinuous data. We consider the Dirichlet problem for a class of quasilinear elliptic systems in domain with irregular boundary. The principal part satisfies componentwise coercivity condition and the nonlinear terms are Carath 'eodory maps having Morrey regularity in x and verifying controlled growth conditions with respect to the other variables. We have obtained boundedness of the weak solution to the problem that permits to apply an iteration procedure in order to find optimal Morrey regularity of its gradient.


Thermodynamically consistent physics informed neural networks for hyperbolic systems. Physics informed neural network architectures have emerged as a powerful tool for developing flexible PDE solvers which easily assimilate data, but face challenges related to the PDE discretization underpinning them. By instead adapting a least squares space time control volume scheme, we circumvent issues particularly related to imposition of boundary conditions and conservation while reducing solution regularity requirements. Additionally, connections to classical finite volume methods allows application of biases toward entropy solutions and total variation diminishing properties. For inverse problems, we may impose further thermodynamic biases, allowing us to fit shock hydrodynamics models to molecular simulation of rarefied gases and metals. The resulting data driven equations of state may be incorporated into traditional shock hydrodynamics codes.


Network Flows that Solve Sylvester Matrix Equations. In this paper, we study distributed methods for solving a Sylvester equation in the form of AX XB C for matrices A, B, C in R n times n with X being the unknown variable. The entries of A, B and C (called data) are partitioned into a number of pieces (or sometimes we permit these pieces to overlap). Then a network with a given structure is assigned, whose number of nodes is consistent with the partition. Each node has access to the corresponding set of data and holds a dynamic state. Nodes share their states among their neighbors defined from the network structure, and we aim to design flows that can asymptotically converge to a solution of this equation. The decentralized data partitions may be resulted directly from networks consisting of physically isolated subsystems, or indirectly from artificial and strategic design for processing large data sets. Natural partial row column partitions, full row column partitions and clustering block partitions of the data A, B and C are assisted by the use of the vectorized matrix equation. We show that the existing "consensus projection" flow and the "local conservation global consensus" flow for distributed linear algebraic equations can be used to drive distributed flows that solve this kind of equations. A "consensus projection symmetrization" flow is also developed for equations with symmetry constraints on the solution matrices. We reveal some fundamental convergence rate limitations for such flows regardless of the choices of node interaction strengths and network structures. For a special case with B A T , where the equation mentioned is reduced to a classical Lyapunov equation, we demonstrate that by exploiting the symmetry of data, we can obtain flows with lower complexity for certain partitions.


Knowledge Based Legal Document Assembly. This paper proposes a knowledge based legal document assembly method that uses a machine readable representation of knowledge of legal professionals. This knowledgebase has two components the formal knowledge of legal norms represented as a rule base and the tacit knowledge represented by a document template. A document assembly system is developed as a proof of concept. It collects input data in the form of an interactive interview, performs legal reasoning over input data, and generates the output document. The system also creates an argument graph as an explanation of the reasoning process providing the user with an interpretation of how the input data and the rule base influence the content of the output document. The system also semantically marks up data in the output document, facilitating its further processing and providing support to the interoperability of information systems in the legal domain.


SimPose Effectively Learning DensePose and Surface Normals of People from Simulated Data. With a proliferation of generic domain adaptation approaches, we report a simple yet effective technique for learning difficult per pixel 2.5D and 3D regression representations of articulated people. We obtained strong sim to real domain generalization for the 2.5D DensePose estimation task and the 3D human surface normal estimation task. On the multi person DensePose MSCOCO benchmark, our approach outperforms the state of the art methods which are trained on real images that are densely labelled. This is an important result since obtaining human manifold's intrinsic uv coordinates on real images is time consuming and prone to labeling noise. Additionally, we present our model's 3D surface normal predictions on the MSCOCO dataset that lacks any real 3D surface normal labels. The key to our approach is to mitigate the "Inter domain Covariate Shift" with a carefully selected training batch from a mixture of domain samples, a deep batch normalized residual network, and a modified multi task learning objective. Our approach is complementary to existing domain adaptation techniques and can be applied to other dense per pixel pose estimation problems.


A two stage algorithm for aircraft conflict resolution with trajectory recovery. As air traffic volume is continuously increasing, it has become a priority to improve traffic control algorithms to handle future air travel demand and improve airspace capacity. We address the conflict resolution problem in air traffic control using a novel approach for aircraft collision avoidance with trajectory recovery. We present a two stage algorithm that first solves all initial conflicts by adjusting aircraft headings and speeds, before identifying the optimal time for aircraft to recover towards their target destination. The collision avoidance stage extends an existing mixed integer programming formulation to heading control. For the trajectory recovery stage, we introduce a novel exact mixed integer programming formulation as well as a greedy heuristic algorithm. The proposed two stage approach guarantees that all trajectories during both the collision avoidance and recovery stages are conflict free. Numerical results on benchmark problems show that the proposed heuristic for trajectory recovery is competitive while also emphasizing the difficulty of this optimization problem. The proposed approach can be used as a decision support tool for introducing automation in air traffic control.


A COVID 19 Epidemiological Model for Community and Policy Maker Use. An epidemiological model for COVID 19 was developed and implemented in MATLAB GNU Octave for use by public health practitioners, policy makers and the general public. The model distinguishes four stages in the disease infected, sick, seriously sick, and better. The model was preliminarily parameterized based on observations of the spread of the disease. The model is consistent with a mortality rate of 1.5 . Preliminary simulations with the model indicate that concepts such as "herd immunity" and "flattening the curve" are highly misleading in the context of this virus. Public policies based on these concepts are inadequate to protect the population. Only reducing the R0 of the virus below 1 is an effective strategy for maintaining the death burden of COVID 19 within the normal range of seasonal flu. As R0 values estimated with the model range from 2.82 worldwide outside of China and 3.83 in the Western world in late February early March 2020, this means social distancing with effectiveness greater than 65 (worldwide) or 75 (Western world) are needed to combat the virus successfully.


Reproducible Workflow on a Public Cloud for Computational Fluid Dynamics. In a new effort to make our research transparent and reproducible by others, we developed a workflow to run and share computational studies on the public cloud Microsoft Azure. It uses Docker containers to create an image of the application software stack. We also adopt several tools that facilitate creating and managing virtual machines on compute nodes and submitting jobs to these nodes. The configuration files for these tools are part of an expanded "reproducibility package" that includes workflow definitions for cloud computing, in addition to input files and instructions. This facilitates re creating the cloud environment to re run the computations under the same conditions. Although cloud providers have improved their offerings, many researchers using high performance computing (HPC) are still skeptical about cloud computing. Thus, we ran benchmarks for tightly coupled applications to confirm that the latest HPC nodes of Microsoft Azure are indeed a viable alternative to traditional on site HPC clusters. We also show that cloud offerings are now adequate to complete computational fluid dynamics studies with in house research software that uses parallel computing with GPUs. Finally, we share with the community what we have learned from nearly two years of using Azure cloud to enhance transparency and reproducibility in our computational simulations.


VOC Sensing Using Batch fabricated Temperature Compensated Self Leveling Microstructures. We present the design, fabrication, and response of a low power, polymer based VOC sensor based on the self leveling of mechanically leveraged structures. The device utilizes folded polymer coated microcantilevers to achieve passive temperature compensation without the need for additional compensating sensors or electronics. We demonstrate that a self leveling vapor sensor provides the same gas response as a simple microcantilever geometry, showing 20 change in device capacitance when subjected to 35 85 RH change while showing nearly zero baseline drift due to changes in ambient temperature when the temperature is increased from 23 72 deg C which is 52 fold better than a simple microcantilever geometry. The response of the VOC sensor was measured using three polymers (Polyimide, Polyurethane, and PDMS) against five different analytes (Ethanol, Acetone, Benzene, Hexane, and Water) and an SVM based model was used to show target specificity. The sensor also showed an absorption response time ( tau 90) of 138s. We propose that the self leveling vapor sensor geometry is a significant improvement to a simple microcantilever vapor sensor as it offers the same performance but shows near complete elimination of temperature induced baseline drift.


Bergman projection on the symmetrized bidisk. We apply the Bekoll 'e Bonami estimate for the (positive) Bergman projection on the weighted L p spaces on the unit disk. We then obtain the boundedness of the Bergman projection on the weighted Sobolev space on the symmetrized bidisk, by the reduction to the (positive) Bergman projection on the weighted L p space on the unit disk,. We also improve the boundedness result of the Bergman projection on the unweighted L p space on the symmetrized bidisk in cite CKY .


Exploring Sequence to Sequence Models for SPARQL Pattern Composition. A booming amount of information is continuously added to the Internet as structured and unstructured data, feeding knowledge bases such as DBpedia and Wikidata with billions of statements describing millions of entities. The aim of Question Answering systems is to allow lay users to access such data using natural language without needing to write formal queries. However, users often submit questions that are complex and require a certain level of abstraction and reasoning to decompose them into basic graph patterns. In this short paper, we explore the use of architectures based on Neural Machine Translation called Neural SPARQL Machines to learn pattern compositions. We show that sequence to sequence models are a viable and promising option to transform long utterances into complex SPARQL queries.


Tethered Aerial Visual Assistance. In this paper, an autonomous tethered Unmanned Aerial Vehicle (UAV) is developed into a visual assistant in a marsupial co robots team, collaborating with a tele operated Unmanned Ground Vehicle (UGV) for robot operations in unstructured or confined environments. These environments pose extreme challenges to the remote tele operator due to the lack of sufficient situational awareness, mostly caused by the unstructuredness and confinement, stationary and limited field of view and lack of depth perception from the robot's onboard cameras. To overcome these problems, a secondary tele operated robot is used in current practices, who acts as a visual assistant and provides external viewpoints to overcome the perceptual limitations of the primary robot's onboard sensors. However, a second tele operated robot requires extra manpower and teamwork demand between primary and secondary operators. The manually chosen viewpoints tend to be subjective and sub optimal. Considering these intricacies, we develop an autonomous tethered aerial visual assistant in place of the secondary tele operated robot and operator, to reduce human robot ratio from 2 2 to 1 2. Using a fundamental viewpoint quality theory, a formal risk reasoning framework, and a newly developed tethered motion suite, our visual assistant is able to autonomously navigate to good quality viewpoints in a risk aware manner through unstructured or confined spaces with a tether. The developed marsupial co robots team could improve tele operation efficiency in nuclear operations, bomb squad, disaster robots, and other domains with novel tasks or highly occluded environments, by reducing manpower and teamwork demand, and achieving better visual assistance quality with trustworthy risk aware motion.


Uniform Error Estimates for the Lanczos Method. The Lanczos method is one of the most powerful and fundamental techniques for solving an extremal symmetric eigenvalue problem. Convergence based error estimates depend heavily on the eigenvalue gap. In practice, this gap is often relatively small, resulting in significant overestimates of error. One way to avoid this issue is through the use of uniform error estimates, namely, bounds that depend only on the dimension of the matrix and the number of iterations. In this work, we prove explicit upper and lower uniform error estimates for the Lanczos method. These lower bounds, paired with numerical results, imply that the maximum error of m iterations of the Lanczos method over all n times n symmetric matrices does indeed depend on the dimension n in practice. The improved bounds for extremal eigenvalues translate immediately to error estimates for the condition number of a symmetric positive definite matrix. In addition, we prove more specific results for matrices that possess some level of eigenvalue regularity or whose eigenvalues converge to some limiting empirical spectral distribution. Through numerical experiments, we show that the theoretical estimates of this paper do apply to practical computations for reasonably sized matrices.


Using double Weil sums in finding the c Boomerang Connectivity Table for monomial functions on finite fields. In this paper we characterize the c Boomerang Connectivity Table (BCT), c neq 0 (thus, including the classical c 1 case), for all monomial function x d in terms of characters and Weil sums on the finite field F p n . We further simplify these expressions for the Gold functions x p k 1 for all 1 leq k n , and p odd. It is the first such attempt for a complete description for the classical BCT and its relative c BCT, for all parameters involved.


Gradient based block coordinate descent algorithms for joint approximate diagonalization of matrices. In this paper, we propose a gradient based block coordinate descent (BCD G) framework to solve the joint approximate diagonalization of matrices defined on the product of the complex Stiefel manifold and the special linear group. Instead of the cyclic fashion, we choose the block for optimization in a way based on the Riemannian gradient. To update the first block variable in the complex Stiefel manifold, we use the well known line search descent method. To update the second block variable in the special linear group, based on four different kinds of elementary rotations, we construct two classes Jacobi GLU and Jacobi GLQ, and then get two BCD G algorithms BCD GLU and BCD GLQ. We establish the weak convergence and global convergence of these two algorithms using the L ojasiewicz gradient inequality under the assumption that the iterates are bounded. In particular, the problem we focus on in this paper includes as special cases the well known joint approximate diagonalization of Hermitian (or complex symmetric) matrices by invertible transformations in blind source separation, and our algorithms specialize as the gradient based Jacobi type algorithms. All the algorithms and convergence results in this paper also apply to the real case.


Toral CW complexes and bifurcation control in Eulerian flows with multiple Hopf singularities. We are concerned with bifurcation analysis and control of nonlinear Eulerian flows with non resonant n tuple Hopf singularity. The analysis is involved with CW complex bifurcations of flow invariant Clifford hypertori, where we refer to these toral manifolds by toral CW complexes. We observe from primary to tertiary flow invariant toral CW complex bifurcations for one parametric systems associated with two most generic cases. In a particular case, a tertiary toral CW complex bifurcates from and resides outside a secondary toral CW complex. When the parameter varies, the secondary internal toral CW complex collapses with the origin. However, the tertiary external toral CW complex continues to live even after the secondary internal toral manifold disappears. Our analysis starts with a flow invariant primary cell decomposition of the state space. Each open cell admits a secondary cell decomposition via a smooth flow invariant foliation. Each leaf of the foliations is a minimal flow invariant realization of the state space configuration for all Eulerian flows with n tuple Hopf singularities. Permissible leaf vector field preserving transformations are introduced via a Lie algebra structure for nonlinear vector fields on the leaf manifold. Complete parametric leaf normal form classification is provided for singular leaf flows. Leaf bifurcation analysis of leaf normal forms are performed for three most leaf generic cases associated with one to three unfolding bifurcation parameters. Leaf bifurcation varieties are derived. Leaf bifurcations provides a venue for cell bifurcation control of invariant toral CW complexes. The results are implemented and verified using Maple for practical bifurcation control of such parametric nonlinear oscillators.


Cohesive fracture in 1D quasi static evolution and derivation from static phase field models. In this paper we propose a notion of irreversibility for the evolution of cracks in presence of cohesive forces, which allows for different responses in the loading and unloading processes, motivated by a variational approximation with damage models. We investigate its applicability to the construction of a quasi static evolution in a simple one dimensional model. The cohesive fracture model arises naturally via Gamma convergence from a phase field model of the generalized Ambrosio Tortorelli type, which may be used as regularization for numerical simulations.


On the Border of the Amyloidogenic Sequences Prefix Analysis of the Parallel Beta Sheets in the PDB Amyloid Collection. The Protein Data Bank (PDB) today contains more than 153,000 entries with the 3 dimensional structures of biological macromolecules. Using the rich resources of this repository, it is possible identifying subsets with specific, interesting properties for different applications. Our research group prepared an automatically updated list of amyloid and probably amyloidogenic molecules, the PDB Amyloid collection, which is freely available at the address url http pitgroup.org amyloid . This resource applies exclusively the geometric properties of the steric structures for identifying amyloids. In the present contribution, we analyze the starting (i.e., prefix) subsequences of the characteristic, parallel beta sheets of the structures in the PDB Amyloid collection, and identify further appearances of these length 5 prefix subsequences in the whole PDB data set. We have identified this way numerous proteins, whose normal or irregular functions involve amyloid formation, structural misfolding, or anti coagulant properties, simply by containing these prefixes including the T cell receptor (TCR), bound with the major histocompatibility complexes MHC 1 and MHC 2 the p53 tumor suppressor protein a mycobacterial RNA polymerase transcription initialization complex the human bridging integrator protein BIN 1 and the tick anti coagulant peptide TAP.


Optimal polynomial time compression for Boolean Max CSP. In the Boolean maximum constraint satisfaction problem Max CSP ( Gamma) one is given a collection of weighted applications of constraints from a finite constraint language Gamma , over a common set of variables, and the goal is to assign Boolean values to the variables so that the total weight of satisfied constraints is maximized. There exists an elegant dichotomy theorem providing a criterion on Gamma for the problem to be polynomial time solvable and stating that otherwise it becomes NP hard. We study the NP hard cases through the lens of kernelization and provide a complete characterization of Max CSP ( Gamma) with respect to the optimal compression size. Namely, we prove that Max CSP ( Gamma) parameterized by the number of variables n is either polynomial time solvable, or there exists an integer d ge 2 depending on Gamma , such that 1. An instance of textsc Max CSP ( Gamma) can be compressed into an equivalent instance with O(n d log n) bits in polynomial time, 2. Max CSP ( Gamma) does not admit such a compression to O(n d epsilon ) bits unless text NP subseteq text co NP text poly . Our reductions are based on interpreting constraints as multilinear polynomials combined with the framework of constraint implementations. As another application of our reductions, we reveal tight connections between optimal running times for solving Max CSP ( Gamma) . More precisely, we show that obtaining a running time of the form O(2 (1 epsilon)n ) for particular classes of Max CSPs is as hard as breaching this barrier for Max d SAT for some d .


Glauber Exclusion dynamics rapid mixing regime. We show that for any attractive Glauber Exclusion process on the one dimensional lattice of size N with periodic boundary condition, if the corresponding hydrodynamic limit equation has a reaction term with a strictly convex potential, then the total variation mixing time is of order O( log N) . In particular, the result covers the full high temperature regime in the original model introduced by De Masi, Ferrari and Lebowitz (1985).


A Counterexample to the phi Dimension Conjecture. In 2005, the second author and Todorov introduced an upper bound on the finitistic dimension of an Artin algebra, now known as the phi dimension. The phi dimension conjecture states that this upper bound is always finite, a fact that would imply the finitistic dimension conjecture. In this paper, we present a counterexample to the phi dimension conjecture and explain where it comes from. We also discuss implications for further research and the finitistic dimension conjecture.


Rupture of liquid bridges on porous tips Competing mechanisms of spontaneous imbibition and stretching. Liquid bridges are commonly encountered in nature and the liquid transfer induced by their rupture are widely used in various industrial applications. In this work, with the focus on the porous tip, we studied the impacts of capillary effects on the liquid transfer induced by the rupture through numerical simulations. To depict the capillary effects of a porous tip, a time scale ratio, R T, is proposed to compare the competing mechanisms of spontaneous imbibitiona and external drag. In terms of R T, we then develop a theoretical model for estimating the liquid retention ratio considering the geometry, porosity and wettability of tips. The mecahnism presented in this work provides a possible approach to control the liquid transfer with better accuracy in microfluidics or microfabrications.


A New Software Framework for Traffic Engineering Path Cardinality and the Effect of Multipath on Residual Capacity. In this paper, we present a new traffic engineering (TE) software framework to analyze, configure, and optimize (with the aid of a linear programming solver) a network for service provisioning. The developed software tool is based on our new data driven traffic engineering approach that analyzes a large volume of network configuration data generated given the user input. By analyzing the data, one can then make efficient decisions later when designing a traffic engineering solution. We focus on three well known traffic engineering objective functions minimum cost routing (MCR), load balancing (LB), and average delay (AD). With this new tool, one can answer numerous traffic engineering questions. For example, what are the differences among the three objective functions? What is the impact of an objective function on link utilization? How many candidate paths are enough to achieve optimality or near optimality with respect to a specific objective. This new software tool allows us to conveniently perform various experiments and visualize the results for performance analysis. As case studies, this paper presents examples that answer the questions for two traffic engineering problems (1) how many paths are required to obtain a solution that is within a few percent from the optimal solution and whether that number is fixed for any network size? (2) how the choice of single path multi path routing affects the load in the network? For the first problem, it turns out that the number of paths needed to achieve optimality increases as the number of links in the network increases.


A Tracking System For Baseball Game Reconstruction. The baseball game is often seen as many contests that are performed between individuals. The duel between the pitcher and the batter, for example, is considered the engine that drives the sport. The pitchers use a variety of strategies to gain competitive advantage against the batter, who does his best to figure out the ball trajectory and react in time for a hit. In this work, we propose a system that captures the movements of the pitcher, the batter, and the ball in a high level of detail, and discuss several ways how this information may be processed to compute interesting statistics. We demonstrate on a large database of videos that our methods achieve comparable results as previous systems, while operating solely on video material. In addition, state of the art AI techniques are incorporated to augment the amount of information that is made available for players, coaches, teams, and fans.


Quantitative estimates for almost constant mean curvature hypersurfaces. Alexandrov's soap bubble theorem asserts that spheres are the only connected closed embedded hypersurfaces in the Euclidean space with constant mean curvature. The theorem can be extended to space forms and it holds for more general functions of the principal curvatures. In this short review, we discuss quantitative stability results regarding Alexandrov's theorem which have been obtained by the author in recent years. In particular, we consider hypersurfaces having mean curvature close to a constant and we quantitatively describe the proximity to a single sphere or to a collection of tangent spheres in terms of the oscillation of the mean curvature. Moreover, we also consider the problem in a non local setting, and we show that the non local effect gives a stronger rigidity to the problem and prevents the appearance of bubbling.


Relative Drone Ground Vehicle Localization using LiDAR and Fisheye Cameras through Direct and Indirect Observations. Estimating the pose of an unmanned aerial vehicle (UAV) or drone is a challenging task. It is useful for many applications such as navigation, surveillance, tracking objects on the ground, and 3D reconstruction. In this work, we present a LiDAR camera based relative pose estimation method between a drone and a ground vehicle, using a LiDAR sensor and a fisheye camera on the vehicle's roof and another fisheye camera mounted under the drone. The LiDAR sensor directly observes the drone and measures its position, and the two cameras estimate the relative orientation using indirect observation of the surrounding objects. We propose a dynamically adaptive kernel based method for drone detection and tracking using the LiDAR. We detect vanishing points in both cameras and find their correspondences to estimate the relative orientation. Additionally, we propose a rotation correction technique by relying on the observed motion of the drone through the LiDAR. In our experiments, we were able to achieve very fast initial detection and real time tracking of the drone. Our method is fully automatic.


Linear conjugacy of chemical kinetic systems. Two networks are said to be linearly conjugate if the solution of their dynamic equations can be transformed into each other by a positive linear transformation. The study on dynamical equivalence in chemical kinetic systems was initiated by Craciun and Pantea in 2008 and eventually led to the Johnston Siegel Criterion for linear conjugacy (JSC). Several studies have applied Mixed Integer Linear Programming (MILP) approach to generate linear conjugates of MAK (mass action kinetic) systems, Bio CRNs (which is a subset of hill type kinetic systems when the network is restricted to digraphs), and PL RDK (complex factorizable power law kinetic) systems. In this study, we present a general computational solution to construct linear conjugates of any "rate constant interaction function decomposable" (RID) chemical kinetic systems, wherein each of its rate function is the product of a rate constant and an interaction function. We generate an extension of the JSC to the complex factorizable (CF) subset of RID kinetic systems and show that any non complex factorizable (NF) RID kinetic system can be dynamically equivalent to a CF system via transformation. We show that linear conjugacy can be generated for any RID kinetic systems by applying the JSC to any NF kinetic system that are transformed to CF kinetic system.


Active Learning for Event Detection in Support of Disaster Analysis Applications. Disaster analysis in social media content is one of the interesting research domains having abundance of data. However, there is a lack of labeled data that can be used to train machine learning models for disaster analysis applications. Active learning is one of the possible solutions to such problem. To this aim, in this paper we propose and assess the efficacy of an active learning based framework for disaster analysis using images shared on social media outlets. Specifically, we analyze the performance of different active learning techniques employing several sampling and disagreement strategies. Moreover, we collect a large scale dataset covering images from eight common types of natural disasters. The experimental results show that the use of active learning techniques for disaster analysis using images results in a performance comparable to that obtained using human annotated images, and could be used in frameworks for disaster analysis in images without tedious job of manual annotation.


ITSELF Iterative Saliency Estimation fLexible Framework. Saliency object detection estimates the objects that most stand out in an image. The available unsupervised saliency estimators rely on a pre determined set of assumptions of how humans perceive saliency to create discriminating features. By fixing the pre selected assumptions as an integral part of their models, these methods cannot be easily extended for specific settings and different image domains. We then propose a superpixel based ITerative Saliency Estimation fLexible Framework (ITSELF) that allows any user defined assumptions to be added to the model when required. Thanks to recent advancements in superpixel segmentation algorithms, saliency maps can be used to improve superpixel delineation. By combining a saliency based superpixel algorithm to a superpixel based saliency estimator, we propose a novel saliency superpixel self improving loop to iteratively enhance saliency maps. We compare ITSELF to two state of the art saliency estimators on five metrics and six datasets, four of which are composed of natural images, and two of biomedical images. Experiments show that our approach is more robust than the compared methods, presenting competitive results on natural image datasets and outperforming them on biomedical image datasets.


Reasoning Over Semantic Level Graph for Fact Checking. Fact checking is a challenging task because verifying the truthfulness of a claim requires reasoning about multiple retrievable evidence. In this work, we present a method suitable for reasoning about the semantic level structure of evidence. Unlike most previous works, which typically represent evidence sentences with either string concatenation or fusing the features of isolated evidence sentences, our approach operates on rich semantic structures of evidence obtained by semantic role labeling. We propose two mechanisms to exploit the structure of evidence while leveraging the advances of pre trained models like BERT, GPT or XLNet. Specifically, using XLNet as the backbone, we first utilize the graph structure to re define the relative distances of words, with the intuition that semantically related words should have short distances. Then, we adopt graph convolutional network and graph attention network to propagate and aggregate information from neighboring nodes on the graph. We evaluate our system on FEVER, a benchmark dataset for fact checking, and find that rich structural information is helpful and both our graph based mechanisms improve the accuracy. Our model is the state of the art system in terms of both official evaluation metrics, namely claim verification accuracy and FEVER score.


Illumination Adaptive Person Re identification. Most person re identification (ReID) approaches assume that person images are captured under relatively similar illumination conditions. In reality, long term person retrieval is common, and person images are often captured under different illumination conditions at different times across a day. In this situation, the performances of existing ReID models often degrade dramatically. This paper addresses the ReID problem with illumination variations and names it as em Illumination Adaptive Person Re identification (IA ReID) . We propose an Illumination Identity Disentanglement (IID) network to dispel different scales of illuminations away while preserving individuals' identity information. To demonstrate the illumination issue and to evaluate our model, we construct two large scale simulated datasets with a wide range of illumination variations. Experimental results on the simulated datasets and real world images demonstrate the effectiveness of the proposed framework.


Varieties of Regular Pseudocomplemented de Morgan Algebras. In this paper, we investigate the varieties mathbf M n and mathbf K n of regular pseudocomplemented de Morgan and Kleene algebras of range n , respectively. Priestley duality as it applies to pseudocomplemented de Morgan algebras is used. We characterise the dual spaces of the simple (equivalently, subdirectly irreducible) algebras in mathbf M n and explicitly describe the dual spaces of the simple algebras in mathbf M 1 and mathbf K 1 . We show that the variety mathbf M 1 is locally finite, but this property does not extend to mathbf M n or even mathbf K n for n geq 2 . We also show that the lattice of subvarieties of mathbf K 1 is an omega 1 chain and the cardinality of the lattice of subvarieties of either mathbf K 2 or mathbf M 1 is 2 omega . A description of the lattice of subvarieties of mathbf M 1 is given.


Categorical primitive forms of Calabi Yau A infty categories with semi simple cohomology. We study categorical primitive forms for Calabi Yau A infty categories with semi simple Hochschild cohomology. We classify these primitive forms in terms of certain grading operators on the Hochschild homology. We use this result to prove that, if the Fukaya category sf Fuk (M) of a symplectic manifold M has semi simple Hochschild cohomology, then its genus zero Gromov Witten invariants may be recovered from the A infty category sf Fuk (M) together with the closed open map. An immediate corollary of this is that in the semi simple case, homological mirror symmetry implies enumerative mirror symmetry.


A Deep Dive on the Impact of COVID 19 in Software Development. Context COVID 19 pandemic has impacted different business sectors around the world. Objective. This study investigates the impact of COVID 19 on software projects and software development professionals. Method We conducted a mining software repository study based on 100 GitHub projects developed in Java using ten different metrics. Next, we surveyed 279 software development professionals for better understanding the impact of COVID 19 on daily activities and wellbeing. Results We identified 12 observations related to productivity, code quality, and wellbeing. Conclusions Our findings highlight that the impact of COVID 19 is not binary (reduce productivity vs. increase productivity) but rather a spectrum. For many of our observations, substantial proportions of respondents have differing opinions from each other. We believe that more research is needed to uncover specific conditions that cause certain outcomes to be more prevalent.


Scattering from rough thin films DDA simulations. We investigate the wave optical light scattering properties of deformed thin circular films of constant thickness using the discrete dipole approximation. Effects on the intensity distribution of the scattered light due to different statistical roughness models, model dependant roughness parameters, and uncorrelated random small scale porosity of the inhomogeneous medium are studied. The usability of discrete dipole approximation to rough surface scattering problems is evaluated by considering thin films as computationally feasible rough surface analogs. The effects due to small scale inhomogeneity of the scattering medium are compared with the analytic approximation by Maxwell Garnett and the results are found to agree with the approximation.


Classifying Intrinsically Linked Tournaments by Score Sequence. A tournament on 8 or more vertices may be intrinsically linked as a directed graph. We begin the classification of intrinsically linked tournaments by examining their score sequences. While many distinct tournaments may have the same score sequence, there exist score sequences S such that any tournament with score sequence S has an embedding with no nonsplit consistently oriented link. We call such score sequences textit linkless , and we show that the vast majority of score sequences for 8 vertex tournaments are linkless. We also extend these results to n vertex tournaments and are able to classify many longer score sequences as well. We show that for any n , there exist at least O(n) linkless score sequences, but we conjecture that the fraction of score sequences of length n that are linkless goes to 0 as n becomes large.


A Review on Brain Mechanisms for Language Acquisition and Comprehension. This paper reviews the main perspectives of language acquisition and language comprehension. In language acquisition, we have reviewed the different types of language acquisitions like first language acquisition, second language acquisition, sign language acquisition and skill acquisition. The experimental techniques for neurolinguistic acquisition detection is also discussed. The findings of experiments for acquisition detection is also discussed, it includes the region of brain activated after acquisition. Findings shows that the different types of acquisition involve different regions of the brain. In language comprehension, native language comprehension and bilingual's comprehension has been considered. Comprehension involve different brain regions for different sentence or word comprehension depending upon their semantic and syntax. The different fMRIEEG analysis techniques (statistical or graph theoretical) are also discoursed in our review. Tools for neurolinguistics computations are also discussed.


The Decidability of Verification under Promising 2.0. In PLDI'20, Lee et al. introduced the emph promising semantics PS 2.0 of the C concurrency that captures most of the common program transformations while satisfying the DRF guarantee. The reachability problem for finite state programs under PS 2.0 with only release acquire accesses is already known to be undecidable. Therefore, we address, in this paper, the reachability problem for programs running under PS 2.0 with relaxed accesses together with promises. We show that this problem is undecidable even in the case where the input program has finite state. Given this undecidability result, we consider the fragment of PS 2.0 with only relaxed accesses allowing bounded number of promises. We show that under this restriction, the reachability is decidable, albeit very expensive it is non primitive recursive. Given this high complexity with bounded number of promises and the undecidability result for the RA fragment of PS 2.0, we consider a bounded version of the reachability problem. To this end, we bound both the number of promises and the "view switches", i.e, the number of times the processes may switch their local views of the global memory. We provide a code to code translation from an input program under PS 2.0, with relaxed and release acquire memory accesses along with promises, to a program under SC. This leads to a reduction of the bounded reachability problem under PS 2.0 to the bounded context switching problem under SC. We have implemented a prototype tool and tested it on a set of benchmarks, demonstrating that many bugs in programs can be found using a small bound.


DART A Lightweight Quality Suggestive Data to Text Annotation Tool. We present a lightweight annotation tool, the Data AnnotatoR Tool (DART), for the general task of labeling structured data with textual descriptions. The tool is implemented as an interactive application that reduces human efforts in annotating large quantities of structured data, e.g. in the format of a table or tree structure. By using a backend sequence to sequence model, our system iteratively analyzes the annotated labels in order to better sample unlabeled data. In a simulation experiment performed on annotating large quantities of structured data, DART has been shown to reduce the total number of annotations needed with active learning and automatically suggesting relevant labels.


AIBench An Industry Standard Internet Service AI Benchmark Suite. Today's Internet Services are undergoing fundamental changes and shifting to an intelligent computing era where AI is widely employed to augment services. In this context, many innovative AI algorithms, systems, and architectures are proposed, and thus the importance of benchmarking and evaluating them rises. However, modern Internet services adopt a microservice based architecture and consist of various modules. The diversity of these modules and complexity of execution paths, the massive scale and complex hierarchy of datacenter infrastructure, the confidential issues of data sets and workloads pose great challenges to benchmarking. In this paper, we present the first industry standard Internet service AI benchmark suite AIBench with seventeen industry partners, including several top Internet service providers. AIBench provides a highly extensible, configurable, and flexible benchmark framework that contains loosely coupled modules. We identify sixteen prominent AI problem domains like learning to rank, each of which forms an AI component benchmark, from three most important Internet service domains search engine, social network, and e commerce, which is by far the most comprehensive AI benchmarking effort. On the basis of the AIBench framework, abstracting the real world data sets and workloads from one of the top e commerce providers, we design and implement the first end to end Internet service AI benchmark, which contains the primary modules in the critical paths of an industry scale application and is scalable to deploy on different cluster scales. The specifications, source code, and performance numbers are publicly available from the benchmark council web site http www.benchcouncil.org AIBench index.html.


On the Cauchy problem of 3D incompressible Navier Stokes Cahn Hilliard system. In this paper, we are concerned with the well posedness and large time behavior of Cauchy problem for 3D incompressible Navier Stokes Cahn Hilliard equations. First, using Banach fixed point theorem, we establish the local well posedness of solutions. Second, assuming (u 0, nabla phi 0) dot H frac12 is sufficiently small, we obtain the global well posedness of solutions. Moreover, the optimal decay rates of the higher order spatial derivatives of the solution are also obtained.


On potential theory of Markov processes with jump kernels decaying at the boundary. Motivated by some recent potential theoretic results on subordinate killed L 'evy processes in open subsets of the Euclidean space, we study processes in an open set D subset mathbb R d defined via Dirichlet forms with jump kernels of the form J D(x,y) j( x y ) mathcal B (x,y) and critical killing functions. Here j( x y ) is the L 'evy density of an isotropic stable process (or more generally, a pure jump isotropic unimodal L 'evy process) in mathbb R d . The main novelty is that the term mathcal B (x,y) tends to 0 when x or y approach the boundary of D . Under some general assumptions on mathcal B (x,y) , we construct the corresponding process and prove that non negative harmonic functions of the process satisfy the Harnack inequality and Carleson's estimate. We give several examples of boundary terms satisfying those assumptions. The examples depend on four parameters, beta 1, beta 2, beta 3 , beta 4 , roughly governing the decay of the boundary term near the boundary of D . In the second part of this paper, we specialise to the case of the half space D mathbb R d x ( widetilde x ,x d) , x d 0 , the alpha stable kernel j( x y ) x y d alpha and the killing function kappa(x) c x d alpha , alpha in (0,2) , where c is a positive constant. Our main result in this part is a boundary Harnack principle which says that, for any p ( alpha 1) , there are values of the parameters beta 1, beta 2, beta 3 , beta 4 , and the constant c such that non negative harmonic functions of the process must decay at the rate x d p if they vanish near a portion of the boundary. We further show that there are values of the parameters beta 1, beta 2, beta 3 , beta 4 , for which the boundary Harnack principle fails despite the fact that Carleson's estimate is valid.


Moving object detection for visual odometry in a dynamic environment based on occlusion accumulation. Detection of moving objects is an essential capability in dealing with dynamic environments. Most moving object detection algorithms have been designed for color images without depth. For robotic navigation where real time RGB D data is often readily available, utilization of the depth information would be beneficial for obstacle recognition. Here, we propose a simple moving object detection algorithm that uses RGB D images. The proposed algorithm does not require estimating a background model. Instead, it uses an occlusion model which enables us to estimate the camera pose on a background confused with moving objects that dominate the scene. The proposed algorithm allows to separate the moving object detection and visual odometry (VO) so that an arbitrary robust VO method can be employed in a dynamic situation with a combination of moving object detection, whereas other VO algorithms for a dynamic environment are inseparable. In this paper, we use dense visual odometry (DVO) as a VO method with a bi square regression weight. Experimental results show the segmentation accuracy and the performance improvement of DVO in the situations. We validate our algorithm in public datasets and our dataset which also publicly accessible.


cs.AI cs.LG physics.flu-dyn physics.soc-ph

Deep Reinforcement Learning for Real Time Optimization of Pumps in Water Distribution Systems. Real time control of pumps can be an infeasible task in water distribution systems (WDSs) because the calculation to find the optimal pump speeds is resource intensive. The computational need cannot be lowered even with the capabilities of smart water networks when conventional optimization techniques are used. Deep reinforcement learning (DRL) is presented here as a controller of pumps in two WDSs. An agent based on a dueling deep q network is trained to maintain the pump speeds based on instantaneous nodal pressure data. General optimization techniques (e.g., Nelder Mead method, differential evolution) serve as baselines. The total efficiency achieved by the DRL agent compared to the best performing baseline is above 0.98, whereas the speedup is around 2x compared to that. The main contribution of the presented approach is that the agent can run the pumps in real time because it depends only on measurement data. If the WDS is replaced with a hydraulic simulation, the agent still outperforms conventional techniques in search speed.


Hydrostatic approximation of the 2D primitive equations in a thin strip. We prove the global wellposedness of the 2D non rotating primitive equations with no slip boundary conditions in a thin strip of width eps for small data which are analytic in the tangential direction. We also prove that the hydrostatic limit (when eps to 0 ) is a couple of a Prandtl like system for the velocity with a transport diffusion equation for the temperature.


On reconstructing subvarieties from their periods. We give a new practical method for computing subvarieties of projective hypersurfaces. By computing the periods of a given hypersurface X, we find algebraic cohomology cycles on X. On well picked algebraic cycles, we can then recover the equations of subvarieties of X that realize these cycles. In practice, a bulk of the computations involve transcendental numbers and have to be carried out with floating point numbers. However, if X is defined over algebraic numbers then the coefficients of the equations of subvarieties can be reconstructed as algebraic numbers. A symbolic computation then verifies the results. As an illustration of the method, we compute generators of the Picard groups of some quartic surfaces. A highlight of the method is that the Picard group computations are proved to be correct despite the fact that the Picard numbers of our examples are not extremal.


MRI Brain Tumor Segmentation using Random Forests and Fully Convolutional Networks. In this paper, we propose a novel learning based method for automated segmentation of brain tumor in multimodal MRI images, which incorporates two sets of machine learned and hand crafted features. Fully convolutional networks (FCN) forms the machine learned features and texton based features are considered as hand crafted features. Random forest (RF) is used to classify the MRI image voxels into normal brain tissues and different parts of tumors, i.e. edema, necrosis and enhancing tumor. The method was evaluated on BRATS 2017 challenge dataset. The results show that the proposed method provides promising segmentations. The mean Dice overlap measure for automatic brain tumor segmentation against ground truth is 0.86, 0.78 and 0.66 for whole tumor, core and enhancing tumor, respectively.


New Solutions to the Reflection Equation with Braces. The reflection equation of Cherednik is a counterpart to the celebrated Yang Baxter equation, with importance in the theory of integrable systems. We obtain several new solutions of the reflection equation using braces building on the work of Smoktunowicz, Vendramin and Weston. In particular, we show that every brace yields several simple solutions and that a natural class of solutions is a near ring. We find more solutions for factorizable rings. Some of our solutions apply to the original parameter dependent equation.


Balancing Accuracy and Diversity in Recommendations using Matrix Completion Framework. Design of recommender systems aimed at achieving high prediction accuracy is a widely researched area. However, several studies have suggested the need for diversified recommendations, with acceptable level of accuracy, to avoid monotony and improve customers experience. However, increasing diversity comes with an associated reduction in recommendation accuracy thereby necessitating an optimum tradeoff between the two. In this work, we attempt to achieve accuracy vs diversity balance, by exploiting available ratings and item metadata, through a single, joint optimization model built over the matrix completion framework. Most existing works, unlike our formulation, propose a 2 stage model, a heuristic item ranking scheme on top of an existing collaborative filtering technique. Experimental evaluation on a movie recommender system indicates that our model achieves higher diversity for a given drop in accuracy as compared to existing state of the art techniques.


Uniqueness, Comparison and Stability for Scalar BSDEs with Lexp( mu sqrt(2log(1 L))) integrable terminal values and monotonic generators. This paper considers a class of scalar backward stochastic differential equations (BSDEs) with L exp( mu sqrt 2 log(1 L) ) integrable terminal values. We associate these BSDEs with BSDEs with integrable parameters through Girsanov change. Using this technique, we prove uniqueness, comparisons and stability for them under an extended monotonicity condition (more precisely one sided Osgood condition).


T VSE Transformer Based Visual Semantic Embedding. Transformer models have recently achieved impressive performance on NLP tasks, owing to new algorithms for self supervised pre training on very large text corpora. In contrast, recent literature suggests that simple average word models outperform more complicated language models, e.g., RNNs and Transformers, on cross modal image text search tasks on standard benchmarks, like MS COCO. In this paper, we show that dataset scale and training strategy are critical and demonstrate that transformer based cross modal embeddings outperform word average and RNN based embeddings by a large margin, when trained on a large dataset of e commerce product image title pairs.


Enriched functional limit theorems for dynamical systems. We prove functional limit theorems for stochastic processes which have clusters of large values which, when summed and suitably normalised in time and space, get collapsed in a jump of the limiting process observed at the same time point. In order to keep track of the clustering information, which gets lost in the usual Skorohod topologies in the space of c adl ag functions, we introduce a new space specially designed for that purpose which generalises the already more general spaces introduced by Whitt. Our main applications are to hyperbolic and non uniformly expanding dynamical systems with heavy tailed observable functions maximised at dynamically linked maximal sets (such as periodic points). We also study limits of extremal processes and record times point processes for observables not necessarily heavy tailed. The applications studied include hyperbolic systems such as Anosov diffeomorphisms, but also non uniformly expanding maps such as maps with critical points of Benedicks Carleson type or indifferent fixed points such as Manneville Pomeau maps. The main tool is a limit theorem for point processes with decorations derived from a bi infinite sequence which we call the piling process, designed to keep the clustering information. This piling process is based on the tail process defined in the regularly varying context rather than the general stationary stochastic processes we consider here.


Jindalrae and Gaenari numbers and polynomials. The aim of this paper is to study Jindalrae and Gaenari numbers and polynomials in connection with Jindalrae Stirling numbers of the first and second kinds. For this purpose, we first introduce Jindalrae Stirling numbers of the first and second kinds as extensions of the notions of the degenerate Stirling numbers of the first and second kinds, and deduce several relations involving those special numbers. Then we introduce Jindalrae and Gaenari numbers and polynomials and obtain some explicit expressions and identities associated with those numbers and polynomials. In addition, we interpret our results by using umbral calculus.


What the MASK ? Making Sense of Language Specific BERT Models. Recently, Natural Language Processing (NLP) has witnessed an impressive progress in many areas, due to the advent of novel, pretrained contextual representation models. In particular, Devlin et al. (2019) proposed a model, called BERT (Bidirectional Encoder Representations from Transformers), which enables researchers to obtain state of the art performance on numerous NLP tasks by fine tuning the representations on their data set and task, without the need for developing and training highly specific architectures. The authors also released multilingual BERT (mBERT), a model trained on a corpus of 104 languages, which can serve as a universal language model. This model obtained impressive results on a zero shot cross lingual natural inference task. Driven by the potential of BERT models, the NLP community has started to investigate and generate an abundant number of BERT models that are trained on a particular language, and tested on a specific data domain and task. This allows us to evaluate the true potential of mBERT as a universal language model, by comparing it to the performance of these more specific models. This paper presents the current state of the art in language specific BERT models, providing an overall picture with respect to different dimensions (i.e. architectures, data domains, and tasks). Our aim is to provide an immediate and straightforward overview of the commonalities and differences between Language Specific (language specific) BERT models and mBERT. We also provide an interactive and constantly updated website that can be used to explore the information we have collected, at https bertlang.unibocconi.it.


Fast Perception, Planning, and Execution for a Robotic Butler Wheeled Humanoid M Hubo. As the aging population grows at a rapid rate, there is an ever growing need for service robot platforms that can provide daily assistance at practical speed with reliable performance. In order to assist with daily tasks such as fetching a beverage, a service robot must be able to perceive its environment and generate corresponding motion trajectories. This becomes a challenging and computationally complex problem when the environment is unknown and thus the path planner must sample numerous trajectories that often are sub optimal, extending the execution time. To address this issue, we propose a unique strategy of integrating a 3D object detection pipeline with a kinematically optimal manipulation planner to significantly increase speed performance at runtime. In addition, we develop a new robotic butler system for a wheeled humanoid that is capable of fetching requested objects at 24 of the speed a human needs to fulfill the same task. The proposed system was evaluated and demonstrated in a real world environment setup as well as in public exhibition.


Factorization in monoids and rings. Let H times be the group of units of a multiplicatively written monoid H . We say H is acyclic if xyz ne y for all x, y, z in H with x notin H times or z notin H times unit cancellative if yx ne x ne xy for all x, y in H with y notin H times f.g.u. if there is a finite set A subseteq H such that every non unit of H is a finite product of elements of the form uav with u, v in H times and a in A l.f.g.u. if, for each x in H , the smallest divisor closed submonoid of H containing x is f.g.u and atomic if every non unit can be written as a finite product of atoms, where an atom is a non unit that does not factor into a product of two non units. We generalize to l.f.g.u. or acyclic l.f.g.u. monoids a few results so far only known for unit cancellative l.f.g.u. commutative monoids (cancellative monoids are unit cancellative, and a commutative monoid is unit cancellative if and only if it is acyclic). In particular, we prove the following bullet If H is an atomic l.f.g.u. monoid, then every non unit has only finitely many factorizations (into atoms) that are "minimal" and "pairwise non equivalent" (with respect to some naturally defined relations on the free monoid over the "alphabet" of atoms). bullet If H is an acyclic l.f.g.u. monoid, then it is atomic and moreover, each element has only finitely many "pairwise non equivalent" factorizations if we additionally assume H to be commutative.


New Directions in Automated Traffic Analysis. Despite the use of machine learning for many network traffic analysis tasks in security, from application identification to intrusion detection, the aspects of the machine learning pipeline that ultimately determine the performance of the model feature selection and representation, model selection, and parameter tuning remain manual and painstaking. This paper presents a method to automate many aspects of traffic analysis, making it easier to apply machine learning techniques to a wider variety of traffic analysis tasks. We introduce nPrint, a tool that generates a unified packet representation that is amenable for representation learning and model training. We integrate nPrint with automated machine learning (AutoML), resulting in nPrintML, a public system that largely eliminates feature extraction and model tuning for a wide variety of traffic analysis tasks. We have evaluated nPrintML on eight separate traffic analysis tasks and released nPrint and nPrintML to enable future work to extend these methods.


Numerical Methods for Flow in Fractured Porous Media. In this work we present the mathematical models for single phase flow in fractured porous media. An overview of the most common approaches is considered, which includes continuous fracture models and discrete fracture models. For the latter, we discuss strategies that are developed in literature for its numerical solution mainly related to the geometrical relation between the fractures and porous media grids.