| 12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152 |
- #| -*-Scheme-*-
- Copyright (C) 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994,
- 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005,
- 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013 Massachusetts
- Institute of Technology
- This file is part of MIT/GNU Scheme.
- MIT/GNU Scheme is free software; you can redistribute it and/or modify
- it under the terms of the GNU General Public License as published by
- the Free Software Foundation; either version 2 of the License, or (at
- your option) any later version.
- MIT/GNU Scheme is distributed in the hope that it will be useful, but
- WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
- General Public License for more details.
- You should have received a copy of the GNU General Public License
- along with MIT/GNU Scheme; if not, write to the Free Software
- Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301,
- USA.
- |#
- ;;; We define the standard singular-matrix failure continuations as follows.
- (define (barf-on-zero-pivot dismiss)
- (singular-matrix-error))
- (define (allow-zero-pivot dismiss)
- (dismiss))
- ;;; default value
- (define (default-singular-matrix-error)
- (error "Singular matrix - zero pivot"))
- ;;; Rebind this to catch errors
- ;; FBE: make this a parameter here and in 'full-pivot.scm', 'gauss-jordan.scm'
- (define singular-matrix-error (make-parameter default-singular-matrix-error))
- ;; (set! singular-matrix-error
- ;; default-singular-matrix-error)
- (define (with-singular-matrix-handler handler thunk)
- (parameterize ((singular-matrix-error handler))
- (thunk)))
- (define (handle-singularity-errors-with error-handler)
- (singular-matrix-error error-handler))
|
