Hi
I modified Eigenschink's matrix-0.5 LU to be able to use slibs' array
the method I used
are just changing a few lines. I placed it here
http://www.tip9ug.jp/who/fernan/mat.tar.gz
Now I am testing it with a simple finite difference problem from here
http://www.adeptscience.co.uk/products/mathsim/maple/powertools/des/unit18.html
The matrix:plu-factor seems to take about 55% of the processing time.
This is
after putting a few (current-time) between the routines.
This is the ****tion that I am trying to find a way to make use of
subarray and maybe
make the whole thing faster.
You guys can look at it and maybe it would be straight-forward for you
guys.
77 (define lu-start (current-time))
78 (do ((k 0 (1+ k))) ((= k (- rows 1)))
79 (do ((i (1+ k) (1+ i))) ((= i rows))
80 (let ((m_ik (/ (matrix:ref a i k) (matrix:ref a k k))))
81 (matrix:set! a i k m_ik)
82 (do ((j (1+ k) (1+ j))) ((= j cols))
83 (matrix:set! a i j
84 (- (matrix:ref a i j)
85 (* m_ik (matrix:ref a k j))))))))
86 (display "LU-Factor time: ")
87 (display (- (current-time) lu-start))
88 (display " sec(s)\n")
fernan
