#include <Eigen/src/IterativeLinearSolvers/IncompleteLUT.h>
keeps off-diagonal entries; drops diagonal entries