High order structure preserving explicit methods for solving linear-quadratic optimal control problems

[EN] We consider the numerical integration of linear-quadratic optimal control problems. This problem requires the solution of a boundary value problem: a non-autonomous matrix Riccati differential equation (RDE) with final conditions coupled with the state vector equation with initial conditions. The RDE has positive definite matrix solution and to numerically preserve this qualitative property we propose first to integrate this equation backward in time with a sufficiently accurate scheme. Then, this problem turns into an initial value problem, and we analyse splitting and Magnus integrators for the forward time integration which preserve the positive definite matrix solutions for the RDE. Duplicating the time as two new coordinates and using appropriate splitting methods, high order methods preserving the desired property can be obtained. The schemes make sequential computations and do not require the storrage of intermediate results, so the storage requirements are minimal. The proposed methods are also adapted for solving linear-quadratic N-player differential games. The performance of the splitting methods can be considerably improved if the system is a perturbation of an exactly solvable problem and the system is properly split. Some numerical examples illustrate the performance of the proposed methods.The author wishes to thank the University of California San Diego for its hospitality where part of this work was done. He also acknowledges the support of the Ministerio de Ciencia e Innovacion (Spain) under the coordinated project MTM2010-18246-C03. The author also acknowledges the suggestions by the referees to improve the presentation of this work.Blanes Zamora, S. (2015). High order structure preserving explicit methods for solving linear-quadratic optimal control problems. Numerical Algorithms. 69:271-290. https://doi.org/10.1007/s11075-014-9894-0S27129069Abou-Kandil, H., Freiling, G., Ionescy, V., Jank, G.: Matrix Riccati equations in control and systems theory. Basel, Burkhäuser Verlag (2003)Al-Mohy, A.H., Higham, N.J.: Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators. SIAM. J. Sci. Comp. 33, 488–511 (2011)Anderson, B.D.O., Moore, J.B.: Optimal control: linear quadratic methods. Dover, New York (1990)Ascher, U.M., Mattheij, R.M., Russell, R.D.: Numerical solutions of boundary value problems for ordinary differential equations. Prentice-Hall, Englewood Cliffs (1988)Bader, P., Blanes, S., Ponsoda, E.: Structure preserving integrators for solving linear quadratic optimal control problems with applications to describe the flight of a quadrotor. J. Comput. Appl. Math. 262, 223–233 (2014)Basar, T., Olsder, G.J.: Dynamic non cooperative game theory, 2nd Ed, SIAM, Philadelphhia (1999)Blanes, S., Casas, F.: On the necessity of negative coefficients for operator splitting schemes of order higher than two. Appl. Num. Math. 54, 23–37 (2005)Blanes, S., Casas, F., Farrés, A., Laskar, J., Makazaga, J., Murua, A.: New families of symplectic splitting methods for numerical integration in dynamical astronomy. Appl. Numer. Math. 68, 58–72 (2013)Blanes, S., Casas, F., Oteo, J.A., Ros, J.: The Magnus expansion and some of its applications. Phys. Rep. 470, 151–238 (2009)Blanes, S., Casas, F., Ros, J.: High order optimized geometric integrators for linear differential equations. BIT 42, 262–284 (2002)Blanes, S., Diele, F., Marangi, C., Ragni, S.: Splitting and composition methods for explicit time dependence in separable dynamical systems. J. Comput. Appl. Math. 235, 646–659 (2010)Blanes, S., Moan, P.C.: Practical symplectic partitioned Runge-Kutta and Runge-Kutta-Nystrm methods. J. Comput. Appl. Math. 142, 313–330 (2002)Blanes, S., Ponsoda, E.: Magnus integrators for solving linear-quadratic differential games. J. Comput. Appl. Math. 236, 3394–3408 (2012)Brif, C., Chakrabarti, R., Rabitz, H.: Control of quantum phenomena: past, present and future. New J. Phys. 12, 075008(68pp) (2010)Cruz, J.B., Chen, C.I.: Series Nash solution of two person non zero sum linear quadratic games. J. Optim. Theory Appl. 7, 240–257 (1971)Dieci, L., Eirola, T.: Positive definitness in the numerical solution of Riccati differential quations. Numer. Math. 67, 303–313 (1994)Engwerda, J.: LQ dynamic optimization and differential games. Wiley (2005)Hairer, E., Lubich, C., Wanner, G.: Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations (2nd edition). Springer Series in Computational Mathematics, 31. Springer-Verlag (2006)Hochbruck, M., Ostermann, A.: Exponential integrators. Acta Numerica 19, 209–286 (2010)Horn, R.A., Johnson, C.R.: Matrix Analysis. Cambridge University Press, New York (1985)Iserles, A., Munthe-Kaas, H.Z., Nørsett, S.P., Zanna, A.: Lie group methods. Acta Numerica 9, 215–365 (2000)Iserles, A., Nørsett, S.P.: On the solution of linear differential equations in Lie groups. Phil. Trans. R. Soc. Lond. A 357, 983–1019 (1999)Jódar, L., Ponsoda, E.: Non-autonomous Riccati-type matrix differential equations: existence interval, construction of continuous numerical solutions and error bounds. IMA. J. Num. Anal. 15, 61–74 (1995)Jódar, L., Ponsoda, E., Company, R.: Solutions of coupled Riccati equations arising in differential games. Control. Cybern. 24, 117–128 (1995)Kaitala, V, Pohjola, M. In: Carraro, Filar (eds.) : Sustainable international agreement on greenhouse warming. A game theory study. Control and Game Theoretic Models of the Environment, pp 67–87. Birkhauser, Boston (1995)Keller, H.B.: Numerical solution of two point boundary value problems. In: CBMS-NSF Regional Conference Series in Applied Mathematics, Vol. 24. SIAM, Philadelphia (1976)McLachlan, R.I.: Composition methods in the presence of small parameters. BIT 35, 258–268 (1995)McLachlan, R.I., Quispel, R.: Splitting Methods. Acta Numer. 11, 341–434 (2002)Moler, C.B., Van Loan, C.F.: Nineteen Dubious Ways to Compute the Exponential of a Matrix, twenty-five years later. SIAM Rev. 45, 3–49 (2003)Na, T.Y.: Computational methods in engineering boundary value problems. In: Mathematics in Science and Engineering, Vol. 145. Accademic Press, New York (1979)Palao, J.P., Kosloff, R.: Quantum computing by an optimal control algorithm for unitry transformations. Phys. Rev. Lett. 28 (2002)Peirce, A.P., Dahleh, M.A., Rabitz, H.: Optimal control of quantum-mechanical systems: existence, numerical approximation, and applications. Phys. Rev. A 37, 4950–4967 (1988)Reid, W.T.: Riccati Differential Equations. Academic, New York (1972)Sanz-Serna, J.M., Calvo, M.P.: Numerical Hamiltonian Problems. Chapman & Hall, London (1994)Sidje, R.B.: Expokit: a software package for computing matrix exponentials. ACM Trans. Math. Software 24, 130–156 (1998)Speyer, J.L., Jacobson, D.H.: Primer on optimal control theory. SIAM, Philadelphia (2010)Starr, A.W., Ho, Y.C.: Non-zero sum differential games. J. Optim. Theory and Appl 3, 179–197 (1969)Zhu, W., Rabitz, H.: A rapid monotonically convergent iteration algorithm for quantum optimal control ever the expectation value of a positive definite operator. J. Chem. Phys. 109, 385–391 (1998

