|
# SPDX-License-Identifier: GPL-2.0-or-later
# GaussForHomalg: Gauss functionality for the homalg project
#
# Implementations
#
## Matrix operations
####################################
#
# constructor functions and methods:
#
####################################
## <#GAPDoc Label="MatElm">
## <ManSection>
## <Meth Arg="M, r, c, R" Name="MatElm"/>
## <Returns><A>M</A>[<A>r</A>,<A>c</A>]</Returns>
## <Description>
## If the Eval attribute of the homalg matrix <A>M</A> over the &homalg;
## ring <A>R</A> is sparse, this calls the corresponding &Gauss;
## command <C>GetEntry</C>.
## </Description>
## </ManSection>
## <#/GAPDoc>
##
InstallMethod( MatElm,
[ IsHomalgInternalMatrixRep, IsPosInt, IsPosInt, IsHomalgInternalRingRep ],
function( M, i, j, R )
local m;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod();
fi;
return GetEntry( m, i, j ); #calls GetEntry for sparse matrices
end
);
## <#GAPDoc Label="SetMatElm">
## <ManSection>
## <Meth Arg="M, r, c, e, R" Name="SetMatElm"/>
## <Returns>nothing</Returns>
## <Description>
## If the Eval attribute of the homalg matrix <A>M</A> over the &homalg;
## ring <A>R</A> is sparse, this calls the corresponding &Gauss;
## command <C>GetEntry</C>, to achieve <C><A>M</A>[<A>r</A>,<A>c</A>]:=<A>e</A></C>.
## </Description>
## </ManSection>
## <#/GAPDoc>
##
InstallMethod( SetMatElm,
[ IsHomalgInternalMatrixRep and IsMutable, IsPosInt, IsPosInt, IsRingElement, IsHomalgInternalRingRep ],
function( M, i, j, e, R )
local m;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod();
fi;
SetEntry( m, i, j, e ); #calls SetEntry for sparse matrices
end
);
## <#GAPDoc Label="AddToMatElm">
## <ManSection>
## <Meth Arg="M, r, c, e, R" Name="AddToMatElm"/>
## <Returns>nothing</Returns>
## <Description>
## If the Eval attribute of the homalg matrix <A>M</A> over the &homalg;
## ring <A>R</A> is sparse, this calls the corresponding &Gauss;
## command <C>AddToEntry</C>, to achieve <C><A>M</A>[<A>r</A>,<A>c</A>]
## := <A>M</A>[<A>r</A>,<A>c</A>] + <A>e</A></C>.
## </Description>
## </ManSection>
## <#/GAPDoc>
##
InstallMethod( AddToMatElm,
[ IsHomalgInternalMatrixRep and IsMutable, IsPosInt, IsPosInt, IsRingElement, IsHomalgInternalRingRep ],
function( M, i, j, e, R )
local m;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod();
fi;
AddToEntry( m, i, j, e ); #calls AddToEntry for sparse matrices
end
);
##
InstallMethod( GetListOfHomalgMatrixAsString,
"for sparse matrices",
[ IsHomalgInternalMatrixRep, IsHomalgInternalRingRep ],
function( M, R )
local m;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod( );
fi;
m := ConvertSparseMatrixToMatrix( m );
if Characteristic( R ) > 0 then
return String( Concatenation( List( m, r -> List( r, Int ) ) ) );
fi;
return String( Concatenation( m ) );
end );
##
InstallMethod( GetListListOfStringsOfHomalgMatrix,
"for sparse matrices",
[ IsHomalgInternalMatrixRep, IsHomalgInternalRingRep ],
function( M, R )
local m;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod( );
fi;
m := ConvertSparseMatrixToMatrix( m );
return GetListListOfStringsOfMatrix( m, R );
end );
##
InstallMethod( GetSparseListOfHomalgMatrixAsString,
"for sparse matrices",
[ IsHomalgInternalMatrixRep, IsHomalgInternalRingRep ],
function( M, R )
local m, s, c, i, j, e;
m := Eval( M );
if not IsSparseMatrix( m ) then
TryNextMethod( );
fi;
s := [ ];
m := ConvertSparseMatrixToMatrix( m );
if Characteristic( R ) > 0 then
m := List( m, r -> List( r, Int ) );
fi;
c := Length( m[1] );
for i in [ 1 .. Length( m ) ] do
for j in [ 1 .. c ] do
e := m[i][j];
if not IsZero( e ) then
Add( s, [ String( i ), String( j ), String( e ) ] );
fi;
od;
od;
s := JoinStringsWithSeparator( List( s, JoinStringsWithSeparator ), "],[" );
return Concatenation( "[[", s, "]]" );
end );
[ Dauer der Verarbeitung: 0.2 Sekunden
(vorverarbeitet)
]
|
