Design of a freeform two reflector system to collimate and shape a point source distribution. In this paper we propose a method to compute a freeform reflector system for collimating and shaping a beam from a point source. We construct these reflectors such that the radiant intensity of the source is converted into a desired target. An important generalization in our approach compared to previous research is that the output beam can be in an arbitrary direction. The design problem is approached by using a generalized Monge Amp ere equation. This equation is solved using a least squares algorithm for non quadratic cost functions. This algorithm calculates the optical map, from which we can then compute the surfaces. We test our algorithm on two cases. First we consider a uniform source and target distribution. Next, we use the model of a laser diode light source and a ring shaped target distribution.


Computing Diverse Shortest Paths Efficiently A Theoretical and Experimental Study. Finding diverse solutions in combinatorial problems recently has received considerable attention (Baste et al. 2020 Fomin et al. 2020 Hanaka et al. 2021). In this paper we study the following type of problems given an integer k , the problem asks for k solutions such that the sum of pairwise (weighted) Hamming distances between these solutions is maximized. Such solutions are called diverse solutions. We present a polynomial time algorithm for finding diverse shortest st paths in weighted directed graphs. Moreover, we study the diverse version of other classical combinatorial problems such as diverse weighted matroid bases, diverse weighted arborescences, and diverse bipartite matchings. We show that these problems can be solved in polynomial time as well. To evaluate the practical performance of our algorithm for finding diverse shortest st paths, we conduct a computational experiment with synthetic and real world instances.The experiment shows that our algorithm successfully computes diverse solutions within reasonable computational time.


Series reversion in Calder 'on's problem. This work derives explicit series reversions for the solution of Calder 'on's problem. The governing elliptic partial differential equation is nabla cdot(A nabla u) 0 in a bounded Lipschitz domain and with a matrix valued coefficient. The corresponding forward map sends A to a projected version of a local Neumann to Dirichlet operator, allowing for the use of partial boundary data and finitely many measurements. It is first shown that the forward map is analytic, and subsequently reversions of its Taylor series up to specified orders lead to a family of numerical methods for solving the inverse problem with increasing accuracy. The convergence of these methods is shown under conditions that ensure the invertibility of the Fr 'echet derivative of the forward map. The introduced numerical methods are of the same computational complexity as solving the linearised inverse problem. The analogous results are also presented for the smoothened complete electrode model.


