Large deviations and a phase transition in the Block Spin Potts models. We introduce and analyze a generalization of the blocks spin Ising (Curie Weiss) models that were discussed in a number of recent articles. In these block spin models each spin in one of s blocks can take one of a finite number of q ge 3 values, hence the name block spin Potts model. The values a spin can take are called colors. We prove a large deviation principle for the percentage of spins of a certain color in a certain block. These values are represented in an s times q matrix. We show that for uniform block sizes and appropriately chosen interaction strength there is a phase transition. In some regime the only equilibrium is the uniform distribution of all colors in all blocks, while in other parameter regimes there is one predominant color, and this is the same color with the same frequency for all blocks. Finally, we establish log Sobolev type inequalities for the block spin Potts model.


CL IMS DIACR Ita Volente o Nolente BERT does not outperform SGNS on Semantic Change Detection. We present the results of our participation in the DIACR Ita shared task on lexical semantic change detection for Italian. We exploit Average Pairwise Distance of token based BERT embeddings between time points and rank 5 (of 8) in the official ranking with an accuracy of .72 . While we tune parameters on the English data set of SemEval 2020 Task 1 and reach high performance, this does not translate to the Italian DIACR Ita data set. Our results show that we do not manage to find robust ways to exploit BERT embeddings in lexical semantic change detection.


SCADA System Testbed for Cybersecurity Research Using Machine Learning Approach. This paper presents the development of a Supervisory Control and Data Acquisition (SCADA) system testbed used for cybersecurity research. The testbed consists of a water storage tank's control system, which is a stage in the process of water treatment and distribution. Sophisticated cyber attacks were conducted against the testbed. During the attacks, the network traffic was captured, and features were extracted from the traffic to build a dataset for training and testing different machine learning algorithms. Five traditional machine learning algorithms were trained to detect the attacks Random Forest, Decision Tree, Logistic Regression, Naive Bayes and KNN. Then, the trained machine learning models were built and deployed in the network, where new tests were made using online network traffic. The performance obtained during the training and testing of the machine learning models was compared to the performance obtained during the online deployment of these models in the network. The results show the efficiency of the machine learning models in detecting the attacks in real time. The testbed provides a good understanding of the effects and consequences of attacks on real SCADA environments


Automatic Tip Detection of Surgical Instruments in Biportal Endoscopic Spine Surgery. Some endoscopic surgeries require a surgeon to hold the endoscope with one hand and the surgical instruments with the other hand to perform the actual surgery with correct vision. Recent technical advances in deep learning as well as in robotics can introduce robotics to these endoscopic surgeries. This can have numerous advantages by freeing one hand of the surgeon, which will allow the surgeon to use both hands and to use more intricate and sophisticated techniques. Recently, deep learning with convolutional neural network achieves state of the art results in computer vision. Therefore, the aim of this study is to automatically detect the tip of the instrument, localize a point, and evaluate detection accuracy in biportal endoscopic spine surgery. The localized point could be used for the controller's inputs of robotic endoscopy in these types of endoscopic surgeries.


Biophysical characterization of membrane phase transition profiles for the discrimination of Outer Membrane Vesicles (OMVs) from Escherichia coli grown at different temperatures. Dynamic Light Scattering (DLS), Small Angle X ray Scattering (SAXS) and Transmission Electron Microscopy (TEM) are physical techniques widely employed to characterize the morphology and the structure of vesicles such as liposomes or human extracellular vesicles (exosomes). Bacterial extracellular vesicles are similar in size to human exosomes, although their function and membrane properties have not been elucidated in such detail as in the case of exosomes. Here, we applied the above cited techniques, in synergy with the thermodynamic characterization of the vesicles lipid membrane using a turbidimetric technique to the study of vesicles produced by Gram negative bacteria (Outer Membrane Vesicles, OMVs) grown at different temperatures. This study demonstrated that our combined approach is useful to discriminate vesicles of different origin or coming from bacteria cultured under different experimental conditions. We envisage that in a near future the techniques employed in our work will be further implemented to discriminate complex mixtures of bacterial vesicles, thus showing great promises for biomedical or diagnostic applications.


Disentangling Image Distortions in Deep Feature Space. Previous literature suggests that perceptual similarity is an emergent property shared across deep visual representations. Experiments conducted on a dataset of human judged image distortions have proven that deep features outperform classic perceptual metrics. In this work we take a further step in the direction of a broader understanding of such property by analyzing the capability of deep visual representations to intrinsically characterize different types of image distortions. To this end, we firstly generate a number of synthetically distorted images and then we analyze the features extracted by different layers of different Deep Neural Networks. We observe that a dimension reduced representation of the features extracted from a given layer permits to efficiently separate types of distortions in the feature space. Moreover, each network layer exhibits a different ability to separate between different types of distortions, and this ability varies according to the network architecture. Finally, we evaluate the exploitation of features taken from the layer that better separates image distortions for i) reduced reference image quality assessment, and ii) distortion types and severity levels characterization on both single and multiple distortion databases. Results achieved on both tasks suggest that deep visual representations can be unsupervisedly employed to efficiently characterize various image distortions.


acados a modular open source framework for fast embedded optimal control. This paper presents the acados software package, a collection of solvers for fast embedded optimization intended for fast embedded applications. Its interfaces to higher level languages make it useful for quickly designing an optimization based control algorithm by putting together different algorithmic components that can be readily connected and interchanged. Since the core of acados is written on top of a high performance linear algebra library, we do not sacrifice computational performance. Thus, we aim to provide both flexibility and performance through modularity, without the need to rely on automatic code generation, which facilitates maintainability and extensibility. The main features of acados are efficient optimal control algorithms targeting embedded devices implemented in C, linear algebra based on the high performance BLASFEO library, user friendly interfaces to Matlab and Python, and compatibility with the modeling language of CasADi. acados is free and open source software released under the permissive BSD 2 Clause license.


Improved Activity Forecasting for Generating Trajectories. An efficient inverse reinforcement learning for generating trajectories is proposed based of 2D and 3D activity forecasting. We modify reward function with L p norm and propose convolution into value iteration steps, which is called convolutional value iteration. Experimental results with seabird trajectories (43 for training and 10 for test), our method is best in terms of MHD error and performs fastest. Generated trajectories for interpolating missing parts of trajectories look much similar to real seabird trajectories than those by the previous works.


Self adapting confidence estimation for stereo. Estimating the confidence of disparity maps inferred by a stereo algorithm has become a very relevant task in the years, due to the increasing number of applications leveraging such cue. Although self supervised learning has recently spread across many computer vision tasks, it has been barely considered in the field of confidence estimation. In this paper, we propose a flexible and lightweight solution enabling self adapting confidence estimation agnostic to the stereo algorithm or network. Our approach relies on the minimum information available in any stereo setup (i.e., the input stereo pair and the output disparity map) to learn an effective confidence measure. This strategy allows us not only a seamless integration with any stereo system, including consumer and industrial devices equipped with undisclosed stereo perception methods, but also, due to its self adapting capability, for its out of the box deployment in the field. Exhaustive experimental results with different standard datasets support our claims, showing how our solution is the first ever enabling online learning of accurate confidence estimation for any stereo system and without any requirement for the end user.


On Counting (Quantum )Graph Homomorphisms in Finite Fields of Prime Order. We study the problem of counting the number of homomorphisms from an input graph G to a fixed (quantum) graph bar H in any finite field of prime order mathbb Z p . The subproblem with graph H was introduced by Faben and Jerrum ToC'15 and its complexity is still uncharacterised despite active research, e.g. the very recent work of Focke, Goldberg, Roth, and Zivn 'y SODA'21 . Our contribution is threefold. First, we introduce the study of quantum graphs to the study of modular counting homomorphisms. We show that the complexity for a quantum graph bar H collapses to the complexity criteria found at dimension 1 graphs. Second, in order to prove cases of intractability we establish a further reduction to the study of bipartite graphs. Lastly, we establish a dichotomy for all bipartite (K 3,3 backslash e , , domino ) free graphs by a thorough structural study incorporating both local and global arguments. This result subsumes all results on bipartite graphs known for all prime moduli and extends them significantly. Even for the subproblem with p 2 this establishes new results.


A Convolutional Neural Network into graph space. Convolutional neural networks (CNNs), in a few decades, have outperformed the existing state of the art methods in classification context. However, in the way they were formalised, CNNs are bound to operate on euclidean spaces. Indeed, convolution is a signal operation that are defined on euclidean spaces. This has restricted deep learning main use to euclidean defined data such as sound or image. And yet, numerous computer application fields (among which network analysis, computational social science, chemo informatics or computer graphics) induce non euclideanly defined data such as graphs, networks or manifolds. In this paper we propose a new convolution neural network architecture, defined directly into graph space. Convolution and pooling operators are defined in graph domain. We show its usability in a back propagation context. Experimental results show that our model performance is at state of the art level on simple tasks. It shows robustness with respect to graph domain changes and improvement with respect to other euclidean and non euclidean convolutional architectures.


Jackson theorem and modulus of continuity in Hilbert spaces and on homogeneous manifolds. We consider a Hilbert space bf H equipped with a set of strongly continuous bounded semigroups satisfying certain conditions. The conditions allow to define a family of moduli of continuity Omega r (s,f), r in mathbb N , s 0, of vectors in bf H and a family of Paley Wiener subspaces PW sigma parametrized by bandwidth sigma 0 . These subspaces are explored to introduce notion of the best approximation mathcal E ( sigma, f) of a general vector in bf H by Paley Wiener vectors of a certain bandwidth sigma 0 . The main objective of the paper is to prove the so called Jackson type estimate mathcal E ( sigma, f) leq C left( Omega r ( sigma 1 ,f) sigma r f right) for sigma 1 . It was shown in our previous publications that our assumptions are satisfied for a strongly continuous unitary representation of a Lie group G in a Hilbert space bf H . This way we obtain the Jackson type estimates on homogeneous manifolds.


