Model Mining and Efficient Verification of Software Product Lines

Software product line modeling aims at capturing a set of software products in an economic yet meaningful way. We introduce a class of variability models that capture the sharing between the software artifacts forming the products of a software product line (SPL) in a hierarchical fashion, in terms of commonalities and orthogonalities. Such models are useful when analyzing and verifying all products of an SPL, since they provide a scheme for divide-and-conquer-style decomposition of the analysis or verification problem at hand. We define an abstract class of SPLs for which variability models can be constructed that are optimal w.r.t. the chosen representation of sharing. We show how the constructed models can be fed into a previously developed algorithmic technique for compositional verification of control-flow temporal safety properties, so that the properties to be verified are iteratively decomposed into simpler ones over orthogonal parts of the SPL, and are not re-verified over the shared parts. We provide tool support for our technique, and evaluate our tool on a small but realistic SPL of cash desks.

Maximal structuring of acyclic process models

This article addresses the transformation of a process model with an arbitrary topology into an equivalent structured process model. In particular, this article studies the subclass of process models that have no equivalent well-structured representation but which, nevertheless, can be partially structured into their maximally-structured representation. The transformations are performed under a behavioral equivalence notion that preserves the observed concurrency of tasks in equivalent process models. The article gives a full characterization of the subclass of acyclic process models that have no equivalent well-structured representation, but do have an equivalent maximally-structured one, as well as proposes a complete structuring method. Together with our previous results, this article completes the solution of the process model structuring problem for the class of acyclic process models.

Experimental and computer simulation validation of ultrasound phase interference created by lateral inhomogeneity of transit time in replica bone marrow composite models

Purpose: The measurement of broadband ultrasonic attenuation (BUA) in cancellous bone for the assessment of osteoporosis follows a parabolic-type dependence with bone volume fraction; having minima values corresponding to both entire bone and entire marrow. Langton has recently proposed that the primary BUA mechanism may be significant phase interference due to variations in propagation transit time through the test sample as detected over the phase-sensitive surface of the receive ultrasound transducer. This fundamentally simple concept assumes that the propagation of ultrasound through a complex solid : liquid composite sample such as cancellous bone may be considered by an array of parallel ‘sonic rays’. The transit time of each ray is defined by the proportion of bone and marrow propagated, being a minimum (tmin) solely through bone and a maximum (tmax) solely through marrow. A Transit Time Spectrum (TTS), ranging from tmin to tmax, may be defined describing the proportion of sonic rays having a particular transit time, effectively describing lateral inhomogeneity of transit time over the surface of the receive ultrasound transducer. Phase interference may result from interaction of ‘sonic rays’ of differing transit times. The aim of this study was to test the hypothesis that there is a dependence of phase interference upon the lateral inhomogenity of transit time by comparing experimental measurements and computer simulation predictions of ultrasound propagation through a range of relatively simplistic solid:liquid models exhibiting a range of lateral inhomogeneities. Methods: A range of test models was manufactured using acrylic and water as surrogates for bone and marrow respectively. The models varied in thickness in one dimension normal to the direction of propagation, hence exhibiting a range of transit time lateral inhomogeneities, ranging from minimal (single transit time) to maximal (wedge; ultimately the limiting case where each sonic ray has a unique transit time). For the experimental component of the study, two unfocused 1 MHz ¾” broadband diameter transducers were utilized in transmission mode; ultrasound signals were recorded for each of the models. The computer simulation was performed with Matlab, where the transit time and relative amplitude of each sonic ray was calculated. The transit time for each sonic ray was defined as the sum of transit times through acrylic and water components. The relative amplitude considered the reception area for each sonic ray along with absorption in the acrylic. To replicate phase-sensitive detection, all sonic rays were summed and the output signal plotted in comparison with the experimentally derived output signal. Results: From qualtitative and quantitative comparison of the experimental and computer simulation results, there is an extremely high degree of agreement of 94.2% to 99.0% between the two approaches, supporting the concept that propagation of an ultrasound wave, for the models considered, may be approximated by a parallel sonic ray model where the transit time of each ray is defined by the proportion of ‘bone’ and ‘marrow’. Conclusions: This combined experimental and computer simulation study has successfully demonstrated that lateral inhomogeneity of transit time has significant potential for phase interference to occur if a phase-sensitive ultrasound receive transducer is implemented as in most commercial ultrasound bone analysis devices.