Self Supervised Inference in State Space Models. We perform approximate inference in state space models with nonlinear state transitions. Without parameterizing a generative model, we apply Bayesian update formulas using a local linearity approximation parameterized by neural networks. This comes accompanied by a maximum likelihood objective that requires no supervision via uncorrupt observations or ground truth latent states. The optimization backpropagates through a recursion similar to the classical Kalman filter and smoother. Additionally, using an approximate conditional independence, we can perform smoothing without having to parameterize a separate model. In scientific applications, domain knowledge can give a linear approximation of the latent transition maps, which we can easily incorporate into our model. Usage of such domain knowledge is reflected in excellent results (despite our model's simplicity) on the chaotic Lorenz system compared to fully supervised and variational inference methods. Finally, we show competitive results on an audio denoising experiment.


Wave Analysis and Homogenization of Spatiotemporally Modulated Wire Medium. In this paper we develop homogenization theory for spatiotemporally modulated wire medium. We first solve for the modal waves that are supported by this composite medium, we show peculiar properties such as extraordinary waves that propagate at frequencies below the cut off frequency of the corresponding stationary medium. We explain how these unique solutions give rise to an extreme Fresnel drag that exists already with weak and slow spatiotemporal modulation. Next, we turn to derive the effective material permittivity that corresponds to each of the first few supported modes, and write the average fields and Poynting's vector. Nonlocality, nonreciprocity, and anisotropy due to the spatiotemporal modulation direction, are three inherent properties of this medium, and are clearly seen in the effective material parameters. As a figure of merit, we also derive the effective permittivity of a plasma medium with spatiotemporally modulated plasma frequency. This comparison is interesting since the plasma medium can be considered as the effective medium that is obtained by a stationary wire medium. We validate that homogenization and spatiotemporal variation are not necessarily interchangeable operations. And indeed, in certain parameter regimes the homogenization should be performed directly on spatiotemporally modulated composite medium, rather than first homogenize the stationary medium and then phenomenologically introduce the effect of the space time modulation.


Tuna A Static Analysis Approach to Optimizing Deep Neural Networks. We introduce Tuna, a static analysis approach to optimizing deep neural network programs. The optimization of tensor operations such as convolutions and matrix multiplications is the key to improving the performance of deep neural networks. Many deep learning model optimization mechanisms today use dynamic analysis, which relies on experimental execution on a target device to build a data driven cost model of the program. The reliance on dynamic profiling not only requires access to target hardware at compilation time but also incurs significant cost in machine resources. We introduce an approach that profiles the program by constructing features based on the target hardware characteristics in order. We use static analysis of the relative performance of tensor operations to optimize the deep learning program. Experiments show that our approach can achieve up to 11x performance compared to dynamic profiling based methods with the same compilation time.


Active Learning for Deep Neural Networks on Edge Devices. When dealing with deep neural network (DNN) applications on edge devices, continuously updating the model is important. Although updating a model with real incoming data is ideal, using all of them is not always feasible due to limits, such as labeling and communication costs. Thus, it is necessary to filter and select the data to use for training (i.e., active learning) on the device. In this paper, we formalize a practical active learning problem for DNNs on edge devices and propose a general task agnostic framework to tackle this problem, which reduces it to a stream submodular maximization. This framework is light enough to be run with low computational resources, yet provides solutions whose quality is theoretically guaranteed thanks to the submodular property. Through this framework, we can configure data selection criteria flexibly, including using methods proposed in previous active learning studies. We evaluate our approach on both classification and object detection tasks in a practical setting to simulate a real life scenario. The results of our study show that the proposed framework outperforms all other methods in both tasks, while running at a practical speed on real devices.


How Should Pre Trained Language Models Be Fine Tuned Towards Adversarial Robustness?. The fine tuning of pre trained language models has a great success in many NLP fields. Yet, it is strikingly vulnerable to adversarial examples, e.g., word substitution attacks using only synonyms can easily fool a BERT based sentiment analysis model. In this paper, we demonstrate that adversarial training, the prevalent defense technique, does not directly fit a conventional fine tuning scenario, because it suffers severely from catastrophic forgetting failing to retain the generic and robust linguistic features that have already been captured by the pre trained model. In this light, we propose Robust Informative Fine Tuning (RIFT), a novel adversarial fine tuning method from an information theoretical perspective. In particular, RIFT encourages an objective model to retain the features learned from the pre trained model throughout the entire fine tuning process, whereas a conventional one only uses the pre trained weights for initialization. Experimental results show that RIFT consistently outperforms the state of the arts on two popular NLP tasks sentiment analysis and natural language inference, under different attacks across various pre trained language models.


CLAMGen Closed Loop Arm Motion Generation via Multi view Vision Based RL. We propose a vision based reinforcement learning (RL) approach for closed loop trajectory generation in an arm reaching problem. Arm trajectory generation is a fundamental robotics problem which entails finding collision free paths to move the robot's body (e.g. arm) in order to satisfy a goal (e.g. place end effector at a point). While classical methods typically require the model of the environment to solve a planning, search or optimization problem, learning based approaches hold the promise of directly mapping from observations to robot actions. However, learning a collision avoidance policy using RL remains a challenge for various reasons, including, but not limited to, partial observability, poor exploration, low sample efficiency, and learning instabilities. To address these challenges, we present a residual RL method that leverages a greedy goal reaching RL policy as the base to improve exploration, and the base policy is augmented with residual state action values and residual actions learned from images to avoid obstacles. Further more, we introduce novel learning objectives and techniques to improve 3D understanding from multiple image views and sample efficiency of our algorithm. Compared to RL baselines, our method achieves superior performance in terms of success rate.


Knowledge Grounded Conversational Symptom Detection with Graph Memory Networks. In this work, we propose a novel goal oriented dialog task, automatic symptom detection. We build a system that can interact with patients through dialog to detect and collect clinical symptoms automatically, which can save a doctor's time interviewing the patient. Given a set of explicit symptoms provided by the patient to initiate a dialog for diagnosing, the system is trained to collect implicit symptoms by asking questions, in order to collect more information for making an accurate diagnosis. After getting the reply from the patient for each question, the system also decides whether current information is enough for a human doctor to make a diagnosis. To achieve this goal, we propose two neural models and a training pipeline for the multi step reasoning task. We also build a knowledge graph as additional inputs to further improve model performance. Experiments show that our model significantly outperforms the baseline by 4 , discovering 67 of implicit symptoms on average with a limited number of questions.


From Distributed Machine Learning to Federated Learning A Survey. In recent years, data and computing resources are typically distributed in the devices of end users, various regions or organizations. Because of laws or regulations, the distributed data and computing resources cannot be directly shared among different regions or organizations for machine learning tasks. Federated learning emerges as an efficient approach to exploit distributed data and computing resources, so as to collaboratively train machine learning models, while obeying the laws and regulations and ensuring data security and data privacy. In this paper, we provide a comprehensive survey of existing works for federated learning. We propose a functional architecture of federated learning systems and a taxonomy of related techniques. Furthermore, we present the distributed training, data communication, and security of FL systems. Finally, we analyze their limitations and propose future research directions.


Post Triangular Rewiring Method for Shorter RRT Robot Path Planning. This paper proposed the 'Post Triangular Rewiring' method that minimizes the sacrifice of planning time and overcomes the limit of Optimality of sampling based algorithm such as Rapidly exploring Random Tree (RRT) algorithm. The proposed 'Post Triangular Rewiring' method creates a closer to the optimal path than RRT algorithm before application through the triangular inequality principle. The experiments were conducted to verify a performance of the proposed method. When the method proposed in this paper are applied to the RRT algorithm, the Optimality efficiency increase compared to the planning time.


On the Value of ML Models. We argue that, when establishing and benchmarking Machine Learning (ML) models, the research community should favour evaluation metrics that better capture the value delivered by their model in practical applications. For a specific class of use cases selective classification we show that not only can it be simple enough to do, but that it has import consequences and provides insights what to look for in a good'' ML model.


On the absoluteness of aleph 1 freeness. aleph 1 free groups, abelian groups for which every countable subgroup is free, exhibit a number of interesting algebraic and set theoretic properties. In this paper, we give a complete proof that the property of being aleph 1 free is absolute that is, if an abelian group G is aleph 1 free in some transitive model textbf M of ZFC, then it is aleph 1 free in any transitive model of ZFC containing G . The absoluteness of aleph 1 freeness has the following remarkable consequence an abelian group G is aleph 1 free in some transitive model of ZFC if and only if it is (countable and) free in some model extension. This set theoretic characterization will be the starting point for further exploring the relationship between the set theoretic and algebraic properties of aleph 1 free groups. In particular, this paper will demonstrate how proofs may be dramatically simplified using model extensions for aleph 1 free groups.


SuMo A Mutation Testing Strategy for Solidity Smart Contracts. Smart Contracts are software programs that are deployed and executed within a blockchain infrastructure. Due to their immutable nature, directly resulting from the specific characteristics of the deploying infrastructure, smart contracts must be thoroughly tested before their release. Testing is one of the main activities that can help to improve the reliability of a smart contract, so as to possibly prevent considerable loss of valuable assets. It is therefore important to provide the testers with tools that permit them to assess the activity they performed. Mutation testing is a powerful approach for assessing the fault detection capability of a test suite. In this paper, we propose SuMo, a novel mutation testing tool for Ethereum Smart Contracts. SuMo implements a set of 44 mutation operators that were designed starting from the latest Solidity documentation, and from well known mutation testing tools. These allow to simulate a wide variety of faults that can be made by smart contract developers. The set of operators was designed to limit the generation of stillborn mutants, which slow down the mutation testing process and limit the usability of the tool. We report a first evaluation of SuMo on open source projects for which test suites were available. The results we got are encouraging, and they suggest that SuMo can effectively help developers to deliver more reliable smart contracts.


Collective management of environmental commons with multiple usages a guaranteed viability approach. In this paper we address the collective management of environmental commons with multiple usages in the framework of the mathematical viability theory. We consider that the stakeholders can derive from the study of their own socioeconomic problem the variables describing their different usages of the commons and its evolution, and a representation of the desirable states for the commons. We then consider the guaranteed viability kernel, subset of the set of desirable states where it is possible to maintain the state of the commons even when its evolution is represented by several conflicting models. This approach is illustrated on a problem of lake eutrophication.


Partition identities from higher level crystals of A 1 (1) . We study perfect crystals for the standard modules of the affine Lie algebra A 1 (1) at all levels using the theory of multi grounded partitions. We prove a family of partition identities which are reminiscent of the Andrews Gordon identities and companions to the Meurman Primc identities, but with simple difference conditions involving absolute values. We also give simple non specialised character formulas with obviously positive coefficients for the three level 2 standard modules.


Cutting and Pasting in the Torelli subgroup of Out( F n ). Using ideas from 3 manifolds, Hatcher Wahl defined a notion of automorphism groups of free groups with boundary. We study their Torelli subgroups, adapting ideas introduced by Putman for surface mapping class groups. Our main results show that these groups are finitely generated, and also that they satisfy an appropriate version of the Birman exact sequence.


Stacked LSTM Based Deep Recurrent Neural Network with Kalman Smoothing for Blood Glucose Prediction. Blood glucose (BG) management is crucial for type 1 diabetes patients resulting in the necessity of reliable artificial pancreas or insulin infusion systems. In recent years, deep learning techniques have been utilized for a more accurate BG level prediction system. However, continuous glucose monitoring (CGM) readings are susceptible to sensor errors. As a result, inaccurate CGM readings would affect BG prediction and make it unreliable, even if the most optimal machine learning model is used. In this work, we propose a novel approach to predicting blood glucose level with a stacked Long short term memory (LSTM) based deep recurrent neural network (RNN) model considering sensor fault. We use the Kalman smoothing technique for the correction of the inaccurate CGM readings due to sensor error. For the OhioT1DM dataset, containing eight weeks' data from six different patients, we achieve an average RMSE of 6.45 and 17.24 mg dl for 30 minutes and 60 minutes of prediction horizon (PH), respectively. To the best of our knowledge, this is the leading average prediction accuracy for the ohioT1DM dataset. Different physiological information, e.g., Kalman smoothed CGM data, carbohydrates from the meal, bolus insulin, and cumulative step counts in a fixed time interval, are crafted to represent meaningful features used as input to the model. The goal of our approach is to lower the difference between the predicted CGM values and the fingerstick blood glucose readings the ground truth. Our results indicate that the proposed approach is feasible for more reliable BG forecasting that might improve the performance of the artificial pancreas and insulin infusion system for T1D diabetes management.


Heat kernel bounds for parabolic equations with singular (form bounded) vector fields. We consider Kolmogorov operator nabla cdot a cdot nabla b cdot nabla with measurable uniformly elliptic matrix a and prove Gaussian lower and upper bounds on its heat kernel under minimal assumptions on the vector field b and its divergence rm div , b . More precisely, we prove (1) Gaussian lower bound, provided that rm div , b geq 0 , and b is in the class of form bounded vector fields (containing e.g. ,the class L d , the weak L d class, as well as some vector fields that are not even in L rm loc 2 varepsilon , varepsilon 0 ) in these assumptions, the Gaussian upper bound is in general invalid (2) Gaussian upper bound, provided that b is form bounded, and the positive part of rm div , b is in the Kato class in these assumptions, the Gaussian lower bound is in general invalid (3) Gaussian upper and lower bounds, provided that b is form bounded, rm div , b is in the Kato class (4) A priori Gaussian upper and lower bounds, provided that b is in a large class containing the class of form bounded vector fields, rm div , b is in the Kato class.


Mathematics and map drawing in the eighteenth century. We consider the mathematical theory of geographical maps, with an emphasis on the eighteenth century works of Euler, Lagrange and Delisle. This period is characterized by the frequent use of maps that are no more obtained by the stereographic projection or its variations, but by much more general maps from the sphere to the plane. More especially, the characteristics of the desired geographical maps were formulated in terms of an appropriate choice of the images of the parallels and meridians, and the mathematical properties required by the map concern the distortion of the maps restricted to these lines. The paper also contains some notes on the general use of mathematical methods in cartography in Greek Antiquity and in the modern period, and on the mutual influence of the two fields, mathematics and geography. The final version of this paper will appear in Ganita Bharati (Indian mathematics).


High order finite difference Hermite WENO fixed point fast sweeping method for static Hamilton Jacobi equations. In this paper, we combine the nonlinear HWENO reconstruction in cite newhwenozq and the fixed point iteration with Gauss Seidel fast sweeping strategy, to solve the static Hamilton Jacobi equations in a novel HWENO framework recently developed in cite mehweno1 . The proposed HWENO frameworks enjoys several advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function phi . They are now computed from the current value of phi and the previous spatial derivatives of phi . This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO scheme. In addition, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed point fast sweeping method is used to update the numerical approximation. It is an explicit method and does not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some known issues, including that the iterative number is very sensitive to the parameter varepsilon used in the nonlinear weights, as observed in previous studies. Finally, in order to further reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two dimensional problems, which demonstrate the good performance of the proposed fixed point fast sweeping HWENO methods.


Automatic Assembly Planning based on Digital Product Descriptions. This paper proposes a new concept in which a digital twin derived from a digital product description will automatically perform assembly planning and orchestrate the production resources in a manufacturing cell. Thus the manufacturing cell has generic services with minimal assumptions about what kind of product will be assembled, while the digital product description is designed collaboratively between the designer at an OEM and automated services at potential manufacturers. This has several advantages. Firstly, the resulting versatile manufacturing facility can handle a broad variety of products with minimal or no reconfiguration effort, so it can cost effectively offer its services to a large number of OEMs. Secondly, a solution is presented to the problem of performing concurrent product design and assembly planning over the organizational boundary. Thirdly, the product design at the OEM is not constrained to the capabilities of specific manufacturing facilities. The concept is presented in general terms in UML and an implementation is provided in a 3D simulation environment using Automation Markup Language for digital product descriptions. Finally, two case studies are presented and applications in a real industrial context are discussed.


A Critical Comparison on Six Static Analysis Tools Detection, Agreement, and Precision. Background. Developers use Automated Static Analysis Tools (ASATs) to control for potential quality issues in source code, including defects and technical debt. Tool vendors have devised quite a number of tools, which makes it harder for practitioners to select the most suitable one for their needs. To better support developers, researchers have been conducting several studies on ASATs to favor the understanding of their actual capabilities. Aims. Despite the work done so far, there is still a lack of knowledge regarding (1) which source quality problems can actually be detected by static analysis tool warnings, (2) what is their agreement, and (3) what is the precision of their recommendations. We aim at bridging this gap by proposing a large scale comparison of six popular static analysis tools for Java projects Better Code Hub, CheckStyle, Coverity Scan, Findbugs, PMD, and SonarQube. Method. We analyze 47 Java projects and derive a taxonomy of warnings raised by 6 state of the practice ASATs. To assess their agreement, we compared them by manually analyzing at line level whether they identify the same issues. Finally, we manually evaluate the precision of the tools. Results. The key results report a comprehensive taxonomy of ASATs warnings, show little to no agreement among the tools and a low degree of precision. Conclusions. We provide a taxonomy that can be useful to researchers, practitioners, and tool vendors to map the current capabilities of the tools. Furthermore, our study provides the first overview on the agreement among different tools as well as an extensive analysis of their precision.


On sums of arithmetic functions involving the greatest common divisor. Let gcd(d 1 , ldots,d k ) be the greatest common divisor of the positive integers d 1 , ldots,d k , for any integer k geq 2 , and let tau and mu denote the divisor function and the M " o bius function, respectively. For an arbitrary arithmetic function g and for any real number x 5 and any integer k geq 3 , we define the sum S g,k (x) sum n leq x sum d 1 cdots d k n g( gcd(d 1 , ldots,d k )) In this paper, we give asymptotic formulas for S tau,k (x) and S mu,k (x) for k geq 3 .


Museum Painting Retrieval. To retrieve images based on their content is one of the most studied topics in the field of computer vision. Nowadays, this problem can be addressed using modern techniques such as feature extraction using machine learning, but over the years different classical methods have been developed. In this paper, we implement a query by example retrieval system for finding paintings in a museum image collection using classic computer vision techniques. Specifically, we study the performance of the color, texture, text and feature descriptors in datasets with different perturbations in the images noise, overlapping text boxes, color corruption and rotation. We evaluate each of the cases using the Mean Average Precision (MAP) metric, and we obtain results that vary between 0.5 and 1.0 depending on the problem conditions.


A DCT based Tensor Completion Approach for Recovering Color Images and Videos from Highly Undersampled Data. Recovering color images and videos from highly undersampled data is a fundamental and challenging task in face recognition and computer vision. By the multi dimensional nature of color images and videos, in this paper, we propose a novel tensor completion approach, which is able to efficiently explore the sparsity of tensor data under the discrete cosine transform (DCT). Specifically, we introduce two DCT based tensor completion models as well as two implementable algorithms for their solutions. The first one is a DCT based weighted nuclear norm minimization model. The second one is called DCT based p shrinking tensor completion model, which is a nonconvex model utilizing p shrinkage mapping for promoting the low rankness of data. Moreover, we accordingly propose two implementable augmented Lagrangian based algorithms for solving the underlying optimization models. A series of numerical experiments including color and MRI image inpainting and video data recovery demonstrate that our proposed approach performs better than many existing state of the art tensor completion methods, especially for the case when the ratio of missing data is high.


Application of deep reinforcement learning for Indian stock trading automation. In stock trading, feature extraction and trading strategy design are the two important tasks to achieve long term benefits using machine learning techniques. Several methods have been proposed to design trading strategy by acquiring trading signals to maximize the rewards. In the present paper the theory of deep reinforcement learning is applied for stock trading strategy and investment decisions to Indian markets. The experiments are performed systematically with three classical Deep Reinforcement Learning models Deep Q Network, Double Deep Q Network and Dueling Double Deep Q Network on ten Indian stock datasets. The performance of the models are evaluated and comparison is made.


Exploiting Expert guided Symmetry Detection in Markov Decision Processes. Offline estimation of the dynamical model of a Markov Decision Process (MDP) is a non trivial task that greatly depends on the data available to the learning phase. Sometimes the dynamics of the model is invariant with respect to some transformations of the current state and action. Recent works showed that an expert guided pipeline relying on Density Estimation methods as Deep Neural Network based Normalizing Flows effectively detects this structure in deterministic environments, both categorical and continuous valued. The acquired knowledge can be exploited to augment the original data set, leading eventually to a reduction in the distributional shift between the true and the learnt model. In this work we extend the paradigm to also tackle non deterministic MDPs, in particular 1) we propose a detection threshold in categorical environments based on statistical distances, 2) we introduce a benchmark of the distributional shift in continuous environments based on the Wilcoxon signed rank statistical test and 3) we show that the former results lead to a performance improvement when solving the learnt MDP and then applying the optimal policy in the real environment.