Medical students' and facilitators' experiences of an Early Professional Contact course: Active and motivated students, strained facilitators

Background: Today, medical students are introduced to patient contact, communication skills, and clinical examination in the preclinical years of the curriculum with the purpose of gaining clinical experience. These courses are often evaluated from the student perspective. Reports with an additional emphasis on the facilitator perspective are scarce. According to constructive alignment, an influential concept from research in higher education, the learning climate between students and teachers is also of great importance. In this paper, we approach the learning climate by studying both students' and facilitators' course experiences.\ud \ud In 2001, a new "Early Professional Contact" longitudinal strand through term 1–4, was introduced at the Sahlgrenska Academy, University of Gothenburg, Sweden. General practitioners and hospital specialists were facilitators.\ud \ud The aim of this study was to assess and analyse students' and clinical facilitators' experiences of the Early Professional Contact course and to illuminate facilitators' working conditions.\ud \ud Methods: Inspired by a Swedish adaptation of the Course Experience Questionnaire, an Early Professional Contact Questionnaire was constructed. In 2003, on the completion of the first longitudinal strand, a student and facilitator version was distributed to 86 students and 21 facilitators. In the analysis, both Chi-square and the Mann-Whitney tests were used.\ud \ud Results: Sixty students (70%) and 15 facilitators (71%) completed the questionnaire. Both students and facilitators were satisfied with the course. Students reported gaining [inspiration] for their future work as doctors along with increased confidence in meeting patients. They also reported increased motivation for biomedical studies. Differences in attitudes between facilitators and students were found. Facilitators experienced a greater workload, less reasonable demands and less support, than students.\ud \ud Conclusion: In this project, a new Early Professional Contact course was analysed from both student and facilitator perspectives. The students experienced the course as providing them with a valuable introduction to the physician's professional role in clinical practice. In contrast, course facilitators often experienced a heavy workload and lack of support, despite thorough preparatory education. A possible conflict between the clinical facilitator's task as educator and member of the workplace is suggested. More research is needed on how doctors combine their professional tasks with work as facilitators