Recurrent Feature Reasoning for Image Inpainting. Existing inpainting methods have achieved promising performance for recovering regular or small image defects. However, filling in large continuous holes remains difficult due to the lack of constraints for the hole center. In this paper, we devise a Recurrent Feature Reasoning (RFR) network which is mainly constructed by a plug and play Recurrent Feature Reasoning module and a Knowledge Consistent Attention (KCA) module. Analogous to how humans solve puzzles (i.e., first solve the easier parts and then use the results as additional information to solve difficult parts), the RFR module recurrently infers the hole boundaries of the convolutional feature maps and then uses them as clues for further inference. The module progressively strengthens the constraints for the hole center and the results become explicit. To capture information from distant places in the feature map for RFR, we further develop KCA and incorporate it in RFR. Empirically, we first compare the proposed RFR Net with existing backbones, demonstrating that RFR Net is more efficient (e.g., a 4 SSIM improvement for the same model size). We then place the network in the context of the current state of the art, where it exhibits improved performance. The corresponding source code is available at https github.com jingyuanli001 RFR Inpainting


Terahertz field induced near cutoff even order harmonics in femtosecond laser. High order harmonic generation by femtosecond laser pulse in the presence of a moderately strong terahertz (THz) field is studied under the strong field approximation, showing a simple proportionality of near cutoff even order harmonic (NCEH) amplitude to the THz electric field. The formation of the THz induced NCEHs is analytically shown for both continuous wave and Gaussian pulse. The perturbation analysis with regard to the frequency ratio of the THz field to the femtosecond pulse shows the THz induced NCEHs originates from its first order correction, and the available parametric conditions for the phenomenon is also clarified. As the complete characterization of the time domain waveform of broadband THz field is essential for a wide variety of applications, the work provides an alternative time resolved field detection technique, allowing for a robust broadband characterization of pulses in THz spectral range.


The Monogeneity of Kummer Extensions and Radical Extensions. We give necessary and sufficient conditions for the Kummer extension K mathbb Q left( zeta n, sqrt n alpha right) to be monogenic over mathbb Q ( zeta n) with sqrt n alpha as a generator, i.e., for mathcal O K mathbb Z left zeta n right left sqrt n alpha right . We generalize these ideas to radical extensions of an arbitrary number field L and provide necessary and sufficient conditions for sqrt n alpha to generate a power mathcal O L basis for mathcal O L left( sqrt n alpha right) . We also give sufficient conditions for K to be non monogenic over mathbb Q and establish a general criterion relating ramification and relative monogeneity. Using this criterion, we find a necessary and sufficient condition for a relative cyclotomic extension of degree phi(n) to have zeta n as a monogenic generator.


A Streaming Machine Learning Framework for Online Aggression Detection on Twitter. The rise of online aggression on social media is evolving into a major point of concern. Several machine and deep learning approaches have been proposed recently for detecting various types of aggressive behavior. However, social media are fast paced, generating an increasing amount of content, while aggressive behavior evolves over time. In this work, we introduce the first, practical, real time framework for detecting aggression on Twitter via embracing the streaming machine learning paradigm. Our method adapts its ML classifiers in an incremental fashion as it receives new annotated examples and is able to achieve the same (or even higher) performance as batch based ML models, with over 90 accuracy, precision, and recall. At the same time, our experimental analysis on real Twitter data reveals how our framework can easily scale to accommodate the entire Twitter Firehose (of 778 million tweets per day) with only 3 commodity machines. Finally, we show that our framework is general enough to detect other related behaviors such as sarcasm, racism, and sexism in real time.


On extremal problems concerning the traces of sets. Given two non negative integers n and s , define m(n,s) to be the maximal number such that in every hypergraph mathcal H on n vertices and with at most m(n,s) edges there is a vertex x such that mathcal H x geq E( mathcal H ) s , where mathcal H x H setminus x H in E( mathcal H ) . This problem has been posed by F "uredi and Pach and by Frankl and Tokushige. While the first results were only for specific small values of s , Frankl determined m(n,2 d 1 1) for all d in mathbb N with d mid n . Subsequently, the goal became to determine m(n,2 d 1 c) for larger c . Frankl and Watanabe determined m(n,2 d 1 c) for c in 0,2 . Other general results were not known so far. Our main result sheds light on what happens further away from powers of two We prove that m(n,2 d 1 c) frac n d (2 d c) for d geq 4c and d mid n and give an example showing that this equality does not hold for c d . The other line of research on this problem is to determine m(n,s) for small values of s . In this line, our second result determines m(n,2 d 1 c) for c in 3,4 . This solves more instances of the problem for small s and in particular solves a conjecture by Frankl and Watanabe.


Exciton Localization on Ru based Photosensitizers Induced by Binding to Lipid Membranes. The characterization of electronic properties of metal complexes embedded in membrane environments is of paramount importance to develop efficient photosensitizers in optogenetic applications. Molecular dynamics and QM MM simulations together with quantitative wavefunction analysis reveal an unforeseen difference between the exciton formed upon excitation of Ru(bpy)2(bpy C17) 2 embedded in a lipid bilayer and in water, despite the media does not influence the charge transfer character and the excitation energy of the absorption spectra. In water the exciton is mainly delocalized over two ligands and the presence of non polar substituents induces directionality of the electronic transition. Instead, when the photosensitizer is embedded into a lipid membrane, the exciton size diminishes and is more localized at one bypyridyl site due to electrostatic interactions with the positively charged surface of the bilayer. These differences show that the electronic structure of metal complexes can be controlled through the binding to external species, underscoring the crucial role of the environment in directing the electronic flow upon excitation and thus helping rational tuning of optogenetic agents.


Beyond I.I.D. Three Levels of Generalization for Question Answering on Knowledge Bases. Existing studies on question answering on knowledge bases (KBQA) mainly operate with the standard i.i.d assumption, i.e., training distribution over questions is the same as the test distribution. However, i.i.d may be neither reasonably achievable nor desirable on large scale KBs because 1) true user distribution is hard to capture and 2) randomly sample training examples from the enormous space would be highly data inefficient. Instead, we suggest that KBQA models should have three levels of built in generalization i.i.d, compositional, and zero shot. To facilitate the development of KBQA models with stronger generalization, we construct and release a new large scale, high quality dataset with 64,331 questions, GrailQA, and provide evaluation settings for all three levels of generalization. In addition, we propose a novel BERT based KBQA model. The combination of our dataset and model enables us to thoroughly examine and demonstrate, for the first time, the key role of pre trained contextual embeddings like BERT in the generalization of KBQA.


Model Reduction of Semistable Distributed Parameter Systems. The model reduction problem for semistable infinite dimensional control systems is studied in this paper. In relation to these systems, we study an object we call the semistability Gramian, which serves as a generalization of the ordinary controllability Gramian valid for semistable systems. This Gramian is then given geometric as well as algebraic characterization via a Lyapunov equation. We then proceed to show that under a commutativity assumption relating the original and reduced systems, and as long as the semistability is preserved, we may derive a priori error formulas in mathcal H 2 norm in terms of the trace of this Gramian.


Reliable Group Communication Protocol for Internet of Things. In this paper, we propose RECOUP, a reliable group communication routing protocol for IoT networks. RECOUP efficiently uses a low overhead cluster based multicast routing technique on top of the RPL protocol. RECOUP increases the probability of message delivery to the intended destination(s), irrespective of the network size and faults (such as broken links or non responsive nodes), and in the presence of misbehaving nodes.We show that the cluster based routing mechanism of RECOUP remains robust in presence of various topology (i.e., rank and sybil) and data communication (i.e., blackhole, wormhole, and jamming) attacks targeting the IoT networking infrastructure.An implementation of RECOUP is realized in Contiki. Our results show the effectiveness of RECOUP over state of art protocols concerning packet delivery ratio to 25 , end to end delay down to 100 ms, low radio transmissions required for per packet delivery to 6 mJ, and most importantly, it improves the robustness and scalability of data communication process in thewhole network.


Continual learning using hash routed convolutional neural networks. Continual learning could shift the machine learning paradigm from data centric to model centric. A continual learning model needs to scale efficiently to handle semantically different datasets, while avoiding unnecessary growth. We introduce hash routed convolutional neural networks a group of convolutional units where data flows dynamically. Feature maps are compared using feature hashing and similar data is routed to the same units. A hash routed network provides excellent plasticity thanks to its routed nature, while generating stable features through the use of orthogonal feature hashing. Each unit evolves separately and new units can be added (to be used only when necessary). Hash routed networks achieve excellent performance across a variety of typical continual learning benchmarks without storing raw data and train using only gradient descent. Besides providing a continual learning framework for supervised tasks with encouraging results, our model can be used for unsupervised or reinforcement learning.