Search for Maximal Flavor Violating Scalars in Same-Charge Lepton Pairs in pp̅ Collisions at √s=1.96 TeV

Models of Maximal Flavor Violation (MxFV) in elementary particle physics may contain at least one new scalar SU$(2)$ doublet field $\Phi_{FV} = (\eta^0,\eta^+)$ that couples the first and third generation quarks ($q_1,q_3$) via a Lagrangian term $\mathcal{L}_{FV} = \xi_{13} \Phi_{FV} q_1 q_3$. These models have a distinctive signature of same-charge top-quark pairs and evade flavor-changing limits from meson mixing measurements. Data corresponding to 2 fb$^{-1}$ collected by the CDF II detector in $p\bar{p}$ collisions at $\sqrt{s} = 1.96$ TeV are analyzed for evidence of the MxFV signature. For a neutral scalar $\eta^0$ with $m_{\eta^0} = 200$ GeV/$c^2$ and coupling $\xi_{13}=1$, $\sim$ 11 signal events are expected over a background of $2.1 \pm 1.8$ events. Three events are observed in the data, consistent with background expectations, and limits are set on the coupling $\xi_{13}$ for $m_{\eta^0} = 180-300$ GeV/$c^2$.

Efficient learning of decomposable models with a bounded clique size

The learning of probability distributions from data is a ubiquitous problem in the fields of Statistics and Artificial Intelligence. During the last decades several learning algorithms have been proposed to learn probability distributions based on decomposable models due to their advantageous theoretical properties. Some of these algorithms can be used to search for a maximum likelihood decomposable model with a given maximum clique size, k, which controls the complexity of the model. Unfortunately, the problem of learning a maximum likelihood decomposable model given a maximum clique size is NP-hard for k > 2. In this work, we propose a family of algorithms which approximates this problem with a computational complexity of O(k · n^2 log n) in the worst case, where n is the number of implied random variables. The structures of the decomposable models that solve the maximum likelihood problem are called maximal k-order decomposable graphs. Our proposals, called fractal trees, construct a sequence of maximal i-order decomposable graphs, for i = 2, ..., k, in k − 1 steps. At each step, the algorithms follow a divide-and-conquer strategy based on the particular features of this type of structures. Additionally, we propose a prune-and-graft procedure which transforms a maximal k-order decomposable graph into another one, increasing its likelihood. We have implemented two particular fractal tree algorithms called parallel fractal tree and sequential fractal tree. These algorithms can be considered a natural extension of Chow and Liu’s algorithm, from k = 2 to arbitrary values of k. Both algorithms have been compared against other efficient approaches in artificial and real domains, and they have shown a competitive behavior to deal with the maximum likelihood problem. Due to their low computational complexity they are especially recommended to deal with high dimensional domains.

Structural dynamic eutrophication models

This article discusses problems of modelling the seasonal succession of algal species in lakes and reservoirs, and the adaptive selection of certain groups of algae in response to changes in the inputs and relative concentrations of nutrients and other environmental variables. A new generation of quantitative models is being developed which attempts to translate some important biological properties of species (survival, variation, inheritance, reproductive rates and population growth) into predictions about the survival of the fittest, where ”fitness” is measured or estimated in thermodynamic terms. The concept of ”exergy” and its calculation is explored to examine maximal exergy as a measure of fitness in ecosystems, and its use for calculating changes in species composition by means of structural dynamic models. These models accomodate short-term changes in parameters that affect the adaptive responses (species selection) of algae.

Stocking models of Chinese mitten crab (Eriocheir japonica sinensis) in Yangtze lakes

