;; -*-theme-d-*-

;; Copyright (C) 2014  Tommi Höynälänmaa
;; Distributed under GNU General Public License version 3,
;; see file doc/GPL-3.

;; Expected results: translation and running OK


(define-proper-program (tests test259)


  (import (standard-library core)
	  (standard-library complex)
	  (standard-library matrix)
	  (standard-library console-io))


  (define-param-method display-matrix (%number)
		       (((mx (:matrix %number)))
			<none>
			(nonpure))
    (let ((rows (field-ref mx 'rows))
	  (columns (field-ref mx 'columns)))
      (do ((i1 <integer> 0 (+ i1 1))) ((>= i1 rows))
	(begin
	  (do ((i2 <integer> 0 (+ i2 1))) ((>= i2 columns))
	    (console-display (matrix-ref mx i1 i2))
	    (if (< i2 (- columns 1))
		(console-display-string " ")))
	  (console-newline)))))


  (define-param-method display-matrix (%number)
		       (((mx (:diagonal-matrix %number)))
			<none>
			(nonpure))
    (let ((rows (field-ref mx 'len))
	  (columns (field-ref mx 'len)))
      (do ((i1 <integer> 0 (+ i1 1))) ((>= i1 rows))
	(begin
	  (do ((i2 <integer> 0 (+ i2 1))) ((>= i2 columns))
	    (console-display (matrix-ref mx i1 i2))
	    (if (< i2 (- columns 1))
		(console-display-string " ")))
	  (console-newline)))))


  (define-param-method display-matrix-nl (%number)
		       (((mx (:matrix %number)))
			<none>
			(nonpure))
    (display-matrix mx)
    (console-newline))


  (define main
    (lambda (() <none> nonpure)
      (let ((mx1 (matrix
		  (list
		   (list (complex 0.0 0.0) (complex 1.0 -4.0) (complex 1.5 2.5))
		   (list (complex 1.1 2.1) (complex 3.5 -1.0) (complex 0.0 0.0))
		   (list (complex 1.0 0.0) (complex 0.0 -2.5) (complex 1.0 1.0)))))
	    (mx2 (diagonal-matrix
		  (list (complex -4.5 0.0) (complex -1.0 1.0) (complex 2.5 1.5))))
	    (cv (column-vector
		 (list (complex 5.0 1.0) (complex 0.0 1.0) (complex 2.5 1.0)))))
	(display-matrix-nl (* mx1 mx2))
	(display-matrix-nl (* mx1 cv))
	(display-matrix-nl (* mx2 cv))))))