Spatio Temporal Bayesian Learning for Mobile Edge Computing Resource Planning in Smart Cities. A smart city improves operational efficiency and comfort of living by harnessing techniques such as the Internet of Things (IoT) to collect and process data for decision making. To better support smart cities, data collected by IoT should be stored and processed appropriately. However, IoT devices are often task specialized and resource constrained, and thus, they heavily rely on online resources in terms of computing and storage to accomplish various tasks. Moreover, these cloud based solutions often centralize the resources and are far away from the end IoTs and cannot respond to users in time due to network congestion when massive numbers of tasks offload through the core network. Therefore, by decentralizing resources spatially close to IoT devices, mobile edge computing (MEC) can reduce latency and improve service quality for a smart city, where service requests can be fulfilled in proximity. As the service demands exhibit spatial temporal features, deploying MEC servers at optimal locations and allocating MEC resources play an essential role in efficiently meeting service requirements in a smart city. In this regard, it is essential to learn the distribution of resource demands in time and space. In this work, we first propose a spatio temporal Bayesian hierarchical learning approach to learn and predict the distribution of MEC resource demand over space and time to facilitate MEC deployment and resource management. Second, the proposed model is trained and tested on real world data, and the results demonstrate that the proposed method can achieve very high accuracy. Third, we demonstrate an application of the proposed method by simulating task offloading. Finally, the simulated results show that resources allocated based upon our models' predictions are exploited more efficiently than the resources are equally divided into all servers in unobserved areas.


Drazin invertibility of linear operators on quaternionic Banach spaces. Let A be a right linear operator on a two sided quaternionic Banach space X . The paper studies the Drazin inverse for right linear operators on a quaternionic Banach space. It is shown that if A is Drazin invertible then the Drazin inverse of A is given by f(A) where f is 0 in an axially symmetric neighborhood of 0 and f(q) q 1 in an axially symmetric neighborhood of the nonzero spherical spectrum of A . Some results analogous to the ones concerning the Drazin inverse of operators on complex Banach spaces are proved in the quaternionic context.


Programming Knowledge Tracing A Comprehensive Dataset and A New Model. In this paper, we study knowledge tracing in the domain of programming education and make two important contributions. First, we harvest and publish so far the most comprehensive dataset, namely BePKT, which covers various online behaviors in an OJ system, including programming text problems, knowledge annotations, user submitted code and system logged events. Second, we propose a new model PDKT to exploit the enriched context for accurate student behavior prediction. More specifically, we construct a bipartite graph for programming problem embedding, and design an improved pre training model PLCodeBERT for code embedding, as well as a double sequence RNN model with exponential decay attention for effective feature fusion. Experimental results on the new dataset BePKT show that our proposed model establishes state of the art performance in programming knowledge tracing. In addition, we verify that our code embedding strategy based on PLCodeBERT is complementary to existing knowledge tracing models to further enhance their accuracy. As a side product, PLCodeBERT also results in better performance in other programming related tasks such as code clone detection.


One dimensional game theoretic differential equations. We provide a very brief introduction to typical paths and the corresponding It o type integration. Relying on this robust It o integration, we prove an existence and uniqueness result for one dimensional differential equations driven by typical paths with non Lipschitz continuous coefficients in the spirit of Yamada Watanabe as well as an approximation result in the spirit of Doss Sussmann.


Functional and spatial rewiring jointly generate convergent divergent units in self organizing networks. Self organization through adaptive rewiring of random neural networks generates brain like topologies comprising modular small world structures with rich club effects, merely as the product of optimizing the network topology. In the nervous system, spatial organization is optimized no less by rewiring, through minimizing wiring distance and maximizing spatially aligned wiring layouts. We show that such spatial organization principles interact constructively with adaptive rewiring, contributing to establish the networks' connectedness and modular structures. We use an evolving neural network model with weighted and directed connections, in which neural traffic flow is based on consensus and advection dynamics, to show that wiring cost minimization supports adaptive rewiring in creating convergent divergent unit structures. Convergent divergent units consist of a convergent input hub, connected to a divergent output hub via subnetworks of intermediate nodes, which may function as the computational core of the unit. The prominence of minimizing wiring distance in the dynamic evolution of the network determines the extent to which the core is encapsulated from the rest of the network, i.e., the context sensitivity of its computations. This corresponds to the central role convergent divergent units play in establishing context sensitivity in neuronal information processing.


