tebdol

Simulation of ultracold atoms in optical lattices
git clone https://miroslavurbanek.com/git/tebdol.git
Log | Files | Refs | README

commit cbb1537f28a49b20d17d09180288240049f72b39
parent 040c8fb8c22a1883384f31a4076f9b81e7da64af
Author: Miroslav Urbanek <mu@miroslavurbanek.com>
Date:   Thu, 29 Jun 2017 01:56:46 +0200

Do not mask floating-point errors

It is better to inform the user that an error occurred.

Diffstat:
tebdol/exp.lisp | 5++---
1 file changed, 2 insertions(+), 3 deletions(-)

diff --git a/tebdol/exp.lisp b/tebdol/exp.lisp @@ -39,9 +39,8 @@ (let ((bwork (double-array-alien work)) (brwork (double-array-alien rwork)) (biwork (integer-array-alien iwork))) - (sb-int:with-float-traps-masked (:overflow :invalid :divide-by-zero) - (zheevr "v" "a" "u" n ba n 0d0 0d0 0 0 abstol 0 bw bz n - bisuppz bwork lwork brwork lrwork biwork liwork)))))) + (zheevr "v" "a" "u" n ba n 0d0 0d0 0 0 abstol 0 bw bz n + bisuppz bwork lwork brwork lrwork biwork liwork))))) (m) (setf lwork (floor (realpart (aref work 0)))) (setf work (make-blas-array lwork))