Discrete Richman bidding Scoring Games. We study zero sum (combinatorial) games, within the framework of so called Richman auctions (Lazarus et al. 1996) namely, we modify the alternating play scoring ruleset Cumulative Subtraction (CS) (Cohensius et al. 2019), to a discrete bidding scheme (similar to Develin and Payne 2010). Players bid to move and the player with the highest bid wins the move, and hands over the winning bidding amount to the other player. The new game is dubbed Bidding Cumulative Subtraction (BCS). In so called unitary games, players remove exactly one item out of a single heap of identical items, until the heap is empty, and their actions contribute to a common score, which increases or decreases by one unit depending on whether the maximizing player won the turn or not. We show that there is a unique bidding equilibrium for a much larger class of games, that generalize standard scoring play in the literature. We prove that for all sufficiently large heap sizes, the equilibrium outcomes of unitary BCS are eventually periodic, with period 2, and we show that the periodicity appears at the latest for heaps of sizes quadratic in the total budget.


SPM Tracker Series Parallel Matching for Real Time Visual Object Tracking. The greatest challenge facing visual object tracking is the simultaneous requirements on robustness and discrimination power. In this paper, we propose a SiamFC based tracker, named SPM Tracker, to tackle this challenge. The basic idea is to address the two requirements in two separate matching stages. Robustness is strengthened in the coarse matching (CM) stage through generalized training while discrimination power is enhanced in the fine matching (FM) stage through a distance learning network. The two stages are connected in series as the input proposals of the FM stage are generated by the CM stage. They are also connected in parallel as the matching scores and box location refinements are fused to generate the final results. This innovative series parallel structure takes advantage of both stages and results in superior performance. The proposed SPM Tracker, running at 120fps on GPU, achieves an AUC of 0.687 on OTB 100 and an EAO of 0.434 on VOT 16, exceeding other real time trackers by a notable margin.


An integral theorem for plurisubharmonic functions. In this paper, we prove an integral theorem for Cegrell class mathcal F (f) and use this result to study the mathcal F equivalence relation.


Long time asymptotic behavior for an extended modified Korteweg de Vries equation. We investigate an integrable extended modified Korteweg de Vries equation on the line with the initial value belonging to the Schwartz space. By performing the nonlinear steepest descent analysis of an associated matrix Riemann Hilbert problem, we obtain the explicit leading order asymptotics of the solution of this initial value problem as time t goes to infinity. For a special case alpha 0 , we present the asymptotic formula of the solution to the extended modified Korteweg de Vries equation in region mathcal P (x,t) in bfR 2 0 x leq Mt frac 1 5 ,t geq3 in terms of the solution of a fourth order Painlev 'e II equation.


Birational types of algebraic orbifolds. We introduce a variant of the birational symbols group of Kontsevich, Pestun, and the second author, and use this to define birational invariants of algebraic orbifolds.


SUNDIALS Multiphysics MPIManyVector Performance Testing. In this report we document performance test results on a SUNDIALS based multiphysics demonstration application. We aim to assess the large scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically (a) SUNDIALS' new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources, (b) ARKode's new multirate integration module, MRIStep, allowing high order accurate calculations that subcycle "fast" processes within "slow" ones, (c) SUNDIALS' new flexible linear solver interfaces, that allow streamlined specification of problem specific linear solvers, and (d) SUNDIALS' new N Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS' generic SUNNonlinearSolver interface and accelerator enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).


Extraction of the translational Eucken factor from light scattering by molecular gas. Although the thermal conductivity of molecular gases can be measured straightforwardly and accurately, it is difficult to experimentally determine its separate contributions from the translational and internal motions of gas molecules. Yet this information is critical in rarefied gas dynamics as the rarefaction effects corresponding to these motions are different. In this paper, we propose a novel methodology to extract the translational thermal conductivity (or equivalently, the translational Eucken factor) of molecular gases from the Rayleigh Brillouin scattering (RBS) experimental data. From the numerical simulation of the cite LeiJFM2015 model we find that, in the kinetic regime, in addition to bulk viscosity, the RBS spectrum is sensitive to the translational Eucken factor, even when the total thermal conductivity is fixed. Thus it is not only possible to extract the bulk viscosity, but also the translational Eucken factor of molecular gases from RBS light scattering spectra measurements. Such experiments bear the additional advantage that gas surface interactions do not affect the measurements. For the first time, bulk viscosities (due to the rotational relaxation of gas molecules only) and translational Eucken factors of N2, CO2 and SF6 are simultaneously extracted from RBS experiments.


Representation Formula for Viscosity Solution to a PDE Problem involving Pucci's Extremal Operator. We provide a representation formula for viscosity solutions to an elliptic Dirichlet problem involving Pucci's extremal operators. This is done through a dynamic programming principle derived from Denis, Hu and Peng (2010). The formula can be seen as a nonlinear extension of the Feynman Kac formula.


Force Feedback Effects on Single Molecule Hopping and Pulling Experiments. Single molecule experiments with optical tweezers have become an important tool to study the properties and mechanisms of biological systems, such as cells and nucleic acids. In particular, force unzipping experiments have been used to extract the thermodynamics and kinetics of folding and unfolding reactions. In hopping experiments, a molecule executes transitions between the unfolded and folded states at a preset value of the force (constant force mode CFM under force feedback) or trap position (passive mode PM without feedback) and the force dependent kinetic rates extracted from the lifetime of each state (CFM) and the rupture force distributions (PM) using the Bell Evans model. However, hopping experiments in the CFM are known to overestimate molecular distances and folding free energies for fast transitions compared to the response time of the feedback. In contrast, kinetic rate measurements from pulling experiments have been mostly done in the PM while the CFM is seldom implemented in pulling protocols. Here, we carry out hopping and pulling experiments in a short DNA hairpin in the PM and CFM at three different temperatures (6 circ C, 25 circ C and 45 circ C) exhibiting largely varying kinetic rates. As expected, we find that equilibrium hopping experiments in the CFM and PM perform well at 6 circ C (where kinetics is slow) whereas the CFM overestimates molecular parameters at 45 circ C (where kinetics is fast). In contrast, nonequilibrium pulling experiments perform well in both modes at all temperatures. This demonstrates that the same kind of feedback algorithm in the CFM leads to more reliable determination of the folding reaction parameters in irreversible pulling experiments.


Learning State Dependent Losses for Inverse Dynamics Learning. Being able to quickly adapt to changes in dynamics is paramount in model based control for object manipulation tasks. In order to influence fast adaptation of the inverse dynamics model's parameters, data efficiency is crucial. Given observed data, a key element to how an optimizer updates model parameters is the loss function. In this work, we propose to apply meta learning to learn structured, state dependent loss functions during a meta training phase. We then replace standard losses with our learned losses during online adaptation tasks. We evaluate our proposed approach on inverse dynamics learning tasks, both in simulation and on real hardware data. In both settings, the structured and state dependent learned losses improve online adaptation speed, when compared to standard, state independent loss functions.


Three dimensional compacted optical waveguide couplers designed by invariant engineerings. Due to the limitations either on the sizes of devices and signal routing channels, the current planar integrated optical waveguide circuits await for the further developments into the three dimensional (3D) integrations, although their designs and fabrications are still challenges. In this paper we demonstrate an analytical method, basing on the invariant engineering, to overcome the complication in the usual method by numerically solving the relevant 3D coupled mode equations for designing various 3D optical waveguide devices such as the typical couplers. Our method is based on the quantum optical analogy, i.e., the Maxwell equation for the electrcomagnetic wave prorogating along the waveguide structure in the spatial domain is formally similar to the Schr "odinger equation for the evolving quantum state in the time domain. We find that the spatial domain invariants can be effectively constructed to solve the 3D coupled mode equations, analogously to solve the dynamical evolutions of quantum systems in the time domain. As a consequence, as long as appropriately set the coupling parameters between the 3D interconnected waveguides, the 3D three waveguide couplers could be designed for various desirably power divisions. As the invariant method is a natural shortcut to the adiabaticity, the compacted devices designed by the invariant based engineerings are robust against the coupling coefficient variations and the coupler lengths.


Global weak solutions for Landau Lifshitz flows and heat flows associated to micromagnetic energy functional. We follow the idea of Wang cite W to show the existence of global weak solutions to the Cauchy problems of Landau Lifshtiz type equations and related heat flows from a n dimensional Euclidean domain Om or a n dimensional closed Riemannian manifold M into a 2 dimensional unit sphere U 2 . Our conclusions extend a series of related results obtained in the previous literature.


The Covid 19 pandemic's effects on poor rural dwellers in sub Saharan Africa A case study of access to basic clean water, sanitary systems and hand washing facilities. The fear of the invisible but prevalent Coronavirus (COVID 19), disease cannot be overemphasized since there is the potential possibility of it wiping out the entire world population within a few months if adequate and quick steps are not taken to curb this menace, and the sub Saharan African (SSAn) region is no exception. It is evident that water, as an essential daily commodity, has long been in a state of emergency in SSAn nations, which is largely attributed to decades of neglect by the successive governments, because it has not been possible to separate the existing bond between water, health, livelihood and the economy. The laudable Millennium Development Goals (MDGs) proposed by the United Nations had yet to achieve the stated objective of improving the standards of living and health conditions of the rural communities in the SSAn region before the COVID 19 pandemic outbreak. This failure has been masked by a sort of delusion in which the people of this region are subjected to the hardship of searching for clean and healthy water in their own ponds, rivers, streams and shallow hand dug local wells on a continuous basis. Less than 17 of the rural population in all the SSAn communities can access basic hand washing facilities and sanitation systems. The total water productivity, as measured by the Gross Domestic Product (GDP) per cubic meter of total freshwater withdrawn, for the people was less than 5 GDP.


A caricature of dilation theory. We present a set theoretic version of some basic dilation results of operator theory. The results we have considered are Wold decomposition, Halmos dilation, Sz. Nagy dilation, inter twining lifting, commuting and non commuting dilations, BCL theorem etc. We point out some natural generalizations and variations.


