Subversion Repositories pavr

// Project
//    pAVR (pipelined AVR) is an 8 bit RISC controller, compatible with Atmel's
//    AVR core, but about 3x faster in terms of both clock frequency and MIPS.
//    The increase in speed comes from a relatively deep pipeline. The original
//    AVR core has only two pipeline stages (fetch and execute), while pAVR has
//    6 pipeline stages:
//       1. PM    (read Program Memory)
//       2. INSTR (load Instruction)
//       3. RFRD  (decode Instruction and read Register File)
//       4. OPS   (load Operands)
//       5. ALU   (execute ALU opcode or access Unified Memory)
//       6. RFWR  (write Register File)
// Version
//    0.32
// Date
//    2002 August 07
// Author
//    Doru Cuturela, doruu@yahoo.com
// License
//    This program is free software; you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation; either version 2 of the License, or
//    (at your option) any later version.
//    This program is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//    You should have received a copy of the GNU General Public License
//    along with this program; if not, write to the Free Software
//    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
 
 
// About this file...
//    This is a pAVR test that simulates waves on the surface of a liquid.
//    This uses floating point numbers.
//    Checking the result is done by simulating this program on pAVR, and
//       comparing the final mesh (read directly from Data Memory) against
//       the reference final mesh obtained by running this program on a PC.
//       To switch between the two situations, #define or not REF (see below).
 
 
 
//#include <stdio.h>
//#include <stdlib.h>
//#include <string.h>
//#include <sfr2313.h>
//#include <sfr8515.h>
#include <float.h>
 
#define NR_OX 5
#define NR_OY 5
#define NR_IT 5
 
// Define this when running reference test on host PC.
// Don't define this when running pAVR test simulation.
#define REF
 
 
#ifdef REF
   #include <stdio.h>
   #include <stdlib.h>
   #include <string.h>
#endif
 
void main() {
   #ifdef REF
   FILE *f1;
   char txt[100];
   #endif
 
   static float zData[NR_OX][NR_OY], vzData[NR_OX][NR_OY];
   static float tmp1, tmp2, tmp3, tmp4;
   static char cData[NR_OX][NR_OY];
   int ix, iy, it;
 
   for (ix=0; ix<NR_OX; ix++) {
      for (iy=0; iy<NR_OY; iy++) {
         zData[ix][iy] = 0.0;
         vzData[ix][iy] = 0.0;
      }
   }
   zData[1][1]=0.9;
 
//*
   for (it=0; it<NR_IT; it++) {
      for (ix=0; ix<NR_OX; ix++) {
         for (iy=0; iy<NR_OY; iy++) {
            if ((ix==0)||(iy==0)||(ix==NR_OX-1)||(iy==NR_OY-1)) {
               // Boundary conditions
               zData[ix][iy] = 0.0;
               vzData[ix][iy] = 0.0;
            }
            else {
               vzData[ix][iy] = 0.90*vzData[ix][iy] - 0.0333*(6.0*zData[ix][iy]
                                                              - (zData[ix-1][iy] + zData[ix+1][iy] + zData[ix][iy-1] + zData[ix][iy+1])
                                                              - 0.5 * (zData[ix-1][iy-1] + zData[ix-1][iy+1] + zData[ix+1][iy-1] + zData[ix+1][iy+1]));
               zData[ix][iy] = zData[ix][iy] + vzData[ix][iy];
 
            }
         }
      }
   }
//*/
 
   for (ix=0; ix<NR_OX; ix++) {
      for (iy=0; iy<NR_OY; iy++) {
         cData[ix][iy] = (char) (zData[ix][iy]*128.0);
//         if (zData[ix][iy]>=1) cData[ix][iy] =  127;
//         if (zData[ix][iy]<-1) cData[ix][iy] = -128;
      }
   }
 
 
   #ifdef REF
   // Write results
   f1 = fopen("results.dat", "wb+");
   if (f1!=NULL) {
      for (ix=0; ix<NR_OX; ix++) {
         for (iy=0; iy<NR_OY; iy++) {
            itoa(cData[ix][iy], txt, 10);
            strcat(txt, "\n");
            fwrite(txt, strlen(txt), 1, f1);
            fflush(f1);
         }
      }
   }
   else {
      printf("Could not write file.\n");
   }
   fclose(f1);
   #endif
 
 
 
 
 
//*/
   #ifndef REF
   while (1) {
   }
   #endif
}