Author Archives: rjl

Latest on Elsevier boycott

To date more than 11,000 people have signed the boycott at thecostofknowledge.com, including more than 2000 mathematicians.  Please consider joining us!

Elsevier recently sent two letters to the mathematics community, responding to the boycott and listing various steps they are taking.  However, for the most part these are not as significant as they make them out to be.  See this blog post of Tim Gowers commenting in detail on the letters.

An excellent article “Mathematicians Take a Stand” was recently
posted on the arXiv by Doug Arnold and Henry Cohn, discussing in more
detail the reasons for the boycott and for targeting Elsevier in particular.  This article will appear in the Notices of the American Mathematical
Society in the June/July issue.

See also the Statement of Purpose  for the boycott, which was signed by 38 mathematicians and gives additional explanation.  Many other links can be found on the PolyMath Journal Publishing Reform page.

At least one major university has responded to the boycott by cancelling all Elsevier subscriptions: see this announcement from TU Munich and Gowers’ comment.

Paper posted on Tohoku tsunami source comparison

I spent the last week of March visiting the Institute of Seismology and Volcanology at the University of Hokkaido in Sapporo, finishing off a paper that was submitted to a special issue of the Bulletin of the Seismological Society of America (BSSA) on the Tohoku-oki tsunami of 11 March 2011. It is now posted here.

This work started last summer when Bre MacInnes was still at UW working with me on an NSF RAPID grant to perform simulations of this tsunami using our GeoClaw code. She’s now a postdoc in Sapporo and together with Adit Gusman (another postdoc there) and Y. Tanioka (director of the Institute), we compared 10 proposed sources for the tsunami by simulating the resulting wave using GeoClaw and comparing with observations. Each source is a model of the seafloor deformation that caused the tsunami. Dozens of similar (but far from identical) source models have been developed by different research groups, by performing some form of inversion based on some combination of seismic, tsunami, and GPS data. Some sources model the slip along the fault plane, 10s of kilometers below the seafloor, which we converted to seafloor deformation using the Okada model (based on a Green’s function solution to the equations of elastic deformation in a half space). Other groups directly modeled the seafloor deformation based primarily on tsunami data.

For each earthquake source, we compared results at 4 DART buoys (pressure gauges on the ocean floor) near Japan and in four locations along the coast, including the Sendai Plain and 3 communities further north on the Sanriku Coast. We found that the sources developed using tsunami data typically performed the best for tsunami modeling, perhaps not surprisingly. Many of the simulated results agreed quite well with observations. However, none of the sources gave good agreement everywhere and it remains to sort out what source is best for tsunami modeling. Of course the numerical model has its own limitations, in particular that the shallow water equations are solved, a good approximation in many cases but not perfect by any means and it would be interesting to study some of the differences between observation and simulation in more detail.

For some information on how the GeoClaw code has been validated on a number of test problems, and other papers on tsunami modeling and the software more generally, see the GeoClaw webpage.

In addition to the figures in the paper, we also created an electronic supplement of additional figures along with most of the raw data used in the simulations and comparisons. Until a few years ago, I had never worked with data from the real world (like most mathematicians) and it has been an interesting experience learning to work with large data sets, and thinking about the best way to archive and curate this data. Reproducibility in computational science research is one of my big interests these days, and this paper was a good case study for this. The journal BSSA encourages the submission of electronic supplements, and I found that the webpage format they provided was easy to work with. We’re also collecting all the simulation and data analysis codes in a Github repository.

Visit to Vietnam

I spent three weeks in March in Vietnam, visiting several universities and getting to know a number of mathematicians. The first week was spent at an International Conference on High Performance Scientific Computing in Hanoi. This was the sixth conference in a series that started in 2000 and has been held every three years since. Organized in collaboration with the University of Heidelberg, it attracted many researchers from Germany, but also from many other countries and featured a number of excellent talks by local researchers as well as many high profile invited speakers. I spoke in a minisymposium organized by Rolf Jeltsch, an old friend and former colleague at ETH-Zurich.

Ha Long Bay

The weekend after the conference there was an organized weekend trip to Ha Long Bay, a 3-hour bus ride from Hanoi. The bay is full of magnificent karst formations and the islands have extensive limestone caves that are also amazing.

 

With Hieu (red jacket) and his class

The second week I spent at Vietnam National University in Hanoi where, in addition to a seminar talk on tsunami modeling, I also gave about 13 hours of guest lectures in a course on Finite Difference Methods for Ordinary and Partial Differential being taught by Nguyen Trung Hieu. Somehow the students survived and were even smiling at the end.