Holomorphic immersions of bi disks into 9 dimensional real hypersurfaces with Levi signature (2, 2) . Inspired by an article of R. Bryant on holomorphic immersions of unit disks into Lorentzian CR manifolds, we discuss the application of Cartan's method to the question of the existence of bi disk mathbb D 2 in a smooth 9 dimensional real analytic real hypersurface M 9 subset mathbb C 5 with Levi signature (2,2) passing through a fixed point. The result is that the lift to M 9 times U(2) of the image of the bi disk in M 9 must lie in the zero set of two complex valued functions in M 9 times U(2) . We then provide an example where one of the functions does not identically vanish, thus obstructing holomorphic immersions.


Interacting Langevin Diffusions Gradient Structure And Ensemble Kalman Sampler. Solving inverse problems without the use of derivatives or adjoints of the forward model is highly desirable in many applications arising in science and engineering. In this paper, we propose a new version of such a methodology, a framework for its analysis, and numerical evidence of the practicality of the method proposed. Our starting point is an ensemble of over damped Langevin diffusions which interact through a single preconditioner computed as the empirical ensemble covariance. We demonstrate that the nonlinear Fokker Planck equation arising from the mean field limit of the associated stochastic differential equation (SDE) has a novel gradient flow structure, built on the Wasserstein metric and the covariance matrix of the noisy flow. Using this structure, we investigate large time properties of the Fokker Planck equation, showing that its invariant measure coincides with that of a single Langevin diffusion, and demonstrating exponential convergence to the invariant measure in a number of settings. We introduce a new noisy variant on ensemble Kalman inversion (EKI) algorithms found from the original SDE by replacing exact gradients with ensemble differences this defines the ensemble Kalman sampler (EKS). Numerical results are presented which demonstrate its efficacy as a derivative free approximate sampler for the Bayesian posterior arising from inverse problems.


Analyzing Cyber Physical Systems from the Perspective of Artificial Intelligence. Principles of modern cyber physical system (CPS) analysis are based on analytical methods that depend on whether safety or liveness requirements are considered. Complexity is abstracted through different techniques, ranging from stochastic modelling to contracts. However, both distributed heuristics and Artificial Intelligence (AI) based approaches as well as the user perspective or unpredictable effects, such as accidents or the weather, introduce enough uncertainty to warrant reinforcement learning based approaches. This paper compares traditional approaches in the domain of CPS modelling and analysis with the AI researcher perspective to exploring unknown complex systems.


A Note on Matrix Valued Truncated Toeplitz Operators. We prove a spatial Isomorphism between two spaces of matrix valued truncated Toeplitz operators.


Refinement invariance of intersection (co)homologies. We study the refinement invariance of several intersection (co)homologies existing in the literature. These (co)homologies have been introduced in order to establish the Poincar 'e Duality in variousl contexts. We found the classical topological invariance of the intersection homology and also various refinement invariance results already proved.


Batalin Vilkovisky formalism in the p adic Dwork theory. The goal of this article is to develop BV (Batalin Vilkovisky) formalism in the p adic Dwork theory. Based on this formalism, we explicitly construct a p adic dGBV algebra (differential Gerstenhaber Batalin Vilkovisky algebra) for a smooth projective complete intersection variety X over a finite field, whose cohomology gives the p adic Dwork cohomology of X , and its cochain endomorphism (the p adic Dwork Frobenius operator) which encodes the information of the zeta function X . As a consequence, we give a modern deformation theoretic interpretation of Dwork's theory of the zeta function of X and derive a formula for the p adic Dwork Frobenius operator in terms of homotopy Lie morphisms and the Bell polynomials.


Scalability of High Performance PDE Solvers. Performance tests and analyses are critical to effective HPC software development and are central components in the design and implementation of computational algorithms for achieving faster simulations on existing and future computing architectures for large scale application problems. In this paper, we explore performance and space time trade offs for important compute intensive kernels of large scale numerical solvers for PDEs that govern a wide range of physical applications. We consider a sequence of PDE motivated bake off problems designed to establish best practices for efficient high order simulations across a variety of codes and platforms. We measure peak performance (degrees of freedom per second) on a fixed number of nodes and identify effective code optimization strategies for each architecture. In addition to peak performance, we identify the minimum time to solution at 80 parallel efficiency. The performance analysis is based on spectral and p type finite elements but is equally applicable to a broad spectrum of numerical PDE discretizations, including finite difference, finite volume, and h type finite elements.


Goal oriented error estimation and adaptivity in MsFEM computations. We introduce a goal oriented strategy for multiscale computations performed using the Multiscale Finite Element Method (MsFEM). In a previous work, we have shown how to use, in the MsFEM framework, the concept of Constitutive Relation Error (CRE) to obtain a guaranteed and fully computable a posteriori error estimate in the energy norm (as well as error indicators on various error sources). Here, the CRE concept is coupled with the solution of an adjoint problem to control the error and drive an adaptive procedure with respect to a given output of interest. Furthermore, a local and non intrusive enrichment technique is proposed to enhance the accuracy of error bounds. The overall strategy, which is fully automatic and robust, enables to reach an appropriate trade off between certified reliability and computational cost in the MsFEM context. The performances of the proposed method are investigated on several illustrative numerical test cases. In particular, the error estimation is observed to be very accurate, yielding a very efficient adaptive procedure.


Gradient Descent Averaging and Primal dual Averaging for Strongly Convex Optimization. Averaging scheme has attracted extensive attention in deep learning as well as traditional machine learning. It achieves theoretically optimal convergence and also improves the empirical model performance. However, there is still a lack of sufficient convergence analysis for strongly convex optimization. Typically, the convergence about the last iterate of gradient descent methods, which is referred to as individual convergence, fails to attain its optimality due to the existence of logarithmic factor. In order to remove this factor, we first develop gradient descent averaging (GDA), which is a general projection based dual averaging algorithm in the strongly convex setting. We further present primal dual averaging for strongly convex cases (SC PDA), where primal and dual averaging schemes are simultaneously utilized. We prove that GDA yields the optimal convergence rate in terms of output averaging, while SC PDA derives the optimal individual convergence. Several experiments on SVMs and deep learning models validate the correctness of theoretical analysis and effectiveness of algorithms.


Sharing Attention Weights for Fast Transformer. Recently, the Transformer machine translation system has shown strong results by stacking attention layers on both the source and target language sides. But the inference of this model is slow due to the heavy use of dot product attention in auto regressive decoding. In this paper we speed up Transformer via a fast and lightweight attention model. More specifically, we share attention weights in adjacent layers and enable the efficient re use of hidden states in a vertical manner. Moreover, the sharing policy can be jointly learned with the MT model. We test our approach on ten WMT and NIST OpenMT tasks. Experimental results show that it yields an average of 1.3X speed up (with almost no decrease in BLEU) on top of a state of the art implementation that has already adopted a cache for fast inference. Also, our approach obtains a 1.8X speed up when it works with the textsc Aan model. This is even 16 times faster than the baseline with no use of the attention cache.


Classification of Open and Closed Convex Codes on Five Neurons. Neural codes, represented as collections of binary strings, encode neural activity and show relationships among stimuli. Certain neurons, called place cells, have been shown experimentally to fire in convex regions in space. A natural question to ask is Which neural codes can arise as intersection patterns of convex sets? While past research has established several criteria, complete conditions for convexity are not yet known for codes with more than four neurons. We classify all neural codes with five neurons as convex non convex codes. Furthermore, we investigate which of these codes can be represented by open versus closed convex sets. Interestingly, we find a code which is an open but not closed convex code and demonstrate a minimal example for this phenomenon.


Multi Point Synchronization for Fog Controlled Internet of Things. This paper presents a fog resident controller architecture for synchronizing the operations of large collections of Internet of Things (IoT) such as drones, Internet of Vehicles, etc. Synchronization in IoT is grouped into different classes, use cases identified and multi point synchronous scheduling algorithms are developed to schedule tasks with varying timing requirements strict (synchronous) and relaxed (asynchronous and local) onto a bunch of worker nodes that are coordinated by a fog resident controller in the presence of disconnections and worker failures. The algorithms use time based or component based redundancy to cope with failures and embed a publish subscribe message update scheme to reduce the message overhead at the controller as the number of workers increase. The performance of the algorithms are evaluated using trace driven experiments and practicability is shown by implementing the time based redundancy synchronous scheduling algorithm in JAMScript a polyglot programming platform for Cloud of Things and report initial findings.


A Case for Reversible Coherence Protocol. We propose the first Reversible Coherence Protocol (RCP), a new protocol designed from ground up that enables invisible speculative load. RCP takes a bold approach by including the speculative loads and merge purge operation in the interface between processor and cache coherence, and allowing them to participate in the coherence protocol. It means, speculative load, ordinary load store, and merge purge can all affect the state of a given cache line. RCP is the first coherence protocol that enables the commit and squash of the speculative load among distributed cache components in a general memory hierarchy. RCP incurs an average slowdown of (3.0 ,8.3 ,7.4 ) on (SPEC2006,SPEC2017,PARSEC), which is lower compared to (26.5 ,12 ,18.3 ) in InvisiSpec and (3.2 ,9.4 ,24.2 ) in CleanupSpec. The coherence traffic overhead is on average 46 , compared to 40 and 27 of InvisiSpec and CleanupSpec, respectively. Even with higher traffic overhead ( 46 ), the performance overhead of RCP is lower than InvisiSpec and comparable to CleanupSpec. It reveals a key advantage of RCP the coherence actions triggered by the merge and purge operations are not in the critical path of the execution and can be performed in the cache hierarchy concurrently with processor execution


