A pipeline structure for the block QR update in digital signal processing

[EN] There exist problems in the field of digital signal processing, such as filtering of acoustic signals that require processing a large amount of data in real time. The beamforming algorithm, for instance, is a process that can be modeled by a rectangular matrix built on the input signals of an acoustic system and, thus, changes in real time. To obtain the output signals, it is required to compute its QR factorization. In this paper, we propose to organize the concurrent computational resources of a given multicore computer in a pipeline structure to perform this factorization as fast as possible. The pipeline has been implemented using both the application programming interface OpenMP and GrPPI, a library interface to design parallel applications based on parallel patterns. We tackle not only the performance challenge but also the programmability of our idea using parallel programming frameworks.This work was supported by the Spanish Ministry of Economy and Competitiveness under MINECO and FEDER projects TIN2014-53495-R and TEC2015-67387-C4-1-R.Dolz, MF.; Alventosa, FJ.; Alonso-Jordá, P.; Vidal Maciá, AM. (2019). A pipeline structure for the block QR update in digital signal processing. The Journal of Supercomputing. 75(3):1470-1482. https://doi.org/10.1007/s11227-018-2666-1S14701482753Huang Y, Benesty J, Chen J (2006) Acoustic MIMO signal processing (signals and communication technology). Springer, BerlinRamiro C, Vidal AM, González A (2015) MIMOPack: a high performance computing library for MIMO communication systems. J Supercomput 71:751–760Alventosa FJ, Alonso P, Piñero G, Vidal AM (2016) Implementation of the Beamformer algorithm for the NVIDIA Jetson. In: Actas de la Conferencia, Granada, Spain, pp 201–211. ISBN 978-3-319-49955-0Alventosa FJ, Alonso P, Vidal AM, Piñero G, Quintana-Ortí ES (2018) Fast block QR update in digital signal processing. J Supercomput. https://doi.org/10.1007/s11227-018-2298-5del Rio D, Dolz MF, Fernández J, García JD (2017) A generic parallel pattern interface for stream and data processing. Concurr Comput Pract Exp 29(24):e4175Benesty J, Chen J, Huang Y, Dmochowski J (2007) On microphone-array Beamforming from a MIMO acoustic signal processing perspective. IEEE Trans Audio Speech Lang Process 15(3):1053–1065Lorente J, Piñero G, Vidal AM, Belloch JA, González A (2011) Parallel implementations of Beamforming design and filtering for microphone array applications. In: 19th European Signal Processing Conference (EUSIPCO), Barcelona, Spain, pp 501–505Belloch JA, Ferrer M, González A, Martínez-Zaldívar FJ, Vidal AM (2013) Headphone-based virtual spatialization of sound with a GPU accelerator. J Audio Eng Soc 61:546–561Belloch JA, González A, Martínez-Zaldívar FJ, Vidal AM (2011) Real-time massive convolution for audio applications on GPU. J Supercomput 58(3):449–457Golub GH, Van Loan CF (2013) Matrix computations. Johns Hopkins studies in the mathematical sciences. Johns Hopkins University Press, BaltimoreGunter BC, van de Geijn RA (2005) Parallel out-of-core computation and updating the QR factorization. ACM Trans Math Softw 31(1):60–78Buttari A, Langou J, Kurzak J, Dongarra J (2009) A class of parallel tiled linear algebra algorithms for multicore architectures. Parallel Comput 35(1):38–53Dolz MF, Alventosa FJ, Alonso-Jordá P, Vidal AM (2018) A pipeline for the QR update in digital signal processing. In: Proceedings of the 18th International Conference on Computational and Mathematical Methods in Science and Engineering (CMMSE 2018), Rota, Cádiz, Spain, pp 1–5Quintana-Ortí G, Quintana-Ortí ES, Van De Geijn RA, Van Zee FG, Chan E (2009) Programming matrix algorithms-by-blocks for thread-level parallelism. ACM Trans Math Softw 36(3):14:1–14:2

Comparative Skull Morphology of Uropeltid Snakes (Alethinophidia: Uropeltidae) with Special Reference to Disarticulated Elements and Variation