The University of Washington has close ties with this university thanks to a grant from the Vietnamese government supporting interchanges, and I’m the fourth member of my department to visit the university in Hanoi. About 10 members of the Mathematics faculty in Hanoi have also visited UW for several months at a time, to sit in on classes and engage in research. In fact two weeks after leaving Hanoi I ran into two professors and a graduate student from this department on my flight from Seoul to Seattle.

After two weeks in Hanoi and a brief visit to Hue, I spent another week in Ho Chi Minh City and gave a talk at the University of Science, another branch of the extensive Vietnam National University. I also spent some time with Mai Duc Thanh from the International University, a student of Philippe LeFloch’s who works on hyperbolic problems, including shallow water equations over bathymetry.

Thanh

In the photo he’s sitting by the Saigon River, which we were watching flow with a strong current upstream due to the tide, even 60 km from its outlet in the South China Sea. The entire Mekong Delta region is so flat that flooding is a major problem and storm surge modeling is of great interest. Tsunamis are also a potential hazard, although there hasn’t been a major one recently on this coast.

Reproduciblity: Recent 60 Minutes segment on the Duke trials

A segment on 60 Minutes on February 12, 2012 concerns the cancer clinical trials at Duke University that were recently found to be based on manipulated  and/or badly analyzed data.  This case is frequently evoked as evidence for the need for higher standards for data sharing and reproducibility in scientific work.

For more details of the statistical work uncovering the problems, see this talk by Keith Baggerly: The Importance of Reproducible Research in High-Throughput Biology: Case Studies in Forensic Bioinformatics

For another video on the importance of reproducibility and the Duke trials, see the talk Freedom (to reproduce) by John Wilbanks, who runs the Science Commons project at Creative Commons.  This was the lead-off talk at the workshop Reproducible Research: Tools and Strategies for Scientific Computing that I organized with Victoria Stodden and Ian Mitchell last summer, see this link for all the talks, abstracts, and slides.  There are several good ones here!

Update 29 Feb 2012: There are some interesting discussions of the Duke case and the 60 Minutes segment at http://groups.google.com/group/reproducible-research.

[HPC]^3 workshop summary

The 2012 [HPC]^3 workshop at KAUST just ended.  It’s been a busy week with lots of exciting developments — this was a workshop with an emphasis on work.  Each morning there were a few talks, but the afternoons were devoted to working in groups developing new software capabilities, much of it related to Clawpack and in particular the PyClaw suite of software.  This was the second such workshop at KAUST and some of the projects were continuations of things started at the 2011 workshop.

The title [HPC]^3 refers to 3 different interpretations of the acronym HPC.  The full title is High Performance Computing and  Hybrid Programming Concepts for Hyperbolic PDE Codes.

The workshop program is on the web and eventually videos and slides from the talks will also appear, so I won’t say too much about these here.  Instead I’ll mention a few of the highlights and accomplishments of the working groups.  Eventually a final report from each group will be posted.  For now you can find out more about what went on in some groups on the workshop wiki.

Some highlights

  • The  AMR group (Donna Calhoun, Tobias Weinzierl, Carsten Burstedde, Kristof Unterweger, Amal Alghamdi, Qi Tang, Kyle Mandli, and Marsha Berger via Skype) made progress on two different approaches to incorporating tree-based adaptive mesh refinement into PyClaw, one based on p4est and the other on Peano.
  • The Manycore group (George Turkiyyah, Nathan Bell, Andy Terell, Garune Ohannessian, Rio Yokota) made it possible to call CUDAClaw from Python and started re-engineering it to simplify other manycore implementations,  including OpenCL, TBB, and ISPC.
  • The implicit group (Jed Brown, Matteo Parsani, Lulu Liu) did some interesting work on using preconditioners for hyperbolic systems that equalize the wave speeds and on an implementation of downwind WENO for Runge-Kutta stages with negative coefficients.
  • The Discontinuous Galerkin group (James Rossmanith and Scott Moe) wrapped DoGPack in Python and got it working with VisClaw. They started planning how to more fully incorporate DG into a DoGClaw code.
  • Matt Emmett started working on applying WENO on mapped grids, and has a working code in 1d.
  • The geosciences group (Chris Kees, Marc Hesse, Robert Weiss) worked on Clawpack implementations of two problems: two-phase flow in porous media with an interesting nonconvex and multimodal flux function, and sediment transport with an eye towards tsunami deposits.  The final presentation also had a nice demo of the new IPython Notebook.
  • The visualization and steering group (Madhu Srinivasan, Chris Knox, Atanas Atanasov, Bruce D’Amora) made progress on several fronts: improvement in the HDF5 PyClaw output routines and conversion of this output to xdmf form, starting to add VisIt and Paraview capabilities into VisClaw using this form,  and coupling on-the-fly visualization into PyClaw so that the solution can be plotted during computation rather than as postprocessing. The ultimate goal is to use this to help steer the computation.
  • A documentation sprint one evening got us started working on using Sphinx to better document each example and produce a better set of galleries of sample results.  Yiannis Hadjimichael and Amal Alghamdi fixed many documentation pages to reflect recent changes to the code, and got the doctests working.

