faasm: add includes to make kernels run with WAMR

Author: csegarragonz

MPI1/Nstream/nstream.c

@@ -1,32 +1,32 @@
 /*
 Copyright (c) 2013, Intel Corporation
 
-Redistribution and use in source and binary forms, with or without 
-modification, are permitted provided that the following conditions 
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions
 are met:
 
-* Redistributions of source code must retain the above copyright 
+* Redistributions of source code must retain the above copyright
       notice, this list of conditions and the following disclaimer.
-* Redistributions in binary form must reproduce the above 
-      copyright notice, this list of conditions and the following 
-      disclaimer in the documentation and/or other materials provided 
+* Redistributions in binary form must reproduce the above
+      copyright notice, this list of conditions and the following
+      disclaimer in the documentation and/or other materials provided
       with the distribution.
-* Neither the name of Intel Corporation nor the names of its 
-      contributors may be used to endorse or promote products 
-      derived from this software without specific prior written 
+* Neither the name of Intel Corporation nor the names of its
+      contributors may be used to endorse or promote products
+      derived from this software without specific prior written
       permission.
 
-THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 
-"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 
-LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 
-FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE 
-COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, 
-INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, 
-BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; 
-LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER 
-CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 
-LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN 
-ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
+FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
+COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
+INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
+BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
+LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
+LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
+ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 POSSIBILITY OF SUCH DAMAGE.
 */
 
@@ -67,44 +67,44 @@ POSSIBILITY OF SUCH DAMAGE.
 
 
 /**********************************************************************
- 
+
 NAME:      nstream
- 
+
 PURPOSE:   To compute memory bandwidth when adding a vector of a given
-           number of double precision values to the scalar multiple of 
+           number of double precision values to the scalar multiple of
            another vector of the same length, and storing the result in
-           a third vector. 
- 
+           a third vector.
+
 USAGE:     The program takes as input the number of iterations to loop
-           over the triad vectors, the length of the vectors, and the 
+           over the triad vectors, the length of the vectors, and the
            offset between vectors.
- 
+
            <progname> <# iterations> <vector length> <offset>
- 
-           The output consists of diagnostics to make sure the 
+
+           The output consists of diagnostics to make sure the
            algorithm worked, and of timing statistics.
- 
+
 FUNCTIONS CALLED:
- 
-           Other than MPI or standard C functions, the following 
+
+           Other than MPI or standard C functions, the following
            external functions are used in this program:
- 
+
            wtime()
            bail_out()
            checkTRIADresults()
- 
-NOTES:     Bandwidth is determined as the number of words read, plus the 
-           number of words written, times the size of the words, divided 
-           by the execution time. For a vector length of N, the total 
+
+NOTES:     Bandwidth is determined as the number of words read, plus the
+           number of words written, times the size of the words, divided
+           by the execution time. For a vector length of N, the total
            number of words read and written is 4*N*sizeof(double).
- 
+
 HISTORY:   This code is loosely based on the Stream benchmark by John
            McCalpin, but does not follow all the Stream rules. Hence,
            reported results should not be associated with Stream in
            external publications
-REVISION:  Modified by Rob Van der Wijngaart, December 2005, to     
+REVISION:  Modified by Rob Van der Wijngaart, December 2005, to
            parameterize vector size and offsets through compiler flags.
-           Also removed all Stream cases except TRIAD.                 
+           Also removed all Stream cases except TRIAD.
 REVISION:  Modified by Rob Van der Wijngaart, March 2006, to handle MPI.
 REVISION:  Modified by Rob Van der Wijngaart, May 2006, to introduce
            dependence between successive triad operations. This is
@@ -115,16 +115,23 @@ REVISION:  Modified by Rob Van der Wijngaart, November 2014, replaced
            between iterations (must now be included in timing) with
            accumulation: a[] += b[] + scalar*c[]
 **********************************************************************/
- 
+
 #include <par-res-kern_general.h>
 #include <par-res-kern_mpi.h>
- 
+
+#ifdef __faasm
+#include <faasm/faasm.h>
+#endif
+
 #define SCALAR  3.0