SARS-CoV-2 variants, spike mutations and immune escape.

Although most mutations in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome are expected to be either deleterious and swiftly purged or relatively neutral, a small proportion will affect functional properties and may alter infectivity, disease severity or interactions with host immunity. The emergence of SARS-CoV-2 in late 2019 was followed by a period of relative evolutionary stasis lasting about 11 months. Since late 2020, however, SARS-CoV-2 evolution has been characterized by the emergence of sets of mutations, in the context of 'variants of concern', that impact virus characteristics, including transmissibility and antigenicity, probably in response to the changing immune profile of the human population. There is emerging evidence of reduced neutralization of some SARS-CoV-2 variants by postvaccination serum; however, a greater understanding of correlates of protection is required to evaluate how this may impact vaccine effectiveness. Nonetheless, manufacturers are preparing platforms for a possible update of vaccine sequences, and it is crucial that surveillance of genetic and antigenic changes in the global virus population is done alongside experiments to elucidate the phenotypic impacts of mutations. In this Review, we summarize the literature on mutations of the SARS-CoV-2 spike protein, the primary antigen, focusing on their impacts on antigenicity and contextualizing them in the protein structure, and discuss them in the context of observed mutation frequencies in global sequence datasets

Vitalism in contemporary chiropractic: a help or a hinderance?

