{"id":1229,"date":"2015-07-13T12:00:33","date_gmt":"2015-07-13T17:00:33","guid":{"rendered":"https:\/\/blogs.mathworks.com\/cleve\/?p=1229"},"modified":"2015-07-12T18:29:19","modified_gmt":"2015-07-12T23:29:19","slug":"trip-report-naconf-2015-and-sparse-days-iii","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/cleve\/2015\/07\/13\/trip-report-naconf-2015-and-sparse-days-iii\/","title":{"rendered":"Trip Report: NACONF 2015 and Sparse Days III"},"content":{"rendered":"<div class=\"content\"><!--introduction--><p>I have just returned from two meetings in Europe, the 26th Biennial Conference on Numerical Analysis at the University of Strathclyde in Glasgow, Scotland, and Sparse Days III in Saint-Girons, France.<\/p><!--\/introduction--><h3>Contents<\/h3><div><ul><li><a href=\"#23be48a9-3344-4d3a-b521-97a86ff8854c\">26th Biennial Conference on Numerical Analysis<\/a><\/li><li><a href=\"#2a0636b0-2fea-4f0c-a145-8e4962dee8b8\">University of Strathclyde<\/a><\/li><li><a href=\"#ff7c50fa-0eb6-4ab7-8ccc-bcd63f8e1f00\">The Conference<\/a><\/li><li><a href=\"#f537bd90-5be3-48a9-954c-b489f657476d\">A. R. Mitchell Lecture<\/a><\/li><li><a href=\"#1d80cdb3-680c-4629-84a9-eae6976ebf98\">Fletcher-Powell Lecture<\/a><\/li><li><a href=\"#bce4fbb8-ccb0-468e-99e6-24f02f6b68b5\">Sparse Days in St Girons III<\/a><\/li><li><a href=\"#d2300a3f-ea1b-43b6-813d-04ae1a690a98\">St Girons<\/a><\/li><li><a href=\"#29152f11-96bd-491b-a325-bbef05f97611\">The Conference<\/a><\/li><li><a href=\"#518e0c3c-2362-457f-b73c-fcd265535569\">Interesting Talks<\/a><\/li><li><a href=\"#65eddf8b-36a3-46d0-8cca-bee2e8aa7879\">Another Cleve<\/a><\/li><\/ul><\/div><h4>26th Biennial Conference on Numerical Analysis<a name=\"23be48a9-3344-4d3a-b521-97a86ff8854c\"><\/a><\/h4><p>The <a href=\"http:\/\/numericalanalysisconference.org.uk\/\">Biennial Conference on Numerical Analysis<\/a> in Scotland has a long history. The first two meetings were held at the University of St. Andrews in 1965 and 1967.  The meetings moved to the University of Dundee in 1969. They were held there every two years, under the leadership of Ron Mitchell, until his retirement in 2007.  They moved to the University of Strathclyde in 2009, under the leadership of Alison Ramage.  This meeting was the 26th and the 50th anniversary.  Although this is a premier conference series in numerical analysis, it was the first time I participated.<\/p><h4>University of Strathclyde<a name=\"2a0636b0-2fea-4f0c-a145-8e4962dee8b8\"><\/a><\/h4><p>The University of Strathclyde is a public research university located in Glasgow, Scotland.  The university's name comes from the \"Valley of the River Clyde\".  The school is Scotland's foremost institution of science and engineering.<\/p><h4>The Conference<a name=\"ff7c50fa-0eb6-4ab7-8ccc-bcd63f8e1f00\"><\/a><\/h4><p>The conference was held June 23rd through 26th.  I was one of a dozen invited speakers -- five from the US, four from continental Europe, and three from the UK.<\/p><p>Almost all of the participants gave half-hour talks.  There were 168 talks organized in seven parallel sessions.  Many of the talks were in one of twelve minisymposia.  Sample minisymposia titles include \"Stable and accurate discretisations for convection-dominated problems\", \"City analytics\", \"Chebfun: new developments, cool applications and on the horizon\", and \"Numerical linear algebra for optimsation and data assimilation\".<\/p><h4>A. R. Mitchell Lecture<a name=\"f537bd90-5be3-48a9-954c-b489f657476d\"><\/a><\/h4><p>Two of the invited lectures are given special emphasis to honor British numerical analysts.  The A. R. Mitchell Lecture honors the University of Dundee's Ron Mitchell who was the dominant force in this conference for forty years.<\/p><p>This year the lecture was given by Mike Giles of the University of Oxford on \"Multilevel Monte Carlo Methods\".  I knew nothing about this subject before his talk and I learned a great deal.  His abstract provides a link to a web page featuring a survey paper and MATLAB codes, <a href=\"http:\/\/people.maths.ox.ac.uk\/gilesm\/acta\/\">&lt;http:\/\/people.maths.ox.ac.uk\/gilesm\/acta<\/a>&gt;<\/p><h4>Fletcher-Powell Lecture<a name=\"1d80cdb3-680c-4629-84a9-eae6976ebf98\"><\/a><\/h4><p>The Fletcher-Powell Lecture honors two mathematicians known for their work in optimization algorithms, including the Davidon-Fletcher-Powell, DFP, formula.  Mike Powell, from Cambridge University, had passed away in April.  Roger Fletcher, from the University of Dundee, attended the lecture.<\/p><p>This year the lecture was given by my long-time friend Mike Saunders from Stanford.  He talked about \"Experiments with linear and nonlinear optimization using Quad precision.\"  He has used the GFortran compiler to build a version of his MINOS optimization software with the REAL*16 floating point datatype.  He tackled flux balance analysis models of metabolic networks that involve coefficients ranging over 15 or 16 orders of magnitude.  Ordinary double precision, that is Fortran's REAL*8, cannot do the job.<\/p><h4>Sparse Days in St Girons III<a name=\"bce4fbb8-ccb0-468e-99e6-24f02f6b68b5\"><\/a><\/h4><p>The history of the <a href=\"\">Sparse Days in St Girons conference<\/a> is itself sparse. There have been just two previous conferences, one in 1994 and one in 2003. The conference is part of a program on \"High Performance Linear and Nonlinear Methods for Large Scale Applications\" organized by CIMI, the French Centre International de Mathematique et d'Informatique in Toulouse.<\/p><p>Organizers of Sparse Days in St Girons III included my buddies Iain Duff, who splits his time between Rutherford Appleton Laboratory in the UK and the Parallel Algorithms Group at CERFACS in Toulouse; Jack Dongarra, from the University of Tennessee; and Jim Demmel, from the University of California, Berkeley.<\/p><h4>St Girons<a name=\"d2300a3f-ea1b-43b6-813d-04ae1a690a98\"><\/a><\/h4><p>St Girons is a picturesque town in the French Pyrenees, not far from the border with Spain.  The Tour de France sometimes passes through here.  We met in the town's only cinema, which is air conditioned. This turned out to be fortunate planning because a heat wave hit that week, June 28th through July 2nd, with temperatures in the mid 30s C, which is mid 90s F.<\/p><h4>The Conference<a name=\"29152f11-96bd-491b-a325-bbef05f97611\"><\/a><\/h4><p>The format for this conference was quite different from the one the previous week in Scotland.  There were only about half as many attendees, around 100.  There were no invited talks and no parallel sessions.  Anyone who wanted to give a talk gave one.  There were 42 talks, most of them half an hour.<\/p><h4>Interesting Talks<a name=\"518e0c3c-2362-457f-b73c-fcd265535569\"><\/a><\/h4><p>Tim Davis, MathWorks consultant on sparse matrices, who recently moved to Texas A&amp;M, talked about \"Sparse SVD, and a GPU-accelerated sparse QR\".<\/p><p>Bora Ucar of Ecoles Normales Superieures in Lyon talked about \"Two approximation algorithms for bipartite matching on multicore architectures.\" The analysis of his algorithm happened to involve my old friend, the <a href=\"https:\/\/blogs.mathworks.com\/cleve\/2013\/09\/02\/the-lambert-w-function\/\">Lambert W function<\/a>.<\/p><p>Joost Rommes, from Mentor Graphics in the Netherlands, talked about \"Challenges in numerical simulation of electrical networks.\"<\/p><p>Alex Pothen, from Purdue, talked about \"The Virtual Scalpel: real-time matrix computations and finite elements for surgery.\"<\/p><p>Dan Sorensen, from Rice, talked about \"A DEIM induced CUR factorization.\" This title means that C is some columns of a matrix A, R is some rows of A, and DEIM is an algorithm that finds U so that the product CUR is a good low approximation to A.<\/p><p>Martin Gander, from Universite de Geneve, talked about \"Five decades of time parallel integration.\"  I was particularly interested in his talk because I knew something about the subject many years ago, but was not familiar with recent developments.  At first, it might seem impossible to parallelize a computation that involves stepping along in time.  But Martin gave an overview of four different approaches to how this can be done.<\/p><h4>Another Cleve<a name=\"65eddf8b-36a3-46d0-8cca-bee2e8aa7879\"><\/a><\/h4><p>There is another Cleve in this business.  Cleve Ashcraft of LSTC, Livermore Software Technology Corporation, is actually a \"grandstudent\".  His PhD advisor at Yale was Stan Eisenstat, who was my PhD advisee when I was a visiting professor at Stanford.  At Sparse Days, Cleve talked about \"Separability, partitions and coverings.\".  His talk didn't mention the applications he is deeply involved in at LSTC, which are the developers, among other things, of the crash analysis software LS-DYNA.  If you want to see a hint of some serious use of dynamic finite element calculations, take at the <a href=\"http:\/\/www.lstc.com\/\">LSTC web pages<\/a>.<\/p><p>This is the group photo from Sparse Days in St Girons III.  There are two Cleves.  Cleve Ashcraft is the guy with the green shirt and grey beard in the front row.  Can you find the other one? (Thanks to Pierre-Henri Cros for the photo, and for handling local arrangements.)<\/p><p><img decoding=\"async\" vspace=\"5\" hspace=\"5\" src=\"https:\/\/blogs.mathworks.com\/images\/cleve\/sparse_days.jpg\" alt=\"\"> <\/p><script language=\"JavaScript\"> <!-- \r\n    function grabCode_8c72cc5a5eef45b1843ceb74d036b5fd() {\r\n        \/\/ Remember the title so we can use it in the new page\r\n        title = document.title;\r\n\r\n        \/\/ Break up these strings so that their presence\r\n        \/\/ in the Javascript doesn't mess up the search for\r\n        \/\/ the MATLAB code.\r\n        t1='8c72cc5a5eef45b1843ceb74d036b5fd ' + '##### ' + 'SOURCE BEGIN' + ' #####';\r\n        t2='##### ' + 'SOURCE END' + ' #####' + ' 8c72cc5a5eef45b1843ceb74d036b5fd';\r\n    \r\n        b=document.getElementsByTagName('body')[0];\r\n        i1=b.innerHTML.indexOf(t1)+t1.length;\r\n        i2=b.innerHTML.indexOf(t2);\r\n \r\n        code_string = b.innerHTML.substring(i1, i2);\r\n        code_string = code_string.replace(\/REPLACE_WITH_DASH_DASH\/g,'--');\r\n\r\n        \/\/ Use \/x3C\/g instead of the less-than character to avoid errors \r\n        \/\/ in the XML parser.\r\n        \/\/ Use '\\x26#60;' instead of '<' so that the XML parser\r\n        \/\/ doesn't go ahead and substitute the less-than character. \r\n        code_string = code_string.replace(\/\\x3C\/g, '\\x26#60;');\r\n\r\n        copyright = 'Copyright 2015 The MathWorks, Inc.';\r\n\r\n        w = window.open();\r\n        d = w.document;\r\n        d.write('<pre>\\n');\r\n        d.write(code_string);\r\n\r\n        \/\/ Add copyright line at the bottom if specified.\r\n        if (copyright.length > 0) {\r\n            d.writeln('');\r\n            d.writeln('%%');\r\n            if (copyright.length > 0) {\r\n                d.writeln('% _' + copyright + '_');\r\n            }\r\n        }\r\n\r\n        d.write('<\/pre>\\n');\r\n\r\n        d.title = title + ' (MATLAB code)';\r\n        d.close();\r\n    }   \r\n     --> <\/script><p style=\"text-align: right; font-size: xx-small; font-weight:lighter;   font-style: italic; color: gray\"><br><a href=\"javascript:grabCode_8c72cc5a5eef45b1843ceb74d036b5fd()\"><span style=\"font-size: x-small;        font-style: italic;\">Get \r\n      the MATLAB code <noscript>(requires JavaScript)<\/noscript><\/span><\/a><br><br>\r\n      Published with MATLAB&reg; R2015a<br><\/p><\/div><!--\r\n8c72cc5a5eef45b1843ceb74d036b5fd ##### SOURCE BEGIN #####\r\n%% Trip Report: NACONF 2015 and Sparse Days III\r\n% I have just returned from two meetings in Europe, the 26th Biennial\r\n% Conference on Numerical Analysis at the University of Strathclyde\r\n% in Glasgow, Scotland, and Sparse Days III in Saint-Girons, France.\r\n\r\n%% 26th Biennial Conference on Numerical Analysis\r\n% The <http:\/\/numericalanalysisconference.org.uk\/\r\n% Biennial Conference on Numerical Analysis> in Scotland has a long history.\r\n% The first two meetings were held at the University of St. Andrews in\r\n% 1965 and 1967.  The meetings moved to the University of Dundee in 1969.\r\n% They were held there every two years, under the leadership of Ron Mitchell,\r\n% until his retirement in 2007.  They moved to the University of\r\n% Strathclyde in 2009, under the leadership of Alison Ramage.  This meeting\r\n% was the 26th and the 50th anniversary.  Although this is a premier conference\r\n% series in numerical analysis, it was the first time I participated.\r\n\r\n%% University of Strathclyde\r\n% The University of Strathclyde is a public research university located in\r\n% Glasgow, Scotland.  The university's name comes from the \"Valley of the River\r\n% Clyde\".  The school is Scotland's foremost institution of science and\r\n% engineering. \r\n\r\n%% The Conference\r\n% The conference was held June 23rd through 26th.  I was one of a dozen\r\n% invited speakers REPLACE_WITH_DASH_DASH five from the US, four from continental Europe,\r\n% and three from the UK.  \r\n\r\n%%\r\n% Almost all of the participants gave half-hour talks.  There were 168 talks\r\n% organized in seven parallel sessions.  Many of the talks were in one of\r\n% twelve minisymposia.  Sample minisymposia titles include \"Stable and\r\n% accurate discretisations for convection-dominated problems\", \"City\r\n% analytics\", \"Chebfun: new developments, cool applications and on the\r\n% horizon\", and \"Numerical linear algebra for optimsation and data\r\n% assimilation\".\r\n\r\n%% A. R. Mitchell Lecture\r\n% Two of the invited lectures are given special emphasis to honor\r\n% British numerical analysts.  The A. R. Mitchell Lecture honors\r\n% the University of Dundee's Ron Mitchell who was the dominant force in\r\n% this conference for forty years.\r\n\r\n%%\r\n% This year the lecture was given by Mike Giles of the University of Oxford\r\n% on \"Multilevel Monte Carlo Methods\".  I knew nothing about this subject\r\n% before his talk and I learned a great deal.  His abstract provides a link\r\n% to a web page featuring a survey paper and MATLAB codes,\r\n% <http:\/\/people.maths.ox.ac.uk\/gilesm\/acta\/\r\n% http:\/\/people.maths.ox.ac.uk\/gilesm\/acta>\r\n\r\n%% Fletcher-Powell Lecture\r\n% The Fletcher-Powell Lecture honors two mathematicians known for their\r\n% work in optimization algorithms, including the Davidon-Fletcher-Powell,\r\n% DFP, formula.  Mike Powell, from Cambridge University, had passed away\r\n% in April.  Roger Fletcher, from the University of Dundee, attended the\r\n% lecture.\r\n\r\n%%\r\n% This year the lecture was given by my long-time friend Mike Saunders\r\n% from Stanford.  He talked about \"Experiments with linear and nonlinear\r\n% optimization using Quad precision.\"  He has used the GFortran compiler\r\n% to build a version of his MINOS optimization software with the REAL*16\r\n% floating point datatype.  He tackled flux balance analysis models of\r\n% metabolic networks that involve coefficients ranging over 15 or 16 orders\r\n% of magnitude.  Ordinary double precision, that is Fortran's REAL*8,\r\n% cannot do the job.\r\n\r\n%% Sparse Days in St Girons III\r\n% The history of the <\r\n% Sparse Days in St Girons conference> is itself sparse.\r\n% There have been just two previous conferences, one in 1994 and one in 2003.\r\n% The conference is part of a program on \"High Performance Linear and\r\n% Nonlinear Methods for Large Scale Applications\" organized by CIMI, the\r\n% French Centre International de Mathematique et d'Informatique in Toulouse.\r\n\r\n%%\r\n% Organizers of Sparse Days in St Girons III included my buddies Iain Duff,\r\n% who splits his time between Rutherford Appleton Laboratory in the UK and the\r\n% Parallel Algorithms Group at CERFACS in Toulouse; Jack Dongarra, from the\r\n% University of Tennessee; and Jim Demmel, from the University of California,\r\n% Berkeley.\r\n\r\n%% St Girons\r\n% St Girons is a picturesque town in the French Pyrenees, not far\r\n% from the border with Spain.  The Tour de France sometimes passes through\r\n% here.  We met in the town's only cinema, which is air conditioned.\r\n% This turned out to be fortunate planning because a heat wave hit that\r\n% week, June 28th through July 2nd, with temperatures in the mid 30s C,\r\n% which is mid 90s F.  \r\n\r\n%% The Conference\r\n% The format for this conference was quite different from the one the\r\n% previous week in Scotland.  There were only about half as many attendees,\r\n% around 100.  There were no invited talks and no parallel sessions.  Anyone\r\n% who wanted to give a talk gave one.  There were 42 talks, most of them half\r\n% an hour.\r\n\r\n%% Interesting Talks\r\n% Tim Davis, MathWorks consultant on sparse matrices, who recently moved to\r\n% Texas A&M, talked about \"Sparse SVD, and a GPU-accelerated sparse QR\".\r\n\r\n%%\r\n% Bora Ucar of Ecoles Normales Superieures in Lyon talked about \"Two\r\n% approximation algorithms for bipartite matching on multicore architectures.\"\r\n% The analysis of his algorithm happened to involve my old friend, the\r\n% <https:\/\/blogs.mathworks.com\/cleve\/2013\/09\/02\/the-lambert-w-function\/\r\n% Lambert W function>.\r\n\r\n%%\r\n% Joost Rommes, from Mentor Graphics in the Netherlands, talked about\r\n% \"Challenges in numerical simulation of electrical networks.\"\r\n\r\n%%\r\n% Alex Pothen, from Purdue, talked about \"The Virtual Scalpel: real-time\r\n% matrix computations and finite elements for surgery.\"\r\n\r\n%%\r\n% Dan Sorensen, from Rice, talked about \"A DEIM induced CUR factorization.\"\r\n% This title means that C is some columns of a matrix A, R is some rows of A,\r\n% and DEIM is an algorithm that finds U so that the product CUR is a good low\r\n% approximation to A.\r\n\r\n%%\r\n% Martin Gander, from Universite de Geneve, talked about \"Five decades of\r\n% time parallel integration.\"  I was particularly interested in his talk \r\n% because I knew something about the subject many years ago, but was not\r\n% familiar with recent developments.  At first, it might seem impossible\r\n% to parallelize a computation that involves stepping along in time.  But\r\n% Martin gave an overview of four different approaches to how this can be done.\r\n\r\n%% Another Cleve\r\n% There is another Cleve in this business.  Cleve Ashcraft of LSTC, Livermore\r\n% Software Technology Corporation, is actually a \"grandstudent\".  His PhD\r\n% advisor at Yale was Stan Eisenstat, who was my PhD advisee when I was a\r\n% visiting professor at Stanford.  At Sparse Days, Cleve talked about\r\n% \"Separability, partitions and coverings.\".  His talk didn't mention the\r\n% applications he is deeply involved in at LSTC, which are the developers,\r\n% among other things, of the crash analysis software LS-DYNA.  If you want\r\n% to see a hint of some serious use of dynamic finite element calculations, \r\n% take at the <http:\/\/www.lstc.com\/ LSTC web pages>.\r\n\r\n%%\r\n% This is the group photo from Sparse Days in St Girons III.  There are two\r\n% Cleves.  Cleve Ashcraft is the guy with the green shirt and grey beard\r\n% in the front row.  Can you find the other one?\r\n% (Thanks to Pierre-Henri Cros for the photo, and for handling local\r\n% arrangements.) \r\n%\r\n% <<sparse_days.jpg>>\r\n%\r\n\r\n##### SOURCE END ##### 8c72cc5a5eef45b1843ceb74d036b5fd\r\n-->","protected":false},"excerpt":{"rendered":"<div class=\"overview-image\"><img decoding=\"async\"  class=\"img-responsive\" src=\"https:\/\/blogs.mathworks.com\/images\/cleve\/sparse_days.jpg\" onError=\"this.style.display ='none';\" \/><\/div><!--introduction--><p>I have just returned from two meetings in Europe, the 26th Biennial Conference on Numerical Analysis at the University of Strathclyde in Glasgow, Scotland, and Sparse Days III in Saint-Girons, France.... <a class=\"read-more\" href=\"https:\/\/blogs.mathworks.com\/cleve\/2015\/07\/13\/trip-report-naconf-2015-and-sparse-days-iii\/\">read more >><\/a><\/p>","protected":false},"author":78,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[4,16,8,10],"tags":[],"_links":{"self":[{"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/posts\/1229"}],"collection":[{"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/users\/78"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/comments?post=1229"}],"version-history":[{"count":1,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/posts\/1229\/revisions"}],"predecessor-version":[{"id":1230,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/posts\/1229\/revisions\/1230"}],"wp:attachment":[{"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/media?parent=1229"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/categories?post=1229"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/cleve\/wp-json\/wp\/v2\/tags?post=1229"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}