View difference between Paste ID: KwNF1F25 and iJpUBh4w

 #|;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-HW4 #5 A
+HW4 #6 A
-ScalarProductDoLoop.rkt by Ryan Liszewski
+DisplayMatrixRowCol.rkt by Ryan Liszewski
-Calculates the scalar product of the given 2 vectors using a do loop
+Reads matrix from a file and display a specified row or column
-produces an error message when the vectors have different lengths or is empty
+
 #lang racket
 ;read matrix from file
-(define scalar-product
+(define (read-matrix filename)
-  (lambda (x y)
+  (let* ((inport (open-input-file filename))
-    (cond((= (vector-length x) (vector-length y))
+         (nrow (read inport))
-           (let ((z 0))
+         (ncol (read inport))
-                  (do ((i 0 (add1 i)))
+         (mat (make-vector nrow)))
-                    ((>= i (vector-length x)) (display z))
+    (do ((i 0 (add1 i)))
-                    (set! z
+      ((>= i nrow) (close-input-port inport) mat)
-                          (+ z (* (vector-ref x i) (vector-ref y i)))))))
+      (let ((row (make-vector ncol)))
-          (else (display "Error: different sizes of vectors!")))))
+        (do ((j 0 (add1 j)))
           ((>= j ncol) (vector-set! mat i row))
           (vector-set! row j (read inport)))))))
 
-HW4 #5 B
+;display vector
-ScalarProductRecursive.rkt by Ryan Liszewski
+(define (display-vector v)
-Calculates the scalar product of the given 2 vectors using recursion
+  (do ((i 0 (add1 i)))
-produces an error message when the vectors have different lengths
+    ((>= i (vector-length v)) (display ""))
     (display (vector-ref v i)) (display " ")))
 
-(define scalar-product
+;return i-th row
-  (lambda (x y)
+(define (ro filename i)
-    ;If they equal the same size
+  (define mat (read-matrix filename))
-    (cond ((= (vector-length x) (vector-length y))
+  (vector-ref mat i))
-            ;FInd scalar product
+;display i-th row
-           (recur (vector->list x) (vector->list y)))
+(define (row filename i)
-           ;else the size are not the same
+  (display-vector (ro filename (- i 1))))
-          (else (display "Error: different sizes of vectors!")))))
+
-;Calls function recursively
+;return j-th col
-(define recur
+(define (co filename j)
-    (lambda (list1 list2)
+  (define mat (read-matrix filename))
-      (cond ((null? (cdr list1)) (* (car list1) (car list2)))
+  (define nrow (vector-length mat))
-            (else (+ (* (car list1) (car list2))
+  (define column (make-vector nrow))
-                     (recur (cdr list1) (cdr list2)))))))
+  (do ((i 0 (add1 i)))
     ((>= i nrow) column)
     (vector-set! column i (vector-ref (vector-ref mat i) j))))
 
 ;display j-th col
 (define (col filename j)
   (display-vector (co filename (- j 1))))
 
 
 
 
 
 
 
 #|;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 HW4 #6 B
 MatrixMultiplication.rkt by Ryan Liszewski
 
 Multiplies two matrices from two input files and produces the result in an output file
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;|#
 #lang racket
 ;read matrix from file
 (define (read-matrix filename)
   (let* ((inport (open-input-file filename))
          (nrow (read inport))
          (ncol (read inport))
          (mat (make-vector nrow)))
     (do ((i 0 (add1 i)))
       ((>= i nrow) (close-input-port inport) mat)
       (let ((row (make-vector ncol)))
         (do ((j 0 (add1 j)))
           ((>= j ncol) (vector-set! mat i row))
           (vector-set! row j (read inport)))))))
 
 ;display components of the given vector
 (define (display-vector v)
   (do ((i 0 (add1 i)))
     ((>= i (vector-length v)) (display ""))
     (display (vector-ref v i)) (display " ")))
 
 ;return i-th row
 (define (ro filename i)
   (define mat (read-matrix filename))
   (vector-ref mat i))
 
 ;display i-th row
 (define (row filename i)
   (display-vector (ro filename (- i 1))))
 
 ;return j-th col
 (define (co filename j)
   (define mat (read-matrix filename))
   (define nrow (vector-length mat))
   (define column (make-vector nrow))
   (do ((i 0 (add1 i)))
     ((>= i nrow) column)
     (vector-set! column i (vector-ref (vector-ref mat i) j))))
 
 ;display j-th col
 (define (col filename j)
   (display-vector (co filename (- j 1))))
 
 ;return dot product
 (define (dot-product v1 v2)
   (do ((i 0 (add1 i)) (sum 0 (+ sum (* (vector-ref v1 i)
                                        (vector-ref v2 i)))))
     ((>= i (vector-length v1)) sum)))
 
 ;matrix multiplication
 (define (mmul f1 f2 f3)
   (define m1 (read-matrix f1))
   (define m2 (read-matrix f2))
   (define nrow (vector-length m1))
   (define ncol (vector-length m2))
   (define outport (open-output-file f3))
   (display nrow outport) (display " " outport)
   (display ncol outport) (newline outport)
   (do ((i 0 (add1 i)))
     ((>= i nrow) (close-output-port outport) (display ""))
     (let ((row (make-vector ncol)))
       (do ((j 0 (add1 j)))
         ((>= j ncol) (display-vector row) (newline) (newline outport))
         (vector-set! row j (dot-product (ro f1 i) (co f2 j)))
         (display (vector-ref row j) outport) (display " \n" outport)))))