NVIDIA/AMGX

Right-hand-side of Poisson examples should be changed

开放

#30 创建于 2018年5月31日

 (4 条评论) (0 个反应) (0 位负责人)Cuda (186 个派生)auto 404
enhancementgood first issue

仓库指标

星标
 (679 个星标)
PR 合并指标
 (PR 指标待抓取)

描述

AMGX ships with many nice sample programs, including simple implementation of Poisson-like equations. Currently, the charge function (e.g. in A x = b, the array b) is set to 1.d0 everywhere.

This is a poor choice, since integrating a constant over larger and large volumes gives you an exponentially growing solution. As a consequence, the sample cases have terrible convergence rates and do not scale as the problem size is increased.

I suggest a better choice: set b = sin( 2 * pi * x / Lx ) * sin( 2 * pi * y / Ly ) * sin( 2 * pi * z / Lz ), which satisfies either periodic or Dirichlet boundary conditions. This is more representative of the use of a Poisson solver.

贡献者指南