Multi-dimensional filtering: Reducing the dimension through rotation

Julia Docampo, Jennifer K. Ryan, Mahsa Mirzargar, Robert M. Kirby

Research output: Contribution to journalArticlepeer-review

8 Downloads (Pure)


Over the past few decades there has been a strong effort towards the development of Smoothness-Increasing Accuracy-Conserving (SIAC) filters for Discontinuous Galerkin (DG) methods, designed to increase the smoothness and improve the convergence rate of the DG solution through this post-processor. These advantages can be exploited during flow visualization, for example by applying the SIAC filter to the DG data before streamline computations [Steffan et al., IEEE-TVCG 14(3): 680-692]. However, introducing these filters in engineering applications can be challenging since a tensor product filter grows in support size as the field dimension increases, becoming computationally expensive. As an alternative, [Walfisch et al., JOMP 38(2);164-184] proposed a univariate filter implemented along the streamline curves. Until now, this technique remained a numerical experiment. In this paper we introduce the line SIAC filter and explore how the orientation, structure and filter size affect the order of accuracy and global errors. We present theoretical error estimates showing how line filtering preserves the properties of traditional tensor product filtering, including smoothness and improvement in the convergence rate. Furthermore, numerical experiments are included, exhibiting how these filters achieve the same accuracy at significantly lower computational costs, becoming an attractive tool for the scientific visualization
Original languageEnglish
Pages (from-to)A2179–A2200
JournalSIAM Journal on Scientific Computing
Issue number5
Publication statusPublished - 27 Sep 2017


  • Discontinuous Galerkin
  • post-processing
  • SIAC filtering
  • accuracy enhancement
  • error reduction

Cite this