WiFiNet WiFi based indoor localisation using CNNs. Different technologies have been proposed to provide indoor localisation magnetic field, bluetooth , WiFi, etc. Among them, WiFi is the one with the highest availability and highest accuracy. This fact allows for an ubiquitous accurate localisation available for almost any environment and any device. However, WiFi based localisation is still an open problem. In this article, we propose a new WiFi based indoor localisation system that takes advantage of the great ability of Convolutional Neural Networks in classification problems. Three different approaches were used to achieve this goal a custom architecture called WiFiNet designed and trained specifically to solve this problem and the most popular pre trained networks using both transfer learning and feature extraction. Results indicate that WiFiNet is as a great approach for indoor localisation in a medium sized environment (30 positions and 113 access points) as it reduces the mean localisation error (33 ) and the processing time when compared with state of the art WiFi indoor localisation algorithms such as SVM.


Invertibility of the 3 core of Erdos Renyi Graphs with Growing Degree. Let A in mathbb R n times n be the adjacency matrix of an Erd H o s R 'enyi graph G(n, d n) for d omega(1) and d leq 3 log(n) . We show that as n goes to infinity, with probability that goes to 1 , the adjacency matrix of the 3 core of G(n, d n) is invertible.


Erdos Hajnal for graphs with no 5 hole. The Erdos Hajnal conjecture says that for every graph H there exists c 0 such that every graph G not containing H as an induced subgraph has a clique or stable set of cardinality at least G c. We prove that this is true when H is a cycle of length five. We also prove several further results for instance, that if C is a cycle and H is the complement of a forest, there exists c 0 such that every graph G containing neither of C,H as an induced subgraph has a clique or stable set of cardinality at least G c.


Rethink, Revisit, Revise A Spiral Reinforced Self Revised Network for Zero Shot Learning. Current approaches to Zero Shot Learning (ZSL) struggle to learn generalizable semantic knowledge capable of capturing complex correlations. Inspired by emph Spiral Curriculum , which enhances learning processes by revisiting knowledge, we propose a form of spiral learning which revisits visual representations based on a sequence of attribute groups (e.g., a combined group of emph color and emph shape ). Spiral learning aims to learn generalized local correlations, enabling models to gradually enhance global learning and thus understand complex correlations. Our implementation is based on a 2 stage emph Reinforced Self Revised (RSR) framework emph preview and emph review . RSR first previews visual information to construct diverse attribute groups in a weakly supervised manner. Then, it spirally learns refined localities based on attribute groups and uses localities to revise global semantic correlations. Our framework outperforms state of the art algorithms on four benchmark datasets in both zero shot and generalized zero shot settings, which demonstrates the effectiveness of spiral learning in learning generalizable and complex correlations. We also conduct extensive analysis to show that attribute groups and reinforced decision processes can capture complementary semantic information to improve predictions and aid explainability.


ARShoe Real Time Augmented Reality Shoe Try on System on Smartphones. Virtual try on technology enables users to try various fashion items using augmented reality and provides a convenient online shopping experience. However, most previous works focus on the virtual try on for clothes while neglecting that for shoes, which is also a promising task. To this concern, this work proposes a real time augmented reality virtual shoe try on system for smartphones, namely ARShoe. Specifically, ARShoe adopts a novel multi branch network to realize pose estimation and segmentation simultaneously. A solution to generate realistic 3D shoe model occlusion during the try on process is presented. To achieve a smooth and stable try on effect, this work further develop a novel stabilization method. Moreover, for training and evaluation, we construct the very first large scale foot benchmark with multiple virtual shoe try on task related labels annotated. Exhaustive experiments on our newly constructed benchmark demonstrate the satisfying performance of ARShoe. Practical tests on common smartphones validate the real time performance and stabilization of the proposed approach.


Unsupervised Domain Adaptation in LiDAR Semantic Segmentation with Self Supervision and Gated Adapters. In this paper, we focus on a less explored, but more realistic and complex problem of domain adaptation in LiDAR semantic segmentation. There is a significant drop in performance of an existing segmentation model when training (source domain) and testing (target domain) data originate from different LiDAR sensors. To overcome this shortcoming, we propose an unsupervised domain adaptation framework that leverages unlabeled target domain data for self supervision, coupled with an unpaired mask transfer strategy to mitigate the impact of domain shifts. Furthermore, we introduce the gated adapter module with a small number of parameters into the network to account for target domain specific information. Experiments adapting from both real to real and synthetic to real LiDAR semantic segmentation benchmarks demonstrate the significant improvement over prior arts.


Hosting exceptional points in 1D photonic bandgap waveguide for mode engineering. We report a planar 1D photonic bandgap waveguide exhibiting four second order exceptional points (EP2s). The interactions between the selective pairs of supported quasi guided TE modes are modulated by spatial distribution of transverse in homogeneous gain loss profile.


Learning from Temporal Gradient for Semi supervised Action Recognition. Semi supervised video action recognition tends to enable deep neural networks to achieve remarkable performance even with very limited labeled data. However, existing methods are mainly transferred from current image based methods (e.g., FixMatch). Without specifically utilizing the temporal dynamics and inherent multimodal attributes, their results could be suboptimal. To better leverage the encoded temporal information in videos, we introduce temporal gradient as an additional modality for more attentive feature extraction in this paper. To be specific, our method explicitly distills the fine grained motion representations from temporal gradient (TG) and imposes consistency across different modalities (i.e., RGB and TG). The performance of semi supervised action recognition is significantly improved without additional computation or parameters during inference. Our method achieves the state of the art performance on three video action recognition benchmarks (i.e., Kinetics 400, UCF 101, and HMDB 51) under several typical semi supervised settings (i.e., different ratios of labeled data).


HTN Planning Domain for Deployment of Cloud Applications. Cloud providers are facing a complex problem in configuring software applications ready for deployment on their infrastructures. Hierarchical Task Network (HTN) planning can provide effective means to solve such deployment problems. We present an HTN planning domain that models deployment problems as found in realistic Cloud environments.


Feudal Reinforcement Learning by Reading Manuals. Reading to act is a prevalent but challenging task which requires the ability to reason from a concise instruction. However, previous works face the semantic mismatch between the low level actions and the high level language descriptions and require the human designed curriculum to work properly. In this paper, we present a Feudal Reinforcement Learning (FRL) model consisting of a manager agent and a worker agent. The manager agent is a multi hop plan generator dealing with high level abstract information and generating a series of sub goals in a backward manner. The worker agent deals with the low level perceptions and actions to achieve the sub goals one by one. In comparison, our FRL model effectively alleviate the mismatching between text level inference and low level perceptions and actions and is general to various forms of environments, instructions and manuals and our multi hop plan generator can significantly boost for challenging tasks where multi step reasoning form the texts is critical to resolve the instructed goals. We showcase our approach achieves competitive performance on two challenging tasks, Read to Fight Monsters (RTFM) and Messenger, without human designed curriculum learning.


Variance bounds for disc polygons. We prove asymptotic lower bounds on the variance of the number of vertices and missed area of random disc polygons in convex discs whose boundary is C 2 smooth. The established lower bounds are of the same order as the upper bounds proved previously by Fodor and V ' i gh (2018).


Forecast Analysis of the COVID 19 Incidence in Lebanon Prediction of Future Epidemiological Trends to Plan More Effective Control Programs. Ever since the COVID 19 pandemic started, all the governments have been trying to limit its effects on their citizens and countries. This pandemic was harsh on different levels for almost all populations worldwide and this is what drove researchers and scientists to get involved and work on several kinds of simulations to get a better insight into this virus and be able to stop it the earliest possible. In this study, we simulate the spread of COVID 19 in Lebanon using an Agent Based Model where people are modeled as agents that have specific characteristics and behaviors determined from statistical distributions using Monte Carlo Algorithm. These agents can go into the world, interact with each other, and thus, infect each other. This is how the virus spreads. During the simulation, we can introduce different Non Pharmaceutical Interventions or more commonly NPIs that aim to limit the spread of the virus (wearing a mask, closing locations, etc). Our Simulator was first validated on concepts (e.g. Flattening the Curve and Second Wave scenario), and then it was applied on the case of Lebanon. We studied the effect of opening schools and universities on the pandemic situation in the country since the Lebanese Ministry of Education is planning to do so progressively, starting from 21 April 2021. Based on the results we obtained, we conclude that it would be better to delay the school openings while the vaccination campaign is still slow in the country.


Arbitrarily high order energy conserving methods for Poisson problems. In this paper we are concerned with energy conserving methods for Poisson problems, which are effectively solved by defining a suitable generalization of HBVMs, a class of energy conserving methods for Hamiltonian problems. The actual implementation of the methods is fully discussed, with a particular emphasis on the conservation of Casimirs. Some numerical tests are reported, in order to assess the theoretical findings.


Sharp pointwise in time error estimate of L1 scheme for nonlinear subdiffusion equations. An essential feature of the subdiffusion equations with the alpha order time fractional derivative is the weak singularity at the initial time. The weak regularity of the solution is usually characterized by a regularity parameter sigma in (0,1) cup(1,2) . Under this general regularity assumption, we here obtain the pointwise in time error estimate of the widely used L1 scheme for nonlinear subdiffusion equations. To the end, we present a refined discrete fractional type Gr "onwall inequality and a rigorous analysis for the truncation errors. Numerical experiments are provided to demonstrate the effectiveness of our theoretical analysis.


