Changeset 5667 in project


Ignore:
Timestamp:
08/24/07 05:18:11 (12 years ago)
Author:
Ivan Raikov
Message:

Added Econst field to the nordsieck record definition.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • nordsieck/trunk/nordsieck.scm

    r5643 r5667  
    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 d=~S d=~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 m=~S n=~S)"
    5757           (if (= blas:RowMajor (nordsieck-def-order x))
    5858               "row" "col")
     
    6262           (nordsieck-def-C x)
    6363           (nordsieck-def-V x)
     64           (nordsieck-def-Econst x)
    6465           (nordsieck-def-d x)
    6566           (nordsieck-def-dhist x)
     
    197198
    198199
    199 (define (nordsieck-form order a b . rest)
     200(define (nordsieck-form order Econst a b . rest)
    200201  (let-optionals rest ((va #f) (vb #f) (dhistlen 3) (m 1))
    201202    (let* ((implicit? (and va vb))
     
    226227                    (print "C = " C)
    227228                    (print "V = " V)
    228                     (make-nordsieck-def order Ainv A B C V d (list) dhistlen p1 p2 m n))))))))))
     229                    (make-nordsieck-def order Ainv A B C V Econst d (list) dhistlen p1 p2 m n))))))))))
    229230 
    230231(define (nordsieck-step def prev h frhs xrhs)
     
    265266                           (blas:dgemv! order blas:NoTrans n n 1.0 DKV gamma 0.0 v))))
    266267                (match def
    267                        (($ nordsieck-def order Ainv A B C V d dhist dhistlen p1 p2 m n)
     268                       (($ nordsieck-def order Ainv A B C V Econst d dhist dhistlen p1 p2 m n)
    268269                        (let ((dhist  (if (fx>= (length dhist) dhistlen)
    269270                                          (take dhist (fx- dhistlen 1)) dhist)))
    270                           (make-nordsieck-def order Ainv A B C V d1 (cons d dhist) dhistlen p1 p2 m n))))))))))))
     271                          (make-nordsieck-def order Ainv A B C V Econst d1 (cons d dhist) dhistlen p1 p2 m n))))))))))))
    271272       
Note: See TracChangeset for help on using the changeset viewer.