CBC Annual Report 2015

Total Page:16

File Type:pdf, Size:1020Kb

CBC Annual Report 2015 Center for Biomedical Computing Annual Report 2015 Last year’s annual report emphasized that CBC is approach- work. In addition to CBC researchers being core member of the ing the end of its funding period, and this fact is of no less impor- FEniCS development community, we continue to fund key tech- tance this year. The focus has remained on reaching milestones nical personnel responsible for software maintenance and less and deliverables laid out in the CBC research plan, and on posi- research-oriented development tasks. As pointed out in several tioning the research team for securing continued funding. Both previous annual reports, this kind of work is essential for bring- focus areas have seen considerable success through 2015. Sev- ing advanced software tools out to a larger community, and it is eral new research projects have been granted from the Research normally very difficult to fund with smaller and more short-term Council of Norway and various international funding agencies, research grants. The substantial long-term funding of a Centerof which serve to reduce the negative impact of the Centre of Ex- Excellence has therefore been essential in providing these ser- cellence grant ending in 2017. However, none of the granted vices to the FEniCS community, and this contribution may end projects have the broad, long-term scope of the CBC, and the up as one of the most important long-term results of the CBC loss of a large and stable funding source will impact the research funding. In May 2016, CBC will host the annual FEniCS work- environment. The group therefore continues to seek new oppor- shop, again manifesting its role as a key partner in the FEniCS tunitiesforcontinuedfundingof thekey activities, and we expect consortium. this situtation to be the normal in years to come. With respect Since the very beginning, CBC has had the aim of both pro- to completing the planned research for CBC milestones, a sub- viding generally applicable computational software tools, and stantial number of the original milestones have been reached at address specific biomedical applications. This split focus has this point, and we are making goodprogresstowardsreaching the continued through 2015, and the continuation and expansion of major software releases and milestones scheduled for the end of our biomedical research is an important part of the Centre’s exit the CBC funding period. strategy. The integration of advanced ICT tools in health care A particularly welcome eventforCBCin 2015wasthe award is a strategic focus area for both National and European author- of the Wilkinson Prize for Numerical Software to Patrick Farrell, ities and funding agencies. Moving towards a stronger focus on Simon Funke, David Ham, and Marie Rognes, for the develop- biomedical applications, and expanding our target areas beyond ment of dolfin-adjoint. The Wilkinson Prize is awarded every the current focus on biofluids and cardiovascular research, are fourth year at the International Congress on Industrial and Ap- important steps for increasing the potentialfor continued funding plied Mathematics, and is one of the most prestigious awards in and for high impact research. Over the last year, computational the fields of numerical mathematics and computationalsoftware. neuroscience has grown from a niche activity on cerebrospinal The dolfin-adjoint software automatically derives and solves ad- fluid flow to become a major application area for CBC research. joint and tangent linear equations from high-level mathematical Many of the tools and techniques developed for cardiovascular specifications of finite element discretizations of partial differen- research can be adapted for use in neuroscience, and the field tial equations. Adjoint mathematical problems have a wide range has a large potentialfor high-impact research. We are pleasedto of applicability in engineering, in particular for optimization and see a number of small-scale research collaborations in this area, sensitivityanalysis,andwhilederivingtheadjointofalinear equa- as well as funding initiatives on various scales. tion is relatively straightforward, performing the same task for a In conclusion, we are pleased with the current state of CBC non-linear or time dependent problem may be extremely difficult. and of computational science in general at our host institution. This task is automated by dolfin-adjoint, and thereby making the The long-term, stable funding combined with the visibility of the powerful analysis based on adjoints available to a much wider Center of Excellence brand has helped us create a vibrant re- group of researchers and engineers. search environment that attracts excellent research talent. At Although dolfin-adjoint has been the most celebrated soft- the time of writing CBC has just entered its last year of funding wareprojectthrough2015,workhascontinuedsteadilyonadding and CoE status, but the future still looks bright for the research featuresand improving other partsof the FEniCS software frame- environment that has been created. 2 WORDS FROM THE DIRECTORS 2 RUNNING THE CENTER 4 Project overview 4 The Scientific Advisory Board (SAB) 4 People, Recruitment and Gender Diversity 4 Gender diversity 5 SCIENTIFIC ACTIVITIES 6 Computational Middleware 6 Robust Solvers 7 Biomedical Flows and Structures 8 CBC at NTNU 9 Cardiac Computations 10 FEATURED RESEARCH 2015: 11 Generalized geometry handling using Multimesh 11 EDUCATION AND OUTREACH 13 The Simula School of Research and Innovation 13 University Teaching 13 APPENDIX 14 Staff 14 Accounting and Budget 18 Publications 20 Conferences, Workshops and Seminars 28 Other Activities 30 Refereeing Activities 30 Prizes and Recognition 30 Committee Work 30 Organization of Minisymposia and Workshops at Conferences 30 Editorial Boards 31 Invited talks 31 Collaboration partners 33 List of International Guests in 2015 36 3 Project overview The SAB convened at CBC in February 2015, and was pre- sented with research, education and innovation activities within All research activities in CBC are organized in large projects with the scope of the Center’s activity. The generalfeedback from the an appointed project leader and a well-defined project plan. board was that they were pleased with the obtained results and 1) Computational Middleware. This project is devoted to devel- supported the overall choice of research topics and strategy. On oping generic high-performance software componentsforbuilding a more detailed level, the SAB made numerous useful comments the simulation programs needed in the Center. The results also and suggestions, in particular for increased collaboration with serve the global computational science community with a new medical industry. Although most of these collaboration leads are generation of widely applicable computational software. At the of a long-term nature, they may be extremely valuable for in- core of the developmentis the finite elementsoftwaresuite FEn- creasing visibility and relevance of CBC research, and potentially iCS. for securing valuable funding. Several of the SAB’s recommen- dations regarding sustainability and funding opportunities have 2) Robust Solvers. The Robust Solvers project focuses on ef- already been addressed and acted upon, as they will have impact ficient and stable numerical methods with error and uncertainty on our exit strategy. estimation, as well as implementation of such methods for prob- lems arising in the two application projects Cardiac Computations (CC) and Biomedical Flows and Structures (BFS). People, Recruitment and 3) Cardiac Computations. This project performs research on med- Gender Diversity ical problems involving models of heart electrophysiology and Recruitmentto the Center. We recruitedseveral newpeopleto mechanics. The current activity is closely aligned with the goals the Center’s activity during 2015, both through the RCN funded ofthe Center for CardiologicalInnovation(CCI), whichisaCenter part of the Center and through research collaboration and exter- for Research Based Innovation funded by the Research Council nallyfundedactivities. All in all we recruited9 new people to the of Norway. Center: 1 postdoctoral fellow, 6 PhDs, and 2 Research trainees. The Cardiac Computations staff was expanded with a Re- 4) Biomedical Flows and Structures. Research in the Biomedi- search Trainee (Kariline Horgmo Jæger), funded by our host in- cal Flows and Structures project centers around biomedical flow stitution. Dr. Benjamin Ragan-Kelley (funded through the Moore and tissue interaction problems of high clinical importance. In Foundation) joined the Computational Middleware project to con- a short to medium time frame, the applications to be targeted tinue hiswork on IPythonand Jupyter. The BiomedicalFlows and are aneurysm formation and rupture in the Circle of Willis, the Structures project has been strengthened with 3 PhD students: relation between Chiari I malformation and cyst formation in the Giulia Pizzichelli (IIT), Karl Erik Holter (SUURPh), Lars Magnus spinal cord, modeling of large cardiovascular networks, and fluid- Valnes (UiO), and one Research Trainee (Aslak Bergersen). In structure interactions in cardiovascular biomechanical systems. addition the PhD student Jacob Sturdy (NTNU) has joined the The research on the lattertwo topics is headedby our partnersat CBC@NTNU node. Finally, the Robust Solvers project strength- the Biomechanics Division at the Department of Structural Engi- ened its work on high performance computing by collaboration neering at the Norwegian University
Recommended publications
  • CV & Publication List
    +44 1865 270743 [email protected] https://pefarrell.org Patrick E. Farrell orcid:0000-0002-1241-7060 Employment history 2016–date Associate Professor, Mathematical Institute, University of Oxford. 2016–date Tutorial Fellow in Applied Mathematics, Oriel College, University of Oxford. 2013–2018 EPSRC Early Career Research Fellow, Mathematical Institute, University of Oxford. 2013–2016 Postdoctoral Research Fellow, Christ Church College, University of Oxford. 2012–2017 Adjunct Research Scientist, Simula Research Laboratory, Oslo. 2010–2013 Postdoctoral Research Associate, Earth Science & Engineering, Imperial College London. Academic history 2006–2009 PhD in Computational Physics, Imperial College London. Thesis Galerkin projection of discrete fields via supermesh construction. Viva: 27 Nov 2009. Prizes Association of Computational Mechanics in Engineering award, 2010; Finalist, European Community on Computational Methods in Applied Sciences award, 2010; Imperial College Research Excellence Award, 2010; Janet Watson award, Imperial College London, 2009. 2002–2006 BSc (Hons) in Mathematics, National University of Ireland, Galway. Thesis Cryptographic applications of polycyclic groups. Prizes Hamilton Prize, Royal Irish Academy, 2006; Blayney Exhibition, National University of Ireland, 2006. Notable prizes 2021 Whitehead Prize, awarded by the London Mathematical Society. 2021 Broyden Prize in Optimization, for the best paper published in Optimization Methods and Software in 2020, with M. Croci and T. M. Surowiec. 2015 Wilkinson Prize for Numerical Software, for the development of dolfin-adjoint, with D. A. Ham, S. W. Funke and M. E. Rognes. 2015 Leslie Fox Prize in Numerical Analysis, second place, for the development of deflation techniques for partial differential equations. Research funding 2020–2021 Gen X: ExCALIBUR working group on exascale continuum mechanics through code generation, EP/V001493/1, £174,252, Co-I.
    [Show full text]
  • The Leslie Fox Prize
    The leslie Fox Prize As many people will know, Professor Leslie Fox retired successful doctoral candidates plus twelve MSc students. in September 1983, from the Chair of Numerical Analy­ Several of these are themselves now lecturing in univer­ sis at Oxford University and in the previous year from sities and have supervised postgraduates-giving Leslie the Directorship of the Oxford University Computing many "grand-children" and even "great-grand­ Laboratory. Leslie has inspired many talented scholars children," in whom he takes justifiable pride. A "family to enter the field of numerical analysis and has done tree" is being constructed but is so far very incomplete, much to develop educational and research opportunities and I would welcome any information from Leslie's for students of computational mathematics. To show our former students, especially those living abroad. (See appreciation of his devotion and to continue his work of Table 1.) The records available show that Leslie's own encouragement to others, a capital fund has been es­ students are now employed in fourteen different coun­ tablished in his honour fot the award of an annual prize tries and his students' students in (at least) a further seven to a promising numerical analyst. This year the Leslie countries. Fox Prize is to be awarded for the first time. Besides the many postgraduates from the Computing Laboratory, Leslie has touched a much wider audience through his undergraduate teaching and especially through his work with the Open University. He has continued to encourage the introduction of computa­ tional mathematics into the undergraduate syllabus of universities and polytechnics and also into the school curriculum.
    [Show full text]
  • The 1988 Leslie Fox Prize Meeting
    The 1988 Leslie Fox Prize Meeting THE third Leslie Fox Prize Meeting was held at Imperial Jackson (Cambridge) on "An order of convergence for College, London on Monday, March 28th, 1988. Five radieU basis functions." After lunch T. Tang (Leeds) finalists had been selected to present papers from an spoke on "The spline collocation methods for non­ entry of seventeen, which included entries from the UK, standard Volterra integro-differential equations," Europe, China and America. The finalists were chosen followed by the final speaker of the day, N. Higham by the adjudicating committee for the originality and (Manchester) on "Analysis of the Cholesky decomposi­ quality of their papers, together with the suitability of the tion of a semi-definite matrix." The adjudicators then material for presentation to a general audience of ' met to consider the awarding of prizes. numerical analysts. The five papers chosen for A first prize deservedly went to N. Higham, the only presentation provided a varied programme, covering person to date to have been a finalist at more than one widely differing areas of numerical analysis. Leslie Fox Prize Meeting. The talk given by N. Higham The meeting was divided into a morning and an not only contained new and interesting results but was afternoon session. The morning session opened with a clearly and carefully presented. Each of the other talk by T. Hagstrom (SUNY) on "Asymptotic boundary finalists was presented with a second prize. conditions for computational studies of wave propaga­ The whole day at Imperial College was, I believe, tion." This was followed by P.
    [Show full text]
  • An Interview with IAIN S. DUFF Conducted by Thomas Haigh on 30
    An interview with IAIN S. DUFF Conducted by Thomas Haigh On 30 & 31st August, 2005 In East Hagbourne and Chilton, United Kingdom Interview conducted by the Society for Industrial and Applied Mathematics, as part of grant # DE-FG02-01ER25547 awarded by the US Department of Energy. Transcript and original tapes donated to the Computer History Museum by the Society for Industrial and Applied Mathematics © Computer History Museum Mountain View, California ABSTRACT Iain S. Duff explores the whole of his career to date. Duff grew up in Glasgow, attending The High School before earning a first class degree in Mathematics and Natural Philosophy from the University of Glasgow in 1969. He discusses the mathematical curriculum there, and recounts his early computing experiences as a summer fellow with IBM. Duff won a Carnegie Fellowship for postgraduate work at Oxford, culminating in a D.Phil. conferred in 1972. At Oxford he worked with the numerical analysis group and wrote his thesis on the analysis of sparse systems under the direction of Leslie Fox. While working toward his degree, Duff was already collaborating with Alan Curtis, Michael Powell, and John Reid at the nearby AERE Harwell nuclear power research center. His interest in sparse matrices was cemented by a 1971 visit to IBM’s Yorktown Heights laboratory for a conference, where he first met colleagues such as Alan George, Ralph Willoughby and Fred Gustavson. In 1972, Duff went to the United States on a Harkness postdoctoral fellowship, spending time at SUNY Stony Brook with Reggie Tewarson and Stanford, where he got to know Gene Golub.
    [Show full text]
  • 2016 EMI International Conference Engineering Mechanics Institute
    10434Ce qui S8MS2 ² 2016 EMI International Conference Engineering Mechanics Institute PROGRAM 2016 OCTOBER 2 5 – 2 7 , 2016, M E T Z F R A N C E Conference area layout 1 3 Welcome from the conference Chair Professor M. Hattab Université de Lorraine, Metz, France On behalf of the Steering Committee, the Local Organizing Committee and the International Scientific Committee, I would like to welcome you to the Université de Lorraine, at the Campus of Ile du Saulcy-Metz, site of the second International Conference, after the 2015 conference at the Hong Kong Polytechnic University that the Engineering Mechanics Institute (EMI) of the American Society of Civil Engineers (ASCE) organized outside the United States. This important event will be taking place for the first tim e in Europe. 2016 EMI International Conference feature about 320 oral presentations in the 24 Mini-Symposia organized in parallel sessions. The participants come from 36 countries in the six continents , with impressive participation of young researchers and PhD students who will have a great opportunity to exchange on an international platform. The Plenary Keynote lectures by six distinguished research leaders come from International Universities, the contributions will cover the most challenging fields i n the context of Engineering Mechanics and Materials Science. In addition, several large mini-symposia will start with keynote presentations by leading researchers in the mini-symposia topics. Large areas of research will be discussed during the conference in the different Mini-symposia, including: Porous Media, Severe Plastic Deformation, Damage and Failure, Biomechanics, Biomaterials, Granular Mechanics, Uncertainty Quantification, Non-Linear Dynamics, Instability and Bifurcations, Multiscale Modeling, Co mposite Materials and Structures, Creep and Relaxation Mechanics, Cementitious Materials, Geomechanics, Confined Structures, Laterites Clayey Soils, Model Order Reduction, and Modeling and Testing of Infrastructure.
    [Show full text]
  • An Introduction to the Quality of Computed Solutions1
    An Introduction to the Quality of Computed Solutions1 Sven Hammarling The Numerical Algorithms Group Ltd Wilkinson House Jordan Hill Road Oxford, OX2 8DR, UK [email protected] October 12, 2005 1Based on [Hammarling, 2005] CONTENTS 1 Contents 1 Introduction 2 2 Floating Point Numbers and IEEE Arithmetic 2 3 Why Worry about Computed Solutions? 5 4 Condition, Stability and Error Analysis 9 4.1 Condition ....................................... 9 4.2 Stability....................................... 14 4.3 ErrorAnalysis................................... 19 5 Floating Point Error Analysis 23 6 Posing the Mathematical Problem 29 7 Error Bounds and Software 30 8 Other Approaches 34 9 Summary 34 Bibliography 35 1 INTRODUCTION 2 1 Introduction This report is concerned with the quality of the computed numerical solutions of mathematical problems. For example, suppose we wish to solve the system of linear equations Ax = b using a numerical software package. The package will return a computed solution, say x˜, and we wish to judge whether or not x˜ is a reasonable solution to the equations. Sadly, all too often software packages return poor, or even incorrect, numerical results and give the user no means by which to judge the quality of the numerical results. In 1971, Leslie Fox commented [Fox, 1971, p. 296] “I have little doubt that about 80 per cent. of all the results printed from the computer are in error to a much greater extent than the user would believe, ...” More than thirty years on that paper is still very relevant and worth reading. Another very readable article is Forsythe [1970]. The quality of computed solutions is concerned with assessing how good a computed solution is in some appropriate measure.
    [Show full text]
  • IMA-Fox-Prize-1993.Pdf
    Symposium in Honour of Leslie Fox and the Sixth Leslie Fox Prize Competition I. S. DUFF, C.Math., FIMA Atlas Centre, Rut herford Appleton Laborato ry, Didcot, Oxon THE Sixth Leslie Fox Prize Competition was the first to be held since Leslie's death and was combined with a Symposium in his honour. Both were held in the Nuclear Physics Lecture Theatre in Oxford, the Prize meeting on Thursday June 24th, 1993 with the Symposium starting that evening and continuing on the Friday. The attendance on both days of around 80 was very good given the fact that many British Universities were still involved with examiners' meetings and other end-of-session delights. Even more gratifying was the number of participants from overseas, several of whom were visitors to Oxford when Leslie was there. The audience was treated to two excellent days of numerical analysis, the first in the main serious, the second with many reminiscences mixed into some very entertaining and thought provoking presentations. Speakers, chairmen, and attendees all combined to make Clemency Fox presenting First Prize to Yuying Li a very fitting tribute to Leslie Fox. The Adjudicating Committee for the Fox Prize consisted of Nancy Nichols (Chairman from University of Reading), Charlie Elliott (University of Sussex) and Christopher Baker (University of Manchester). They had the difficult task, first of selecting finalists from nineteen very high quality entrants, then the even more unenviable task of choosing prizewinners from the finalists. The finalists and their subjects are given below and, since I do not feel I could do justice to their excellent papers and presentations by summarising them here, I leave it to readers to enquire further directly with the authors.
    [Show full text]
  • Brian Ford SIAM Oral History, with Steve Hague and Sven Hammarling; 102746792
    An interview with BRIAN FORD, STEVE HAGUE & SVEN HAMMARLING Conducted by Thomas Haigh On 29th and 30th June, 2004 In Oxford, United Kingdom Interview conducted by the Society for Industrial and Applied Mathematics, as part of grant # DE-FG02-01ER25547 awarded by the US Department of Energy. Transcript and original tapes donated to the Computer History Museum by the Society for Industrial and Applied Mathematics © Computer History Museum Mountain View, California ABSTRACT Brian Ford discusses his entire career. Born in Nottingham in 1940, Ford studied mathematics at Imperial College, London. He worked as a school teacher for several years before earning a M.Sc. and Ph.D. from Nottingham University. His work on a QR algorithm brought him into contact with Jim Wilkinson and Leslie Fox. Ford worked with the university’s KDF-9 computer center, building a user advice service and numerical software library. In 1970, as the university planned for the arrival of an ICL 1906A, Ford took part in the launch of a collaborative effort with other universities to create a mathematical library for the new machine. This became NAG, initially the Nottingham Algorithms Group but later the Numerical Algorithms Group. Facing lack of institutional support at Nottingham, NAG moved to Oxford in 1973. NAG and the American IMSL were the main vendors of commercial numerical libraries from the mid-1970s on. The library was initially offered in both FORTRAN and Algol versions. In 1976 it was chartered as a not-for-profit company, led by Ford, and it has been financially self sustaining since 1980. Ford discusses the changing institutional form of NAG, the organization and development of its library, its products and marketing, the tools it created to automate aspects of its work, the sources of its code and relationship to its contributors, and its relationship to its customer base.
    [Show full text]
  • An Introduction to the Quality of Computed Solutions
    An Introduction to the Quality of Computed Solutions Hammarling, Sven 2005 MIMS EPrint: 2005.29 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available from: http://eprints.maths.manchester.ac.uk/ And by contacting: The MIMS Secretary School of Mathematics The University of Manchester Manchester, M13 9PL, UK ISSN 1749-9097 An Introduction to the Quality of Computed Solutions1 Sven Hammarling The Numerical Algorithms Group Ltd Wilkinson House Jordan Hill Road Oxford, OX2 8DR, UK [email protected] December 5, 2005 1Based on [Hammarling, 2005], Chapter 4 of [Einarsson, 2005], which is available for pur- chase from SIAM at http://ec-securehost.com/SIAM/SE18.html. The royalties go to the SIAM Student Travel Fund. CONTENTS 1 Contents 1 Introduction 2 2 Floating Point Numbers and IEEE Arithmetic 2 3 Why Worry about Computed Solutions? 5 4 Condition, Stability and Error Analysis 9 4.1 Condition ....................................... 9 4.2 Stability . 14 4.3 ErrorAnalysis..................................... 19 5 Floating Point Error Analysis 23 6 Posing the Mathematical Problem 29 7 Error Bounds and Software 30 8 Other Approaches 34 9 Summary 34 Bibliography 35 1 INTRODUCTION 2 1 Introduction This report is concerned with the quality of the computed numerical solutions of mathematical problems. For example, suppose we wish to solve the system of linear equations Ax = b using a numerical software package. The package will return a computed solution, say x˜, and we wish to judge whether or not x˜ is a reasonable solution to the equations. Sadly, all too often software packages return poor, or even incorrect, numerical results and give the user no means by which to judge the quality of the numerical results.
    [Show full text]