Fortifying Vehicular Security Through Low Overhead Physically Unclonable Functions. Within vehicles, the Controller Area Network (CAN) allows efficient communication between the electronic control units (ECUs) responsible for controlling the various subsystems. The CAN protocol was not designed to include much support for secure communication. The fact that so many critical systems can be accessed through an insecure communication network presents a major security concern. Adding security features to CAN is difficult due to the limited resources available to the individual ECUs and the costs that would be associated with adding the necessary hardware to support any additional security operations without overly degrading the performance of standard communication. Replacing the protocol is another option, but it is subject to many of the same problems. The lack of security becomes even more concerning as vehicles continue to adopt smart features. Smart vehicles have a multitude of communication interfaces would an attacker could exploit to gain access to the networks. In this work we propose a security framework that is based on physically unclonable functions (PUFs) and lightweight cryptography (LWC). The framework does not require any modification to the standard CAN protocol while also minimizing the amount of additional message overhead required for its operation. The improvements in our proposed framework results in major reduction in the number of CAN frames that must be sent during operation. For a system with 20 ECUs for example, our proposed framework only requires 6.5 of the number of CAN frames that is required by the existing approach to successfully authenticate every ECU.


Learning adaptive coarse spaces of BDDC algorithms for stochastic elliptic problems with oscillatory and high contrast coefficients. In this paper, we consider the balancing domain decomposition by constraints (BDDC) algorithm with adaptive coarse spaces for a class of stochastic elliptic problems. The key ingredient in the construction of the coarse space is the solutions of local spectral problems, which depend on the coefficient of the PDE. This poses a significant challenge for stochastic coefficients as it is computationally expensive to solve the local spectral problems for every realisation of the coefficient. To tackle this computational burden, we propose a machine learning approach. Our method is based on the use of a deep neural network (DNN) to approximate the relation between the stochastic coefficients and the coarse spaces. For the input of the DNN, we apply the Karhunen Lo eve expansion and use the first few dominant terms in the expansion. The output of the DNN is the resulting coarse space, which is then applied with the standard adaptive BDDC algorithm. We will present some numerical results with oscillatory and high contrast coefficients to show the efficiency and robustness of the proposed scheme.


The Biot Stokes coupling using total pressure formulation, analysis and application to interfacial flow in the eye. We consider a multiphysics model for the flow of Newtonian fluid coupled with Biot consolidation equations through an interface, and incorporating total pressure as an unknown in the poroelastic region. A new mixed primal finite element scheme is proposed solving for the pairs fluid velocity pressure and displacement total poroelastic pressure using Stokes stable elements, and where the formulation does not require Lagrange multipliers to set up the usual transmission conditions on the interface. The stability and well posedness of the continuous and semi discrete problems are analysed in detail. Our numerical study is framed in the context of applicative problems pertaining to heterogeneous geophysical flows and to eye poromechanics. For the latter, we investigate different interfacial flow regimes in Cartesian and axisymmetric coordinates that could eventually help describe early morphologic changes associated with glaucoma development in canine species.


Advances in adversarial attacks and defenses in computer vision A survey. Deep Learning (DL) is the most widely used tool in the contemporary field of computer vision. Its ability to accurately solve complex problems is employed in vision research to learn deep neural models for a variety of tasks, including security critical applications. However, it is now known that DL is vulnerable to adversarial attacks that can manipulate its predictions by introducing visually imperceptible perturbations in images and videos. Since the discovery of this phenomenon in 2013 1 , it has attracted significant attention of researchers from multiple sub fields of machine intelligence. In 2 , we reviewed the contributions made by the computer vision community in adversarial attacks on deep learning (and their defenses) until the advent of year 2018. Many of those contributions have inspired new directions in this area, which has matured significantly since witnessing the first generation methods. Hence, as a legacy sequel of 2 , this literature review focuses on the advances in this area since 2018. To ensure authenticity, we mainly consider peer reviewed contributions published in the prestigious sources of computer vision and machine learning research. Besides a comprehensive literature review, the article also provides concise definitions of technical terminologies for non experts in this domain. Finally, this article discusses challenges and future outlook of this direction based on the literature reviewed herein and 2 .


Moment functions and exponential monomials on commutative hypergroups. The purpose of this paper is to prove that if on a commutative hypergroup an exponential monomial has the property that the linear subspace of all sine functions in its variety is one dimensional, then this exponential monomial is a linear combination of generalized moment functions.


Dual Graph Convolutional Networks with Transformer and Curriculum Learning for Image Captioning. Existing image captioning methods just focus on understanding the relationship between objects or instances in a single image, without exploring the contextual correlation existed among contextual image. In this paper, we propose Dual Graph Convolutional Networks (Dual GCN) with transformer and curriculum learning for image captioning. In particular, we not only use an object level GCN to capture the object to object spatial relation within a single image, but also adopt an image level GCN to capture the feature information provided by similar images. With the well designed Dual GCN, we can make the linguistic transformer better understand the relationship between different objects in a single image and make full use of similar images as auxiliary information to generate a reasonable caption description for a single image. Meanwhile, with a cross review strategy introduced to determine difficulty levels, we adopt curriculum learning as the training strategy to increase the robustness and generalization of our proposed model. We conduct extensive experiments on the large scale MS COCO dataset, and the experimental results powerfully demonstrate that our proposed method outperforms recent state of the art approaches. It achieves a BLEU 1 score of 82.2 and a BLEU 2 score of 67.6. Our source code is available at em color magenta url https github.com Unbear430 DGCN for image captioning .


Prediction of the facial growth direction with Machine Learning methods. First attempts of prediction of the facial growth (FG) direction were made over half of a century ago. Despite numerous attempts and elapsed time, a satisfactory method has not been established yet and the problem still poses a challenge for medical experts. To our knowledge, this paper is the first Machine Learning approach to the prediction of FG direction. Conducted data analysis reveals the inherent complexity of the problem and explains the reasons of difficulty in FG direction prediction based on 2D X ray images. To perform growth forecasting, we employ a wide range of algorithms, from logistic regression, through tree ensembles to neural networks and consider three, slightly different, problem formulations. The resulting classification accuracy varies between 71 and 75 .


Expectation based Minimalist Grammars. Expectation based Minimalist Grammars (e MGs) are simplified versions of the (Conflated) Minimalist Grammars, (C)MGs, formalized by Stabler (Stabler, 2011, 2013, 1997) and Phase based Minimalist Grammars, PMGs (Chesi, 2005, 2007 Stabler, 2011). The crucial simplification consists of driving structure building only by relying on lexically encoded categorial top down expectations. The commitment on a top down derivation (as in e MGs and PMGs, as opposed to (C)MGs, Chomsky, 1995 Stabler, 2011) allows us to define a core derivation that should be the same in both parsing and generation (Momma Phillips, 2018).


M3DeTR Multi representation, Multi scale, Mutual relation 3D Object Detection with Transformers. We present a novel architecture for 3D object detection, M3DeTR, which combines different point cloud representations (raw, voxels, bird eye view) with different feature scales based on multi scale feature pyramids. M3DeTR is the first approach that unifies multiple point cloud representations, feature scales, as well as models mutual relationships between point clouds simultaneously using transformers. We perform extensive ablation experiments that highlight the benefits of fusing representation and scale, and modeling the relationships. Our method achieves state of the art performance on the KITTI 3D object detection dataset and Waymo Open Dataset. Results show that M3DeTR improves the baseline significantly by 1.48 mAP for all classes on Waymo Open Dataset. In particular, our approach ranks 1st on the well known KITTI 3D Detection Benchmark for both car and cyclist classes, and ranks 1st on Waymo Open Dataset with single frame point cloud input. Our code is available at https github.com rayguan97 M3DETR.


Off policy Reinforcement Learning with Optimistic Exploration and Distribution Correction. Improving sample efficiency of reinforcement learning algorithms requires effective exploration. Following the principle of textit optimism in the face of uncertainty , we train a separate exploration policy to maximize an approximate upper confidence bound of the critics in an off policy actor critic framework. However, this introduces extra differences between the replay buffer and the target policy in terms of their stationary state action distributions. To mitigate the off policy ness, we adapt the recently introduced DICE framework to learn a distribution correction ratio for off policy actor critic training. In particular, we correct the training distribution for both policies and critics. Empirically, we evaluate our proposed method in several challenging continuous control tasks and show superior performance compared to state of the art methods. We also conduct extensive ablation studies to demonstrate the effectiveness and the rationality of the proposed method.


Stable rank 3 vector bundles on mathbb P 3 with c 1 0 , c 2 3 . We clarify the undecided case c 2 3 of a theorem of Ein, Hartshorne and Vogelaar Math. Ann. 259 (1982), 541 569 about the restriction of a stable rank 3 vector bundle with c 1 0 on the projective 3 space to a general plane. It turns out that there are more exceptions to the stable restriction property than those conjectured by the three authors. One of them is a Schwarzenberger bundle (twisted by 1 ) it has c 3 6 . There are also some exceptions with c 3 2 (plus, of course, their duals). We also prove, for completeness, the basic properties of the corresponding moduli spaces they are all nonsingular and connected, of dimension 28.


