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.

貢獻者指南