LaplacePDD: public Laplacian¶
start(const FieldPerp &b, PDD_data &data)¶
Laplacian inversion using Parallel Diagonal Dominant (PDD) method
July 2008: Adapted from serial version to run in parallel (split in X) for tridiagonal system i.e. no 4th order inversion yet.
- This code stores intermediate results and takes significantly more memory than the serial version. This can be balanced against communication time i.e. faster communications can allow less memory use.
b: RHS values (Ax = b)
data: Internal data used for multiple calls in parallel mode
Data structure for PDD algorithm.