A Semantic Indexing Structure for Image Retrieval. In large scale image retrieval, many indexing methods have been proposed to narrow down the searching scope of retrieval. The features extracted from images usually are of high dimensions or unfixed sizes due to the existence of key points. Most of existing index structures suffer from the dimension curse, the unfixed feature size and or the loss of semantic similarity. In this paper a new classification based indexing structure, called Semantic Indexing Structure (SIS), is proposed, in which we utilize the semantic categories rather than clustering centers to create database partitions, such that the proposed index SIS can be combined with feature extractors without the restriction of dimensions. Besides, it is observed that the size of each semantic partition is positively correlated with the semantic distribution of database. Along this way, we found that when the partition number is normalized to five, the proposed algorithm performed very well in all the tests. Compared with state of the art models, SIS achieves outstanding performance.


Optimal convergence rate for periodic homogenization of convex Hamilton Jacobi equations. We study the rate of convergence in periodic homogenization of convex Hamilton Jacobi equations. We show that the convergence rate is always O( varepsilon) , which is optimal. Our proof uses the optimal control formula, which is similar to that in 7 , and a beautiful topological lemma on decompositions of continuous curves, Lemma 2.1, by Burago 1 .


uTHCD A New Benchmarking for Tamil Handwritten OCR. Handwritten character recognition is a challenging research in the field of document image analysis over many decades due to numerous reasons such as large writing styles variation, inherent noise in data, expansive applications it offers, non availability of benchmark databases etc. There has been considerable work reported in literature about creation of the database for several Indic scripts but the Tamil script is still in its infancy as it has been reported only in one database 5 . In this paper, we present the work done in the creation of an exhaustive and large unconstrained Tamil Handwritten Character Database (uTHCD). Database consists of around 91000 samples with nearly 600 samples in each of 156 classes. The database is a unified collection of both online and offline samples. Offline samples were collected by asking volunteers to write samples on a form inside a specified grid. For online samples, we made the volunteers write in a similar grid using a digital writing pad. The samples collected encompass a vast variety of writing styles, inherent distortions arising from offline scanning process viz stroke discontinuity, variable thickness of stroke, distortion etc. Algorithms which are resilient to such data can be practically deployed for real time applications. The samples were generated from around 650 native Tamil volunteers including school going kids, homemakers, university students and faculty. The isolated character database will be made publicly available as raw images and Hierarchical Data File (HDF) compressed file. With this database, we expect to set a new benchmark in Tamil handwritten character recognition and serve as a launchpad for many avenues in document image analysis domain. Paper also presents an ideal experimental set up using the database on convolutional neural networks (CNN) with a baseline accuracy of 88 on test data.


Optimal Transmit Strategy for Multi user MIMO WPT Systems With Non linear Energy Harvesters. In this paper, we study multi user multi antenna wireless power transfer (WPT) systems, where each antenna at the energy harvesting (EH) nodes is connected to a dedicated non linear rectifier. We propose an optimal transmit strategy which maximizes a weighted sum of the average harvested powers at the EH nodes under a constraint on the power budget of the transmitter. First, for multiple input single output (MISO) WPT systems, we show that it is optimal to transmit scalar symbols with an arbitrary phase and an amplitude, whose probability density function (pdf) has at most two mass points, using maximum ratio transmission (MRT) beamforming. Then, we prove that for single input multiple output (SIMO) WPT systems, the optimal transmit symbol amplitudes are discrete random variables, whose pdf also has no more than two mass points. For general multi user MIMO WPT systems, we show that the optimal transmit strategy involves scalar unit norm symbols with arbitrary phase and at most two beamforming vectors. In order to determine these vectors, we formulate a non convex optimization problem and obtain an optimal solution based on monotonic optimization. Since the computational complexity of the optimal solution is high, we also propose a low complexity iterative algorithm to obtain a suboptimal solution, which achieves nearly optimal performance. Our simulation results reveal that the proposed transmit strategy for multi user MIMO WPT systems outperforms two baseline schemes, which are based on a linear EH model and a single beamforming vector, respectively. For a given transmit power budget, we show that the harvested power saturates when increasing the number of transmit antennas. Finally, we observe that the harvested power region spanned by multiple EH nodes is convex and the power harvested at one EH node can be traded for a higher harvested power at the other nodes.


Conjectural invariance with respect to the fusion system of an almost source algebra. We show that, given an almost source algebra A of a p block of a finite group G , then the unit group of A contains a basis stabilized by the left and right multiplicative action of the defect group if and only if, in a sense to be made precise, certain relative multiplicities of local pointed groups are invariant with respect to the fusion system. We also show that, when G is p solvable, those two equivalent conditions hold for some almost source algebra of the given p block. One motive lies in the fact that, by a theorem of Linckelmann, if the two equivalent conditions hold for A , then any stable basis for A is semicharacteristic for the fusion system.


The S weak global dimension of commutative rings. In this paper, we introduce and study the S weak global dimension S w.gl.dim (R) of a commutative ring R for some multiplicative subset S of R . Moreover, commutative rings with S weak global dimension at most 1 are studied. Finally, we investigated the S weak global dimension of factor rings and polynomial rings.


Bayesian Poroelastic Aquifer Characterization from InSAR Surface Deformation Data Part II Quantifying the Uncertainty. Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. Here we develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE) based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one dimensional problem and demonstrate it numerically for a three dimensional aquifer model. We design a generalized preconditioned Crank Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low rank based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.


Generators for K theoretic Hall algebras of quivers with potential. K theoretic Hall algebras (KHAs) of quivers with potential (Q,W) are a generalization of preprojective KHAs of quivers, which are conjecturally positive parts of the Okounkov Smironov quantum affine algebras. In particular, preprojective KHAs are expected to satisfy a Poincar 'e Birkhoff Witt theorem. We construct semi orthogonal decompositions of categorical Hall algebras using techniques developed by Halpern Leistner, Ballard Favero Katzarkov, and v S penko Van den Bergh. For a quotient of text KHA (Q,W) mathbb Q , we refine these decompositions and prove a PBW type theorem for it. The spaces of generators of text KHA (Q,0) mathbb Q are given by (a version of) intersection K theory of coarse moduli spaces of representations of Q .


Erratum Leveraging Flexible Tree Matching to Repair Broken Locators in Web Automation Scripts. Web applications are constantly evolving to integrate new features and fix reported bugs. Even an imperceptible change can sometimes entail significant modifications of the Document Object Model (DOM), which is the underlying model used by browsers to render all the elements included in a web application. Scripts that interact with web applications (e.g. web test scripts, crawlers, or robotic process automation) rely on this continuously evolving DOM which means they are often particularly fragile. More precisely, the major cause of breakages observed in automation scripts are element locators, which are identifiers used by automation scripts to navigate across the DOM. When the DOM evolves, these identifiers tend to break, thus causing the related scripts to no longer locate the intended target elements. For this reason, several contributions explored the idea of automatically repairing broken locators on a page. These works attempt to repair a given broken locator by scanning all elements in the new DOM to find the most similar one. Unfortunately, this approach fails to scale when the complexity of web pages grows, leading either to long computation times or incorrect element repairs. This article, therefore, adopts a different perspective on this problem by introducing a new locator repair solution that leverages tree matching algorithms to relocate broken locators. This solution, named Erratum, implements a holistic approach to reduce the element search space, which greatly eases the locator repair task and drastically improves repair accuracy. We compare the robustness of Erratum on a large scale benchmark composed of realistic and synthetic mutations applied to popular web applications currently deployed in production. Our empirical results demonstrate that Erratum outperforms the accuracy of WATER, a state of the art solution, by 67 .


Modeling Interactions of Autonomous Vehicles and Pedestrians with Deep Multi Agent Reinforcement Learning for Collision Avoidance. Reliable pedestrian crash avoidance mitigation (PCAM) systems are crucial components of safe autonomous vehicles (AVs). The nature of the vehicle pedestrian interaction where decisions of one agent directly affect the other agent's optimal behavior, and vice versa, is a challenging yet often neglected aspect of such systems. We address this issue by defining a Markov decision process (MDP) for a simulated driving scenario in which an AV driving along an urban street faces a pedestrian trying to cross an unmarked crosswalk. The AV's PCAM decision policy is learned through deep reinforcement learning (DRL). Since modeling pedestrians realistically is challenging, we compare two levels of intelligent pedestrian behavior. While the baseline model follows a predefined strategy, our advanced pedestrian model is defined as a second DRL agent. This model captures continuous learning and the uncertainty inherent in human behavior, making the vehicle pedestrian interaction a deep multi agent reinforcement learning (DMARL) problem. We benchmark the developed PCAM systems according to the agents' collision rate and the resulting traffic flow efficiency with a focus on the influence of observation uncertainty or noise on the decision making of the agents. The results show that the AV is able to completely mitigate collisions under the majority of the investigated conditions, and that the DRL pedestrian model indeed learns a more intelligent crossing behavior.


Estimating Instance dependent Label noise Transition Matrix using DNNs. In label noise learning, estimating the transition matrix is a hot topic as the matrix plays an important role in building statistically consistent classifiers. Traditionally, the transition from clean distribution to noisy distribution (i.e., clean label transition matrix) has been widely exploited to learn a clean label classifier by employing the noisy data. Motivated by that classifiers mostly output Bayes optimal labels for prediction, in this paper, we study to directly model the transition from Bayes optimal distribution to noisy distribution (i.e., Bayes label transition matrix) and learn a Bayes optimal label classifier. Note that given only noisy data, it is ill posed to estimate either the clean label transition matrix or the Bayes label transition matrix. But favorably, Bayes optimal labels are less uncertain compared with the clean labels, i.e., the class posteriors of Bayes optimal labels are one hot vectors while those of clean labels are not. This enables two advantages to estimate the Bayes label transition matrix, i.e., (a) we could theoretically recover a set of Bayes optimal labels under mild conditions (b) the feasible solution space is much smaller. By exploiting the advantages, we estimate the Bayes label transition matrix by employing a deep neural network in a parameterized way, leading to better generalization and superior classification performance.