- 
+
 static int checkTRIADresults(int, long int, double *);
- 
-int main(int argc, char **argv) 
+
+int main(int argc, char **argv)
 {
+  const char* message = "Hello Faasm!";
+  faasmSetOutput((uint8_t*)message, 12);
+
   long int j, iter;       /* dummies                                     */
   double   scalar;        /* constant used in Triad operation            */
   int      iterations;    /* number of times vector loop gets repeated   */
@@ -134,7 +141,7 @@ int main(int argc, char **argv)
   double   bytes;         /* memory IO size                              */
   size_t   space;         /* memory used for a single vector             */
   double   local_nstream_time,/* timing parameters                       */
-           nstream_time, 
+           nstream_time,
            avgtime;
   int      Num_procs,     /* number of ranks                             */
            my_ID,         /* rank                                        */
@@ -143,11 +150,11 @@ int main(int argc, char **argv)
   double * RESTRICT a;    /* main vector                                 */
   double * RESTRICT b;    /* main vector                                 */
   double * RESTRICT c;    /* main vector                                 */
- 
+
 /**********************************************************************************
-* process and test input parameters    
+* process and test input parameters
 ***********************************************************************************/
- 
+
   MPI_Init(&argc,&argv);
   MPI_Comm_size(MPI_COMM_WORLD,&Num_procs);
   MPI_Comm_rank(MPI_COMM_WORLD,&my_ID);
@@ -205,7 +212,7 @@ int main(int argc, char **argv)
   c = b + length + offset;
 
   bytes   = 4.0 * sizeof(double) * length * Num_procs;
- 
+
   if (my_ID == root) {
     printf("Number of ranks      = %d\n", Num_procs);
     printf("Vector length        = %ld\n", total_length);
@@ -218,31 +225,31 @@ int main(int argc, char **argv)
     b[j] = 2.0;
     c[j] = 2.0;
   }
-    
+
   /* --- MAIN LOOP --- repeat Triad iterations times --- */
- 
+
   scalar = SCALAR;
- 
+
   for (iter=0; iter<=iterations; iter++) {
- 
+
     /* start timer after a warmup iteration */
-    if (iter == 1) { 
+    if (iter == 1) {
       MPI_Barrier(MPI_COMM_WORLD);
       local_nstream_time = wtime();
     }
 
     for (j=0; j<length; j++) a[j] += b[j]+scalar*c[j];
 
   } /* end iterations */
- 
+
   /*********************************************************************
   ** Analyze and output results.
   *********************************************************************/
 
   local_nstream_time = wtime() - local_nstream_time;
   MPI_Reduce(&local_nstream_time, &nstream_time, 1, MPI_DOUBLE, MPI_MAX, root,
              MPI_COMM_WORLD);
-  
+
   if (my_ID == root) {
     if (checkTRIADresults(iterations, length, a)) {
       avgtime = nstream_time/iterations;
@@ -254,34 +261,34 @@ int main(int argc, char **argv)
   bail_out(error);
   MPI_Finalize();
 }
- 
- 
+
+
 int checkTRIADresults (int iterations, long int length, double *a) {
   double aj, bj, cj, scalar, asum;
   double epsilon = 1.e-8;
   long int j;
   int iter;
- 
+
   /* reproduce initialization */
   aj = 0.0;
   bj = 2.0;
   cj = 2.0;
- 
+
   /* now execute timing loop */
   scalar = SCALAR;
   for (iter=0; iter<=iterations; iter++) aj += bj+scalar*cj;
- 
+
   aj = aj * (double) (length);
- 
+
   asum = 0.0;
   for (j=0; j<length; j++) asum += a[j];
- 
+
 #if VERBOSE
   printf ("Results Comparison: \n");
   printf ("        Expected checksum: %f\n",aj);
   printf ("        Observed checksum: %f\n",asum);
 #endif
- 
+
   if (ABS(aj-asum)/asum > epsilon) {
     printf ("Failed Validation on output array\n");
 #if !VERBOSE
@@ -295,4 +302,4 @@ int checkTRIADresults (int iterations, long int length, double *a) {
     return (1);
   }
 }
- 
+