The fanning of Chinese mitten crab, a quality aquatic product in China and neighbouring Asian countries, has been developing rapidly in China since last decade. It reached a total yield of 3.4 X 10(5) tonnes in 2002. Due to the successive over-stocking year after year, many lakes in the mid-lower Yangtze Basin, the main farming area, are under deterioration, leading to a reduction of crab yield and quality, and, subsequently, a loss of fanning profits. Aiming at a normal development of crab culture and the sustainable use of lakes, an annual investigation dealing with lake environmental factors in relation to stocked crab populations was carried out at 20 farms in 4 lakes. The results show that the submersed macrophyte biomass (B-Mac) is the key factor affecting annual crab yield (CY). Using the ratio of Secchi depth to mean depth (Z(SD)/Z(M)), an easily measured parameter closely correlated to BMac, as driving variable, 10 regression models of maximal crab yields were generated (r(2) ranging 0.49-0.81). Based on the theory of MSY (Maximum Sustainable Yield), in combination with body-weight (BW) and recapture rate (RR) of adult crabs, a general optimal stocking model was eventually formulated. All models are simple and easy to operate. Comments on their applications and prospects are given in brief. (c) 2006 Elsevier B.V. All rights reserved.

Search for maximal flavor violating scalars in same-charge lepton pairs in pp̄ collisions at s=1.96TeV

Models of maximal flavor violation (MxFV) in elementary particle physics may contain at least one new scalar SU(2) doublet field ΦFV=(η0,η+) that couples the first and third generation quarks (q1, q3) via a Lagrangian term LFV=ξ13ΦFVq1q3. These models have a distinctive signature of same-charge top-quark pairs and evade flavor-changing limits from meson mixing measurements. Data corresponding to 2fb-1 collected by the Collider Dectector at Fermilab II detector in pp̄ collisions at s=1.96TeV are analyzed for evidence of the MxFV signature. For a neutral scalar η0 with mη0=200GeV/c2 and coupling ξ13=1, ∼11 signal events are expected over a background of 2.1±1.8 events. Three events are observed in the data, consistent with background expectations, and limits are set on the coupling ξ13 for mη0=180-300GeV/c2. © 2009 The American Physical Society.

Comparaison des réponses physiologiques lors d’un exercice incrémental maximal sur vélo immergé et sur terrain sec : aspects biomécaniques, cardiopulmonaires et hémodynamiques.

L’exercice en immersion dans l'eau peut générer des réponses hémodynamiques et cardiorespiratoires différentes à celles de l’exercice sur terraine sec. Cependant, aucune étude n’a comparé ces réponses sur vélo aquatique (VA) à celles sur vélo sur terrain sec (VS) à une même puissance mécanique externe (Pext). À cet égard, le premier travail de cette thèse visait, d’abord, à trouver les équivalences de Pext lors du pédalage sur VA en immersion à la poitrine par rapport au VS au laboratoire, en considérant que cela restait non déterminé à ce jour. Une équation de mécanique des fluides fut utilisée pour calculer la force déployée pour le système de pédalage (pales, leviers, pédales) et des jambes à chaque tour de pédale. Ensuite, cette force totale a été multipliée par la vitesse de pédalage pour estimer la Pext sur VA. Ayant trouvé les équivalences de Pext sur VA et VS, nous nous sommes fixés comme objectif dans la deuxième étude de comparer les réponses hémodynamiques et cardiorespiratoires lors d'un exercice maximal progressif sur VS par rapport au VA à une même Pext. Les résultats ont montré que le VO2 (p<0.0001) et la différence artério-veineuse (C(a-v)O2) (p<0.0001) étaient diminués lors de l’exercice sur VA comparativement à celui sur VS. Parmi les variables hémodynamiques, le volume d’éjection systolique (VES) (p˂0.05) et le débit cardiaque (Qc) (p˂0.05) étaient plus élevés sur VA. En plus, on nota une diminution significative de la fréquence cardiaque (FC) (p˂0.05). Étant donné qu’à une même Pext les réponses physiologiques sont différentes sur VA par rapport à celles sur VS, nous avons effectué une troisième étude pour établir la relation entre les différentes expressions de l'intensité relative de l'exercice (% du VO2max,% de la FCmax,% du VO2 de réserve (% de VO2R) et % de la FC réserve (% FCR)). Les résultats ont démontré que la relation % FCR vs % VO2R était la plus corrélée (régression linéaire) et la plus proche de la ligne d’identité. Ces résultats pourraient aider à mieux prescrire et contrôler l’intensité de l'exercice sur VA pour des sujets sains. Finalement, une dernière étude comparant la réactivation parasympathique après un exercice maximal incrémental effectué sur VA et VS en immersion au niveau de la poitrine a montré que la réactivation parasympathique à court terme était plus prédominante sur VA (i,e. t, delta 10 à delta 60 et T30, p<0.05). Cela suggérait, qu’après un exercice maximal sur VA, la réactivation parasympathique à court terme était accélérée par rapport à celle après l'effort maximal sur VS chez de jeunes sujets sains. En conclusion, nous proposons une méthode de calcul de la puissance mécanique externe sur VA en fonction de la cadence de pédalage. Nous avons démontré que pendant l’exercice sur VA les réponses hémodynamiques et cardiorespiratoires sont différentes de celles sur VS à une même Pext et nous proposons des équations pour le calcul du VO2 dans l’eau ainsi qu’une méthode pour la prescription et le contrôle de l’exercice sur VA. Finalement, la réactivation parasympathique à court terme s’est trouvée accélérée après un effort maximal incrémental sur VA comparativement à celle sur VS.