Named Entity Recognition with Small Strongly Labeled and Large Weakly Labeled Data. Weak supervision has shown promising results in many natural language processing tasks, such as Named Entity Recognition (NER). Existing work mainly focuses on learning deep NER models only with weak supervision, i.e., without any human annotation, and shows that by merely using weakly labeled data, one can achieve good performance, though still underperforms fully supervised NER with manually strongly labeled data. In this paper, we consider a more practical scenario, where we have both a small amount of strongly labeled data and a large amount of weakly labeled data. Unfortunately, we observe that weakly labeled data does not necessarily improve, or even deteriorate the model performance (due to the extensive noise in the weak labels) when we train deep NER models over a simple or weighted combination of the strongly labeled and weakly labeled data. To address this issue, we propose a new multi stage computational framework NEEDLE with three essential ingredients (1) weak label completion, (2) noise aware loss function, and (3) final fine tuning over the strongly labeled data. Through experiments on E commerce query NER and Biomedical NER, we demonstrate that NEEDLE can effectively suppress the noise of the weak labels and outperforms existing methods. In particular, we achieve new SOTA F1 scores on 3 Biomedical NER datasets BC5CDR chem 93.74, BC5CDR disease 90.69, NCBI disease 92.28.


BPPChecker An SMT based Model Checker on Basic Parallel Processes(Full Version). Program verification on concurrent programs is a big challenge due to general undecidable results. Petri nets and its extensions are used in most works. However, existing verifiers based on Petri nets are difficult to be complete and efficient. Basic Parallel Process (BPP), as a subclass of Petri nets, can be used as a model for describing and verifying concurrent programs with lower complexity. We propose and implement BPPChecker, the first model checker for verifying a subclass of CTL on BPP. We propose constraint based algorithms for the problem of model checking on BPPs and handle formulas by SMT solver Z3. For EF operator, we reduce the model checking of EF formulas to the satisfiability problem of existential Presburger formula. For EG operator, we provide a k step bounded semantics and reduce the model checking of EG formulas to the satisfiability problem of linear integer arithmetic. Besides, we give Actor Communicating System (ACS) the over approximation BPP based semantics and evaluate BPPChecker on ACSs generated from real Erlang programs. Experimental results show that BPPChecker performs more efficiently than the existing tools for a series of branching time property verification problems of Erlang programs.


Sentence Embeddings by Ensemble Distillation. This paper contributes a new State Of The Art (SOTA) for Semantic Textual Similarity (STS). We compare and combine a number of recently proposed sentence embedding methods for STS, and propose a novel and simple ensemble knowledge distillation scheme that improves on previous approaches. Our experiments demonstrate that a model trained to learn the average embedding space from multiple ensemble students outperforms all the other individual models with high robustness. Utilizing our distillation method in combination with previous methods, we significantly improve on the SOTA unsupervised STS, and by proper hyperparameter tuning of previous methods we improve the supervised SOTA scores.


Refinement of polygonal grids using Convolutional Neural Networks with applications to polygonal Discontinuous Galerkin and Virtual Element methods. We propose new strategies to handle polygonal grids refinement based on Convolutional Neural Networks (CNNs). We show that CNNs can be successfully employed to identify correctly the "shape" of a polygonal element so as to design suitable refinement criteria to be possibly employed within adaptive refinement strategies. We propose two refinement strategies that exploit the use of CNNs to classify elements' shape, at a low computational cost. We test the proposed idea considering two families of finite element methods that support arbitrarily shaped polygonal elements, namely Polygonal Discontinuous Galerkin (PolyDG) methods and Virtual Element Methods (VEMs). We demonstrate that the proposed algorithms can greatly improve the performance of the discretization schemes both in terms of accuracy and quality of the underlying grids. Moreover, since the training phase is performed off line and is independent of the differential model the overall computational costs are kept low.


Ensemble inversion for brain tumor growth models with mass effect. We propose a method for extracting physics based biomarkers from a single multiparametric Magnetic Resonance Imaging (mpMRI) scan bearing a glioma tumor. We account for mass effect, the deformation of brain parenchyma due to the growing tumor, which on its own is an important radiographic feature but its automatic quantification remains an open problem. In particular, we calibrate a partial differential equation (PDE) tumor growth model that captures mass effect, parameterized by a single scalar parameter, tumor proliferation, migration, while localizing the tumor initiation site. The single scan calibration problem is severely ill posed because the precancerous, healthy, brain anatomy is unknown. To address the ill posedness, we introduce an ensemble inversion scheme that uses a number of normal subject brain templates as proxies for the healthy precancer subject anatomy. We verify our solver on a synthetic dataset and perform a retrospective analysis on a clinical dataset of 216 glioblastoma (GBM) patients. We analyze the reconstructions using our calibrated biophysical model and demonstrate that our solver provides both global and local quantitative measures of tumor biophysics and mass effect. We further highlight the improved performance in model calibration through the inclusion of mass effect in tumor growth models including mass effect in the model leads to 10 increase in average dice coefficients for patients with significant mass effect. We further evaluate our model by introducing novel biophysics based features and using them for survival analysis. Our preliminary analysis suggests that including such features can improve patient stratification and survival prediction.


Condition Sensing for Electricity Infrastructure in Disasters by Mining Public Topics from Social Media. Timely and reliable sensing of infrastructure conditions is critical in disaster management for planning effective infrastructure restorations. Social media, a near real time information source, has been widely used in disasters for forming timely situational awareness. Yet, using social media to sense electricity infrastructure conditions has not been explored. This study aims to address the research gap through mining public topics from social media. To achieve this purpose, we proposed a systematic and customized approach wherein (1) electricity related social media data is extracted by the classifier developed based on Bidirectional Encoder Representations from Transformers (BERT) and (2) public topics are modeled with unigrams, bigrams, and trigrams to incorporate the formulaic expressions of infrastructure conditions in social media. Electricity infrastructures in Florida impacted by Hurricane Irma are studied for illustration and demonstration. Results show that the proposed approach is capable of sensing the temporal evolutions and geographic differences of electricity infrastructure conditions.


Biomedical Interpretable Entity Representations. Pre trained language models induce dense entity representations that offer strong performance on entity centric NLP tasks, but such representations are not immediately interpretable. This can be a barrier to model uptake in important domains such as biomedicine. There has been recent work on general interpretable representation learning (Onoe and Durrett, 2020), but these domain agnostic representations do not readily transfer to the important domain of biomedicine. In this paper, we create a new entity type system and training set from a large corpus of biomedical texts by mapping entities to concepts in a medical ontology, and from these to Wikipedia pages whose categories are our types. From this mapping we derive Biomedical Interpretable Entity Representations(BIERs), in which dimensions correspond to fine grained entity types, and values are predicted probabilities that a given entity is of the corresponding type. We propose a novel method that exploits BIER's final sparse and intermediate dense representations to facilitate model and entity type debugging. We show that BIERs achieve strong performance in biomedical tasks including named entity disambiguation and entity label classification, and we provide error analysis to highlight the utility of their interpretability, particularly in low supervision settings. Finally, we provide our induced 68K biomedical type system, the corresponding 37 million triples of derived data used to train BIER models and our best performing model.


FauxWard A Graph Neural Network Approach to Fauxtography Detection Using Social Media Comments. Online social media has been a popular source for people to consume and share news content. More recently, the spread of misinformation online has caused widespread concerns. In this work, we focus on a critical task of detecting fauxtography on social media where the image and associated text together convey misleading information. Many efforts have been made to mitigate misinformation online, but we found that the fauxtography problem has not been fully addressed by existing work. Solutions focusing on detecting fake images or misinformed texts alone on social media often fail to identify the misinformation delivered together by the image and the associated text of a fauxtography post. In this paper, we develop FauxWard, a novel graph convolutional neural network framework that explicitly explores the complex information extracted from a user comment network of a social media post to effectively identify fauxtography. FauxWard is content free in the sense that it does not analyze the visual or textual contents of the post itself, which makes it robust against sophisticated fauxtography uploaders who intentionally craft image centric posts by editing either the text or image content. We evaluate FauxWard on two real world datasets collected from mainstream social media platforms (i.e., Reddit and Twitter). The results show that FauxWard is both effective and efficient in identifying fauxtography posts on social media.


Homo Entropicus, the emotional agent and societies of Neural Networks. A neural network with a learning algorithm optimized by information theory entropic dynamics is used to build an agent dubbed Homo Entropicus. The algorithm can be described at a macroscopic level in terms of aggregate variables interpretable as quantitative markers of proto emotions. We use systems of such interacting neural networks to construct a framework for modeling societies that show complex emergent behavior. A few applications are presented to investigate the role the interactions of opinions about multidimensional issues and trust on the information source play on the state of the agent society. These include the case of a class of N agents learning from a fixed teacher two dynamical agents panels of three agents modeling the interactions that occur in decisions of the US Court of Appeals, where we quantify how politically biased are the agents, how trustful of other agents judges of other parties, how much the agents follow a common understanding of the law. Finally we address under which conditions ideological polarization follows or precedes affective polarization in large societies and how simpler versions of the learning algorithm may change these relations.