Strong Progress for Session Typed Processes in a Linear Metalogic with Circular Proofs. We introduce an infinitary first order linear logic with least and greatest fixed points. To ensure cut elimination, we impose a validity condition on infinite derivations. Our calculus is designed to reason about rich signatures of mutually defined inductive and coinductive linear predicates. In a major case study we use it to prove the strong progress property for binary session typed processes under an asynchronous communication semantics. As far as we are aware, this is the first proof of this property.


POSITION PAPER Credibility of In Silico Trial Technologies A Theoretical Framing. Different research communities have developed various approaches to assess the credibility of predictive models. Each approach usually works well for a specific type of model, and under some epistemic conditions that are normally satisfied within that specific research domain. Some regulatory agencies recently started to consider evidences of safety and efficacy on new medical products obtained using computer modelling and simulation (which is referred to as In Silico Trials) this has raised the attention in the computational medicine research community on the regulatory science aspects of this emerging discipline. But this poses a foundational problem in the domain of biomedical research the use of computer modelling is relatively recent, without a widely accepted epistemic framing for problem of model credibility. Also, because of the inherent complexity of living organisms, biomedical modellers tend to use a variety of modelling methods, sometimes mixing them in the solution of a single problem. In such context merely adopting credibility approaches developed within other research community might not be appropriate. In this position paper we propose a theoretical framing for the problem of assessing the credibility of a predictive models for In Silico Trials, which accounts for the epistemic specificity of this research field and is general enough to be used for different type of models.


On the morphology of two dimensional laminar vortex streets behind triangles. The two dimensional laminar vortex streets behind a triangle have two morphologically distinct structures depending on the Reynolds number and the aspect ratio of the triangle. These two structures are the conventional structure and the separated rows structure, where the latter is characterized by a thin layer of irrotational fluid between two vortex rows. In this paper, by means of numerical simulation, we find that the separated rows structure occurs when the thickness of boundary layers is less than 25 of their separation distance. We also show from the linear stability analysis that the criterion is related to the coupling of two boundary layers in producing unstable modes.


On Pseudo Umbilical Spacelike Submanifolds in Indefinite Space Form Mn p p (c). In the present note, first we derive an intrinsic inequality for Pseudo umbilical spacelike submanifold in an indefinite space form. We use this inequality to show that such submanifold is totally geodesic. In the rest part of this paper, using a result of Aiyama 1 , we prove that Pseudo umbilical spacelike subamnifold is totally umbilical.


Unfolding cubes nets, packings, partitions, chords. We show that every ridge unfolding of an n cube is without self overlap, yielding a valid net. The results are obtained by developing machinery that translates cube unfolding into combinatorial frameworks. Moreover, the geometry of the bounding boxes of these cube nets are classified using integer partitions, as well as the combinatorics of path unfoldings seen through the lens of chord diagrams.


A Survey on Semantic Parsing from the perspective of Compositionality. Different from previous surveys in semantic parsing (Kamath and Das, 2018) and knowledge base question answering(KBQA)(Chakraborty et al., 2019 Zhu et al., 2019 Hoffner et al., 2017) we try to takes a different perspective on the study of semantic parsing. Specifically, we will focus on (a)meaning composition from syntactical structure(Partee, 1975), and (b) the ability of semantic parsers to handle lexical variation given the context of a knowledge base (KB). In the following section after an introduction of the field of semantic parsing and its uses in KBQA, we will describe meaning representation using grammar formalism CCG (Steedman, 1996). We will discuss semantic composition using formal languages in Section 2. In section 3 we will consider systems that uses formal languages e.g. lambda calculus (Steedman, 1996), lambda DCS (Liang, 2013). Section 4 and 5 consider semantic parser using structured language for logical form. Section 6 is on different benchmark datasets ComplexQuestions (Bao et al.,2016) and GraphQuestions (Su et al., 2016) that can be used to evaluate semantic parser on their ability to answer complex questions that are highly compositional in nature.


Free2Shard Adaptive adversary resistant sharding via Dynamic Self Allocation. Propelled by the growth of large scale blockchain deployments, much recent progress has been made in designing sharding protocols that achieve throughput scaling linearly in the number of nodes. However, existing protocols are not robust to an adversary adaptively corrupting a fixed fraction of nodes. In this paper, we propose Free2Shard a new architecture that achieves near linear scaling while being secure against a fully adaptive adversary. The focal point of this architecture is a dynamic self allocation algorithm that lets users allocate themselves to shards in response to adversarial action, without requiring a central or cryptographic proof. This architecture has several attractive features unusual for sharding protocols, including (a) the ability to handle the regime of large number of shards (relative to the number of nodes) (b) heterogeneous shard demands (c) requiring only a small minority to follow the self allocation (d) asynchronous shard rotation (e) operation in a purely identity free proof of work setting. The key technical contribution is a deep mathematical connection to the classical work of Blackwell in dynamic game theory.


Jordan blocks of nilpotent elements in some irreducible representations of classical groups in good characteristic. Let G be a classical group with natural module V and Lie algebra mathfrak g over an algebraically closed field K of good characteristic. For rational irreducible representations f G rightarrow operatorname GL (W) occurring as composition factors of V otimes V , wedge 2(V) , and S 2(V) , we describe the Jordan normal form of mathrm d f(e) for all nilpotent elements e in mathfrak g . The description is given in terms of the Jordan block sizes of the action of e on V otimes V , wedge 2(V) , and S 2(V) , for which recursive formulae are known. Our results are in analogue to earlier work (Proc. Amer. Math. Soc., 147 (2019) 4205 4219), where we considered these same representations and described the Jordan normal form of f(u) for every unipotent element u in G .


Digitalization of exceptional simple Lie algebras into matrices over complex numbers. We give the images of the adjoint representations of exceptional simple Lie algebras by matrices over complex numbers. Next, we digitalize these matrices by the use of Maxima, which is a computer algebra system. These digitalized matrices are provided by using Maxima's function. We prove that these digitalized matrices are closed by Lie bracket operations and make up simple Lie algebras. Moreover, to prove the type of the exceptional simple Lie algebra, we calculate the root system using Maxima for Lie bracket operations as matrix calculations. We show some examples of classical Lie subalgebras of these digitalized matrices.


Quadtree Generating Networks Efficient Hierarchical Scene Parsing with Sparse Convolutions. Semantic segmentation with Convolutional Neural Networks is a memory intensive task due to the high spatial resolution of feature maps and output predictions. In this paper, we present Quadtree Generating Networks (QGNs), a novel approach able to drastically reduce the memory footprint of modern semantic segmentation networks. The key idea is to use quadtrees to represent the predictions and target segmentation masks instead of dense pixel grids. Our quadtree representation enables hierarchical processing of an input image, with the most computationally demanding layers only being used at regions in the image containing boundaries between classes. In addition, given a trained model, our representation enables flexible inference schemes to trade off accuracy and computational cost, allowing the network to adapt in constrained situations such as embedded devices. We demonstrate the benefits of our approach on the Cityscapes, SUN RGBD and ADE20k datasets. On Cityscapes, we obtain an relative 3 mIoU improvement compared to a dilated network with similar memory consumption and only receive a 3 relative mIoU drop compared to a large dilated network, while reducing memory consumption by over 4 times .


BiPart A Parallel and Deterministic Multilevel Hypergraph Partitioner. Hypergraph partitioning is used in many problem domains including VLSI design, linear algebra, Boolean satisfiability, and data mining. Most versions of this problem are NP complete or NP hard, so practical hypergraph partitioners generate approximate partitioning solutions for all but the smallest inputs. One way to speed up hypergraph partitioners is to exploit parallelism. However, existing parallel hypergraph partitioners are not deterministic, which is considered unacceptable in domains like VLSI design where the same partitions must be produced every time a given hypergraph is partitioned. In this paper, we describe BiPart, the first deterministic, parallel hypergraph partitioner. Experimental results show that BiPart outperforms state of the art hypergraph partitioners in runtime and partition quality while generating partitions deterministically.


Underwater Fish Detection with Weak Multi Domain Supervision. Given a sufficiently large training dataset, it is relatively easy to train a modern convolution neural network (CNN) as a required image classifier. However, for the task of fish classification and or fish detection, if a CNN was trained to detect or classify particular fish species in particular background habitats, the same CNN exhibits much lower accuracy when applied to new unseen fish species and or fish habitats. Therefore, in practice, the CNN needs to be continuously fine tuned to improve its classification accuracy to handle new project specific fish species or habitats. In this work we present a labelling efficient method of training a CNN based fish detector (the Xception CNN was used as the base) on relatively small numbers (4,000) of project domain underwater fish no fish images from 20 different habitats. Additionally, 17,000 of known negative (that is, missing fish) general domain (VOC2012) above water images were used. Two publicly available fish domain datasets supplied additional 27,000 of above water and underwater positive fish images. By using this multi domain collection of images, the trained Xception based binary (fish not fish) classifier achieved 0.17 false positives and 0.61 false negatives on the project's 20,000 negative and 16,000 positive holdout test images, respectively. The area under the ROC curve (AUC) was 99.94 .


