Ifpack2: Jacobi tests fail with 2 OpenMP threads, pass with 1
Created by: mhoemmen
@trilinos/ifpack2 @trilinos/tpetra
I was warming up by working on #237 (closed) (easy fix), and found that some of Ifpack2's Jacobi tests fail with OpenMP, when OMP_NUM_THREADS=2, but pass when OMP_NUM_THREADS=1. The difference in iteration counts was significant enough to be a real bug. My fix for #237 (closed) didn't affect it.
I'm wondering if there's something wrong with diagonal extraction (Tpetra::CrsMatrix::getLocalDiagCopy), but I'm not sure.