Efficient State space Exploration in Massively Parallel Simulation Based Inference. Simulation based Inference (SBI) is a widely used set of algorithms to learn the parameters of complex scientific simulation models. While primarily run on CPUs in HPC clusters, these algorithms have been shown to scale in performance when developed to be run on massively parallel architectures such as GPUs. While parallelizing existing SBI algorithms provides us with performance gains, this might not be the most efficient way to utilize the achieved parallelism. This work proposes a new algorithm, that builds on an existing SBI method Approximate Bayesian Computation with Sequential Monte Carlo(ABC SMC). This new algorithm is designed to utilize the parallelism not only for performance gain, but also toward qualitative benefits in the learnt parameters. The key idea is to replace the notion of a single 'step size' hyperparameter, which governs how the state space of parameters is explored during learning, with step sizes sampled from a tuned Beta distribution. This allows this new ABC SMC algorithm to more efficiently explore the state space of the parameters being learnt. We test the effectiveness of the proposed algorithm to learn parameters for an epidemiology model running on a Tesla T4 GPU. Compared to the parallelized state of the art SBI algorithm, we get similar quality results in sim 100 x fewer simulations and observe 80x lower run to run variance across 10 independent trials.


HyperTeNet Hypergraph and Transformer based Neural Network for Personalized List Continuation. The personalized list continuation (PLC) task is to curate the next items to user generated lists (ordered sequence of items) in a personalized way. The main challenge in this task is understanding the ternary relationships among the interacting entities (users, items, and lists) that the existing works do not consider. Further, they do not take into account the multi hop relationships among entities of the same type. In addition, capturing the sequential information amongst the items already present in the list also plays a vital role in determining the next relevant items that get curated. In this work, we propose HyperTeNet a self attention hypergraph and Transformer based neural network architecture for the personalized list continuation task to address the challenges mentioned above. We use graph convolutions to learn the multi hop relationship among the entities of the same type and leverage a self attention based hypergraph neural network to learn the ternary relationships among the interacting entities via hyperlink prediction in a 3 uniform hypergraph. Further, the entity embeddings are shared with a Transformer based architecture and are learned through an alternating optimization procedure. As a result, this network also learns the sequential information needed to curate the next items to be added to the list. Experimental results demonstrate that HyperTeNet significantly outperforms the other state of the art models on real world datasets. Our implementation and datasets are available at https github.com mvijaikumar HyperTeNet.


Twisted Conjugacy in Big Mapping Class Groups. Let G be a group and varphi be an automorphism of G . Two elements x, y of G are said to be varphi twisted conjugate if y gx varphi(g) 1 for some g in G . A group G has the R infty property if the number of varphi twisted conjugacy classes is infinite for every automorphism varphi of G . In this paper we prove that the big mapping class group MCG(S) possesses the R infty property under some suitable conditions on the infinite type surface S . As an application we also prove that the big mapping class group possesses the R infty property if and only if it satisfies the S infty property.


Empirical Analysis on Productivity Prediction and Locality for Use Case Points Method. Use Case Points (UCP) method has been around for over two decades. Although, there was a substantial criticism concerning the algebraic construction and factors assessment of UCP, it remains an efficient early size estimation method. Predicting software effort from UCP is still an ever present challenge. The earlier version of UCP method suggested using productivity as a cost driver, where fixed or a few pre defined productivity ratios have been widely agreed. While this approach was successful when no enough historical data is available, it is no longer acceptable because software projects are different in terms of development aspects. Therefore, it is better to understand the relationship between productivity and other UCP variables. This paper examines the impact of data locality approaches on productivity and effort prediction from multiple UCP variables. The environmental factors are used as partitioning factors to produce local homogeneous data either based on their influential levels or using clustering algorithms. Different machine learning methods, including solo and ensemble methods, are used to construct productivity and effort prediction models based on the local data. The results demonstrate that the prediction models that are created based on local data surpass models that use entire data. Also, the results show that conforming the hypothetical assumption between productivity and environmental factors is not necessarily a requirement for success of locality.


Local and uniform moduli of continuity of chi square processes. Let eta eta(t) t in 0,1 be a mean zero continuous Gaussian process with covariance U U(s,t),s,t in 0,1 , with U(0,0) 0 . Let eta i i 1, ldots, k be independent copies of eta and set Y k (t) sum i 1 k eta 2 i (t), t in 0,1 . The stochastic process Y k Y k (t),t in 0,1 is referred to as a chi square process of order k with kernel U . Let phi(t) be a positive function on 0, delta for some delta 0 . If limsup t to 0 frac eta(t) eta(0) phi(t) 1 qquad a.s., then for all integers k ge 1 , limsup t to 0 frac Y k (t) Y k (0) phi (t) 2 Y 1 2 k (0) qquad a.s. Set sigma 2(u,v) E( eta(u) eta(v)) 2 quad text and quad widetilde sigma 2(x) sup u v le x sigma 2(u,v). Assume that inf t in 0,1 U(t,t) 0 and, lim x to0 widetilde sigma 2(x) log 1 x 0. Let varphi(t) be a positive function on 0,1 . Then if lim h to 0 sup stackrel u v le h u,v in Delta frac eta(u) eta(v) varphi( u v ) 1 qquad a.s. for all intervals Delta subset 0,1 , it follows that for all intervals Delta subset 0,1 and all integers k ge 1 , lim h to 0 sup stackrel u v le h u,v in Delta frac Y k (u) Y k (v) varphi ( u v ) 2 sup u in Delta Y k 1 2 (u), hspace .2 in a.s.


The Card Guessing Game A generating function approach. Consider a card guessing game with complete feedback in which a deck of n cards ordered 1, dots, n is riffle shuffled once. With the goal to maximize the number of correct guesses, a player guesses cards from the top of the deck one at a time under the optimal strategy until no cards remain. We provide an expression for the expected number of correct guesses with arbitrary number of terms, an accuracy improvement over the results of Liu (2021). In addition, using generating functions, we give a unified framework for systematically calculating higher order moments. Although the extension of the framework to k geq2 shuffles is not immediately straightforward, we are able to settle a long standing conjectured optimal strategy of Bayer and Diaconis (1992) by showing that the optimal guessing strategy for k 1 riffle shuffle does not necessarily apply to k geq2 shuffles.


Study of similarity rules for electromagnetic process in partially ionized plasmas. As proved by previous study, the similarity of electromagnetic processes in plasmas will be violated by Coulomb collisions between electron and ions. Therefore, there is no similarity in highly ionized collisional plasma. But the situation will be completely different for collisional plasmas with low ionization degree. The main collision type will change from electron ion Coulomb collision to electron molecule collision, and a new variable (the number density of neutral molecules) will be introduced into the similarity constraint, which can increase the degree of freedom.Thus, in this condition, the similarity restriction caused by the collision process does not conflict with the other restrictions. Therefore, the similarity for electromagnetic process in collisional plasmas can be valid for partially ionized plasmas. In this paper, we propose the similarity in partially ionized plasmas, and prove it via particle in cell Monte Carlo (PIC MCC) simulation. Our research has a wide range of engineering applications.


Time Masking for Temporal Language Models. Our world is constantly evolving, and so is the content on the web. Consequently, our languages, often said to mirror the world, are dynamic in nature. However, most current contextual language models are static and cannot adapt to changes over time. In this work, we propose a temporal contextual language model called TempoBERT, which uses time as an additional context of texts. Our technique is based on modifying texts with temporal information and performing time masking specific masking for the supplementary time information. We leverage our approach for the tasks of semantic change detection and sentence time prediction, experimenting on diverse datasets in terms of time, size, genre, and language. Our extensive evaluation shows that both tasks benefit from exploiting time masking.


Response of the GAGG(Ce) scintillator to charged particles compared with the CsI(Tl) scintillator. GAGG(Ce) is a novel scintillator with a fast response and high light output without a hygroscopic nature. It is expected to be a useful detector for charged particles at high counting rates. However, the response of the GAGG(Ce) scintillator to charged particles has not been fully examined. In the present work, the light output and energy resolution of the GAGG(Ce) scintillator were measured for protons and alpha particles at E p 5 68 MeV and E alpha 8 54 MeV as well as gamma rays at E gamma 662 keV from a 137 Cs source. The results were compared with those of the CsI(Tl) scintillator. The scintillation efficiencies dL dE of the GAGG(Ce) and CsI(Tl) scintillators were obtained and parametrized as a function of linear energy transfer dE dx .


A Sequential Convex Programming Approach to Solving Quadratic Programs and Optimal Control Problems with Linear Complementarity Constraints. Mathematical programs with complementarity constraints are notoriously difficult to solve due to their nonconvexity and lack of constraint qualifications in every feasible point. This work focuses on the subclass of quadratic programs with linear complementarity constraints. A novel approach to solving a penalty reformulation using sequential convex programming and a homotopy on the penalty parameter is introduced. Linearizing the necessarily nonconvex penalty function yields convex quadratic subproblems, which have a constant Hessian matrix throughout all iterates. This allows solution computation with a single KKT matrix factorization. Furthermore, a globalization scheme is introduced in which the underlying merit function is minimized analytically, and guarantee of descent is provided at each iterate. The algorithmic features and possible computational speedups are illustrated in a numerical experiment.