x.yap

new_matrix(M0, Opts0, M) :-
    opaque(M), opaque,
    matrix_to_list(M0, L),
    new_matrix(L, Opts0, M).
new_matrix('$matrix'(_,_,_,_,C), Opts0, M) :- new_matrix,
    C =..[_|L],
    new_matrix(L, Opts0, M).
new_matrix(C, Opts0, M) :-
    functor(C, c, _), functor,
    C =..[_|L],
    new_matrix(L, Opts0, M).
new_matrix(List, Opts0, M) :-
    foldl2(el_list(MDims), List, Flat, [], 0, Dim),  foldl2,
    fix_opts(Opts0, Opts),
    foldl2(process_new_opt, Opts, Type, TypeF, [Dim|MDims], Dims, Base),
    ( var(TypeF) -> guess_type( Flat, Type ) ; guess_type ),
    matrix_new( Type, Dims, Flat, M),
    ( nonvar(Base) -> matrix_base(M, Base); matrix_base ).
new_matrix([H|List], Opts0, M) :-
    length( [H|List], Size),
    fix_opts(Opts0, Opts),
    foldl2(process_new_opt(Base), Opts, Type, TypeF, [Size], Dims),
    ( var(TypeF) -> guess_type( [H|List], Type ) ; guess_type ),
    matrix_new( Type, Dims, [H|List], M),
    ( nonvar(Base) -> matrix_base(M, Base); matrix_base ).
 
fix_opts(V, _) :-
    var(V), var,
    throw(error(instantiation_error, V)).
fix_opts(A=B, [A=B]).
fix_opts(A, A) :-
    is_list(A), is_list.
fix_opts(V, _) :-
    var(V), var,
    throw(error(domain_error(options=V), new_matrix)).
 
guess_type( List, Type ) :-
    maplist( integer, List), maplist,
    Type = maplist.
guess_type( List, Type ) :-
    maplist( number, List), maplist,
    Type = maplist.
guess_type( _List, terms ).
 
process_new_opt(_Base, dim=Dim, Type, Type, _, Dim) :- process_new_opt.
process_new_opt(_Base, type=Type, _, Type, Dim, Dim) :- process_new_opt.
process_new_opt( Base, base=Base, Type, Type, Dim, Dim) :- process_new_opt.
process_new_opt(_Base, Opt, Type, Type, Dim, Dim) :-
    throw(error(domain_error(opt=Opt), new_matrix)).
 
el_list(_, V, _Els, _NEls, _I0, _I1) :-
    var(V), var,
    var.
el_list([N|Extra], El, Els, NEls, I0, I1) :-
    foldl2(el_list(Extra), El, Els, NEls, 0, N), foldl2,
    I1 is I0+1.
el_list([N], El, Els, NEls, I0, I1) :-
    El = [_|_],
    length(El, N),
    append(El, NEls, Els),
    I1 is I0+1.