Nowhere to Hide Cross modal Identity Leakage between Biometrics and Devices. Along with the benefits of Internet of Things (IoT) come potential privacy risks, since billions of the connected devices are granted permission to track information about their users and communicate it to other parties over the Internet. Of particular interest to the adversary is the user identity which constantly plays an important role in launching attacks. While the exposure of a certain type of physical biometrics or device identity is extensively studied, the compound effect of leakage from both sides remains unknown in multi modal sensing environments. In this work, we explore the feasibility of the compound identity leakage across cyber physical spaces and unveil that co located smart device IDs (e.g., smartphone MAC addresses) and physical biometrics (e.g., facial vocal samples) are side channels to each other. It is demonstrated that our method is robust to various observation noise in the wild and an attacker can comprehensively profile victims in multi dimension with nearly zero analysis effort. Two real world experiments on different biometrics and device IDs show that the presented approach can compromise more than 70 of device IDs and harvests multiple biometric clusters with 94 purity at the same time.


Stochastic Proximal Gradient Methods for Nonconvex Problems in Hilbert Spaces. For finite dimensional problems, stochastic approximation methods have long been used to solve stochastic optimization problems. Their application to infinite dimensional problems is less understood, particularly for nonconvex objectives. This paper presents convergence results for the stochastic proximal gradient method applied to Hilbert spaces, motivated by optimization problems with partial differential equation (PDE) constraints with random inputs and coefficients. We study stochastic algorithms for nonconvex and nonsmooth problems, where the nonsmooth part is convex and the nonconvex part is the expectation, which is assumed to have a Lipschitz continuous gradient. The optimization variable is an element of a Hilbert space. We show almost sure convergence of strong limit points of the random sequence generated by the algorithm to stationary points. We demonstrate the stochastic proximal gradient algorithm on a tracking type functional with a L 1 penalty term constrained by a semilinear PDE and box constraints, where input terms and coefficients are subject to uncertainty. We verify conditions for ensuring convergence of the algorithm and show a simulation.


Root subgroups on affine spherical varieties. Given a connected reductive algebraic group G and a Borel subgroup B subseteq G , we study B normalized one parameter additive group actions on affine spherical G varieties. We establish basic properties of such actions and their weights and discuss many examples exhibiting various features. We propose a construction of such actions that generalizes the well known construction of normalized one parameter additive group actions on affine toric varieties. Using this construction, for every affine horospherical G variety X we obtain a complete description of all G normalized one parameter additive group actions on X and show that the open G orbit in X can be connected with every G stable prime divisor via a suitable choice of a B normalized one parameter additive group action. Finally, when G is of semisimple rank 1 , we obtain a complete description of all B normalized one parameter additive group actions on affine spherical G varieties having an open orbit of a maximal torus T subseteq B .


Equivariant splitting of the Hodge de Rham exact sequence. Let X be an algebraic curve with an action of a finite group G over a field k . We show that if the Hodge de Rham short exact sequence of X splits G equivariantly then the action of G on X is weakly ramified. In particular, this generalizes the result of K " o ck and Tait for hyperelliptic curves. We discuss also converse statements and tie this problem to lifting coverings of curves to the ring of Witt vectors of length 2 .


On limitations of learning algorithms in competitive environments. We discuss conceptual limitations of generic learning algorithms pursuing adversarial goals in competitive environments, and prove that they are subject to limitations that are analogous to the constraints on knowledge imposed by the famous theorems of G "odel and Turing. These limitations are shown to be related to intransitivity, which is commonly present in competitive environments.


LP Based Upper and Lower Bounds for Ces aro and Abel Limits of the Optimal Values in Problems of Control of Stochastic Discrete Time Systems. In this paper, we study asymptotic properties of problems of control of stochastic discrete time systems (also known as Markov decision processes) with time averaging and time discounting optimality criteria, and we establish that the Ces aro and Abel limits of the optimal values in such problems can be evaluated with the help of a certain infinite dimensional linear programming problem and its dual.


Higher order Comparisons of Sentence Encoder Representations. Representational Similarity Analysis (RSA) is a technique developed by neuroscientists for comparing activity patterns of different measurement modalities (e.g., fMRI, electrophysiology, behavior). As a framework, RSA has several advantages over existing approaches to interpretation of language encoders based on probing or diagnostic classification namely, it does not require large training samples, is not prone to overfitting, and it enables a more transparent comparison between the representational geometries of different models and modalities. We demonstrate the utility of RSA by establishing a previously unknown correspondence between widely employed pretrained language encoders and human processing difficulty via eye tracking data, showcasing its potential in the interpretability toolbox for neural models


An Improved and Parallel Version of a Scalable Algorithm for Analyzing Time Series Data. Today, very large amounts of data are produced and stored in all branches of society including science. Mining these data meaningfully has become a considerable challenge and is of the broadest possible interest. The size, both in numbers of observations and dimensionality thereof, requires data mining algorithms to possess time complexities with both variables that are linear or nearly linear. One such algorithm, see Comput. Phys. Commun. 184, 2446 2453 (2013), arranges observations into a sequence called the progress index. The progress index steps through distinct regions of high sampling density sequentially. By means of suitable annotations, it allows a compact representation of the behavior of complex systems, which is encoded in the original data set. The only essential parameter is a notion of distance between observations. Here, we present the shared memory parallelization of the key step in constructing the progress index, which is the calculation of an approximation of the minimum spanning tree of the complete graph of observations. We demonstrate that excellent parallel efficiencies are obtained for up to 72 logical (CPU) cores. In addition, we introduce three conceptual advances to the algorithm that improve its controllability and the interpretability of the progress index itself.


Boundary Conditions for Continuum Simulations of Wall bounded Kinetic Plasmas. Continuum kinetic simulations of plasmas, where the distribution function of the species is directly discretized in phase space, permits fully kinetic simulations without the statistical noise of particle in cell methods. Recent advances in numerical algorithms have made continuum kinetic simulations computationally competitive. This work presents the first continuum kinetic description of high fidelity wall boundary conditions that utilize the readily available particle distribution function. The boundary condition is realized through a reflection function that can capture a wide range of cases from simple specular reflection to more involved first principles models. Examples with detailed discontinuous Galerkin implementation are provided for secondary electron emission using phenomenological and first principles quantum mechanical models. Results presented in this work demonstrate the effect of secondary electron emission on a classical plasma sheath.


A Hybrid Algorithm Based Robust Big Data Clustering for Solving Unhealthy Initialization, Dynamic Centroid Selection and Empty clustering Problems with Analysis. Big Data is a massive volume of both structured and unstructured data that is too large and it also difficult to process using traditional techniques. Clustering algorithms have developed as a powerful learning tool that can exactly analyze the volume of data that produced by modern applications. Clustering in data mining is the grouping of a particular set of objects based on their characteristics. The main aim of clustering is to classified data into clusters such that objects are grouped in the same clusters when they are corresponding according to similarities and features mainly. Till now, K MEANS is the best utilized calculation connected in a wide scope of zones to recognize gatherings where cluster separations are a lot than between gathering separations. Our developed algorithm works with K MEANS for high quality clustering during clustering from big data. Our proposed algorithm EG K MEANS Extended Generation K MEANS solves mainly three issues of K MEANS unhealthy initialization, dynamic centroid selection and empty clustering. It ensures the best way of preventing unhealthy initialization, dynamic centroid selection and empty clustering problems for getting high quality clustering.


TAPA MVS Textureless Aware PAtchMatch Multi View Stereo. One of the most successful approaches in Multi View Stereo estimates a depth map and a normal map for each view via PatchMatch based optimization and fuses them into a consistent 3D points cloud. This approach relies on photo consistency to evaluate the goodness of a depth estimate. It generally produces very accurate results however, the reconstructed model often lacks completeness, especially in correspondence of broad untextured areas where the photo consistency metrics are unreliable. Assuming the untextured areas piecewise planar, in this paper we generate novel PatchMatch hypotheses so to expand reliable depth estimates in neighboring untextured regions. At the same time, we modify the photo consistency measure such to favor standard or novel PatchMatch depth hypotheses depending on the textureness of the considered area. We also propose a depth refinement step to filter wrong estimates and to fill the gaps on both the depth maps and normal maps while preserving the discontinuities. The effectiveness of our new methods has been tested against several state of the art algorithms in the publicly available ETH3D dataset containing a wide variety of high and low resolution images.


Constructing cocyclic Hadamard matrices of order 4p. Cocyclic Hadamard matrices (CHMs) were introduced by de Launey and Horadam as a class of Hadamard matrices with interesting algebraic properties. 'O Cath 'ain and R "oder described a classification algorithm for CHMs of order 4n based on relative difference sets in groups of order 8n this led to the classification of all CHMs of order at most 36. Based on work of de Launey and Flannery, we describe a classification algorithm for CHMs of order 4p with p a prime we prove refined structure results and provide a classification for p leqslant 13 . Our analysis shows that every CHM of order 4p with p equiv 1 bmod 4 is equivalent to a Hadamard matrix with one of five distinct block structures, including Williamson type and (transposed) Ito matrices. If p equiv 3 bmod 4 , then every CHM of order 4p is equivalent to a Williamson type or (transposed) Ito matrix.