Background: Chiropractic emerged in 1895 and was promoted as a viable health care substitute in direct competition with the medical profession. This was an era when there was a belief that one cause and one cure for all disease would be discovered. The chiropractic version was a theory that most diseases were caused by subluxated (slightly displaced) vertebrae interfering with “nerve vibrations” (a supernatural, vital force) and could be cured by adjusting (repositioning) vertebrae, thereby removing the interference with the body’s inherent capacity to heal. DD Palmer, the originator of chiropractic, established chiropractic based on vitalistic principles. Anecdotally, the authors have observed that many chiropractors who overtly claim to be “vitalists” cannot define the term. Therefore, we sought the origins of vitalism and to examine its effects on chiropractic today. Discussion: Vitalism arose out of human curiosity around the biggest questions: Where do we come from? What is life? For some, life was derived from an unknown and unknowable vital force. For others, a vital force was a placeholder, a piece of knowledge not yet grasped but attainable. Developments in science have demonstrated there is no longer a need to invoke vitalistic entities as either explanations or hypotheses for biological phenomena. Nevertheless, vitalism remains within chiropractic. In this examination of vitalism within chiropractic we explore the history of vitalism, vitalism within chiropractic and whether a vitalistic ideology is compatible with the legal and ethical requirements for registered health care professionals such as chiropractors. Conclusion: Vitalism has had many meanings throughout the centuries of recorded history. Though only vaguely defined by chiropractors, vitalism, as a representation of supernatural force and therefore an untestable hypothesis, sits at the heart of the divisions within chiropractic and acts as an impediment to chiropractic legitimacy, cultural authority and integration into mainstream health care

Mating dances and the evolution of language: What’s the next step?

The Darwinian protolanguage hypothesis is one of the most popular theories of the evolution of human language. According to this hypothesis, language evolved through a three stage process involving general increases in intelligence, the emergence of grammatical structure as a result of sexual selection on protomusical songs, and finally the attachment of meaning to the components of those songs. The strongest evidence for the second stage of this process has been considered to be birdsong, and as a result researchers have investigated the existence of various forms of grammar in the production and comprehension of songs by birds. Here, we argue that mating dances are another relevant source of sexually-selected complexity that has until now been largely overlooked by proponents of Darwinian protolanguage, focusing especially on the dances of long-tailed manakins. We end by sketching several lines of research that should be pursued to determine the relevance of mating dances to the evolution of language

Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function

Correction Volume: 10, Article Number: 2068 DOI: 10.1038/s41467-019-10160-w WOS:000466339700001General cognitive function is a prominent and relatively stable human trait that is associated with many important life outcomes. We combine cognitive and genetic data from the CHARGE and COGENT consortia, and UK Biobank (total N = 300,486; age 16-102) and find 148 genome-wide significant independent loci (P <5 x 10(-8)) associated with general cognitive function. Within the novel genetic loci are variants associated with neurodegenerative and neurodevelopmental disorders, physical and psychiatric illnesses, and brain structure. Gene-based analyses find 709 genes associated with general cognitive function. Expression levels across the cortex are associated with general cognitive function. Using polygenic scores, up to 4.3% of variance in general cognitive function is predicted in independent samples. We detect significant genetic overlap between general cognitive function, reaction time, and many health variables including eyesight, hypertension, and longevity. In conclusion we identify novel genetic loci and pathways contributing to the heritability of general cognitive function.Peer reviewe

Mitochondrial physiology

As the knowledge base and importance of mitochondrial physiology to evolution, health and disease expands, the necessity for harmonizing the terminology concerning mitochondrial respiratory states and rates has become increasingly apparent. The chemiosmotic theory establishes the mechanism of energy transformation and coupling in oxidative phosphorylation. The unifying concept of the protonmotive force provides the framework for developing a consistent theoretical foundation of mitochondrial physiology and bioenergetics. We follow the latest SI guidelines and those of the International Union of Pure and Applied Chemistry (IUPAC) on terminology in physical chemistry, extended by considerations of open systems and thermodynamics of irreversible processes. The concept-driven constructive terminology incorporates the meaning of each quantity and aligns concepts and symbols with the nomenclature of classical bioenergetics. We endeavour to provide a balanced view of mitochondrial respiratory control and a critical discussion on reporting data of mitochondrial respiration in terms of metabolic flows and fluxes. Uniform standards for evaluation of respiratory states and rates will ultimately contribute to reproducibility between laboratories and thus support the development of data repositories of mitochondrial respiratory function in species, tissues, and cells. Clarity of concept and consistency of nomenclature facilitate effective transdisciplinary communication, education, and ultimately further discovery

