For JACOBI preconditioner, this%diag_lapl should not be allocated to np but 1 in case of constant weights.