Changeset 5833 in project


Ignore:
Timestamp:
08/31/07 08:58:48 (12 years ago)
Author:
Ivan Raikov
Message:

Removed m from method definition structure.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • nordsieck/trunk/nordsieck.scm

    r5831 r5833  
    5151(define-matrix-map f64)
    5252
    53 (define-record nordsieck-def order Ainv A B C V Econst d dhist dhistlen p1 p2 m n)
     53(define-record nordsieck-def order Ainv A B C V Econst d dhist dhistlen p1 p2 n)
    5454
    5555(define-record-printer (nordsieck-def x out)
    56   (fprintf out "#(nordsieck order=~A Ainv=~S A=~S B=~S C=~S V=~S Econst=~S d=~S dhist=~S p1=~S p2=~S m=~S n=~S)"
     56  (fprintf out "#(nordsieck order=~A Ainv=~S A=~S B=~S C=~S V=~S Econst=~S d=~S dhist=~S p1=~S p2=~S n=~S)"
    5757           (if (= blas:RowMajor (nordsieck-def-order x))
    5858               "row" "col")
     
    199199
    200200(define (nordsieck-form order Econst a b . rest)
    201   (let-optionals rest ((va #f) (vb #f) (dhistlen 3) (m 1))
     201  (let-optionals rest ((va #f) (vb #f) (dhistlen 3))
    202202    (let* ((implicit? (and va vb))
    203203           (p1  (- (f64vector-length a) 1))
     
    227227                    (print "C = " C)
    228228                    (print "V = " V)
    229                     (make-nordsieck-def order Ainv A B C V Econst d (list) dhistlen p1 p2 m n))))))))))
     229                    (make-nordsieck-def order Ainv A B C V Econst d (list) dhistlen p1 p2 n))))))))))
    230230
    231231 
    232232(define (nordsieck-step def h prev frhs)
    233233  (let ((order (nordsieck-def-order def))
    234         (m     (nordsieck-def-m def))
     234        (p2     (nordsieck-def-p2 def))
    235235        (n     (nordsieck-def-n def))
    236236        (C     (nordsieck-def-C def))
     
    238238  (let ((y0  (let ((v (matrix-zeros n 1)))
    239239               (blas:dgemv! blas:NoTrans n n 1.0 C prev 0.0 v))))
    240     (let loop ((i 0) (m m) (ym (list y0)))
    241       (if (fx<= m 1)  ym
     240    (let loop ((i 0) (m p2) (ym (list y0)))
     241      (if (fx<= m 0)  ym
    242242          (let ((yi  (let ((fy  (fp- (fp* h (frhs (f64vector-ref (car ym) 0)))
    243243                                     (fp* h (f64vector-ref (car ym) 1)))))
Note: See TracChangeset for help on using the changeset viewer.