Thanks to everyone who participated and worked so hard, often late into the night.  Thanks to my co-organizers and the other plenary speakers who moved between groups and helped out with many issues.  And a special thanks to [David K]^2 (Ketcheson and Keyes) for helping create the wonderful environment at KAUST for computational science and multi-cultural interaction (including elliptic/hyperbolic and FEM/FV!), and for the financial support of the workshop and participants.

 

A more formal statement about mathematical publishing

Recently Tim Gowers brought together a group of mathematicians to discuss the Elsevier boycott he started. The group prepared a statement that explains some of the background to, and reasons for, the boycott, along with a more general discussion of mathematical publishing.

This statement has just been published in  a new post on Gowers’s Weblog.

We plan to advertise this statement widely among science writers and journalists and hope that the story will be picked up to a greater extent in the media.  We also hope that word will spread among researchers, referees, and editors in other fields.  Please help spread the word!

Update on 13 February 2012:

This statement and the boycott movement have attracted quite a bit of media attention, for example an article in the New York Times today.

A long list of media articles and editorials has been compiled by Terry Tao.

Boycotting Elsevier

 

I’m signing the pledge at http://thecostofknowledge.com/ and thought I’d say a bit here about why I’m doing so.

In case you haven’t heard, this was set up in response to a  blog post of Tim Gowers proposing such a site, and has gained several thousand signatories in the first few weeks. See also his follow-up post and the  PolyMath journal publishing reform page, where a large number of links are posted, many to articles, editorials and blog posts that have rapidly appeared in the two weeks since Gowers’s original post.

I have been mostly boycotting Elsevier for several years. In particular I have not submitted a paper to the  Journal of Computational Physics since they were acquired by Elsevier, and rarely referee for this journal any more, in spite of the fact that I was once on the editorial board and was a big fan of the journal. But Elsevier has had a reputation for decades for being the worst sort of the commercial publisher from the viewpoint of scientists. Of course some point out they are a business and they can run it as they see fit, but we are equally free to take our business elsewhere. Except that we’re not quite, which is part of the problem. They now bundle many journals in such a way that libraries cannot buy the ones essential to research without buying a bunch of other low-quality and/or un-needed journals at inflated prices. Elsevier made  profits of more than a billion dollars last year. Meanwhile, the University of Washington libraries (a major academic library system) had to cancel several thousand journal subscriptions in the past few years due to severe budget cuts. The requirement to spend huge sums on journals that nobody on campus needs or wants is impairing our ability to access the journals we do need, or to use the funds for other worthy research or educational purposes.

