Page: 23 prev next



Sample Code for Computing the Mean



Original Fortran Code:
CMP+ PRECISION LEVEL 500
CMP+ MIXED MODE SAFE
      program reg
      parameter (mdat=100, mpar=10)
CMP+  IMPLICIT MULTIP REAL (A-H, O-Z)
             .
             .
             .
      read (10,*) ndat
      do 10 i=1, ndat
         read (10,*) x(i)
         mean = mean +x(i)
 10   continue
      mean = mean /ndat

Code From Multiple-Precision Translator
      read (10,*) ndat
      do 10 i=1, ndat
         CALL MPINP (10, x(1,i), MPA1)
         CALL MPADD (mean, x(1,i), MPM1)
         CALL MPEQ (MPM1, mean)
 10   continue
      MPD1 = ndat
      CALL MPDIVD (mean, MPD1, 0, MPM1)
      CALL MPEQ (MPM1, mean)