matrix/0000755000076500007650000000000010560207605013601 5ustar ptolemyptolemy00000000000000matrix/._.DS_Store0000400000076500000000000000012210560207607015055 0ustar ptolemywheel00000000000000Mac OS X  2 R@matrix/.DS_Store0000644000076500007650000001400410560207607015265 0ustar ptolemyptolemy00000000000000Bud1%  @€ @€ @€ @ E%DSDB`€ @€ @€ @matrix/matrices.m0000644000076500007650000000142510557745114015600 0ustar ptolemyptolemy00000000000000:- module matrices. % Some common matrix instance types. :- interface. :- import_module int, float, rational, matrix. :- instance number(int). :- instance number(float). :- instance number(rational). :- implementation. :- instance number(int) where [ func(matrix.(+)/2) is int.(+), func(matrix.(-)/2) is int.(-), func(matrix.(*)/2) is int.(*), pred(matrix.(<)/2) is int.(<) ]. :- instance number(float) where [ func(matrix.(+)/2) is float.(+), func(matrix.(-)/2) is float.(-), func(matrix.(*)/2) is float.(*), pred(matrix.(<)/2) is float.(<) ]. :- instance number(rational) where [ func(matrix.(+)/2) is rational.(+), func(matrix.(-)/2) is rational.(-), func(matrix.(*)/2) is rational.(*), pred(matrix.(<)/2) is rational.(<) ]. matrix/matrix.m0000644000076500007650000000723710557776302015307 0ustar ptolemyptolemy00000000000000:- module matrix. % A module that implements fast and comprehensive operations for large % matrices. :- interface. :- import_module int, io, array2d. :- typeclass number(T) where [ func T * T = T, func T + T = T, func T - T = T, pred (T::in) < (T::in) is semidet ]. :- type matrix(T). :- func as_array2d(matrix(T)) = array2d(T). :- func from_array2d(array2d(T)) = matrix(T) <= number(T). :- func fill(T, int, int) = matrix(T) <= number(T). :- func id(T, int) = matrix(T) <= number(T). % dot product (matrix x scalar; scalar x matrix) % :- func o(matrix(T), T) = matrix(T). :- func o(T, matrix(T)) = matrix(T) <= number(T). % :- mode o(in, array2d_di) = out(array2d) is det. :- func maxin(matrix(T)) = T <= number(T). :- pred print_matrix(matrix(T)::in, io::di, io::uo) is det. :- implementation. :- import_module array2d. :- type matrix(T) == array2d(T). as_array2d(M) = M. from_array2d(M) = M. fill(Val, Y, X) = new(X, Y, Val). id(Id, Side) = Matrix :- M0 = new(Side, Side, Id - Id), diag_set(0, Side, Id, M0, Matrix). o(T, M0) = Matrix :- bounds(M0, Y, X), destructive_copy_of(M0, M1), scale_rows(0, Y, X, T, M1, Matrix). maxin(M) = Ans :- bounds(M, Y, X), maxin_rows(M, 0, Y, X, M ^ elem(0, 0), Ans). print_matrix(Matrix) --> { bounds(Matrix, Y, X) }, print_rows(Matrix, 0, Y, X). :- pred maxin_rows(matrix(T), int, int, int, T, T) <= number(T). :- mode maxin_rows(in, in, in, in, in, out) is det. maxin_rows(M, Start, Rows, Cols, Guess, Ans) :- Start = Rows -> Guess = Ans ; maxin_row(M, Start, 0, Cols, Guess, NewGuess), maxin_rows(M, Start + 1, Rows, Cols, NewGuess, Ans). :- pred maxin_row(matrix(T), int, int, int, T, T) <= number(T). :- mode maxin_row(in, in, in, in, in, out) is det. maxin_row(M, Row, Start, Cols, Guess, Ans) :- Start = Cols -> Guess = Ans ; NewGuess = M ^ elem(Row, Start), (NewGuess < Guess -> New = Guess ; New = NewGuess), maxin_row(M, Row, Start + 1, Cols, New, Ans). :- pred destructive_copy_of(matrix(T), matrix(T)) <= number(T). :- mode destructive_copy_of(in, out(array2d)) is det. destructive_copy_of(M, M). :- pred print_rows(matrix(T), int, int, int, io, io). :- mode print_rows(in, in, in, in, di, uo) is det. print_rows(Matrix, Start, Finish, Length) --> { Start = Finish } -> [] ; print("| "), print_each_row(Matrix, Start, 0, Length), print("|\n"), print_rows(Matrix, Start + 1, Finish, Length). :- pred print_each_row(matrix(T), int, int, int, io, io). :- mode print_each_row(in, in, in, in, di, uo) is det. print_each_row(Matrix, Row, Start, Finish) --> { Start = Finish } -> [] ; print(Matrix ^ elem(Row, Start)), print(" "), print_each_row(Matrix, Row, Start + 1, Finish). :- pred diag_set(int, int, T, matrix(T), matrix(T)). :- mode diag_set(in, in, in, array2d_di, out(array2d)) is det. diag_set(Start, Finish, Id, M0, Matrix) :- Start = Finish -> M0 = Matrix ; M1 = M0 ^ elem(Start, Start) := Id, diag_set(Start + 1, Finish, Id, M1, Matrix). :- pred scale_rows(int, int, int, T, matrix(T), matrix(T)) <= number(T). :- mode scale_rows(in, in, in, in, array2d_di, out(array2d)) is det. scale_rows(Start, Finish, Length, Value) --> { Start = Finish } -> [] ; scale_row_elems(Start, 0, Length, Value), scale_rows(Start + 1, Finish, Length, Value). :- pred scale_row_elems(int, int, int, T, matrix(T), matrix(T)) <= number(T). :- mode scale_row_elems(in, in, in, in, array2d_di, out(array2d)) is det. scale_row_elems(Row, Start, Finish, Value, M0, Matrix) :- Start = Finish -> M0 = Matrix ; Elem = M0 ^ elem(Row, Start), M1 = M0 ^ elem(Row, Start) := Elem * Value, scale_row_elems(Row, Start + 1, Finish, Value, M1, Matrix). matrix/test_matrix.m0000644000076500007650000000101710557776600016335 0ustar ptolemyptolemy00000000000000:- module test_matrix. :- interface. :- import_module io. :- pred main(io::di, io::uo) is det. :- implementation. :- import_module matrix, matrices, array2d, list. main --> { M = id(1.0, 5), N = fill(7.3, 3, 6), O = o(0.5, M), P = o(3.0, N), Q = from_array2d(array2d([[3,1,4], [1,5,9], [2,6,5]])), Maxy = maxin(Q) }, % We know this an int; does the compiler? print_matrix(O), nl, nl, nl, print_matrix(P), print("The max element in:\n"), print_matrix(Q), print("\nis "), print(Maxy), nl.