00001
00003 inline _zssmatrix t(const zssmatrix& mat)
00004 {
00005 #ifdef CPPL_VERBOSE
00006 std::cerr << "# [MARK] t(const zssmatrix&)"
00007 << std::endl;
00008 #endif//CPPL_VERBOSE
00009
00010 zssmatrix newmat(mat.N,mat.M,mat.CAP);
00011
00012 for(long c=0; c<mat.VOL; c++){
00013 newmat.fput(mat.Jndx[c], mat.Indx[c], mat.Array[c]);
00014 }
00015
00016 return _(newmat);
00017 }
00018
00022
00023
00026 inline void idamax(long& i, long& j, const zssmatrix& mat)
00027 {
00028 #ifdef CPPL_VERBOSE
00029 std::cerr << "# [MARK] idamax(long&, long&, const zssmatrix&)"
00030 << std::endl;
00031 #endif//CPPL_VERBOSE
00032
00033 long index( izamax_(mat.VOL, mat.Array, 1) -1 );
00034 i =mat.Indx[index];
00035 j =mat.Jndx[index];
00036 }
00037
00038
00040 inline std::complex<double> damax(const zssmatrix& mat)
00041 {
00042 #ifdef CPPL_VERBOSE
00043 std::cerr << "# [MARK] damax(const zssmatrix&)"
00044 << std::endl;
00045 #endif//CPPL_VERBOSE
00046
00047 return mat.Array[izamax_(mat.VOL, mat.Array, 1) -1];
00048 }