Uropeltids form a diverse clade of highly derived, fossorial snakes that, because of their phylogenetic position among other alethinophidian lineages, may play a key role in understanding the early evolution of cranial morphology in snakes. We include detailed osteological descriptions of crania and mandibles for eight uropeltid species from three nominal genera (Uropeltis, Rhinophis, and Brachyophidium) and emphasize disarticulated elements and the impact of intraspecific variation on previously proposed morphological characters used for phylogenetic analysis. Preliminary analysis of phylogenetic relationships strongly supports a clade composed exclusively of species of Plectrurus, Uropeltis, and Rhinophis. However, monophyly of each of those genera and Melanophidium is not upheld. There is moderate support that Sri Lankan species (e.g., Rhinophis and Uropeltis melanogaster) are monophyletic with respect to Indian uropeltids. Previously proposed characters that are phylogenetically informative include the shape of the nasals, length of the occipital condyle, level of development of the posteroventral process of the dentary, and participation of the parietal in the optic foramen. Additionally, thirty new features that may be systematically informative are identified and described, but were not verified for their utility. Such verification must await availability of additional disarticulated cranial material from a larger sample of taxa. All characters require further testing through increased focus on sources and patterns of intraspecific variation, inclusion of broader taxonomic samples in comparative studies, and exploration of skeletal development, sexual dimorphism, and biogeographic patterns. Additionally, trends in the relative enlargement of the sensory capsules, reduction in cranial ossification and dentition, fusion of elements, and the appearance of novel morphological conditions, such as the structure and location of the suspensorium, may be related to fossoriality and miniaturization in some uropeltid taxa, and may complicate analysis of relationships within Uropeltidae and among alethinophidian snakes

Comparative cytogenetics of three species of Dichotomius (Coleoptera, Scarabaeidae)

Meiotic and mitotic chromosomes of Dichotomius nisus, D. semisquamosus and D. sericeus were analyzed after conventional staining, C-banding and silver nitrate staining. In addition, Dichotomius nisus and D. semisquamosus chromosomes were also analyzed after fluorescent in situ hybridization (FISH) with an rDNA probe. The species analyzed had an asymmetrical karyotype with 2n = 18 and meta-submetacentric chromosomes. The sex determination mechanism was of the Xyp type in D. nisus and D. semisquamosus and of the Xy r type in D. sericeus. C-banding revealed the presence of pericentromeric blocks of constitutive heterochromatin (CH) in all the chromosomes of the three species. After silver staining, the nucleolar organizer regions (NORs) were located in autosomes of D. semisquamosus and D. sericeus and in the sexual bivalent of D. nisus. FISH with an rDNA probe confirmed NORs location in D. semisquamosus and in D. nisus. Our results suggest that chromosome inversions and fusions occurred during the evolution of the group

Reduced Bone Mass and Muscle Strength in Male 5α-Reductase Type 1 Inactivated Mice

Androgens are important regulators of bone mass but the relative importance of testosterone (T) versus dihydrotestosterone (DHT) for the activation of the androgen receptor (AR) in bone is unknown. 5α-reductase is responsible for the irreversible conversion of T to the more potent AR activator DHT. There are two well established isoenzymes of 5α-reductase (type 1 and type 2), encoded by separate genes (Srd5a1 and Srd5a2). 5α-reductase type 2 is predominantly expressed in male reproductive tissues whereas 5α-reductase type 1 is highly expressed in liver and moderately expressed in several other tissues including bone. The aim of the present study was to investigate the role of 5α-reductase type 1 for bone mass using Srd5a1−/− mice. Four-month-old male Srd5a1−/− mice had reduced trabecular bone mineral density (−36%, p<0.05) and cortical bone mineral content (−15%, p<0.05) but unchanged serum androgen levels compared with wild type (WT) mice. The cortical bone dimensions were reduced in the male Srd5a1−/− mice as a result of a reduced cortical periosteal circumference compared with WT mice. T treatment increased the cortical periosteal circumference (p<0.05) in orchidectomized WT mice but not in orchidectomized Srd5a1−/− mice. Male Srd5a1−/− mice demonstrated a reduced forelimb muscle grip strength compared with WT mice (p<0.05). Female Srd5a1−/− mice had slightly increased cortical bone mass associated with elevated circulating levels of androgens. In conclusion, 5α-reductase type 1 inactivated male mice have reduced bone mass and forelimb muscle grip strength and we propose that these effects are due to lack of 5α-reductase type 1 expression in bone and muscle. In contrast, the increased cortical bone mass in female Srd5a1−/− mice, is an indirect effect mediated by elevated circulating androgen levels

A search for the decay modes B+/- to h+/- tau l

We present a search for the lepton flavor violating decay modes B+/- to h+/- tau l (h= K,pi; l= e,mu) using the BaBar data sample, which corresponds to 472 million BBbar pairs. The search uses events where one B meson is fully reconstructed in one of several hadronic final states. Using the momenta of the reconstructed B, h, and l candidates, we are able to fully determine the tau four-momentum. The resulting tau candidate mass is our main discriminant against combinatorial background. We see no evidence for B+/- to h+/- tau l decays and set a 90% confidence level upper limit on each branching fraction at the level of a few times 10^-5.Comment: 15 pages, 7 figures, submitted to Phys. Rev.

