00001
00002
00003 inline void dgematrix::transpose()
00004 {
00005 _dgematrix newmat(N,M);
00006
00007 for(long i=0; i<newmat.M; i++){ for(long j=0; j<newmat.N; j++){
00008 newmat(i,j)=operator()(j,i);
00009 }}
00010
00011 shallow_copy(newmat);
00012 }
00013
00014
00015
00016 inline _dgematrix t(const dgematrix& mat)
00017 {
00018 _dgematrix newmat(mat.N,mat.M);
00019
00020 for(long i=0; i<newmat.M; i++){ for(long j=0; j<newmat.N; j++){
00021 newmat(i,j) =mat(j,i);
00022 }}
00023
00024 return newmat;
00025 }
00026