Simple and optimal methods for stochastic variational inequalities, II Markovian noise and policy evaluation in reinforcement learning. The focus of this paper is on stochastic variational inequalities (VI) under Markovian noise. A prominent application of our algorithmic developments is the stochastic policy evaluation problem in reinforcement learning. Prior investigations in the literature focused on temporal difference (TD) learning by employing nonsmooth finite time analysis motivated by stochastic subgradient descent leading to certain limitations. These encompass the requirement of analyzing a modified TD algorithm that involves projection to an a priori defined Euclidean ball, achieving a non optimal convergence rate and no clear way of deriving the beneficial effects of parallel implementation. Our approach remedies these shortcomings in the broader context of stochastic VIs and in particular when it comes to stochastic policy evaluation. We developed a variety of simple TD learning type algorithms motivated by its original version that maintain its simplicity, while offering distinct advantages from a non asymptotic analysis point of view. We first provide an improved analysis of the standard TD algorithm that can benefit from parallel implementation. Then we present versions of a conditional TD algorithm (CTD), that involves periodic updates of the stochastic iterates, which reduce the bias and therefore exhibit improved iteration complexity. This brings us to the fast TD (FTD) algorithm which combines elements of CTD and the stochastic operator extrapolation method of the companion paper. For a novel index resetting policy FTD exhibits the best known convergence rate. We also devised a robust version of the algorithm that is particularly suitable for discounting factors close to 1.


An Attention Based System for Damage Assessment Using Satellite Imagery. When disaster strikes, accurate situational information and a fast, effective response are critical to save lives. Widely available, high resolution satellite images enable emergency responders to estimate locations, causes, and severity of damage. Quickly and accurately analyzing the extensive amount of satellite imagery available, though, requires an automatic approach. In this paper, we present Siam U Net Attn model a multi class deep learning model with an attention mechanism to assess damage levels of buildings given a pair of satellite images depicting a scene before and after a disaster. We evaluate the proposed method on xView2, a large scale building damage assessment dataset, and demonstrate that the proposed approach achieves accurate damage scale classification and building segmentation results simultaneously.


A note on long powers of paths in tournaments. A square of a path on k vertices is a directed path x 1 ldots x k , where x i is directed to x i 2 , for every i in 1, ldots, k 1 . Recently, Yuster showed that any tournament on n vertices contains a square of a path of length at least n 0.295 . In this short note, we improve this bound. More precisely, we show that for every varepsilon 0 , there exists c varepsilon 0 such that any tournament on n vertices contains a square of a path on at least c varepsilon n 1 varepsilon vertices.


The Scalability, Efficiency and Complexity of Universities and Colleges A New Lens for Assessing the Higher Educational System. The growing need for affordable and accessible higher education is a major global challenge for the 21st century. Consequently, there is a need to develop a deeper understanding of the functionality and taxonomy of universities and colleges and, in particular, how their various characteristics change with size. Scaling has been a powerful tool for revealing systematic regularities in systems across a range of topics from physics and biology to cities, and for understanding the underlying principles of their organization and growth. Here, we apply this framework to institutions of higher learning in the United States and show that, like organisms, ecosystems and cities, they scale in a surprisingly systematic fashion following simple power law behavior. We analyze the entire spectrum encompassing 5,802 institutions ranging from large research universities to small professional schools, organized in seven commonly used sectors, which reveal distinct regimes of institutional scaling behavior. Metrics include variation in expenditures, revenues, graduation rates and estimated economic added value, expressed as functions of total enrollment, our fundamental measure of size. Our results quantify how each regime of institution leverages specific economies of scale to address distinct priorities. Taken together, the scaling of features within a sector and shifts in scaling across sectors implies that there are generic mechanisms and constraints shared by all sectors which lead to tradeoffs between their different societal functions and roles. We particularly highlight the strong complementarity between public and private research universities, and community and state colleges, four sectors that display superlinear returns to scale.


Exact solution of stochastic gene expression models with bursting, cell cycle and replication dynamics. The bulk of stochastic gene expression models in the literature do not have an explicit description of the age of a cell within a generation and hence they cannot capture events such as cell division and DNA replication. Instead, many models incorporate cell cycle implicitly by assuming that dilution due to cell division can be described by an effective decay reaction with first order kinetics. If it is further assumed that protein production occurs in bursts then the stationary protein distribution is a negative binomial. Here we seek to understand how accurate these implicit models are when compared with more detailed models of stochastic gene expression. We derive the exact stationary solution of the chemical master equation describing bursty protein dynamics, binomial partitioning at mitosis, age dependent transcription dynamics including replication, and random interdivision times sampled from Erlang or more general distributions the solution is different for single lineage and population snapshot settings. We show that protein distributions are well approximated by the solution of implicit models (a negative binomial) when the mean number of mRNAs produced per cycle is low and the cell cycle length variability is large. When these conditions are not met, the distributions are either almost bimodal or else display very flat regions near the mode and cannot be described by implicit models. We also show that for genes with low transcription rates, the size of protein noise has a strong dependence on the replication time, it is almost independent of cell cycle variability for lineage measurements and increases with cell cycle variability for population snapshot measurements. In contrast for large transcription rates, the size of protein noise is independent of replication time and increases with cell cycle variability for both lineage and population measurements.


Smooth branch of travelling waves for the Gross Pitaevskii equation in mathbb R 2 for small speed. We construct a smooth branch of travelling wave solutions for the 2 dimensional Gross Pitaevskii equations for small speed. These travelling waves exhibit two vortices far away from each other. We also compute the leading order term of the derivatives with respect to the speed. We construct these solutions by an implicit function type argument.


Diamond Message Set Groupcasting From an Inner Bound for the DM Broadcast Channel to the Capacity Region of the Combination Network. Multiple groupcasting over the broadcast channel (BC) is studied. In particular, an inner bound is obtained for the K receiver discrete memoryless (DM) BC for the diamond message set which consists of four groupcast messages one desired by all receivers, one by all but two receivers, and two more desired by all but each one of those two receivers. The inner bound is based on rate splitting and superposition coding and is given in explicit form herein as a union over coding distributions of four dimensional polytopes. When specialized to the so called combination network, which is a class of three layer (two hop) broadcast networks parameterized by 2 K 1 finite and arbitrary capacity noiseless links from the source node in the first layer to as many nodes of the second layer, our top down approach from the DM BC to the combination network yields an explicit inner bound as a single polytope via the identification of a single coding distribution. This inner bound consists of inequalities which are then identified to be within the class of a plethora of (indeed, infinitely many) generalized cut set outer bounds recently obtained by Salimi et al for broadcast networks. We hence establish the capacity region of the general K user combination network for the diamond message set, and do so in explicit form. Such a result implies a certain strength of our inner bound for the DM BC in that it (a) produces a hitherto unknown capacity region when specialized to the combination network and (b) may capture many combinatorial aspects of the capacity region of the K receiver DM BC itself (for the diamond message set). Moreover, we further extend that inner bound by adding binning to it and providing that inner bound also in explicit form as a union over coding distributions of four dimensional polytopes in the message rates.


A microscopic description of acid base equilibrium. Acid base reactions are ubiquitous in nature. Understanding their mechanisms is crucial in many fields, from biochemistry to industrial catalysis. Unfortunately, experiments only give limited information without much insight into the molecular behaviour. Atomistic simulations could complement experiments and shed precious light on microscopic mechanisms. The large free energy barriers connected to proton dissociation however make the use of enhanced sampling methods mandatory. Here we perform an ab initio molecular dynamics (MD) simulation and enhance sampling with the help of methadynamics. This has been made possible by the introduction of novel descriptors or collective variables (CVs) that are based on a conceptually new outlook on acid base equilibria. We test successfully our approach on three different aqueous solutions of acetic acid, ammonia, and bicarbonate. These are representative of acid, basic, and amphoteric behaviour.


A Graph Based Modeling Abstraction for Optimization Concepts and Implementation in Plasmo.jl. We present a general graph based modeling abstraction for optimization that we call an OptiGraph. Under this abstraction, any optimization problem is treated as a hierarchical hypergraph in which nodes represent optimization subproblems and edges represent connectivity between such subproblems. The abstraction enables the modular construction of highly complex models in an intuitive manner, facilitates the use of graph analysis tools (to perform partitioning, aggregation, and visualization tasks), and facilitates communication of structures to decomposition algorithms. We provide an open source implementation of the abstraction in the Julia based package Plasmo.jl. We provide tutorial examples and large application case studies to illustrate the capabilities.


On the maximum modulus of integers in Kummer extensions. We study the extension of a result of Loxton (1972) on representation of algebraic integers as sums of roots of unity to Kummer extensions.


Nb3Sn multicell cavity coating system at JLAB. SRF niobium cavities are the building blocks of modern accelerators for scientific applications. Lower surface resistance, higher fields, and high operating temperatures advance the reach of the future accelerators for scientific discovery as well as potentially enabling cost effective industrial solutions. We describe the design and performance of an Nb 3 Sn coating system that converts the inner surface of niobium cavities to Nb 3 Sn film. Niobium surface, heated by radiation from the niobium retort, is exposed to Sn and SnCl 2 vapor during the heat cycle, which results in about 2 mu m Nb 3 Sn film on the niobium surface. Film composition and structure as well as RF properties with 1 cell R D cavities and 5 cell practical accelerator cavities are presented.


Between reduced powers and ultrapowers, II. We prove that, consistently with ZFC, no ultraproduct of countably infinite (or separable metric, non compact) structures is isomorphic to a reduced product of countable (or separable metric) structures associated to the Fr 'echet filter. Since such structures are countably saturated, the Continuum Hypothesis implies that they are isomorphic when elementarily equivalent.


Street to Cloud Improving Flood Maps With Crowdsourcing and Semantic Segmentation. To address the mounting destruction caused by floods in climate vulnerable regions, we propose Street to Cloud, a machine learning pipeline for incorporating crowdsourced ground truth data into the segmentation of satellite imagery of floods. We propose this approach as a solution to the labor intensive task of generating high quality, hand labeled training data, and demonstrate successes and failures of different plausible crowdsourcing approaches in our model. Street to Cloud leverages community reporting and machine learning to generate novel, near real time insights into the extent of floods to be used for emergency response.