Here are some of my thoughts when deciding to sign.  For many more discussions (and opposing viewpoints) you should see the links above.

  • Elsevier has a reputation for publishing more bogus papers (and entire journals) than any other publisher I know.  For a discussion of some specific cases, see this post by Doug Arnold explaining why he is signing on to the boycott.  Doug looked into these issues extensively as part of a joint working group on rankings of mathematical journals for the IMU and ICIAM, and while investigating plagiarism as President of SIAM (see also his article on this topic in SIAM News).
  • Commercial publishers add relatively little value to articles these days, beyond giving a stamp of approval via the prestige of the journal and the hard work of the editorial board and referees (who are generally unpaid volunteers). Once upon a time publishers were necessary for typesetting and copy editing, but these days every mathematician I know does all their own typesetting using LaTeX. This is hardly new; I started using TeX at Stanford in 1977 and within a few years it was available everywhere. Some commercial publishers (particularly Elsevier, I’ve heard) don’t even provide much copy editing these days.
  • Rather than volunteering our time to write papers, typeset them, referee them, edit them, etc., only to be forced to buy them back at high prices from companies that add little value, we should be smart enough to figure out better ways. Until we do, commercial publishers have a stranglehold on libraries. People have been bemoaning this for years and I’m happy to see some action being taken that might get their attention.
  • Elsevier is probably not the only commercial publisher who should be targeted. My last blog post was a draft of an article I’m writing for the Encyclopedia of Applied and Computational Mathematics to be published by Springer. When I searched for it to add a link I was rather shocked to see the list price will be “approximately $5000″. As far as I know authors are being paid nothing (beyond a 33% discount on Springer books and online access to the encyclopedia, which I may have anyway if our library feels compelled to buy a copy). However, as Gowers says in his original blog post, we have to start somewhere and Elsevier is the obvious first choice for many people I know. For the time being some researchers must rely on journals from commercial publishers and we can’t expect people to boycott all publishers.  But I do try to favor non-profit journals in general, particularly when prioritizing accepting the many invitations I receive to referee papers.  I often point editors to my webpage on refereeing.
  • I’m not opposed to selling journals at reasonable prices. Many professional societies publish journals and make money on them, including my main society SIAM (full disclosure: I’m chair of the SIAM Journals Subcommittee). But there’s an order of magnitude difference in the prices being charged, and SIAM does copy editing and still makes enough of a profit on the journals that they help support many other aspects of the society. This is another advantage of nonprofit society journals: whatever profits there are go back to the members one way or another.
  • I would not necessarily encourage young (un-tenured) people to sign this pledge, particularly in fields where Elsevier journals are prevalent. And I have some qualms about doing so myself since I often write papers with students and younger colleagues whose careers depend on publishing in the recognized top journals. But for the most part I’ve managed to avoid publishing in Elsevier journals for several years now and haven’t missed it. One exception was a recent paper for a special issue of Advances in Water Resources on software, but in the future I can find plenty of other alternative venues for publishing in my field. (On the flip side, some comments on Gowers blog suggest that soon it may be seen as a negative for young people to publish with Elsevier. Probably not soon, but perhaps some day.)
  • I’m often surprised that so many mathematicians and scientists are oblivious to the issue of journal pricing and I hope this campaign will raise awareness if nothing else.  Personally I was introduced to the issue very early in my career.  My father was Executive Director of the American Mathematical Society in the 1980′s when the Society was threatened with a lawsuit by Gordon & Breach for publishing a survey of journal pricing and value.  Other societies (and authors) were sued.  The suits went nowhere and I understand that the name of Gordon & Breach was so tarnished in academic circles that it eventually disappeared (acquired by Taylor & Francis).  Elsevier take note!

Draft of EACM Article on Reproducibility Methods

I have prepared a draft of an article on reproducible research methods that I was invited to contribute to the Encyclopedia of Applied and Computational Mathematics to be published by Springer. [Update: I just checked with Springer and the $5000 price quoted is not a typo.  Needless to say I'm not happy about this!  I'm looking into it further.] [Further update on 29 Feb 2012: I've been told that the price on that webpage was a mistake and it's a mystery how it got in the database.  That was last week and it's still there...]

I welcome comments on this draft, either sent by email or posted below. Note that this article was greatly constrained in size (it may be too long already) and so I was very selective in what to include, but if you think I’ve neglected something important I will consider it.

I have avoided mentioning specific tools (other than some version control
systems) since there are a so many out there it would be hard to choose
which to point out, and hard to know which will stand the test of time.

I have included only two references, selected primarily as sources of other links.

New computational science stackexchange question and answer site

David Ketcheson just sent this mail to the  SIAM Computational Science and Engineering Activity Group (SIAG/CSE) mailing list.  I saw an earlier post of his on NA-Digest, but now that I have some time I’ve finally taken a look at it.  I think it will be a valuable forum for discussing questions of interest in our community!

David’s message:

There is a new question and answer site for computational science that uses the same engine as the well-known stackoverflow site for programming.  Take a look at
Some examples of the kind of questions being asked and answered there are:

Gene Golub SIAM Summer School on Simulation and Supercomputing in the Geosciences

This summer I’ll be spending two weeks in Monterey teaching part of the 2012
Gene Golub SIAM Summer School on Simulation and Supercomputing in the Geosciences, taking place July 29 to August 10.  The primary focus will be on techniques for simulating tsunami propagation and inundation, storm surge, and related numerical analysis and high-performance computing topics. Other organizers/lecturers are Michael, Bader, Joern Behrens, and Frank Giraldo.

This summer school is open to graduate students from all countries, and we encourage applicants from under-represented groups.  At least partial funding is available for all participants.

We’ll be selecting about 50 participants. The deadline to apply is February 1, 2012. For more details and instructions on how to apply, see the website http://www.mac.tum.de/g2s3.

This will be the third annual Gene Golub Summer School, funded in part by an endowment from the estate of Gene Golub.  This series is a very appropriate way to honor Gene’s memory, since he always cared greatly about students, not only in his group at Stanford but around the world.   Mike Heath, a fellow former student in the numerical analysis group at Stanford, wrote a nice article remembering Gene in SIAM News.