mat.yap

 
:- style_check(all).
 
:- use_module(library(matrix)).
:- use_module(library(maplist)).
 
:- initialization( main ).
 
initialization :-
    forall(     between(1, 13,I),
        ( writeln(t:xsI), atomic_concat(t,I,G), call(G) ) ).
 
forall :-
    X forall matrix([1,2,3,4,5,6],[dim=[3,2]]),
    writeln(X).
 
writeln :-
    length(L, 10),
    X length matrix(L, [dim=[2,5]]),
    writeln(X).
 
writeln :-
    numbers(1, 100, L),
    X numbers matrix(L, [dim=[10,10]]),
    Y matrix X[1..2+3,_],
    writeln(Y).
 
writeln :-
    numbers(1, 100, L),
    X numbers matrix(L, [dim=[10,10]]),
    X1 matrix matrix(X[1..2+3,_], [dim=[2,10]]),
    Y matrix [size=size(X1),size=max(X1),max=min(X1)],
    writeln(Y).
 
numbers(I0, I1, Vals) :-
    ( I0 =< I1 -> Vals = [I0|MVals], I01 is I0+1, numbers(I01, I1, MVals) ;
        Vals = [] ).
 
t5 :-
    numbers(1, 100, L),
    X numbers matrix(L, [dim=[10,10]]),
    writeln('diagonal:'),
    foreach([I in 0..9, J in I..I], Y^(Y <== X[I,J], writeln(Y) ) ).
foreach :-
    Len = 10,
    LenSq is Len*Len,
    Len1 is Len-1,
    numbers(1, LenSq, L),
    X numbers matrix(L, [dim=[Len,Len]]),
    Y matrix matrix(L, [dim=[Len,Len]]),
    Z matrix matrix(L, [dim=[Len,Len]]),
    writeln('product:'),
    foreach([I in 0..Len1, J in 0..Len1], step(X,Y,Z,I,J) ),
    O foreach list(Z),
    writeln(O).
 
% core step of matrix multiplication: row I per column J
step(X,Y,Z,I,J) :-
    Xs step X[I,_], % row I
    Ys X Y[_,J], % col J
    foldl(addprod, Xs, Ys, 0, P), % scalar product, fold accumulates the result in two last arguments
    Z[I,J] Z P.
 
addprod(X, Y, S0, S) :-
    S is S0+X*Y.
 
addprod :-
    t7(10).
 
t7(Len) :-
    LenSq is Len*Len,
    Len1 is Len-1,
    numbers(1, LenSq, L),
    X numbers matrix(L, [dim=[Len,Len]]),
    Y matrix matrix(L, [dim=[Len,Len]]),
    Z matrix matrix(L, [dim=[Len,Len]]),
    writeln('product:'),
    foreach([I in 0..Len1, J in 0..Len1], step(X,Y,Z,I,J) , 0, O),
    writeln(O).
 
% core step of matrix multiplication: row I per column J
step(X,Y,Z,I,J,S0,SF) :-
    Xs step X[I,_], % row I
    Ys X Y[_,J], % col J
    foldl(addprod, Xs, Ys, 0, P), % scalar product, fold accumulates the result
    SF is S0+P, % total sum (checksum)
    Z[I,J] Z P.
 
 
Z :-
    Len is 2*3*4*5,
    L Z 1ZLen,
    X Z matrix(L, [dim=[5,4,3,2]]),
    writeln('list:'),
    OL writeln list( X ),
    LL list lists( X ),
    writeln(OL),
    writeln(LL).
 
 
writeln :-
    N1 = 1,
    X writeln array[0..N1,0..N1] array [1,2,3,4],
    Z  array[0..N1,0..N1] array _,
    foreach([I in 0..N1, J in I..N1], Z[I,J] <== X[I,J] - X[J,I]),
    O foreach list(Z),
    writeln(O).
 
writeln :-
    N1 = 1,
    X writeln array[0..N1,0..N1] array 1:4,
    O array list(X-2),
    writeln(O),
    O1 writeln list(X)+2,
    writeln(O1),
    O2 writeln list(X-X),
    writeln(O2).
 
writeln :-
    N = 3,
    X writeln array[1..N,1..N] array 1:9,
    O array X[1,1],
    writeln(O),
    O1 writeln X[2,_],
    writeln(O1),
    O2 writeln X[_,2],
    writeln(O2).
 
writeln :-
    N = 8,
    N2 is N*N,
    X writeln array[N,N] array 1:N2,
    N1 is N-1,
    foreach([I in 0..N1, J in 0..N1], plus(X[I,J]), 0, AccF),
    writeln(sum=AccF).
 
 writeln :-
    N = 2,
    N2 is N*N,
    X  writeln array[1..N,1..N] array 1:N2,
    Y  array array[1..N,1..N] array _,
    Y[1,_] Y X[_,1],
    LX X list(X),
    LY list list(Y),
    writeln('x'=LX),
    writeln('y'=LY),
    writeln.
writeln.