Retinal Vascular Fractal Dimension, Childhood IQ, and Cognitive Ability in Old Age: The Lothian Birth Cohort Study 1936

<div><p>Purpose</p><p>Cerebral microvascular disease is associated with dementia. Differences in the topography of the retinal vascular network may be a marker for cerebrovascular disease. The association between cerebral microvascular state and non-pathological cognitive ageing is less clear, particularly because studies are rarely able to adjust for pre-morbid cognitive ability level. We measured retinal vascular fractal dimension (<i>D</i>_<i>f</i>) as a potential marker of cerebral microvascular disease. We examined the extent to which it contributes to differences in non-pathological cognitive ability in old age, after adjusting for childhood mental ability.</p><p>Methods</p><p>Participants from the Lothian Birth Cohort 1936 Study (LBC1936) had cognitive ability assessments and retinal photographs taken of both eyes aged around 73 years (<i>n</i> = 648). IQ scores were available from childhood. Retinal vascular <i>D</i>_<i>f</i> was calculated with monofractal and multifractal analysis, performed on custom-written software. Multiple regression models were applied to determine associations between retinal vascular <i>D</i>_<i>f</i> and general cognitive ability (<i>g</i>), processing speed, and memory.</p><p>Results</p><p>Only three out of 24 comparisons (two eyes × four <i>D</i>_<i>f</i> parameters × three cognitive measures) were found to be significant. This is little more than would be expected by chance. No single association was verified by an equivalent association in the contralateral eye.</p><p>Conclusions</p><p>The results show little evidence that fractal measures of retinal vascular differences are associated with non-pathological cognitive ageing.</p></div

Genetic and lifestyle risk factors for MRI-defined brain infarcts in a population-based setting

OBJECTIVE: To explore genetic and lifestyle risk factors of MRI-defined brain infarcts (BI) in large population-based cohorts.METHODS: We performed meta-analyses of genome-wide association studies (GWAS) and examined associations of vascular risk factors and their genetic risk scores (GRS) with MRI-defined BI and a subset of BI, namely, small subcortical BI (SSBI), in 18 population-based cohorts (n = 20,949) from 5 ethnicities (3,726 with BI, 2,021 with SSBI). Top loci were followed up in 7 population-based cohorts (n = 6,862; 1,483 with BI, 630 with SBBI), and we tested associations with related phenotypes including ischemic stroke and pathologically defined BI.RESULTS: The mean prevalence was 17.7% for BI and 10.5% for SSBI, steeply rising after age 65. Two loci showed genome-wide significant association with BI: FBN2, p = 1.77 × 10-8; and LINC00539/ZDHHC20, p = 5.82 × 10-9. Both have been associated with blood pressure (BP)-related phenotypes, but did not replicate in the smaller follow-up sample or show associations with related phenotypes. Age- and sex-adjusted associations with BI and SSBI were observed for BP traits (p value for BI, p [BI] = 9.38 × 10-25; p [SSBI] = 5.23 × 10-14 for hypertension), smoking (p [BI] = 4.4 × 10-10; p [SSBI] = 1.2 × 10-4), diabetes (p [BI] = 1.7 × 10-8; p [SSBI] = 2.8 × 10-3), previous cardiovascular disease (p [BI] = 1.0 × 10-18; p [SSBI] = 2.3 × 10-7), stroke (p [BI] = 3.9 × 10-69; p [SSBI] = 3.2 × 10-24), and MRI-defined white matter hyperintensity burden (p [BI] = 1.43 × 10-157; p [SSBI] = 3.16 × 10-106), but not with body mass index or cholesterol. GRS of BP traits were associated with BI and SSBI (p ≤ 0.0022), without indication of directional pleiotropy.CONCLUSION: In this multiethnic GWAS meta-analysis, including over 20,000 population-based participants, we identified genetic risk loci for BI requiring validation once additional large datasets become available. High BP, including genetically determined, was the most significant modifiable, causal risk factor for BI.</p