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