refactor solver with FFTs to store N 2D matrices
p-costa opened this issue · comments
The approach that is used to solve N 2D systems of equations when FFTs are employed (by modifying the diagonal using hypre's AddToBoxValues in a loop) does not seem to scale well with many cores.
Probably better, instead, initializing N 2D matrices at the beginning of the calculation.
prototype implemented in branch multi_block_fft_n_matrix_2d_rev