Can preferential concentration of finite size particles in plane Couette turbulence be reproduced with the aid of equilibrium solutions?. This work employs for the first time invariant solutions of the Navier Stokes equations to study the interaction between finite size particles and near wall coherent structures. We consider horizontal plane Couette flow and focus on Nagata's upper branch equilibrium solution (Nagata, 1990) at low Reynolds numbers where this solution is linearly stable. When adding a single heavy particle with a diameter equivalent to 2.5 wall units, we observe that the solution remains stable and is essentially unchanged away from the particle. This result demonstrates that it is technically feasible to utilize exact coherent structures in conjunction with particle resolved DNS. While translating in the streamwise direction, the particle migrates laterally under the action of the quasi streamwise vortices until it reaches the region occupied by the low speed streak, where it attains a periodic state of motion. As a result of the ensuing preferential particle location, the time average streamwise particle velocity differs from the plane average fluid phase velocity at the same wall distance as the particle center, as previously observed in experiments and in numerical data for fully turbulent wall bounded flows. Additional constrained simulations where the particle is maintained at a fixed spanwise position while freely translating in the other two directions reveal the existence of two equilibria located in the low speed and in the high speed streak, respectively, the former being an unstable point. A parametric study with different particle to fluid density ratios is conducted which shows how inertia affects the spanwise fluctuations of the periodic particle motion. Finally, we discuss a number of potential future investigations of solid particle dynamics which can be conducted with the aid of invariant solutions (exact coherent structures) of the Navier Stokes equations.


Highly Confined In plane Exciton Polaritons in Monolayer Semiconductors. 2D materials support unique excitations of quasi particles that consist of a material excitation and photons called polaritons. Especially interesting are in plane propagating polaritons which can be confined to a single monolayer and carry large momentum. In this work, we report the existence of a new type of in plane propagating polariton, supported on monolayer transition metal dicalcogonide (TMD) in the visible spectrum, which has not yet been observed. This 2D in plane exciton polariton (2DEP) is described by the coupling of an electromagnetic light field with the collective oscillations of the excitons supported by monolayer TMDs. We expose the specific experimental conditions required for the excitation of the 2DEP and show that these can be created if the TMD is encapsulated with hexagonal boron nitride (hBN) and cooled to cryogenic temperatures. In addition, we compare the properties of the 2DEPs with those of surface plasmons polaritons (SPPs) at the same spectral range, and find that the 2DEP exhibit over two orders of magnitude larger wavelength confinement. Finally, we propose two configurations for the possible experimental observation of 2DEPs.


Minimal Session Types (Extended Version). Session types are a type based approach to the verification of message passing programs. They have been much studied as type systems for the pi calculus and for languages such as Java. A session type specifies what and when should be exchanged through a channel. Central to session typed languages are constructs in types and processes that specify sequencing in protocols. Here we study minimal session types, session types without sequencing. This is arguably the simplest form of session types. By relying on a core process calculus with sessions and higher order concurrency (abstraction passing), we prove that every process typable with usual (non minimal) session types can be compiled down into a process typed with minimal session types. This means that having sequencing constructs in both processes and session types is redundant only sequentiality in processes is indispensable, as it can precisely codify sequentiality in types. Our developments draw inspiration from work by Parrow on behavior preserving decompositions of untyped processes. By casting Parrow's results in the realm of typed processes, our results reveal a conceptually simple formulation of session types and a principled avenue to the integration of session types into languages without sequencing in types.


Elliptic theory in domains with boundaries of mixed dimension. Take an open domain Omega subset mathbb R n whose boundary may be composed of pieces of different dimensions. For instance, Omega can be a ball on mathbb R 3 , minus one of its diameters D , or Omega subset mathbb R 3 could be a so called saw tooth domain, with a boundary consisting of pieces of 1 dimensional curves intercepted by 2 dimensional spheres. Under appropriate geometric assumptions, such as the existence of doubling measures on Omega and partial Omega with appropriate size conditions, we construct a class of degenerate elliptic operators L adapted to the geometry, and establish key estimates of elliptic theory associated to those operators. This includes boundary Poincar 'e and Harnack inequalities, maximum principle, and H "older continuity of solutions at the boundary. We introduce Hilbert spaces naturally associated to the geometry, construct appropriate trace and extension operators, and use them to define weak solutions to Lu 0 . Then we prove De Giorgi Nash Moser estimates inside Omega and on the boundary, solve the Dirichlet problem and thus construct an elliptic measure omega L associated to L . At last, we introduce Green functions, and use them to prove a comparison principle. Since our theory emphasizes measures, rather than the geometry per se, the results are new even in the classical setting of a half plane mathbb R 2 when the boundary partial mathbb R 2 mathbb R is equipped with a doubling measure mu singular with respect to the Lebesgue measure on mathbb R . Finally, the present paper provides a generalization of the celebrated Caffarelli Sylvestre extension operator from its classical setting of mathbb R n 1 to general open sets, and hence, an extension of the concept of fractional Laplacian to Ahlfors regular boundaries and beyond.


Mabuchi geometry of big cohomology classes with prescribed singularities. Let X be a compact K "ahler unibranch complex analytic space of pure dimension. Fix a big class alpha with smooth representative theta and a model potential phi with positive mass. We define and the study non pluripolar products of quasi plurisubharmonic functions on X . We study the spaces mathcal E p(X, theta phi ) of finite energy K "ahler potentials with prescribed singularities for each p geq 1 . We define a metric d p and show that ( mathcal E p(X, theta phi ),d p) is a complete metric space. This construction generalizes the usual d p metric defined for an ample class.


Domino tilings and flips in dimensions 4 and higher. In this paper we consider domino tilings of bounded regions in dimension n geq 4 . We define the twist of such a tiling, an elements of mathbb Z (2) , and prove it is invariant under flips, a simple local move in the space of tilings. We investigate which regions D are regular, i.e. whenever two tilings t 0 and t 1 of D times 0,N have the same twist then t 0 and t 1 can be joined by a sequence of flips provided some extra vertical space is allowed. We prove that all boxes are regular except D 0,2 3 . Furthermore, given a regular region D , we show that there exists a value M (depending only on D ) such that if t 0 and t 1 are tilings of equal twist of D times 0,N then the corresponding tilings can be joined by a finite sequence of flips in D times 0,N M . As a corollary we deduce that, for regular D and large N , the set of tilings of D times 0,N has two twin giant components under flips, one for each value of the twist.


Uniform bounds for norms of theta series and arithmetic applications. We prove uniform bounds for the Petersson norm of the cuspidal part of the theta series. This gives an improved asymptotic formula for the number of representations by a quadratic form. As an application, we show that every integer n neq 0,4,7 ,( operatorname mod 8) is represented as n x 1 2 x 2 2 x 3 3 for integers x 1,x 2,x 3 such that the product x 1x 2x 3 has at most 72 prime divisors.


A Liouville theorem on asymptotically Calabi spaces. In this paper, we will study harmonic functions on the complete and incomplete spaces with nonnegative Ricci curvature which exhibit inhomogeneous collapsing behaviors at infinity. The main result states that any nonconstant harmonic function on such spaces yields a definite exponential growth rate which depends explicitly on the geometric data at infinity.


Smart IoT Cameras for Crowd Analysis based on augmentation for automatic pedestrian detection, simulation and annotation. Smart video sensors for applications related to surveillance and security are IOT based as they use Internet for various purposes. Such applications include crowd behaviour monitoring and advanced decision support systems operating and transmitting information over internet. The analysis of crowd and pedestrian behaviour is an important task for smart IoT cameras and in particular video processing. In order to provide related behavioural models, simulation and tracking approaches have been considered in the literature. In both cases ground truth is essential to train deep models and provide a meaningful quantitative evaluation. We propose a framework for crowd simulation and automatic data generation and annotation that supports multiple cameras and multiple targets. The proposed approach is based on synthetically generated human agents, augmented frames and compositing techniques combined with path finding and planning methods. A number of popular crowd and pedestrian data sets were used to validate the model, and scenarios related to annotation and simulation were considered.


A PLS blockchain for IoT applications protocols and architecture. This paper proposes an architecture and a protocol suite for a permissioned blockchain for a local IoT network. The architecture is based on a sealed Sequencer and a Fog Server running (post quantum) Guy Fawkes protocols. The blocks of the blockchain are stored in networked Content Addressable Storage alongside any user data and validity proofs. We maintain that a typical IoT device can, despite its resource limitations, use our blockchain protocols directly, without a trusted intermediary. This includes posting and monitoring transactions as well as off chain (post quantum) emergency communications without an explicit public key. Keywords blockchain, Guy Fawkes protocol, post quantum, HORS OTS, LoRa, concurrent transmission


An Iterative Least Squares Method for Proton CT Image Reconstruction. Clinically useful proton Computed Tomography images will rely on algorithms to find the three dimensional proton stopping power distribution that optimally fits the measured proton data. We present a least squares iterative method with many features to put proton imaging into a more quantitative framework. These include the definition of a unique solution that optimally fits the protons, the definition of an iteration vector that takes into account proton measurement uncertainties, the definition of an optimal step size for each iteration individually, the ability to simultaneously optimize the step sizes of many iterations, the ability to divide the proton data into arbitrary numbers of blocks for parallel processing and use of graphical processing units, and the definition of stopping criteria to determine when to stop iterating. We find that it is possible, for any object being imaged, to provide assurance that the image is quantifiably close to an optimal solution, and the optimization of step sizes reduces the total number of iterations required for convergence. We demonstrate the use of these algorithms on real data.