Precise calculations of the existence of multiple AMOC equilibria in coupled climate models Part I: equilibrium states

This study examines criteria for the existence of two stable states of the Atlantic Meridional Overturning Circulation (AMOC) using a combination of theory and simulations from a numerical coupled atmosphere–ocean climate model. By formulating a simple collection of state parameters and their relationships, the authors reconstruct the North Atlantic Deep Water (NADW) OFF state behavior under a varying external salt-flux forcing. This part (Part I) of the paper examines the steady-state solution, which gives insight into the mechanisms that sustain the NADW OFF state in this coupled model; Part II deals with the transient behavior predicted by the evolution equation. The nonlinear behavior of the Antarctic Intermediate Water (AAIW) reverse cell is critical to the OFF state. Higher Atlantic salinity leads both to a reduced AAIW reverse cell and to a greater vertical salinity gradient in the South Atlantic. The former tends to reduce Atlantic salt export to the Southern Ocean, while the latter tends to increases it. These competing effects produce a nonlinear response of Atlantic salinity and salt export to salt forcing, and the existence of maxima in these quantities. Thus the authors obtain a natural and accurate analytical saddle-node condition for the maximal surface salt flux for which a NADW OFF state exists. By contrast, the bistability indicator proposed by De Vries and Weber does not generally work in this model. It is applicable only when the effect of the AAIW reverse cell on the Atlantic salt budget is weak.

Cannabidivarin is anticonvulsant in mouse and rat in vitro and in seizure models

Summary Background and purpose: Phytocannabinoids in Cannabis sativa have diverse pharmacological targets extending beyond cannabinoid receptors and several exert notable anticonvulsant effects. For the first time, we investigated the anticonvulsant profile of the phytocannabinoid cannabidivarin (CBDV) in vitro and in in vivo seizure models. Experimental approach: The effect of CBDV (1-100μM) on epileptiform local field potentials (LFPs) induced in rat hippocampal brain slices by 4-AP application or Mg2+-free conditions was assessed by in vitro multi-electrode array recordings. Additionally, the anticonvulsant profile of CBDV (50-200 mg kg-1) in vivo was investigated in four rodent seizure models: maximal electroshock (mES) and audiogenic seizures in mice, and pentylenetetrazole (PTZ) and pilocarpine-induced seizures in rat. CBDV effects in combination with commonly-used antiepileptic drugs were investigated in rat seizures. Finally, the motor side effect profile of CBDV was investigated using static beam and gripstrength assays. Key results: CDBV significantly attenuated status epilepticus-like epileptiform LFPs induced by 4-AP and Mg2+-free conditions. CBDV had significant anticonvulsant effects in mES (≥100 mg kg-1), audiogenic (≥50 mg kg-1) and PTZ-induced seizures (≥100 mg kg-1). CBDV alone had no effect against pilocarpine-induced seizures, but significantly attenuated these seizures when administered with valproate or phenobarbital at 200 mg kg-1 CBDV. CBDV had no effect on motor function. Conclusions and Implications: These results indicate that CBDV is an effective anticonvulsant across a broad range of seizure models, does not significantly affect normal motor function and therefore merits further investigation in chronic epilepsy models to justify human trials.

Among- and within-individual correlations between basal and maximal metabolic rates in birds

