Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation

Allen J. Toreja, Rizwan Uddin

Research output: ResearchConference contribution

Abstract

An existing implementation of the nodal integral method for the time-dependent convection-diffusion equation is modified to incorporate various PETSc (Portable, Extensible Toolkit for Scientific Computation) solver and preconditioner routines. In the modified implementation, the default iterative Gauss-Seidel solver is replaced with one of the following PETSc iterative linear solver routines: Generalized Minimal Residuals, Stabilized Biconjugate Gradients, or Transpose-Free Quasi-Minimal Residuals. For each solver, a Jacobi or a Successive Over-Relaxation preconditioner is used. Two sample problems, one with a low Peclet number and one with a high Peclet number, are solved using the new implementation. In all the cases tested, the new implementation with the PETSc solver routines outperforms the original Gauss-Seidel implementation. Moreover, the PETSc Stabilized Biconjugate Gradients routine performs the best on the two sample problems leading to CPU times that are less than half the CPU times of the original implementation.

LanguageEnglish (US)
Title of host publicationInternational Conference on Nuclear Engineering, Proceedings, ICONE
Pages883-890
Number of pages8
Volume4
DOIs
StatePublished - 2002
Event10th International Conference on Nuclear Engineering (ICONE 10) - Arlington, VA, United States
Duration: Apr 14 2002Apr 18 2002

Other

Other10th International Conference on Nuclear Engineering (ICONE 10)
CountryUnited States
CityArlington, VA
Period4/14/024/18/02

Fingerprint

Program processors
Convection

ASJC Scopus subject areas

  • Engineering(all)
  • Energy(all)

Cite this

Toreja, A. J., & Uddin, R. (2002). Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation. In International Conference on Nuclear Engineering, Proceedings, ICONE (Vol. 4, pp. 883-890). DOI: 10.1115/ICONE10-22684

Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation. / Toreja, Allen J.; Uddin, Rizwan.

International Conference on Nuclear Engineering, Proceedings, ICONE. Vol. 4 2002. p. 883-890.

Research output: ResearchConference contribution

Toreja, AJ & Uddin, R 2002, Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation. in International Conference on Nuclear Engineering, Proceedings, ICONE. vol. 4, pp. 883-890, 10th International Conference on Nuclear Engineering (ICONE 10), Arlington, VA, United States, 4/14/02. DOI: 10.1115/ICONE10-22684
Toreja AJ, Uddin R. Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation. In International Conference on Nuclear Engineering, Proceedings, ICONE. Vol. 4. 2002. p. 883-890. Available from, DOI: 10.1115/ICONE10-22684
Toreja, Allen J. ; Uddin, Rizwan. / Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation. International Conference on Nuclear Engineering, Proceedings, ICONE. Vol. 4 2002. pp. 883-890
@inbook{4ed62be7facb4a64b7ad08988deb2a3d,
title = "Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation",
abstract = "An existing implementation of the nodal integral method for the time-dependent convection-diffusion equation is modified to incorporate various PETSc (Portable, Extensible Toolkit for Scientific Computation) solver and preconditioner routines. In the modified implementation, the default iterative Gauss-Seidel solver is replaced with one of the following PETSc iterative linear solver routines: Generalized Minimal Residuals, Stabilized Biconjugate Gradients, or Transpose-Free Quasi-Minimal Residuals. For each solver, a Jacobi or a Successive Over-Relaxation preconditioner is used. Two sample problems, one with a low Peclet number and one with a high Peclet number, are solved using the new implementation. In all the cases tested, the new implementation with the PETSc solver routines outperforms the original Gauss-Seidel implementation. Moreover, the PETSc Stabilized Biconjugate Gradients routine performs the best on the two sample problems leading to CPU times that are less than half the CPU times of the original implementation.",
author = "Toreja, {Allen J.} and Rizwan Uddin",
year = "2002",
doi = "10.1115/ICONE10-22684",
volume = "4",
pages = "883--890",
booktitle = "International Conference on Nuclear Engineering, Proceedings, ICONE",

}

TY - CHAP

T1 - Improving the efficiency of the nodal integral method with the portable, extensible toolkit for scientific computation

AU - Toreja,Allen J.

AU - Uddin,Rizwan

PY - 2002

Y1 - 2002

N2 - An existing implementation of the nodal integral method for the time-dependent convection-diffusion equation is modified to incorporate various PETSc (Portable, Extensible Toolkit for Scientific Computation) solver and preconditioner routines. In the modified implementation, the default iterative Gauss-Seidel solver is replaced with one of the following PETSc iterative linear solver routines: Generalized Minimal Residuals, Stabilized Biconjugate Gradients, or Transpose-Free Quasi-Minimal Residuals. For each solver, a Jacobi or a Successive Over-Relaxation preconditioner is used. Two sample problems, one with a low Peclet number and one with a high Peclet number, are solved using the new implementation. In all the cases tested, the new implementation with the PETSc solver routines outperforms the original Gauss-Seidel implementation. Moreover, the PETSc Stabilized Biconjugate Gradients routine performs the best on the two sample problems leading to CPU times that are less than half the CPU times of the original implementation.

AB - An existing implementation of the nodal integral method for the time-dependent convection-diffusion equation is modified to incorporate various PETSc (Portable, Extensible Toolkit for Scientific Computation) solver and preconditioner routines. In the modified implementation, the default iterative Gauss-Seidel solver is replaced with one of the following PETSc iterative linear solver routines: Generalized Minimal Residuals, Stabilized Biconjugate Gradients, or Transpose-Free Quasi-Minimal Residuals. For each solver, a Jacobi or a Successive Over-Relaxation preconditioner is used. Two sample problems, one with a low Peclet number and one with a high Peclet number, are solved using the new implementation. In all the cases tested, the new implementation with the PETSc solver routines outperforms the original Gauss-Seidel implementation. Moreover, the PETSc Stabilized Biconjugate Gradients routine performs the best on the two sample problems leading to CPU times that are less than half the CPU times of the original implementation.

UR - http://www.scopus.com/inward/record.url?scp=0036388103&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0036388103&partnerID=8YFLogxK

U2 - 10.1115/ICONE10-22684

DO - 10.1115/ICONE10-22684

M3 - Conference contribution

VL - 4

SP - 883

EP - 890

BT - International Conference on Nuclear Engineering, Proceedings, ICONE

ER -