Update example/linalg/example_solve_custom.f90

Co-authored-by: Jeremie Vandenplas <[email protected]>

Author: jalvesz

example/linalg/example_solve_custom.f90

@@ -53,7 +53,7 @@ subroutine solve_pcg_custom(A,b,x,di,tol,maxiter,restart,workspace)
         !-------------------------
         ! Jacobi preconditioner factorization
         call diag(A,diagonal)
-        where(abs(diagonal)>epsilon(0.d0)) diagonal = 1._dp/diagonal
+        where(abs(diagonal)>epsilon(0._dp)) diagonal = 1._dp/diagonal
         !-------------------------
         ! main call to the solver
         call solve_pcg_kernel(op,M,b,x,tol_,maxiter_,workspace_)