The aerobic capacity model proposes that endothermy is a by-product of selection favouring high maximal metabolic rates (MMR) and its mechanistic coupling with basal metabolic rate (BMR). Attempts to validate this model in birds are equivocal and restricted to phenotypic correlations (rP), thus failing to distinguish among- and within-individual correlations (rind and re). We examined 300 paired measurements of BMR and MMR from 60 house sparrows before and after two levels of experimental manipulation - testosterone implants and immune challenge. Overall, repeatability was significant in both BMR (R=0.25±0.06) and MMR (R=0.52±0.06). Only the testosterone treatment altered the rP between BMR and MMR, which resulted from contrasting effects on rind and re. While rind was high and significant (0.62±0.22) in sham-implanted birds, re was negative and marginally non-significant (-0.15±0.09) in testosterone-treated birds. Thus, the expected mechanistic link between BMR and MMR was apparent, but only in birds with low testosterone levels.

The highest intensity and the shortest duration permitting attainment of maximal oxygen uptake during cycling: effects of different methods and aerobic fitness level

The aims of this study were: (1) to verify the validity of previous proposed models to estimate the lowest exercise duration (T (LOW)) and the highest intensity (I (HIGH)) at which VO(2)max is reached (2) to test the hypothesis that parameters involved in these models, and hence the validity of these models are affected by aerobic training status. Thirteen cyclists (EC), eleven runners (ER) and ten untrained (U) subjects performed several cycle-ergometer exercise tests to fatigue in order to determine and estimate T (LOW) (ET (LOW)) and I (HIGH) (EI (HIGH)). The relationship between the time to achieved VO(2)max and time to exhaustion (T (lim)) was used to estimate ET (LOW). EI (HIGH) was estimated using the critical power model. I (HIGH) was assumed as the highest intensity at which VO2 was equal or higher than the average of VO(2)max values minus one typical error. T (LOW) was considered T (lim) associated with I (HIGH). No differences were found in T (LOW) between ER (170 +/- 31 s) and U (209 +/- 29 s), however, both showed higher values than EC (117 +/- 29 s). I (HIGH) was similar between U (269 +/- 73 W) and ER (319 +/- 50 W), and both were lower than EC (451 +/- 33 W). EI (HIGH) was similar and significantly correlated with I-HIGH only in U (r = 0.87) and ER (r = 0.62). ET (LOW) and T (LOW) were different only for U and not significantly correlated in all groups. These data suggest that the aerobic training status affects the validity of the proposed models for estimating I (HIGH).

Super-WZNW with reductions to supersymmetric and fermionic integrable models

A systematic construction for an action describing a class of supersymmetric integrable models as well as for pure fermionic theories is discussed in terms of the gauged WZNW model associated to half integer graded affine Kac-Moody algebras. Explicit examples of the N = 1. 2 super-sinh(sine)-Gordon models are discussed in detail. Pure fermionic theories arises for cosets sl(p, 1)/sl(p) circle times u(1) when a maximal kernel condition is fulfilled. The integrability condition for such models is discussed and it is shown that the simplest example when p = 2 (cads to the constrained Bukhvostov-Lipatov, Thirring, scalar massive and pseudo-scalar massless Gross-Neveu models. (C) 2009 Published by Elsevier B.V.

Skin and pulmonary models using coated bentonite particles for the study of the inflammation evoked by Paracoccidioides brasiliensis antigens in previously immunized mice

Bentonite particles coated with polysaccharide antigen or crude soluble antigen of Paracoccidioides brasiliensis were injected intradermally or intravenously in mice. In control animals that were not pre-immunized with P. brasiliensis antigens, coated and uncoated bentonite caused minimal and nonspecific inflammation around the cutaneous injection site or around the bentonite thrombi in small lung vessels after intravenous injection. However, in mice previously immunized with P. brasiliensis antigens, the coated bentonite particles boosted the humoral and cellular immune responses to P. brasiliensis and evoked intense inflammatory reactions. Twelve days after intradermal injection, the inflammatory reaction around the bentonite was rich in neutrophils, macrophages, lymphocytes and plasma cells associated with young granulation tissue. In intravenously injected mice, the pulmonary inflammation was maximal at day 2, and was characterized by a florid neutrophilic and macrophagic cellular infiltration around bentonite thrombi; in some foci, there was incipient organization to mature granuloma. However, in both models, there was no formation of epithelioid granulomata, demonstrating that in paracoccidioidomycosis cellular immunity alone, without the presence of intact micro-organisms, may not be enough for the development of this type of granuloma.