NA Digest, V. 17, # 22

NA Digest Sunday, August 27, 2017 Volume 17 : Issue 22

Today's Editor:

Daniel M. Dunlavy
Sandia National Labs

Today's Topics: Subscribe, unsubscribe, change address, or for na-digest archives:

Submissions for NA Digest:

From: Field G. Van Zee
Date: August 26, 2017
Subject: Dense linear algebra software survey

On behalf of The Science of High-Performance Computing group at The
University of Texas at Austin, we are conducting a brief (5-10 minute)
survey. We are interested in learning about your experiences with
dense linear algebra software.

To participate in the survey, please click on the following link:

It will close at 23:00 UTC on September 11, 2017. Since we all tend
to procrastinate leading to things not getting done, please complete
this survey immediately!

We thank you for your time and look forward to your response. Your
opinion is important to us!

From: Weijian Zhang
Date: August 26, 2017
Subject: New Search Engine for Machine Learning Research

Hello, we are working on a new search engine, called Etymo
(, for data scientists to find interesting research
papers. It visualizes search results to reveal interesting structures
within the research area (for example, it automatically identifies
small clusters within a broad general search). We would love to hear
your thoughts and suggestions!

Here is a blog post:

From: JS Hysing
Date: August 20, 2017
Subject: FEATool v1.7 - Matlab FEM Toolbox with FEniCS Integration

Precise Simulation Ltd. developer of the Matlab and Octave Finite
Element Multiphysics Toolbox FEATool (, has
now released FEATool Version 1.7 which introduces the modeling and
simulation tool concept Multi-Simulation. Multi-Simulation allows for
seamless switching between both built-in and external solvers and
simulation tools. Starting with FEATool 1.7, Multi-Simulation supports
the the FeatFlow high performance computational fluid dynamics CFD and
FEniCS multiphysics solvers, as well as the visualization and
postprocessing tool Plotly. The multi-simulation functionality is
available directly in the FEATool Matlab GUI, and with m-file
scripting as all other FEATool functionality.

FEATool 1.7 also features significant updates and new physics modes
for porous media flows (Darcy's Law and Brinkman equations),
structural mechanics (Euler-Bernoulli Beam modeling), and
electromagnetics (electrostatics and magnetostatics). Moreover, all
applicable physics modes have been updated to support axisymmetric
cylindrical coordinate systems with corresponding GUI support.
Furthermore, FEATool now also supports solving different equations in
different coupled geometrical subdomains, for example allowing for
fluid-structure interaction simulations.

Please visit for more information.

From: Zdenek Dostal
Date: August 21, 2017
Subject: New Book, Scalable Algorithms for Contact Problems

Scalable Algorithms for Contact Problems
by Z. Dostal, T. Kozubek, M. Sadowska, V. Vondrak
Advances in Mechanics and Mathematics 36 / xix+340 pages /
Hardcover- eBook; 978-1-4939-6834-3 / Springer 2016

More information on

This book presents a comprehensive and self-contained treatment of
scalable algorithms for the solutions of multibody contact problems of
linear elasticity. The brand new feature of these algorithms is
theoretically supported numerical scalability and parallel scalability
demonstrated on problems discretized by billions of degrees of
freedom. The theory supports solving multibody frictionless contact
problems, contact problems with possibly orthotropic Tresca's
friction, and transient contact problems. It covers BEM
discretization, jumping coefficients, floating bodies, mortar
non-penetration conditions, etc. The exposition is divided into four
parts, the first of which reviews appropriate facets of linear
algebra, optimization, and analysis. The most important algorithms and
optimality results are presented in the third part of the volume. The
presentation is complete, including continuous formulation,
discretization, decomposition, optimality results, and numerical
experiments. The final part includes extensions to contact shape
optimization, plasticity, and HPC implementation.

From: Patrick Farrell
Date: August 24, 2017
Subject: FEniCS Conference 2018, University of Oxford, Mar 2018

We are pleased to invite you to attend and/or submit an abstract to
the FEniCS Conference 2018 ( The
conference will take place at the University of Oxford on 21-23 March

The FEniCS Project ( is a popular computing
platform for solving partial differential equations. FEniCS enables
users to quickly translate scientific models into efficient finite
element code. The FEniCS'18 conference is an opportunity for all those
interested in the FEniCS Project and related projects to exchange
ideas, communicate their results and interact with the automated
scientific computing community.

We welcome developers, existing and potential users of the FEniCS
ecosystem as well as mathematicians, computer scientists and
application domain specialists interested in numerical methods, their
implementation and applications.

The FEniCS'18 conference will emphasise an open and inclusive
atmosphere, contributed talks and posters from a diverse range of
scientific areas, and dedicated time for discussions and coding.

For further information, please see:

From: Ralf Werner
Date: August 21, 2017
Subject: SIGOPT 2018, Germany, Mar 2018

We would like to draw your attention to the upcoming SIGOPT
conference, SIGOPT 2018 International Conference on Optimization, to
be held on March 21-23, 2018, in Kloster Irsee, Germany. You can find
a preliminary web page with important dates and additional information at

The SIGOPT conference series addresses the full range of theory and
practice in all branches of mathematical optimization, e.g., Linear
and Nonlinear Programming; Discrete and Combinatorial Optimization;
Optimal Control; Stochastic Programming; Multiobjective Programming;
Robust Optimization; Game Theory and Bilevel Optimization; and
Application of Optimization in Real-World Problems.

Young researchers from all areas of optimization are especially
encouraged to participate and to present their work. The conference
will provide a good opportunity for exchanging ideas, especially
across disciplines, and for establishing contacts with colleagues in
the friendly atmosphere of Kloster Irsee.

Plenary Speakers:
- Immanuel Bomze (Universitaet Wien, Austria)
- Jacek Gondzio (The University of Edinburgh, Scotland, UK)
- Britta Peis (RWTH Aachen, Germany)
- Karl Schmedders (Universitaet Zuerich, Switzerland)
- Rico Zenklusen (ETH Zuerich, Switzerland)

From: Rajesh Kumar Sharma
Date: August 23, 2017
Subject: Frontiers in Industrial and Applied Mathematics, India, Apr 2018

It gives us immense pleasure to inform you that Department of
Mathematics, NIT Hamirpur, Hamirpur, Himachal Pradesh, India is
organizing an International conference on Frontiers in Industrial and
Applied Mathematics (FIAM-2018) in association with The Indian Society
of Industrial and Applied Mathematics. Enclosed please find the
conference brochure.
Dates: April 26-27, 2018
Venue: NIT Hamirpur, Hamirpur, Himachal Pradesh, India

You are requested to submit an extended abstract of scholarly
unpublished and original work as per format and instruction given at
the link

Important Deadlines:
Abstract Submission: October 30, 2017
Acceptance of Abstract: November 15, 2017
Full length paper Submission: December 15, 2017
Acceptance of Paper: January 30, 2018

The accepted full length papers will be published in one of the
following publication categories: Category 1: AIP (American Institute
of Physics) Conference Proceedings indexed in Scopus, ISI Web of
Knowledge, MathSciNet, Scirus, Google Scholar, etc and in the past had
been included in SCIE; Category 2: IJETAE Journal indexed in UGC
approved Journal; Category 3: IJSER Journal indexed in UGC approved
Journal, Thomson Reuters etc.

For more information, you may visit the conference website

From: J.J.W. van der Vegt
Date: August 21, 2017
Subject: Assistant Professor Position, Numerical Analysis and Scientific Computing

The University of Twente (UT) in Enschede, The Netherlands, is seeking
researchers with an outstanding record in Mathematics and its
Applications in Science and Engineering for a Faculty position in
Numerical Analysis and Scientific Computing within the Department of
Applied Mathematics.

We are looking for a highly talented research scientist, whose focus
is on the development, analysis and application of finite elements
methods for partial differential equations, which may also include
inverse problems, optimization and optimal control with partial
differential equations. The successful candidate is expected to have
strong affinity with numerical analysis and scientific computing. In
addition, the capability to develop collaborations between mathematics
and the multi-disciplinary research context at the University of
Twente is important.

The position will be at the assistant professor level. Initially, the
contract will be for one year. After a positive evaluation the
contract will be for an indefinite period. Criteria for promotion to
the Associate Professor level comprise a strong track record of
publications, evaluation of teaching capacities and a demonstrated
ability to acquire research funding and to achieve an international
network and reputation.

Applications (including curriculum vitae, list of publications,
research and teaching statement, and three names of potential
references) can be submitted online at
before September 14, 2017 to the University of Twente,
attn. Prof.dr. P.M.G. Apers, Dean of the Faculty Electrical
Engineering, Mathematics and Computer Science. The interviews with
selected candidates are scheduled on September 27, 2017.

For more information please visit or
contact J.J.W. van der Vegt (phone +31 534895628, email

From: Professor Barry Hughes
Date: August 23, 2017
Subject: Lecturer Position, Comp& Appl Math, Univ of Melbourne

From: Jichao Zhao
Date: August 26, 2017
Subject: Postdoc Position, Univ of Auckland, New Zealand

Improving our understanding of human atrial ionic mechanisms of

About the project: Atrial fibrillation (AF), characterized by rapid or
irregular electrical activity in the upper chambers of the heart, is
the most common heart rhythm disturbance. AF is associated with
substantial morbidity and mortality. Anti-arrhythmic drug therapy is
a frontline treatment for AF by targeting specific ionic channels with
variable success rates due to our incomplete understanding of the
basic mechanisms of atrial cellular behaviours in patients under
normal and diseased conditions. Computer activation models provide a
powerful framework for investigating exact mechanisms behind AF and
testing effectiveness of proposed novel treatment strategies by
providing a flexible way to dissect highly intertwined contributing
factors and set up series of control studies.

The potential projects for the Postdoc are
1) To develop biophysics-based human atrial cellular models to
understand the pathophysiology of AF associated with gene mutations
and to identify a potential effective, targeted treatment;
2) To develop biophysics-based human sino-atrial node (SAN) models to
determine the contribution of heart-failure induced structural and
molecular substrates to SAN dysfunction, and potential targeted
3) To develop biophysics-based human atrial cellular models to
understand the precise ionic mechanism of new-onset AF under
diabetic conditions and to identify a potential upstream therapy.

To find out more about this project please contact Dr Jichao Zhao
(, Auckland Bioengieering Institute, The
University of Auckland, New Zealand.

From: Lieven De Lathauwer
Date: August 21, 2017
Subject: PhD or Postdoc Position, Tensor Methods, KULeuven, Belgium

Topic: Higher-order tensors become increasingly important in all
branches of mathematical engineering. We study algebraic aspects such
as the uniqueness of tensor decompositions and develop numerical
algorithms for their computation. We develop advanced tensor-based
methods for signal processing, data analysis, compressed sensing and
machine learning. Code is made available in Tensorlab
( We are looking for new collaborators to conduct
research on relevant topics.

Profile: We are seeking outstanding candidates who have obtained a
master's or doctoral degree in mathematical engineering, computer
science, electrical engineering, mathematics or physics, with a strong
interest in the mathematical aspects of signal processing and a good
knowledge of linear algebra. Experience with numerical mathematics,
optimization, signal processing and/or machine learning is a
plus. Good implementation skills are valued as well. Knowledge of
tensor techniques is not required.

Successful applicants will join a well-equipped and experienced
multidisciplinary research group at a top university. The position is
initially for one year and can be extended. Candidates are requested
to send an application letter and detailed CV to Prof. Lieven De
Lathauwer (URL:

From: Charis Edworthy
Date: August 21, 2017
Subject: Contents, IMA Journal of Numerical Analysis, 37 (3)

IMA Journal of Numerical Analysis
Links to all articles in this issue are available online at:

Banded, stable, skew-symmetric differentiation matrices of high order,
Ernst Hairer, Arieh Iserles,

Collocation methods for third-kind VIEs, Sonia Seyed Allaei, Zhan-Wen
Yang, Hermann Brunner,

Refinable spaces and local approximation estimates for hierarchical
splines, Annalisa Buffa, Eduardo M. Garau,

Fast and accurate computation of Chebyshev coefficients in the complex
plane, Haiyong Wang, Daan Huybrechs,

Electromagnetic wave scattering by random surfaces: uncertainty
quantification via sparse tensor boundary elements, Carlos
Jerez-Hanckes, Christoph Schwab,

Existence of H-matrix approximants to the inverse of BEM matrices: the
hyper-singular integral operator, Markus Faustmann, Jens Markus
Melenk, Dirk Praetorius,

Numerical approximation of fractional powers of regularly accretive
operators, Andrea Bonito, Joseph E. Pasciak,

A stabilized finite element method for the Darcy problem on surfaces,
Peter Hansbo, Mats G. Larson,

L2(H\gamma1) Finite Element Convergence for Degenerate Isotropic
Hamilton-Jacobi-Bellman Equations, Max Jensen,

Conforming and nonconforming virtual element methods for elliptic
problems, Andrea Cangiani, Gianmarco Manzini, Oliver J. Sutton,

On the numerical stability of the second barycentric formula for
trigonometric interpolation in shifted equispaced points, Anthony
P. Austin, Kuan Xu,

Approximation of the invariant law of SPDEs: error analysis using a
Poisson equation for a full-discretization scheme, Charles-Edouard
Brehier, Marie Kopec,

Spectral Galerkin methods for a weakly singular Volterra integral
equation of the second kind, Can Huang, Martin Stynes,

Analysis of a Full Space-Time Discretization of the Navier-Stokes
Equations by a Local Projection Stabilization Method, Naveed Ahmed,
Tomas Chacon Rebollo, Volker John, Samuele Rubino,

Gauss quadrature for quasi-definite linear functionals, Stefano Pozza,
Miroslav S. Pranic, Zdenek Strakos,

Additive domain decomposition operator splittings-convergence analyses
in a dissipative framework, Eskil Hansen, Erik Henningsson,

A convergent adaptive finite element method for electrical impedance
Bangti Jin, Yifeng Xu, Jun Zou,

Gradient discretization of hybrid-dimensional Darcy flow in fractured
porous media with discontinuous pressures at matrix-fracture
interfaces, K. Brenner, J. Hennicker, R. Masson, P. Samier,

Convergence rates of finite difference schemes for the linear
advection and wave equation with rough coefficient, Franziska Weber,

From: Chi-Wang Shu
Date: August 24, 2017
Subject: Contents, Journal of Scientific Computing, 72 (3)

Journal of Scientific Computing
Volume 72, Number 3, September 2017

Superconvergence of Finite Element Approximations for the Fractional
Diffusion-Wave Equation, Jincheng Ren, Xiaonian Long, Shipeng Mao and
Jiwei Zhang

An Adaptive Finite Element Method for the Wave Scattering with
Transparent Boundary Condition, Xue Jiang, Peijun Li, Junliang Lv and
Weiying Zheng

Fast Iterative Method with a Second-Order Implicit Difference Scheme
for Time-Space Fractional Convection-Diffusion Equation, Xian-Ming Gu,
Ting-Zhu Huang, Cui-Cui Ji, Bruno Carpentieri and Anatoly A. Alikhanov

Adaptive WENO Methods Based on Radial Basis Function Reconstruction,
Caterina Bigoni and Jan S. Hesthaven

A Numerical Strategy for Freestream Preservation of the High Order
Weighted Essentially Non-oscillatory Schemes on Stationary Curvilinear
Grids, Yujie Zhu, Zhensheng Sun, Yuxin Ren, Yu Hu and Shiying Zhang

A New Mixed Finite Element Method for Elastodynamics with Weak
Symmetry, Carlos Garcia, Gabriel N. Gatica and Salim Meddahi

Regularization of Singularities in the Weighted Summation of
Dirac-Delta Functions for the Spectral Solution of Hyperbolic
Conservation Laws, Jean-Piero Suarez and Gustaaf B. Jacobs

Unconditional Superconvergence Analysis of a Crank-Nicolson Galerkin
FEM for Nonlinear Schrodinger Equation, Dongyang Shi and Junjun Wang

A Parallel Finite Element Method for 3D Two-Phase Moving Contact Line
Problems in Complex Domains, Li Luo, Qian Zhang, Xiao-Ping Wang and
Xiao-Chuan Cai

A Non-oscillatory Multi-Moment Finite Volume Scheme with Boundary
Gradient Switching, Xi Deng, Ziyao Sun, Bin Xie, Kensuke Yokoi,
Chungang Chen and Feng Xiao

High-Order Accurate Adaptive Kernel Compression Time-Stepping Schemes
for Fractional Differential Equations, Daniel Baffet and Jan
S. Hesthaven

Relaxing the CFL Condition for the Wave Equation on Adaptive Meshes,
Daniel Peterseim and Mira Schedensack

Numerical Analysis of Fully Discretized Crank-Nicolson Scheme for
Fractional-in-Space Allen-Cahn Equations, Tianliang Hou, Tao Tang and
Jiang Yang

Construction of a p-Adaptive Continuous Residual Distribution Scheme,
R. Abgrall, Q. Viville, H. Beaugendre and C. Dobrzynski

Linearized Conservative Finite Element Methods for the
Nernst-Planck-Poisson Equations, Huadong Gao and Dongdong He

Improved Approximation of Phase-Space Densities on Triangulated
Domains Using Discrete Flow Mapping with p-Refinement, Janis Bajars,
David J. Chappell, Timo Hartmann and Gregor Tanner

A Three-Stage Approach for Segmenting Degraded Color Images:
Smoothing, Lifting and Thresholding (SLaT), Xiaohao Cai, Raymond Chan,
Mila Nikolova and Tieyong Zeng

Segmentation of Three-Dimensional Images with Parametric Active
Surfaces and Topology Changes, Heike Benninghoff and Harald Garcke

From: Claude Brezinski
Date: August 27, 2017
Subject: Contents, Numerical Algorithms, 76 (1)

Table of Contents
Numerical Algorithms, Vol. 76, No. 1

Ryo Kobayashi, Takuma Kimura, Shinichi Oishi, A method for verifying
the accuracy of numerical solutions of symmetric saddle point linear

Rachid Ait-Haddou, q-Blossoming and Hermite-Pade approximants to the
q-exponential function

Ky Anh Pham, Ngoc Hai Trinh, Splitting extragradient-like algorithms
for strongly pseudomonotone equilibrium problems

Shi-Liang Wu, Cui-Xia Li, Modified complex-symmetric and
skew-Hermitian splitting iteration method for a class of
complex-symmetric indefinite linear systems

Sascha Woerz, Heinz Bernhardt, A global convergent derivative-free
method for solving a system of non-linear equations

Lukasz Plociniczak, Szymon Sobieszek, Numerical schemes for
integro-differential equations with Erdelyi-Kober fractional operator

Robert A. Van Gorder, On the utility of the homotopy analysis method
for non-analytic and global solutions to nonlinear differential

Ali Baharev, Ferenc Domes, Arnold Neumaier, A robust approach for
finding all well-separated solutions of sparse systems of nonlinear

Benjamin Kehlet, Anders Logg, A posteriori error analysis of round-off
errors in the numerical solution of ordinary differential equations

Michael A. Jandron, Anthony A. Ruffa, James Baglama, An asynchronous
direct solver for banded linear systems

Jiyong Li, Trigonometrically fitted multi-step Runge-Kutta methods for
solving oscillatory initial value problems

Yekini Shehu, Olaniyi S. Iyiola, Strong convergence result for
monotone variational inequalities

Mona Narang, Saurabh Bhatia, Vinay Kanwar, New efficient derivative
free family of seventh-order methods for solving systems of nonlinear

End of Digest