svdsum

求矩阵奇异分解svd算法

	   shift programfiles
    end
    echo "|" >> $output


    echo -n "|  COLUMN DIMENSION OF A  |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
	while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
	   shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          set value = `grep "NO. OF TERMS " $coutputfile | awk '{printf ("%-13s ", $7) }'`
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
		  set value = `grep "NO. OF TERMS " $coutputfile | awk '{printf ("%-13s ", $7) }'`
	   else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
		  set value = `grep "NO. OF TERMS " $coutputfile | awk '{printf ("%-13s ", $7) }'`
	   else
          set value = `grep "NO. OF TERMS " $coutputfile | awk '{printf ("%-13s ", $9) }'`
       endif
       if ($value !~ [0-9]*) then
          set ncolerror = `grep -c 'NCOL MUST NOT BE GREATER THAN NROW' $coutputfile`
          if ($ncolerror == 1) then
             echo -n "|   ERROR 6   " >> $output
             set error6 = 1
             set ncolerror = 0
          else
             echo -n "|      -        " >> $output
          endif
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          grep "NO. OF TERMS " $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
		  grep "NO. OF TERMS " $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
	   else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
		  grep "NO. OF TERMS " $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else
          grep "NO. OF TERMS " $coutputfile | awk '{printf ("| %-13s ", $9) >> "'$output'" }'
       endif
       @ c++
	   shift programfiles
    end
    echo "|" >> $output


    echo -n "| AUXILLARY MEMORY(BYTES) |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
	while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
	   shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          set value = `grep "ALLOCATED MEMORY" $coutputfile | awk '{printf ("%-13s ", $5) }'`
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
		  set value = `grep "MEMORY REQUIRED" $coutputfile | awk '{printf ("%-13s ", $6) }'`
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
		  set value = `grep "MEM REQUIRED" $coutputfile | awk '{printf ("%-13s ", $6) }'`
       else
          set value = `grep "MEMORY NEEDED" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       endif
       if ($value !~ [0-9]*) then
          echo -n "|      -        " >> $output
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          grep "ALLOCATED MEMORY" $coutputfile | awk '{printf ("| %-13s ", $5) >> "'$output'" }'
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
		  grep "MEMORY REQUIRED" $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
		  grep "MEM REQUIRED" $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
       else
          grep "MEMORY NEEDED" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       endif
       @ c++
	   shift programfiles
    end
    echo "|" >> $output


    echo -n "|  WANT S-VECTORS? [T/F]  |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
	while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
       shift programfiles
	   if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
          set value = `grep "WANT S-VECTORS?" $coutputfile | awk '{printf ("%-13s ", $6) }'`
       else
          set value = `grep "WANT S-VECTORS" $coutputfile | awk '{printf ("%-13s ", $6) }'`
       endif
       if ($value !~ [T,F]*) then
          echo -n "|      -        " >> $output
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
          grep "WANT S-VECTORS?" $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
       else
          grep "WANT S-VECTORS" $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
       endif
       @ c++
	   shift programfiles
    end
    echo "|" >> $output



    echo -n "| NO. OF STEPS/ITERATIONS |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
    while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
       shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          set value = `grep "NUMBER OF LANCZOS STEPS" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       else
          set value = `grep "NO. OF ITERATIONS TAKEN" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       endif
       if ($value !~ [0-9]*) then
		  set nw1error = `grep -c '6\*N\+4\*LANMAX\+1 \+ LANMAX\*LANMAX CANNOT EXCEED NW' $coutputfile`
          if ($nw1error == 1) then
             echo -n "|   ERROR 13   " >> $output
             set error13 = 1
             set nw1error = 0
          else
             echo -n "|      -        " >> $output
          endif
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          grep "NUMBER OF LANCZOS STEPS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else
          grep "NO. OF ITERATIONS TAKEN" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       endif
       @ c++
       shift programfiles
    end
    echo "|" >> $output


   echo -n "|        TOLERANCE        |" >> $output
   @ c = 1 + ($count * 5)
   set programfiles = `echo $programlist`
   while ($#programfiles > 0)
      set coutputfile = $programfiles[1].out$id[$c]
      shift programfiles
      if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
         set value = `grep "KAPPA" $coutputfile | awk '{printf ("%12s ", $4) }'`
      else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
         set value = `grep "TOLERANCE" $coutputfile | awk '{printf ("%12s ", $6) }'`
      else
         set value = `grep "TOLERANCE" $coutputfile | awk '{printf ("%-13s ", $4) }'`
      endif
      if ($value !~ [0-9]*) then
	     set toleranceerror = `grep -c 'TOLERANCE IS INVALID' $coutputfile`
		 if ($toleranceerror == 1) then
            echo -n "|   ERROR 7    " >> $output
            set error7 = 1
            set toleranceerror = 0
         else 
            echo -n "|      -        " >> $output
         endif
      else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
         grep "KAPPA" $coutputfile | awk '{printf ("| %-13s ", $4) >> "'$output'" }'
      else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
         grep "TOLERANCE" $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
      else
         grep "TOLERANCE" $coutputfile | awk '{printf ("| %-13s ", $4) >> "'$output'" }'
      endif
      @ c++
      shift programfiles
   end
   echo "|" >> $output


    echo -n "|  USER  CPU TIME (SECS)  |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
    while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
       shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
	      set value = `grep "LASVD EXECUTION TIME" $coutputfile | awk '{printf ("%-13s ", $5) }'`
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
	      set value = `grep "SISVD EXECUTION TIME" $coutputfile | awk '{printf ("%-13s ", $5) }'`
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
	      set value = `grep "TOTAL CPU  TIME" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       else
          set value = `grep "BLSVD EXECUTION TIME" $coutputfile | awk '{printf ("%-13s ", $5) }'`
       endif
       if ($value !~ [0-9]*) then
        echo -n "|      -        " >> $output
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          grep "LASVD EXECUTION TIME" $coutputfile | awk '{printf ("| %-13s ", $5) >> "'$output'" }'
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
          grep "SISVD EXECUTION TIME" $coutputfile | awk '{printf ("| %-13s ", $5) >> "'$output'" }'
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
          grep "TOTAL CPU  TIME" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else
          grep "BLSVD EXECUTION TIME" $coutputfile | awk '{printf ("| %-13s ", $5) >> "'$output'" }'
       endif
       @ c++
       shift programfiles
    end
    echo "|" >> $output


    echo -n "| NO. OF TRIPLETS FOUND   |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
    while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
       shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          set value = `grep "NSIG" $coutputfile | awk '{printf ("%-13s ", $4) }'`
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
	   	  set value = `grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("% -13s ", $7) }'`
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
	   	  set value = `grep "NO. OF APPROX. EIGENPAIRS" $coutputfile | awk '{printf ("% -13s ", $7) }'`
	   else
          set value = `grep "NO. OF TRIPLETS FOUND" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       endif
       if ($value !~ [0-9]*) then
          echo -n "|      -        " >> $output
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          grep "NSIG" $coutputfile | awk '{printf ("| %-13s ", $4) >> "'$output'" }'
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
    	  grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
    	  grep "NO. OF APPROX. EIGENPAIRS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
	   else
          grep "NO. OF TRIPLETS FOUND" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       endif
       @ c++
       shift programfiles
    end
    echo "|" >> $output


    echo -n "| NO. OF TRIPLETS SOUGHT  |" >> $output
    @ c = 1 + ($count * 5)
	set programfiles = `echo $programlist`
    while ($#programfiles > 0)
       set coutputfile = $programfiles[1].out$id[$c]
       shift programfiles
       if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
          set value = `grep "NO. OF EIGENPAIRS" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
          set value = `grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("%-13s ", $7) }'`
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
	   	  set value = `grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("% -13s ", $7) }'`
       else
          set value = `grep "NO. OF TRIPLETS SOUGHT" $coutputfile | awk '{printf ("%-13s ", $7) }'`
	   endif
       if ($value !~ [0-9]*) then
          set singularerror = `grep -c 'NUMBER OF SINGULAR VALUES DESIRED IS INVALID' $coutputfile`
          set initblockerror = `grep -c 'INIT BLOCK SIZE MUST BE LESS THAN NO. OF S-VALUES DESIRED' $coutputfile`
          if ($singularerror == 1) then
		     echo -n "|   ERROR 1    " >> $output
			 set error1 = 1
		     set singularerror = 0
          else if ($initblockerror == 1) then
             echo -n "|   ERROR 3    " >> $output
			 set error3 = 1
			 set initblockerror = 0
          else
             echo -n "|      -        " >> $output
          endif
       else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
	      grep "NO. OF EIGENPAIRS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
	      grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
    	  grep "NO. OF DESIRED EIGENPAIRS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       else
	      grep "NO. OF TRIPLETS SOUGHT" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
       endif
	   @ c++
	   shift programfiles
	end
	echo "|" >> $output


    echo -n "| NO. MULTIPLICATIONS BY A|" >> $output
    @ c = 1 + ($count * 5)
   set programfiles = `echo $programlist`
   while ($#programfiles > 0)
      set coutputfile = $programfiles[1].out$id[$c]
      shift programfiles
      if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
         set value = `grep "NO. MULTIPLICATIONS" $coutputfile | awk '{printf ("%-13s ", $7) }'`
      else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
         set value = `grep "MXV" $coutputfile | awk '{printf ("%-13s ", $4) }'`
      else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
         set value = `grep "MULTIPLICATIONS BY A " $coutputfile | awk '{printf ("%-13s ", $6) }'`
      else
         set value = `grep "NO. MULTIPLICATIONS" $coutputfile | awk '{printf ("%-13s ", $7) }'`
      endif
      if ($value !~ [0-9]*) then
         echo -n "|      -        " >> $output
      else if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
         grep "NO. MULTIPLICATIONS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
      else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
         grep "MXV" $coutputfile | awk '{printf ("| %-13s ", $4) >> "'$output'" }'
      else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
         grep "MULTIPLICATIONS BY A " $coutputfile | awk '{printf ("| %-13s ", $6) >> "'$output'" }'
      else
         grep "NO. MULTIPLICATIONS" $coutputfile | awk '{printf ("| %-13s ", $7) >> "'$output'" }'
      endif
      @ c++
	  shift programfiles
    end
    echo "|" >> $output
 

    echo -n "|NO. MULT. BY TRANSPOSE(A)|" >> $output
    @ c = 1 + ($count * 5)
   set programfiles = `echo $programlist`
   while ($#programfiles > 0)
      set coutputfile = $programfiles[1].out$id[$c]
      shift programfiles
      if ($programfiles[1] =~ las1 || $programfiles[1] =~ las2) then
         set value = `grep "NO. MULT. BY TRANSPOSE" $coutputfile | awk '{printf ("%-13s ", $7) }'`
      else if ($programfiles[1] =~ sis1 || $programfiles[1] =~ sis2) then
         set value = `grep "MXV" $coutputfile | awk '{printf ("%-13s ", $4) }'`
      else if ($programfiles[1] =~ tms1 || $programfiles[1] =~ tms2) then
         set value = `grep "MULTIPLICATIONS BY A^T" $coutputfile | awk '{printf ("%-13s ", $6) }'`
      else
         set value = `grep "NO. MULT. BY TRANSPOSE" $coutputfile | awk '{printf ("%-13s ", $7) }'`
      endif
      if ($value !~ [0-9]*) then