Search For Heavy Pointlike Dirac Monopoles

We have searched for central production of a pair of photons with high transverse energies in $p\bar p$ collisions at $\sqrt{s} = 1.8$ TeV using $70 pb^{-1}$ of data collected with the D\O detector at the Fermilab Tevatron in 1994--1996. If they exist, virtual heavy pointlike Dirac monopoles could rescatter pairs of nearly real photons into this final state via a box diagram. We observe no excess of events above background, and set lower 95% C.L. limits of $610, 870, or 1580 GeV/c^2$ on the mass of a spin 0, 1/2, or 1 Dirac monopole.Comment: 12 pages, 4 figure

The Dijet Mass Spectrum and a Search for Quark Compositeness in bar{p}p Collisions at sqrt{s} = 1.8 TeV

Using the DZero detector at the 1.8 TeV pbarp Fermilab Tevatron collider, we have measured the inclusive dijet mass spectrum in the central pseudorapidity region |eta_jet| < 1.0 for dijet masses greater than 200 Gev/c^2. We have also measured the ratio of spectra sigma(|eta_jet| < 0.5)/sigma(0.5 < |eta_jet| < 1.0). The order alpha_s^3 QCD predictions are in good agreement with the data and we rule out models of quark compositeness with a contact interaction scale < 2.4 TeV at the 95% confidence level.Comment: 11 pages, 4 figures, 2 tables, submitted to Phys. Rev. Let

Search for High Mass Photon Pairs in p-pbar --> gamma-gamma-jet-jet Events at sqrt(s)=1.8 TeV

A search has been carried out for events in the channel p-barp --> gamma gamma jet jet. Such a signature can characterize the production of a non-standard Higgs boson together with a W or Z boson. We refer to this non-standard Higgs, having standard model couplings to vector bosons but no coupling to fermions, as a "bosonic Higgs." With the requirement of two high transverse energy photons and two jets, the diphoton mass (m(gamma gamma)) distribution is consistent with expected background. A 90(95)% C.L. upper limit on the cross section as a function of mass is calculated, ranging from 0.60(0.80) pb for m(gamma gamma) = 65 GeV/c^2 to 0.26(0.34) pb for m(gamma gamma) = 150 GeV/c^2, corresponding to a 95% C.L. lower limit on the mass of a bosonic Higgs of 78.5 GeV/c^2.Comment: 9 pages, 3 figures. Replacement has new H->gamma gamma branching ratios and corresponding new mass limit

Evidence for an excess of B -> D(*) Tau Nu decays

Based on the full BaBar data sample, we report improved measurements of the ratios R(D(*)) = B(B -> D(*) Tau Nu)/B(B -> D(*) l Nu), where l is either e or mu. These ratios are sensitive to new physics contributions in the form of a charged Higgs boson. We measure R(D) = 0.440 +- 0.058 +- 0.042 and R(D*) = 0.332 +- 0.024 +- 0.018, which exceed the Standard Model expectations by 2.0 sigma and 2.7 sigma, respectively. Taken together, our results disagree with these expectations at the 3.4 sigma level. This excess cannot be explained by a charged Higgs boson in the type II two-Higgs-doublet model. We also report the observation of the decay B -> D Tau Nu, with a significance of 6.8 sigma.Comment: Expanded section on systematics, text corrections, improved the format of Figure 2 and included the effect of the change of the Tau polarization due to the charged Higg

Search for the decay modes D^0 → e^+e^-, D^0 → μ^+μ^-, and D^0 → e^±μ∓

We present searches for the rare decay modes D^0→e^+e^-, D^0→μ^+μ^-, and D^0→e^±μ^∓ in continuum e^+e^-→cc events recorded by the BABAR detector in a data sample that corresponds to an integrated luminosity of 468  fb^(-1). These decays are highly Glashow–Iliopoulos–Maiani suppressed but may be enhanced in several extensions of the standard model. Our observed event yields are consistent with the expected backgrounds. An excess is seen in the D^0→μ^+μ^- channel, although the observed yield is consistent with an upward background fluctuation at the 5% level. Using the Feldman–Cousins method, we set the following 90% confidence level intervals on the branching fractions: B(D^0→e^+e^-)<1.7×10^(-7), B(D^0→μ^+μ^-) within [0.6,8.1]×10^(-7), and B(D^0→e^±μ^∓)<3.3×10^(-7)